module Atmosphere Grid CartesianC metirc More...

Functions/Subroutines
subroutine, public	atmos_grid_cartesc_metric_setup
	Setup. More...

subroutine	atmos_grid_cartesc_metric_rotcoef
	Calculate rotation coeffient. More...

Variables
real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_mapf
	map factor More...

real(rp), dimension(:,:,:), allocatable, public	atmos_grid_cartesc_metric_rotc
	rotation coefficient More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_gsqrt
	transformation metrics from Z to Xi, {G}^1/2 More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_j13g
	(1,3) element of Jacobian matrix * {G}^1/2 More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_j23g
	(2,3) element of Jacobian matrix * {G}^1/2 More...

real(rp), public	atmos_grid_cartesc_metric_j33g
	(3,3) element of Jacobian matrix * {G}^1/2 More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_limyz
	flux limiter y-z face More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_limxz
	flux limiter x-z face More...

real(rp), dimension(:,:,:,:), allocatable, public	atmos_grid_cartesc_metric_limxy
	flux limiter x-y face More...

Detailed Description

module Atmosphere Grid CartesianC metirc

Description: Map projection and Terrain-following metrics for the CaresianC grid

Author: Team SCALE

NAMELIST

PARAM_ATMOS_GRID_CARTESC_METRIC

name	type	default value	comment
ATMOS_GRID_CARTESC_METRIC_OUT_BASENAME	character(len=H_LONG)	''	basename of the output file
ATMOS_GRID_CARTESC_METRIC_OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8
ATMOS_GRID_CARTESC_METRIC_TOPO_TYPE	character(len=H_SHORT)	'TERRAINFOLLOWING'	topographical shceme
ATMOS_GRID_CARTESC_METRIC_THINWALL_XDIV	integer	50	number dividing quarter-cell (x)
ATMOS_GRID_CARTESC_METRIC_THINWALL_YDIV	integer	50	number dividing quarter-cell (y)

History Output: No history output

Function/Subroutine Documentation

◆ atmos_grid_cartesc_metric_setup()

subroutine, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_setup

Setup.

Definition at line 72 of file scale_atmos_grid_cartesC_metric.F90.

     use scale_prc, only: &
        prc_abort
     use scale_prc_cartesc, only: &
        prc_twod
     implicit none
  
     namelist / param_atmos_grid_cartesc_metric / &
        atmos_grid_cartesc_metric_out_basename,  &
        atmos_grid_cartesc_metric_out_dtype,     &
        atmos_grid_cartesc_metric_topo_type,     &
        atmos_grid_cartesc_metric_thinwall_xdiv, &
        atmos_grid_cartesc_metric_thinwall_ydiv
  
     integer :: ierr
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("ATMOS_GRID_CARTESC_METRIC_setup",*) 'Setup'
  
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_atmos_grid_cartesc_metric,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("ATMOS_GRID_CARTESC_METRIC_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("ATMOS_GRID_CARTESC_METRIC_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_GRID_CARTESC_METRIC. Check!'
        call prc_abort
     endif
     log_nml(param_atmos_grid_cartesc_metric)
  
     allocate( atmos_grid_cartesc_metric_mapf(ia,ja,2,4) )
  
     allocate( atmos_grid_cartesc_metric_rotc(ia,ja,2) )
  
     if ( prc_twod ) then
        allocate( atmos_grid_cartesc_metric_gsqrt(ka,ia,ja,4) )
        allocate( atmos_grid_cartesc_metric_j13g(ka,ia,ja,4) )
        allocate( atmos_grid_cartesc_metric_j23g(ka,ia,ja,4) )
     else
        allocate( atmos_grid_cartesc_metric_gsqrt(ka,ia,ja,7) )
        allocate( atmos_grid_cartesc_metric_j13g(ka,ia,ja,7) )
        allocate( atmos_grid_cartesc_metric_j23g(ka,ia,ja,7) )
     end if
  
     atmos_grid_cartesc_metric_gsqrt(:,:,:,:) = 1.0_rp
     atmos_grid_cartesc_metric_j13g(:,:,:,:) = 0.0_rp
     atmos_grid_cartesc_metric_j23g(:,:,:,:) = 0.0_rp
     atmos_grid_cartesc_metric_j33g = 1.0_rp
  
     allocate( atmos_grid_cartesc_metric_limyz(ka,ia,ja,7) )
     allocate( atmos_grid_cartesc_metric_limxz(ka,ia,ja,7) )
     allocate( atmos_grid_cartesc_metric_limxy(ka,ia,ja,7) )
     atmos_grid_cartesc_metric_limyz(:,:,:,:) = 1.0_rp
     atmos_grid_cartesc_metric_limxz(:,:,:,:) = 1.0_rp
     atmos_grid_cartesc_metric_limxy(:,:,:,:) = 1.0_rp
  
     ! calc metrics for orthogonal curvelinear coordinate
     call atmos_grid_cartesc_metric_mapfactor
  
     ! calc coeficient for rotaion of velocity vector
     call atmos_grid_cartesc_metric_rotcoef
  
     ! calc metrics for terrain-following,step-mountain,thin-wall coordinate
     log_newline
     log_info("ATMOS_GRID_CARTESC_METRIC_setup",*) 'Terrain coordinate type : ', trim(atmos_grid_cartesc_metric_topo_type)
     select case(atmos_grid_cartesc_metric_topo_type)
     case('TERRAINFOLLOWING')
       log_info_cont(*) '=> Terrain-following method'
       call atmos_grid_cartesc_metric_terrainfollowing
     case('STEPMOUNTAIN')
       log_info_cont(*) '=> Step-mountain method'
       call atmos_grid_cartesc_metric_thin_wall
       call atmos_grid_cartesc_metric_step_mountain
     case('THINWALL')
       log_info_cont(*) '=> Thin-wall approximation method'
       call atmos_grid_cartesc_metric_thin_wall
     case default
        log_error("ATMOS_GRID_CARTESC_METRIC_setup",*) 'Unsupported ATMOS_GRID_CARTESC_METRIC_TOPO_type. STOP'
        call prc_abort
     end select
  
     ! output metrics (for debug)
     call atmos_grid_cartesc_metric_write
  
     return

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_grid_cartesc_metric_rotcoef()

subroutine scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotcoef

Calculate rotation coeffient.

Definition at line 200 of file scale_atmos_grid_cartesC_metric.F90.

     use scale_mapprojection, only: &
        mapprojection_rotcoef
     use scale_atmos_grid_cartesc_real, only: &
        atmos_grid_cartesc_real_lon,  &
        atmos_grid_cartesc_real_lat
     implicit none
     !---------------------------------------------------------------------------
  
     call mapprojection_rotcoef( ia, 1, ia, ja, 1, ja, &
          atmos_grid_cartesc_real_lon(:,:),   & ! [IN]
          atmos_grid_cartesc_real_lat(:,:),   & ! [IN]
          atmos_grid_cartesc_metric_rotc(:,:,1), & ! [OUT]
          atmos_grid_cartesc_metric_rotc(:,:,2)  ) ! [OUT]
  
     return