scale_atmos_grid_cartesc_metric Module Reference

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.

References atmos_grid_cartesc_metric_gsqrt, atmos_grid_cartesc_metric_j13g, atmos_grid_cartesc_metric_j23g, atmos_grid_cartesc_metric_j33g, atmos_grid_cartesc_metric_limxy, atmos_grid_cartesc_metric_limxz, atmos_grid_cartesc_metric_limyz, atmos_grid_cartesc_metric_mapf, atmos_grid_cartesc_metric_rotc, atmos_grid_cartesc_metric_rotcoef(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_calc_areavol(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_latuv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_latuy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_latxv, scale_atmos_grid_cartesc_index::i_uv, scale_atmos_grid_cartesc_index::i_uy, scale_atmos_grid_cartesc_index::i_xv, scale_atmos_grid_cartesc_index::i_xy, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_mapprojection::mapprojection_mapfactor(), and scale_prc::prc_abort().

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

    use scale_prc, only: &
       prc_abort
    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) )

    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) )

    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

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 182 of file scale_atmos_grid_cartesC_metric.F90.

References scale_atmos_grid_cartesc::atmos_grid_cartesc_cx, scale_atmos_grid_cartesc::atmos_grid_cartesc_cy, scale_atmos_grid_cartesc::atmos_grid_cartesc_cz, scale_atmos_grid_cartesc::atmos_grid_cartesc_fx, scale_atmos_grid_cartesc::atmos_grid_cartesc_fy, scale_atmos_grid_cartesc::atmos_grid_cartesc_fz, atmos_grid_cartesc_metric_gsqrt, atmos_grid_cartesc_metric_j13g, atmos_grid_cartesc_metric_j23g, atmos_grid_cartesc_metric_j33g, atmos_grid_cartesc_metric_limxy, atmos_grid_cartesc_metric_limxz, atmos_grid_cartesc_metric_limyz, atmos_grid_cartesc_metric_mapf, atmos_grid_cartesc_metric_rotc, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdy, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_basepoint_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_basepoint_lon, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_czuv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_czuy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_czxv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fzuv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fzuy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fzxv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lon, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdy, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdz, scale_const::const_radius, scale_atmos_grid_cartesc_index::i_fxy, scale_atmos_grid_cartesc_index::i_fxz, scale_atmos_grid_cartesc_index::i_fyz, scale_atmos_grid_cartesc_index::i_uv, scale_atmos_grid_cartesc_index::i_uvz, scale_atmos_grid_cartesc_index::i_uy, scale_atmos_grid_cartesc_index::i_uyw, scale_atmos_grid_cartesc_index::i_uyz, scale_atmos_grid_cartesc_index::i_xv, scale_atmos_grid_cartesc_index::i_xvw, scale_atmos_grid_cartesc_index::i_xvz, scale_atmos_grid_cartesc_index::i_xy, scale_atmos_grid_cartesc_index::i_xyw, scale_atmos_grid_cartesc_index::i_xyz, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_mapprojection::mapprojection_rotcoef(), scale_prc::prc_abort(), scale_topography::topo_zsfc, and scale_vector::vectr_distance().

Referenced by atmos_grid_cartesc_metric_setup().

    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(:,:,:) ) ! [OUT]

    return

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

atmos_grid_cartesc_metric_mapf

real(rp), dimension (:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_mapf

map factor

Definition at line 34 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_setup(), and mod_atmos_vars::atmos_vars_monitor().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_mapf (:,:,:,:)

atmos_grid_cartesc_metric_rotc

real(rp), dimension (:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc

rotation coefficient

Definition at line 35 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), mod_atmos_vars::atmos_vars_get_diagnostic_3d(), scale_comm_cartesc_nest::comm_cartesc_nest_nestdown(), mod_atmos_phy_sf_driver::history_output(), mod_realinput_scale::parentatmosinputscale(), and mod_realinput::realinput_surface().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_rotc (:,:,:)

atmos_grid_cartesc_metric_gsqrt

real(rp), dimension(:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_gsqrt

transformation metrics from Z to Xi, {G}^1/2

Definition at line 37 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by scale_atmos_diagnostic_cartesc::atmos_diagnostic_cartesc_get_vel(), mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), and mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_gsqrt(:,:,:,:)

atmos_grid_cartesc_metric_j13g

real(rp), dimension (:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_j13g

(1,3) element of Jacobian matrix * {G}^1/2

Definition at line 38 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by scale_atmos_diagnostic_cartesc::atmos_diagnostic_cartesc_get_vel(), mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), and mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_j13g (:,:,:,:)

atmos_grid_cartesc_metric_j23g

real(rp), dimension (:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_j23g

(2,3) element of Jacobian matrix * {G}^1/2

Definition at line 39 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by scale_atmos_diagnostic_cartesc::atmos_diagnostic_cartesc_get_vel(), mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), and mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_j23g (:,:,:,:)

atmos_grid_cartesc_metric_j33g

real(rp), public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_j33g

(3,3) element of Jacobian matrix * {G}^1/2

Definition at line 40 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_grid_cartesc_metric_rotcoef(), atmos_grid_cartesc_metric_setup(), and mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency().

  real(RP), public              :: atmos_grid_cartesc_metric_j33g

atmos_grid_cartesc_metric_limyz

real(rp), dimension(:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_limyz

flux limiter y-z face

Definition at line 42 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by atmos_grid_cartesc_metric_rotcoef(), and atmos_grid_cartesc_metric_setup().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_limyz(:,:,:,:)

atmos_grid_cartesc_metric_limxz

real(rp), dimension(:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_limxz

flux limiter x-z face

Definition at line 43 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by atmos_grid_cartesc_metric_rotcoef(), and atmos_grid_cartesc_metric_setup().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_limxz(:,:,:,:)

atmos_grid_cartesc_metric_limxy

real(rp), dimension(:,:,:,:), allocatable, public scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_limxy

flux limiter x-y face

Definition at line 44 of file scale_atmos_grid_cartesC_metric.F90.

Referenced by atmos_grid_cartesc_metric_rotcoef(), and atmos_grid_cartesc_metric_setup().

  real(RP), public, allocatable :: atmos_grid_cartesc_metric_limxy(:,:,:,:)