mod_atmos_mesh_rm Module Reference

module Atmosphere / Mesh More...

Data Types
type	atmosmeshrm
	Derived type to manage a computational mesh of regional atmospheric model. More...

Functions/Subroutines
subroutine	atmosmeshrm_init (this)
	Initialize a object to manage computational mesh.

Detailed Description

module Atmosphere / Mesh

Description

Module for mesh with atmospheric regional model

Author

Yuta kawai, Team SCALE

NAMELIST

• PARAM_ATMOS_MESH

  name	type	default value	comment
  DOM_XMIN	real(RP)	0.0_RP	Minimum x-coordinate value of the computational domain
  DOM_XMAX	real(RP)	100.0E3_RP	Maximum x-coordinate value of the computational domain
  DOM_YMIN	real(RP)	0.0_RP	Minimum y-coordinate value of the computational domain
  DOM_YMAX	real(RP)	100.0E3_RP	Maximum y-coordinate value of the computational domain
  DOM_ZMIN	real(RP)	0.0_RP	Minimum vertical coordinate value of the computational domain
  DOM_ZMAX	real(RP)	10.0E3_RP	Maximum vertical coordinate value of the computational domain
  FZ	real(RP), dimension(FZ_NMAX)		Values of the vertically computational coordinate at the element boundaries
  ISPERIODICX	logical	.true.	Flag whether a periodic boundary condition is applied in the x-direction
  ISPERIODICY	logical	.true.	Flag whether a periodic boundary condition is applied in the y-direction
  ISPERIODICZ	logical	.false.	Flag whether a periodic boundary condition is applied in the vertical direction
  NEX	integer	2	Number of finite element in the x-direction in each MPI process
  NEY	integer	2	Number of finite element in the y-direction in each MPI process
  NEZ	integer	2	Number of finite element in the vertical direction in each MPI process
  POLYORDER_H	integer	2	Polynomial order for the horizontal direction
  POLYORDER_V	integer	2	Polynomial order for the z-direction
  LUMPEDMASSMATFLAG	logical	.false.	Flag whether a mass lumping is applied
  ELEMENT_OPERATION_TYPE	character(len=H_SHORT)	'General'	General or TensorProd3D
  SPMV_STORAGE_FORMAT	character(len=H_SHORT)	'ELL'	CSR or ELL
  NPRCX	integer	1	Number of MPI process in the x-direction
  NPRCY	integer	1	Number of MPI process in the y-direction
  TOPO_IN_BASENAME	character(len=H_LONG)	''	Basename of the input file
  TOPO_IN_VARNAME	character(len=H_MID)	'topo'	Variable name of topography in the input file
  VERTICAL_COORD_NAME	character(len=H_MID)	"TERRAIN_FOLLOWING"	Type of the vertical coordinate
  COMM_USE_MPI_PC	logical	.false.	Flag whether persistent communication is used in MPI

  
  

History Output

No history output

Function/Subroutine Documentation

atmosmeshrm_init()

subroutine mod_atmos_mesh_rm::atmosmeshrm_init

(

class(atmosmeshrm), intent(inout), target

this

)

Initialize a object to manage computational mesh.

Definition at line 82 of file mod_atmos_mesh_rm.F90.

    use scale_file_base_meshfield, only: file_base_meshfield
    use scale_mesh_base2d, only: &
      mftype2d_xy => meshbase2d_dimtypeid_xy
    use scale_meshutil_vcoord, only: &
      meshutil_get_vcoord_typeid
 
    implicit none
    class(AtmosMeshRM), target, intent(inout) :: this
 
    real(RP) :: dom_xmin         = 0.0_rp      
    real(RP) :: dom_xmax         = 100.0e3_rp  
    real(RP) :: dom_ymin         = 0.0_rp      
    real(RP) :: dom_ymax         = 100.0e3_rp  
    real(RP) :: dom_zmin         = 0.0_rp      
    real(RP) :: dom_zmax         = 10.0e3_rp   
    logical  :: isPeriodicX       = .true.     
    logical  :: isPeriodicY       = .true.     
    logical  :: isPeriodicZ       = .false.    
 
    integer  :: NeX               = 2
    integer  :: NeY               = 2
    integer  :: NeZ               = 2
    integer  :: NprcX             = 1
    integer  :: NprcY             = 1
    integer  :: PolyOrder_h       = 2
    integer  :: PolyOrder_v       = 2
    logical  :: LumpedMassMatFlag = .false.  
 
    character(len=H_LONG) :: TOPO_IN_BASENAME    = ''
    character(len=H_MID)  :: TOPO_IN_VARNAME     = 'topo'
    character(len=H_MID)  :: VERTICAL_COORD_NAME = "TERRAIN_FOLLOWING"
 
    logical :: COMM_USE_MPI_PC    = .false.  
 
    character(len=H_SHORT) :: Element_operation_type = 'General'
    character(len=H_SHORT) :: SpMV_storage_format    = 'ELL'
 
    integer, parameter :: ATMOS_MESH_NLocalMeshPerPrc = 1
 
    integer, parameter :: FZ_nmax = 1000
    real(RP) :: FZ(FZ_nmax)
 
    namelist / param_atmos_mesh / &
      dom_xmin, dom_xmax,                          &
      dom_ymin, dom_ymax,                          &
      dom_zmin, dom_zmax,                          &
      fz,                                          &
      isperiodicx, isperiodicy, isperiodicz,       &
      nex, ney, nez,                               &
      polyorder_h, polyorder_v, lumpedmassmatflag, &
      element_operation_type,                      &
      spmv_storage_format,                         &
      nprcx, nprcy,                                &
      topo_in_basename, topo_in_varname,           &
      vertical_coord_name,                         &
      comm_use_mpi_pc
    
    integer :: k
    logical :: is_spec_FZ
 
    integer :: ierr
 
    type(FILE_base_meshfield) :: file_topo
    !-------------------------------------------
 
    log_newline
    log_info("ATMOS_MESH_setup",*) 'Setup'
 
    fz(:) = -1.0_rp
 
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_mesh,iostat=ierr)
    if( ierr < 0 ) then !--- missing
        log_info("ATMOS_MESH_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
        log_error("ATMOS_MESH_setup",*) 'Not appropriate names in namelist PARAM_ATM_MESH. Check!'
        call prc_abort
    endif
    log_nml(param_atmos_mesh)
 
    !----
 
    ! Setup the element
 
    call this%element%Init( polyorder_h, polyorder_v, lumpedmassmatflag )
    call this%element_v1D%Init( polyorder_v, lumpedmassmatflag )
 
    ! Setup the mesh
    
    is_spec_fz = .true.
    do k=1, nez+1
      if (fz(k) < 0.0_rp) then
        is_spec_fz = .false.
      end if
    end do
    if (is_spec_fz) then
      call this%mesh%Init( &
        nprcx*nex, nprcy*ney, nez,                                 &
        dom_xmin, dom_xmax,dom_ymin, dom_ymax, dom_zmin, dom_zmax, &
        isperiodicx, isperiodicy, isperiodicz,                     &
        this%element, atmos_mesh_nlocalmeshperprc, nprcx, nprcy,   &
        fz=fz(1:nez+1)    )
    else
      call this%mesh%Init( &
        nprcx*nex, nprcy*ney, nez,                                 &
        dom_xmin, dom_xmax,dom_ymin, dom_ymax, dom_zmin, dom_zmax, &
        isperiodicx, isperiodicy, isperiodicz,                     &
        this%element, atmos_mesh_nlocalmeshperprc, nprcx, nprcy    )
    end if
    
    call this%mesh%Generate()
    
    !-
    call this%AtmosMesh_Init( this%mesh )
    call this%PrepairElementOperation(  element_operation_type, spmv_storage_format )
 
    !- Set topography & vertical coordinate
 
    if ( topo_in_basename /= '' ) then
      log_info("ATMOS_MESH_setup",*) 'Read topography data'
 
      call file_topo%Init(1, mesh2d=this%mesh%mesh2D )
      call file_topo%Open( topo_in_basename, myrank=prc_myrank )
      call file_topo%Read_Var( mftype2d_xy, topo_in_varname, this%topography%topo )
      call file_topo%Close()
      call file_topo%Final()
    end if
 
    this%vcoord_type_id = meshutil_get_vcoord_typeid( vertical_coord_name )
    call this%Setup_vcoordinate()
 
    !-
    this%comm_use_mpi_pc = comm_use_mpi_pc
 
    return

References mod_atmos_mesh::atm_mesh_max_commnuicator_num, scale_mesh_base2d::meshbase2d_dimtypeid_xy, and scale_meshutil_vcoord::meshutil_get_vcoord_typeid().