virtualgwtmodelmodule Module Reference

Data Types
type	virtualgwtmodeltype

Functions/Subroutines
subroutine, public	add_virtual_gwt_model (model_id, model_name, model)
subroutine	vgwt_create (this, name, id, model)
subroutine	init_virtual_data (this)
subroutine	vgwt_prepare_stage (this, stage)
subroutine	allocate_data (this)
subroutine	deallocate_data (this)
subroutine	vgwt_destroy (this)

Function/Subroutine Documentation

add_virtual_gwt_model()

subroutine, public virtualgwtmodelmodule::add_virtual_gwt_model	(	integer(i4b)	model_id,
		character(len=*)	model_name,
		class(numericalmodeltype), pointer	model
	)

Parameters

model_id	global model id
model_name	model name
model	the actual model (can be null() when remote)

Definition at line 47 of file VirtualGwtModel.f90.

    use virtualdatalistsmodule, only: virtual_model_list
    integer(I4B) :: model_id !< global model id
    character(len=*) :: model_name !< model name
    class(NumericalModelType), pointer :: model !< the actual model (can be null() when remote)
    ! local
    class(VirtualGwtModelType), pointer :: virtual_gwt_model
    class(*), pointer :: obj
 
    allocate (virtual_gwt_model)
    call virtual_gwt_model%create(model_name, model_id, model)
 
    obj => virtual_gwt_model
    call virtual_model_list%Add(obj)
 

Here is the caller graph for this function:

allocate_data()

subroutine virtualgwtmodelmodule::allocate_data ( class(virtualgwtmodeltype)  this )

private

Definition at line 164 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
 
    allocate (this%dsp_idiffc)
    allocate (this%dsp_idisp)
    allocate (this%dsp_diffc)
    allocate (this%dsp_alh)
    allocate (this%dsp_alv)
    allocate (this%dsp_ath1)
    allocate (this%dsp_ath2)
    allocate (this%dsp_atv)
    allocate (this%fmi_igwfspdis)
    allocate (this%fmi_gwfhead)
    allocate (this%fmi_gwfsat)
    allocate (this%fmi_gwfspdis)
    allocate (this%fmi_gwfflowja)
    allocate (this%mst_thetam)
    allocate (this%indsp)
    allocate (this%inmst)
 

deallocate_data()

subroutine virtualgwtmodelmodule::deallocate_data ( class(virtualgwtmodeltype)  this )

private

Definition at line 186 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
 
    deallocate (this%dsp_idiffc)
    deallocate (this%dsp_idisp)
    deallocate (this%dsp_diffc)
    deallocate (this%dsp_alh)
    deallocate (this%dsp_alv)
    deallocate (this%dsp_ath1)
    deallocate (this%dsp_ath2)
    deallocate (this%dsp_atv)
    deallocate (this%fmi_igwfspdis)
    deallocate (this%fmi_gwfhead)
    deallocate (this%fmi_gwfsat)
    deallocate (this%fmi_gwfspdis)
    deallocate (this%fmi_gwfflowja)
    deallocate (this%mst_thetam)
    deallocate (this%indsp)
    deallocate (this%inmst)
 

init_virtual_data()

subroutine virtualgwtmodelmodule::init_virtual_data ( class(virtualgwtmodeltype)  this )

private

Definition at line 79 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
 
    call this%set(this%dsp_idiffc%base(), 'IDIFFC', 'DSP', map_all_type)
    call this%set(this%dsp_idisp%base(), 'IDISP', 'DSP', map_all_type)
    call this%set(this%dsp_diffc%base(), 'DIFFC', 'DSP', map_node_type)
    call this%set(this%dsp_alh%base(), 'ALH', 'DSP', map_node_type)
    call this%set(this%dsp_alv%base(), 'ALV', 'DSP', map_node_type)
    call this%set(this%dsp_ath1%base(), 'ATH1', 'DSP', map_node_type)
    call this%set(this%dsp_ath2%base(), 'ATH2', 'DSP', map_node_type)
    call this%set(this%dsp_atv%base(), 'ATV', 'DSP', map_node_type)
    call this%set(this%fmi_igwfspdis%base(), 'IGWFSPDIS', 'FMI', map_all_type)
    call this%set(this%fmi_gwfhead%base(), 'GWFHEAD', 'FMI', map_node_type)
    call this%set(this%fmi_gwfsat%base(), 'GWFSAT', 'FMI', map_node_type)
    call this%set(this%fmi_gwfspdis%base(), 'GWFSPDIS', 'FMI', map_node_type)
    call this%set(this%fmi_gwfflowja%base(), 'GWFFLOWJA', 'FMI', map_conn_type)
    call this%set(this%mst_thetam%base(), 'THETAM', 'MST', map_node_type)
    call this%set(this%indsp%base(), 'INDSP', '', map_all_type)
    call this%set(this%inmst%base(), 'INMST', '', map_all_type)
 

vgwt_create()

subroutine virtualgwtmodelmodule::vgwt_create	(	class(virtualgwtmodeltype)	this,
		character(len=*)	name,
		integer(i4b)	id,
		class(numericalmodeltype), pointer	model
	)

Definition at line 64 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
    character(len=*) :: name
    integer(I4B) :: id
    class(NumericalModelType), pointer :: model
 
    ! create base
    call this%VirtualModelType%create(name, id, model)
    this%container_type = vdc_gwtmodel_type
 
    call this%allocate_data()
    call this%init_virtual_data()
 

vgwt_destroy()

subroutine virtualgwtmodelmodule::vgwt_destroy ( class(virtualgwtmodeltype)  this )

private

Definition at line 208 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
 
    call this%VirtualModelType%destroy()
    call this%deallocate_data()
 

vgwt_prepare_stage()

subroutine virtualgwtmodelmodule::vgwt_prepare_stage ( class(virtualgwtmodeltype)  this, integer(i4b)  stage  )

private

Definition at line 101 of file VirtualGwtModel.f90.

    class(VirtualGwtModelType) :: this
    integer(I4B) :: stage
    ! local
    integer(I4B) :: nr_nodes, nr_conns
 
    ! prepare base (=numerical) model data items
    call this%VirtualModelType%prepare_stage(stage)
 
    nr_nodes = 0
    nr_conns = 0
 
    if (stage == stg_aft_mdl_df) then
 
      call this%map(this%dsp_idiffc%base(), (/stg_aft_mdl_df/))
      call this%map(this%dsp_idisp%base(), (/stg_aft_mdl_df/))
      call this%map(this%indsp%base(), (/stg_aft_mdl_df/))
      call this%map(this%inmst%base(), (/stg_aft_mdl_df/))
 
    else if (stage == stg_bfr_con_ar) then
 
      nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
      nr_conns = this%element_maps(map_conn_type)%nr_virt_elems
 
      call this%map(this%fmi_igwfspdis%base(), (/stg_bfr_con_ar/))
      call this%map(this%x%base(), nr_nodes, &
                    (/stg_bfr_con_ar, stg_bfr_exg_ad, stg_bfr_exg_cf/))
      call this%map(this%ibound%base(), nr_nodes, (/stg_bfr_con_ar/))
 
      if (this%dsp_idiffc%get() > 0) then
        call this%map(this%dsp_diffc%base(), nr_nodes, (/stg_bfr_con_ar/))
      end if
 
      if (this%dsp_idisp%get() > 0) then
        call this%map(this%dsp_alh%base(), nr_nodes, (/stg_bfr_con_ar/))
        call this%map(this%dsp_alv%base(), nr_nodes, (/stg_bfr_con_ar/))
        call this%map(this%dsp_ath1%base(), nr_nodes, (/stg_bfr_con_ar/))
        call this%map(this%dsp_ath2%base(), nr_nodes, (/stg_bfr_con_ar/))
        call this%map(this%dsp_atv%base(), nr_nodes, (/stg_bfr_con_ar/))
      end if
 
      if (this%indsp%get() > 0 .and. this%inmst%get() > 0) then
        call this%map(this%mst_thetam%base(), nr_nodes, (/stg_aft_con_ar/))
      end if
 
    else if (stage == stg_aft_con_ar) then
 
      nr_nodes = this%element_maps(map_node_type)%nr_virt_elems
      nr_conns = this%element_maps(map_conn_type)%nr_virt_elems
 
      call this%map(this%fmi_gwfhead%base(), nr_nodes, (/stg_bfr_exg_ad/))
      call this%map(this%fmi_gwfsat%base(), nr_nodes, (/stg_bfr_exg_ad/))
      if (this%fmi_igwfspdis%get() > 0) then
        call this%map(this%fmi_gwfspdis%base(), 3, nr_nodes, (/stg_bfr_exg_ad/))
      else
        call this%map(this%fmi_gwfspdis%base(), 3, 0, (/stg_never/))
      end if
      call this%map(this%fmi_gwfflowja%base(), nr_conns, (/stg_bfr_exg_ad/))
 
    end if