scale_atm_dyn_dgm_spongelayer Module Reference

module FElib / Fluid dyn solver / Atmosphere / Common More...

Data Types
type	atmdynspongelayer

Functions/Subroutines
subroutine	atm_dyn_dgm_spongelayer_init (this, mesh3d, dtsec)

Detailed Description

module FElib / Fluid dyn solver / Atmosphere / Common

Description

Sponge layer for Atmospheric dynamical process.

Author

Yuta Kawai, Team SCALE

NAMELIST

• PARAM_ATMOS_DYN_SPONGELAYER

  name	type	default value	comment
  SL_WDAMP_TAU	real(RP)	-1.0_RP	the maximum tau for Rayleigh damping of w [s]
  SL_WDAMP_HEIGHT	real(RP)	-1.0_RP	the height to start apply Rayleigh damping [m]
  SL_WDAMP_LAYER	integer	-1	the vertical number of finite element to start apply Rayleigh damping [num]
  SL_HORIVELDAMP_FLAG	logical	.false.	Is the horizontal velocity damped?

  
  

History Output

No history output

Function/Subroutine Documentation

atm_dyn_dgm_spongelayer_init()

subroutine scale_atm_dyn_dgm_spongelayer::atm_dyn_dgm_spongelayer_init	(	class(atmdynspongelayer), intent(inout), target	this,
		class(meshbase3d), intent(in), target	mesh3d,
		real(rp), intent(in)	dtsec )

Definition at line 57 of file scale_atm_dyn_dgm_spongelayer.F90.

    use scale_prc, only: prc_abort
    use scale_mesh_cubedom3d, only: meshcubedom3d
    use scale_mesh_cubedspheredom3d, only: meshcubedspheredom3d
    implicit none
 
    class(AtmDynSpongeLayer), target, intent(inout) :: this
    class(MeshBase3D), target, intent(in) :: mesh3D
    real(RP), intent(in) :: dtsec
 
    real(RP) :: SL_WDAMP_TAU        = -1.0_rp ! the maximum tau for Rayleigh damping of w [s]
    real(RP) :: SL_WDAMP_HEIGHT     = -1.0_rp ! the height to start apply Rayleigh damping [m]
    integer  :: SL_WDAMP_LAYER      = -1      ! the vertical number of finite element to start apply Rayleigh damping [num]
    logical  :: SL_HORIVELDAMP_FLAG = .false. ! Is the horizontal velocity damped? 
    
    namelist /param_atmos_dyn_spongelayer/ &
      sl_wdamp_tau,                        &                
      sl_wdamp_height,                     &
      sl_wdamp_layer,                      &
      sl_horiveldamp_flag
    
    class(LocalMesh3D), pointer :: lcmesh3D
    class(ElementBase3D), pointer :: elem3D
  
    integer :: NeGZ
    integer :: ierr
    !---------------------------------------------------------------
 
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_dyn_spongelayer,iostat=ierr)
    if( ierr < 0 ) then !--- missing
      log_info("ATMOS_DYN_setup_spongelayer",*) 'Not found namelist. Default used.'
    else if( ierr > 0 ) then !--- fatal error
      log_error("ATMOS_DYN_setup_spongelayer",*) 'Not appropriate names in namelist PARAM_ATMOS_DYN_SPONGELAYER. Check!'
      call prc_abort
    end if
    log_nml(param_atmos_dyn_spongelayer)
 
    this%wdamp_tau    = sl_wdamp_tau 
    this%wdamp_height = sl_wdamp_height
 
    lcmesh3d => mesh3d%lcmesh_list(1)
    elem3d => lcmesh3d%refElem3D
    select type(mesh3d)
    type is (meshcubedom3d)
      negz = mesh3d%NeGZ
    type is (meshcubedspheredom3d)
      negz = mesh3d%NeGZ
    end select
 
    if ( sl_wdamp_layer > negz ) then
      log_error("ATMOS_DYN_setup_spongelayer",*) 'SL_wdamp_layer should be less than total of vertical elements (NeGZ). Check!'
      call prc_abort
    else if( sl_wdamp_layer > 0 ) then
      this%wdamp_height = lcmesh3d%pos_en(1,1+(sl_wdamp_layer-1)*lcmesh3d%NeX*lcmesh3d%NeY,3)
    end if
    if ( this%wdamp_tau < 0.0_rp ) then
      this%wdamp_tau = dtsec * 10.0_rp
    else if ( this%wdamp_tau < dtsec ) then
      log_error("ATMOS_DYN_setup_spongelayer",*) 'SL_wdamp_tau should be larger than TIME_DT (ATMOS_DYN). Check!'
      call prc_abort
    end if
    
    this%hveldamp_flag = sl_horiveldamp_flag
 
    return