117 class(PrtMipType),
intent(inout) :: this
119 character(len=LINELENGTH) :: errmsg
120 type(PrtMipParamFoundType) :: found
121 integer(I4B),
dimension(:),
pointer,
contiguous :: map => null()
124 if (this%dis%nodes < this%dis%nodesuser) map => this%dis%nodeuser
127 call this%allocate_arrays(this%dis%nodes)
130 call mem_set_value(this%porosity,
'POROSITY', this%input_mempath, &
132 call mem_set_value(this%retfactor,
'RETFACTOR', this%input_mempath, &
133 map, found%retfactor)
134 call mem_set_value(this%izone,
'IZONE', this%input_mempath, map, &
138 call mem_set_value(this%zeromethod,
'ZERO_METHOD', this%input_mempath, &
140 if (.not. found%zero_method) this%zeromethod = 1
143 if (.not. found%porosity)
then
144 write (errmsg,
'(a)')
'Error in GRIDDATA block: POROSITY not found'
145 call store_error(errmsg)