97 class(VirtualGwfModelType) :: this
100 integer(I4B) :: nr_nodes
103 call this%VirtualModelType%prepare_stage(stage)
105 if (stage == stg_aft_mdl_df)
then
107 call this%map(this%npf_iangle1%base(), (/stg_aft_mdl_df/))
108 call this%map(this%npf_iangle2%base(), (/stg_aft_mdl_df/))
109 call this%map(this%npf_iangle3%base(), (/stg_aft_mdl_df/))
110 call this%map(this%npf_iwetdry%base(), (/stg_aft_mdl_df/))
111 call this%map(this%inbuy%base(), (/stg_aft_mdl_df/))
113 else if (stage == stg_bfr_con_ar)
then
115 nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
117 call this%map(this%x%base(), nr_nodes, &
118 (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
119 call this%map(this%ibound%base(), nr_nodes, &
120 (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
121 call this%map(this%x_old%base(), nr_nodes, &
122 (/stg_bfr_exg_ad, stg_bfr_exg_cf/))
124 call this%map(this%npf_icelltype%base(), nr_nodes, (/stg_bfr_con_ar/))
125 call this%map(this%npf_k11%base(), nr_nodes, (/stg_bfr_con_ar/))
126 call this%map(this%npf_k22%base(), nr_nodes, (/stg_bfr_con_ar/))
127 call this%map(this%npf_k33%base(), nr_nodes, (/stg_bfr_con_ar/))
129 if (this%npf_iangle1%get() > 0)
then
130 call this%map(this%npf_angle1%base(), nr_nodes, (/stg_bfr_con_ar/))
132 call this%map(this%npf_angle1%base(), 0, (/stg_never/))
134 if (this%npf_iangle2%get() > 0)
then
135 call this%map(this%npf_angle2%base(), nr_nodes, (/stg_bfr_con_ar/))
137 call this%map(this%npf_angle2%base(), 0, (/stg_never/))
139 if (this%npf_iangle3%get() > 0)
then
140 call this%map(this%npf_angle3%base(), nr_nodes, (/stg_bfr_con_ar/))
142 call this%map(this%npf_angle3%base(), 0, (/stg_never/))
144 if (this%npf_iwetdry%get() > 0)
then
145 call this%map(this%npf_wetdry%base(), nr_nodes, (/stg_bfr_con_ar/))
147 call this%map(this%npf_wetdry%base(), 0, (/stg_never/))
150 if (this%inbuy%get() > 0)
then
151 call this%map(this%buy_dense%base(), nr_nodes, (/stg_bfr_exg_cf/))
153 call this%map(this%buy_dense%base(), 0, (/stg_never/))