148 class(GweInterfaceModelType) :: this
150 class(*),
pointer :: disPtr
151 type(TspAdvOptionsType) :: adv_options
152 type(GweCndOptionsType) :: cnd_options
155 adv_options%iAdvScheme = this%iAdvScheme
156 cnd_options%ixt3d = this%ixt3d
160 call this%gridConnection%getDiscretization(castasdisutype(disptr))
161 call this%fmi%fmi_df(this%dis, 0)
163 if (this%inadv > 0)
then
164 call this%adv%adv_df(adv_options)
167 if (this%incnd > 0)
then
168 this%cnd%idisp = this%owner%cnd%idisp
169 call this%cnd%cnd_df(this%dis, cnd_options)
171 if (this%cnd%idisp > 0)
then
172 call mem_reallocate(this%cnd%alh, this%dis%nodes,
'ALH', &
173 trim(this%cnd%memoryPath))
174 call mem_reallocate(this%cnd%alv, this%dis%nodes,
'ALV', &
175 trim(this%cnd%memoryPath))
176 call mem_reallocate(this%cnd%ath1, this%dis%nodes,
'ATH1', &
177 trim(this%cnd%memoryPath))
178 call mem_reallocate(this%cnd%ath2, this%dis%nodes,
'ATH2', &
179 trim(this%cnd%memoryPath))
180 call mem_reallocate(this%cnd%atv, this%dis%nodes,
'ATV', &
181 trim(this%cnd%memoryPath))
182 call mem_reallocate(this%cnd%ktw, this%dis%nodes,
'KTW', &
183 trim(this%cnd%memoryPath))
184 call mem_reallocate(this%cnd%kts, this%dis%nodes,
'KTS', &
185 trim(this%cnd%memoryPath))
188 call mem_allocate(this%est%porosity, this%dis%nodes, &
189 'POROSITY', create_mem_path(this%name,
'EST'))
193 this%neq = this%dis%nodes
194 this%nja = this%dis%nja
195 this%ia => this%dis%con%ia
196 this%ja => this%dis%con%ja
199 call this%gwecommon%gweshared_dat_df(this%neq)
202 call this%allocate_arrays()