mod_atmos_phy_tb_vars Module Reference

module ATMOSPHERE physics / sub-grid scale turbulence More...

Data Types
type	atmosphytbvars

Functions/Subroutines
subroutine	atmosphytbvars_init (this, model_mesh)
subroutine, public	atmosphytbvars_getlocalmeshfields_tend (domid, mesh, tb_tends_list, tb_momx_t, tb_momy_t, tb_momz_t, tb_rhot_t, tb_rhoq_t, lcmesh3d)

Detailed Description

module ATMOSPHERE physics / sub-grid scale turbulence

Description

Container for variables with sub-grid scale turbulence model

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

atmosphytbvars_init()

subroutine mod_atmos_phy_tb_vars::atmosphytbvars_init	(	class(atmosphytbvars), intent(inout), target	this,
		class(modelmeshbase), intent(in), target	model_mesh )

Definition at line 91 of file mod_atmos_phy_tb_vars.F90.

 
    use scale_tracer, only: &
      tracer_name, tracer_desc, tracer_unit    
    use scale_atm_phy_tb_dgm_common, only: &
      atm_phy_tb_dgm_common_setup_variables
    
    implicit none
    class(AtmosPhyTbVars), target, intent(inout) :: this
    class(ModelMeshBase), target, intent(in) :: model_mesh
 
    integer :: iv
    integer :: iq
    logical :: reg_file_hist
 
    class(AtmosMesh), pointer :: atm_mesh
    class(MeshBase2D), pointer :: mesh2D
    class(MeshBase3D), pointer :: mesh3D
    !--------------------------------------------------
 
    log_info('AtmosPhyTbVars_Init',*)
 
    this%TENDS_NUM_TOT = atmos_phy_tb_tends_num1 + qa
 
    !- Initialize auxiliary and diagnostic variables
 
    nullify( atm_mesh )
    select type(model_mesh)
    class is (atmosmesh)
      atm_mesh => model_mesh
    end select
    mesh3d => atm_mesh%ptr_mesh
    
    call mesh3d%GetMesh2D( mesh2d )
 
    !- Get variable information
    ! call atm_phy_tb_dgm_common_get_varinfo( auxvar_info, diagvar_info, tbtend_info ) ! (out)
 
    !- Initialize objects to variables with turbulent model
 
    call this%tends_manager%Init()
    call this%auxvars_manager%Init()
    call this%diagvars_manager%Init()
 
    allocate( this%diagvars(atmos_phy_tb_diag_num) )
    allocate( this%tends(this%TENDS_NUM_TOT) )
    allocate( this%auxvars(atmos_phy_tb_aux_num) )
 
    call atm_phy_tb_dgm_common_setup_variables( &
      this%tends, this%auxvars, this%diagvars,                         & ! (inout)
      this%tends_manager, this%auxvars_manager, this%diagvars_manager, & ! (inout)
      this%TENDS_NUM_TOT, mesh3d )                                       ! (in)
 
    !- Setup communication
 
    call atm_mesh%Create_communicator( &
      atmos_phy_tb_aux_scalar_num, atmos_phy_tb_aux_hvec_num, & ! (in)
      atmos_phy_tb_aux_htensor_num,                           & ! (in)
      this%auxvars_manager,            & ! (inout)
      this%auxvars(:),                 & ! (in)
      this%auxvars_commid              ) ! (out)
 
    return

References scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_setup_variables(), scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_htensor_num, scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_hvec_num, scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_num, scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_scalar_num, scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_num, and scale_atm_phy_tb_dgm_common::atmos_phy_tb_tends_num1.

atmosphytbvars_getlocalmeshfields_tend()

subroutine, public mod_atmos_phy_tb_vars::atmosphytbvars_getlocalmeshfields_tend	(	integer, intent(in)	domid,
		class(meshbase), intent(in)	mesh,
		class(modelvarmanager), intent(inout)	tb_tends_list,
		class(localmeshfieldbase), intent(out), pointer	tb_momx_t,
		class(localmeshfieldbase), intent(out), pointer	tb_momy_t,
		class(localmeshfieldbase), intent(out), pointer	tb_momz_t,
		class(localmeshfieldbase), intent(out), pointer	tb_rhot_t,
		type(localmeshfieldbaselist), dimension(:), intent(out), optional	tb_rhoq_t,
		class(localmesh3d), intent(out), optional, pointer	lcmesh3d )

Definition at line 177 of file mod_atmos_phy_tb_vars.F90.

 
    use scale_mesh_base, only: meshbase
    use scale_meshfield_base, only: meshfieldbase
    implicit none
 
    integer, intent(in) :: domID
    class(MeshBase), intent(in) :: mesh
    class(ModelVarManager), intent(inout) :: tb_tends_list
    class(LocalMeshFieldBase), pointer, intent(out) :: tb_MOMX_t
    class(LocalMeshFieldBase), pointer, intent(out) :: tb_MOMY_t
    class(LocalMeshFieldBase), pointer, intent(out) :: tb_MOMZ_t
    class(LocalMeshFieldBase), pointer, intent(out) :: tb_RHOT_t
    type(LocalMeshFieldBaseList), intent(out), optional :: tb_RHOQ_t(:)
    class(LocalMesh3D), pointer, intent(out), optional :: lcmesh3D
 
    class(MeshFieldBase), pointer :: field   
    class(LocalMeshBase), pointer :: lcmesh
 
    integer :: iq
    !-------------------------------------------------------
 
    !--
    call tb_tends_list%Get(tb_momx_t_vid, field)
    call field%GetLocalMeshField(domid, tb_momx_t)
 
    call tb_tends_list%Get(tb_momy_t_vid, field)
    call field%GetLocalMeshField(domid, tb_momy_t)
 
    call tb_tends_list%Get(tb_momz_t_vid, field)
    call field%GetLocalMeshField(domid, tb_momz_t)
 
    call tb_tends_list%Get(tb_rhot_t_vid, field)
    call field%GetLocalMeshField(domid, tb_rhot_t)
 
    !---
    if ( present(tb_rhoq_t) ) then
      do iq = 1, size(tb_rhoq_t)
        call tb_tends_list%Get(atmos_phy_tb_tends_num1 + iq, field)
        call field%GetLocalMeshField(domid, tb_rhoq_t(iq)%ptr)
      end do    
    end if
 
    if (present(lcmesh3d)) then
      call mesh%GetLocalMesh( domid, lcmesh )
      nullify( lcmesh3d )
 
      select type(lcmesh)
      type is (localmesh3d)
        if (present(lcmesh3d)) lcmesh3d => lcmesh
      end select
    end if
 
    return

Referenced by mod_atmos_phy_tb::atmosphytb_setup(), and mod_atmos_phy_tb_vars::atmosphytbvars::history().