mod_atmos_phy_tb Module Reference

module Atmosphere / Physics / sub-grid scale turbulence process More...

Data Types
type	atmosphytb
	Derived type to manage a component of sub-grid scale turbulent process. More...

Functions/Subroutines
subroutine	atmosphytb_setup (this, model_mesh, tm_parent_comp)
	Setup a component of SGS turbulence process.

Detailed Description

module Atmosphere / Physics / sub-grid scale turbulence process

Description

Module for sub-grid scale turbulence process

Author

Yuta Kawai, Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_TB

  name	type	default value	comment
  TIME_DT	real(DP)	UNDEF8	Timestep for sub-grid scale turbulent process
  TIME_DT_UNIT	character(len=H_SHORT)	'SEC'	Unit of timestep
  TB_TYPE	character(len=H_MID)	'SMAGORINSKY'	Type of sub-grid scale turbulent scheme
  POSTGLFILTER_TYPE	character(len=H_SHORT)	'None'
  POSTGLFILTER_NNODEH1D_RECONST	integer	-1
  POSTGLFILTER_GAUSSINWIDTHFAC	real(RP)	1.5_RP
  POSTMODALFILTER_ALPHA_H	real(RP)	0.0_RP
  POSTMODALFILTER_ETAC_H	real(RP)	0.0_RP
  POSTMODALFILTER_ORDER_H	integer	16
  POSTMODALFILTER_ALPHA_V	real(RP)	0.0_RP
  POSTMODALFILTER_ORDER_V	integer	16

  
  

History Output

No history output

Function/Subroutine Documentation

atmosphytb_setup()

subroutine mod_atmos_phy_tb::atmosphytb_setup	(	class(atmosphytb), intent(inout)	this,
		class(modelmeshbase), intent(in), target	model_mesh,
		class(time_manager_component), intent(inout)	tm_parent_comp )

Setup a component of SGS turbulence process.

Parameters

model_mesh	Object to manage computational mesh of atmospheric model
tm_parent_comp	Object to mange a temporal scheme in a parent component

Definition at line 104 of file mod_atmos_phy_tb.F90.

    use mod_atmos_mesh, only: atmosmesh
    use scale_time_manager, only: time_manager_component
    use scale_element_hexahedral, only: hexahedralelement
 
    implicit none
    class(AtmosPhyTb), intent(inout) :: this
    class(ModelMeshBase), target, intent(in) :: model_mesh
    class(TIME_manager_component), intent(inout) :: tm_parent_comp
 
    real(DP) :: TIME_DT                             = undef8  !< Timestep for sub-grid scale turbulent process
    character(len=H_SHORT) :: TIME_DT_UNIT          = 'SEC'   !< Unit of timestep
 
    character(len=H_MID) :: TB_TYPE = 'SMAGORINSKY'           !< Type of sub-grid scale turbulent scheme
 
    character(len=H_SHORT) :: PostGLFilter_Type = 'None'
    character(len=H_SHORT) :: PostGLFilter_ConvFilterShape = 'GAUSSIAN'
    integer :: PostGLFilter_Nnodeh1D_reconst = -1
    real(RP) :: PostGLFilter_GaussinWidthFac = 1.5_rp
    real(RP) :: PostModalFilter_ALPHA_h = 0.0_rp
    real(RP) :: PostModalFilter_EtaC_h  = 0.0_rp
    integer :: PostModalFilter_ORDER_h  = 16    
    real(RP) :: PostModalFilter_ALPHA_v = 0.0_rp
    integer :: PostModalFilter_ORDER_v = 16
 
    namelist /param_atmos_phy_tb/ &
      time_dt,          & 
      time_dt_unit,     &
      tb_type,          &
      postglfilter_type, &
      postglfilter_nnodeh1d_reconst, &
      postglfilter_gaussinwidthfac, &
      postmodalfilter_alpha_h, & 
      postmodalfilter_etac_h, &
      postmodalfilter_order_h, &
      postmodalfilter_alpha_v, &
      postmodalfilter_order_v      
 
    
    class(AtmosMesh), pointer     :: atm_mesh
    class(MeshBase), pointer      :: ptr_mesh
    real(DP) :: dtsec
 
    integer :: ierr
    class(ElementBase3D), pointer :: elem3D
    type(HexahedralElement) :: elem3D_hex
    !--------------------------------------------------
 
    if (.not. this%IsActivated()) return
 
    log_newline
    log_info("ATMOS_PHY_TB_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_tb,iostat=ierr)
    if( ierr < 0 ) then !--- missing
      log_info("ATMOS_PHY_TB_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
      log_error("ATMOS_PHY_TB_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_TB. Check!'
      call prc_abort
    endif
    log_nml(param_atmos_phy_tb)
 
    !- get mesh --------------------------------------------------
 
    call model_mesh%GetModelMesh( ptr_mesh )
    select type(model_mesh)
    class is (atmosmesh)
      atm_mesh => model_mesh
    end select
 
    !--- Register this component in the time manager
    
    call tm_parent_comp%Regist_process( 'ATMOS_PHY_TB', time_dt, time_dt_unit, & ! (in)
      this%tm_process_id )                                                       ! (out)
 
    dtsec = tm_parent_comp%process_list(this%tm_process_id)%dtsec
 
    !- initialize the variables 
    call this%vars%Init( model_mesh )      
 
    !- Initialize a module for turbulence model
    call this%tb_driver%Init( tb_type, dtsec, atm_mesh )
 
    !--
    this%dyn_bnd => null()
 
 
    !--
    select case( trim(postglfilter_type) )
    case ('ConvolFilter', 'Reconstruction', 'Reconstruction2')
      this%gFilter_flag = .true.
    case ('None', 'ModalFilter')
      this%gFilter_flag = .false.
    case default
      log_info("ATMOS_PHY_TB_setup",*) 'Not appropriate names of PostGLFilter_Type in namelist PARAM_ATMOS_PHY_TB. Check!'
      call prc_abort
    end select
    
    !-
    if ( this%gFilter_flag ) then
      if ( postglfilter_nnodeh1d_reconst < 0 ) then
        elem3d=> atm_mesh%ptr_mesh%lcmesh_list(1)%refElem3D
        postglfilter_nnodeh1d_reconst = elem3d%Nnode_h1D
      end if
      call this%gFilter_phy_tb%Init( postglfilter_type, &
        postglfilter_convfiltershape, postglfilter_gaussinwidthfac, &
        postglfilter_nnodeh1d_reconst,                              &
        2, 0, 0, atm_mesh%ptr_mesh            )
    end if
    if ( postmodalfilter_alpha_h > 0.0_rp .or. postmodalfilter_alpha_v > 0.0_rp ) then
      this%modalFilter_flag = .true.
      elem3d=> atm_mesh%ptr_mesh%lcmesh_list(1)%refElem3D
      call elem3d_hex%Init( elem3d%PolyOrder_h, elem3d%PolyOrder_v, .false. )
      call this%modalfilter_3D%Init( elem3d_hex, postmodalfilter_etac_h, postmodalfilter_alpha_h, postmodalfilter_order_h, 0.0_rp, postmodalfilter_alpha_v, postmodalfilter_order_v )
      call elem3d_hex%Final()
    else
      this%modalFilter_flag = .false.
    end if
 
    return

References scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_kh_id, scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_nu_id, mod_atmos_phy_tb_vars::atmosphytbvars_getlocalmeshfields_tend(), and scale_atm_dyn_dgm_nonhydro3d_common::prgvar_ddens_id.