104 class(VirtualGweModelType) :: this
105 integer(I4B) :: stage
107 integer(I4B) :: nr_nodes, nr_conns
110 call this%VirtualModelType%prepare_stage(stage)
115 if (stage == stg_aft_mdl_df)
then
118 call this%map(this%cnd_idisp%base(), (/stg_aft_mdl_df/))
119 call this%map(this%incnd%base(), (/stg_aft_mdl_df/))
120 call this%map(this%inest%base(), (/stg_aft_mdl_df/))
122 else if (stage == stg_bfr_con_ar)
then
124 nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
125 nr_conns = this%element_maps(map_conn_type)%nr_virt_elems
127 call this%map(this%x%base(), nr_nodes, &
128 (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
129 call this%map(this%ibound%base(), nr_nodes, (/stg_bfr_con_ar/))
135 if (this%cnd_idisp%get() > 0)
then
136 call this%map(this%cnd_alh%base(), nr_nodes, (/stg_bfr_con_ar/))
137 call this%map(this%cnd_alv%base(), nr_nodes, (/stg_bfr_con_ar/))
138 call this%map(this%cnd_ath1%base(), nr_nodes, (/stg_bfr_con_ar/))
139 call this%map(this%cnd_ath2%base(), nr_nodes, (/stg_bfr_con_ar/))
140 call this%map(this%cnd_atv%base(), nr_nodes, (/stg_bfr_con_ar/))
141 call this%map(this%cnd_ktw%base(), nr_nodes, (/stg_bfr_con_ar/))
142 call this%map(this%cnd_kts%base(), nr_nodes, (/stg_bfr_con_ar/))
145 call this%map(this%fmi_gwfhead%base(), nr_nodes, (/stg_bfr_exg_ad/))
146 call this%map(this%fmi_gwfsat%base(), nr_nodes, (/stg_bfr_exg_ad/))
147 call this%map(this%fmi_gwfspdis%base(), 3, nr_nodes, (/stg_bfr_exg_ad/))
148 call this%map(this%fmi_gwfflowja%base(), nr_conns, (/stg_bfr_exg_ad/))
150 if (this%incnd%get() > 0 .and. this%inest%get() > 0)
then
151 call this%map(this%est_porosity%base(), nr_nodes, (/stg_aft_con_ar/))