102 class(VirtualGweModelType) :: this
103 integer(I4B) :: stage
105 integer(I4B) :: nr_nodes, nr_conns
108 call this%VirtualModelType%prepare_stage(stage)
113 if (stage == stg_aft_mdl_df)
then
115 call this%map(this%cnd_idisp%base(), (/stg_aft_mdl_df/))
116 call this%map(this%incnd%base(), (/stg_aft_mdl_df/))
117 call this%map(this%inest%base(), (/stg_aft_mdl_df/))
119 else if (stage == stg_bfr_con_ar)
then
121 nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
122 nr_conns = this%element_maps(map_conn_type)%nr_virt_elems
124 call this%map(this%x%base(), nr_nodes, &
125 (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
126 call this%map(this%ibound%base(), nr_nodes, (/stg_bfr_con_ar/))
128 call this%map(this%fmi_igwfspdis%base(), (/stg_bfr_con_ar/))
130 if (this%cnd_idisp%get() > 0)
then
131 call this%map(this%cnd_alh%base(), nr_nodes, (/stg_bfr_con_ar/))
132 call this%map(this%cnd_alv%base(), nr_nodes, (/stg_bfr_con_ar/))
133 call this%map(this%cnd_ath1%base(), nr_nodes, (/stg_bfr_con_ar/))
134 call this%map(this%cnd_ath2%base(), nr_nodes, (/stg_bfr_con_ar/))
135 call this%map(this%cnd_atv%base(), nr_nodes, (/stg_bfr_con_ar/))
136 call this%map(this%cnd_ktw%base(), nr_nodes, (/stg_bfr_con_ar/))
137 call this%map(this%cnd_kts%base(), nr_nodes, (/stg_bfr_con_ar/))
140 if (this%incnd%get() > 0 .and. this%inest%get() > 0)
then
141 call this%map(this%est_porosity%base(), nr_nodes, (/stg_aft_con_ar/))
144 else if (stage == stg_aft_con_ar)
then
146 nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
147 nr_conns = this%element_maps(map_conn_type)%nr_virt_elems
149 call this%map(this%fmi_gwfhead%base(), nr_nodes, (/stg_bfr_exg_ad/))
150 call this%map(this%fmi_gwfsat%base(), nr_nodes, (/stg_bfr_exg_ad/))
151 if (this%fmi_igwfspdis%get() > 0)
then
152 call this%map(this%fmi_gwfspdis%base(), 3, nr_nodes, (/stg_bfr_exg_ad/))
154 call this%map(this%fmi_gwfspdis%base(), 3, 0, (/stg_never/))
156 call this%map(this%fmi_gwfflowja%base(), nr_conns, (/stg_bfr_exg_ad/))