113 class(PrtMipType),
intent(inout) :: this
115 character(len=LINELENGTH) :: errmsg
116 type(PrtMipParamFoundType) :: found
117 integer(I4B),
dimension(:),
pointer,
contiguous :: map => null()
120 if (this%dis%nodes < this%dis%nodesuser) map => this%dis%nodeuser
123 call this%allocate_arrays(this%dis%nodes)
126 call mem_set_value(this%porosity,
'POROSITY', this%input_mempath, &
128 call mem_set_value(this%retfactor,
'RETFACTOR', this%input_mempath, &
129 map, found%retfactor)
130 call mem_set_value(this%izone,
'IZONE', this%input_mempath, map, &
134 if (.not. found%porosity)
then
135 write (errmsg,
'(a)')
'Error in GRIDDATA block: POROSITY not found'
136 call store_error(errmsg)