scale_atm_dyn_dgm_nonhydro3d_common Module Reference

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

Functions/Subroutines
subroutine, public	atm_dyn_dgm_nonhydro3d_common_init (mesh)
	Initialize a common module for atmospheric nonhydrostatic dynamical core.
subroutine, public	atm_dyn_dgm_nonhydro3d_common_final ()
	Finalize a common module for atmospheric nonhydrostatic dynamical core.
subroutine, public	atm_dyn_dgm_nonhydro3d_common_get_varinfo (prgvar_info, auxvar_info, phytend_info)
	Get variable information for atmospheric nonhydrostatic dynamical core.
subroutine, public	atm_dyn_dgm_nonhydro3d_common_setup_variables (prgvars, qtrcvars, auxvars, phytends, prgvar_manager, qtrcvar_manager, auxvar_manager, phytend_manager, reg_file_hist, do_setup_phytend, phytend_num_tot, mesh3d, prgvar_varinfo)
	Setup variable managers for atmospheric nonhydrostatic dynamical core.
subroutine, public	atm_dyn_dgm_nonhydro3d_common_calc_pressure (pres, dpres, ddens, momx, momy, momz, therm, pres_hyd, dens_hyd, therm_hyd, rtot, cvtot, cptot, mesh3d, entot_conserve_scheme_flag)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_calc_therm_phyd (therm_hyd, pres_hyd, dens_hyd, mesh3d, entot_conserve_scheme_flag)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_drhot2pres (pres, dpres, drhot, pres_hyd, therm_hyd, rtot, cvtot, cptot, lcmesh, elem3d)
	Calculate pressure from the deviation of density-weighted potential temperature (DRHOT)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_drhot2entot (entot, ddens, momx, momy, momz, drhot, dens_hyd, pres_hyd, therm_hyd, rtot, cvtot, cptot, lcmesh, elem3d)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_entot2pres (pres, dpres, ddens, momx, momy, momz, entot, pres_hyd, dens_hyd, rtot, cvtot, lcmesh, elem3d)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_calc_rhot_hyd (rhot_hyd, pres_hyd, lcmesh, elem3d)
subroutine, public	atm_dyn_dgm_nonhydro3d_common_calc_phyd_hgrad_lc (dphyddx, dphyddy, pres_hyd, pres_hyd_ref, element3d_operation, lmesh, elem)
	Calculate horizontal graidient of hydrostatic pressure In this calculation, we assume that PRES_hyd_ref is continuous at element boundaries.

Variables
integer, parameter, public	prgvar_ddens_id = 1
integer, parameter, public	prgvar_therm_id = 2
integer, parameter, public	prgvar_drhot_id = 2
integer, parameter, public	prgvar_etot_id = 2
integer, parameter, public	prgvar_momz_id = 3
integer, parameter, public	prgvar_momx_id = 4
integer, parameter, public	prgvar_momy_id = 5
integer, parameter, public	prgvar_scalar_num = 3
integer, parameter, public	prgvar_hvec_num = 1
integer, parameter, public	prgvar_num = 5
integer, parameter, public	phytend_dens_id = 1
integer, parameter, public	phytend_momx_id = 2
integer, parameter, public	phytend_momy_id = 3
integer, parameter, public	phytend_momz_id = 4
integer, parameter, public	phytend_rhot_id = 5
integer, parameter, public	phytend_rhoh_id = 6
integer, parameter, public	phytend_num = 6
integer, parameter, public	auxvar_preshydro_id = 1
integer, parameter, public	auxvar_denshydro_id = 2
integer, parameter, public	auxvar_thermhydro_id = 3
integer, parameter, public	auxvar_pres_id = 4
integer, parameter, public	auxvar_pt_id = 5
integer, parameter, public	auxvar_rtot_id = 6
integer, parameter, public	auxvar_cvtot_id = 7
integer, parameter, public	auxvar_cptot_id = 8
integer, parameter, public	auxvar_qdry_id = 9
integer, parameter, public	auxvar_preshydro_ref_id = 10
integer, parameter, public	auxvar_num = 10
real(rp), dimension(:,:), allocatable, public	intrpmat_vpordm1

Detailed Description

module FElib / Fluid dyn solver / Atmosphere / Nonhydrostatic model / Common

Description

A common model for atmospheric nonhydrostatic dynamical core

Author

Yuta Kawai, Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
PRES_hyd	hydrostatic part of pressure	Pa	PRES_hyd
DENS_hyd	hydrostatic part of density	kg/m3	DENS_hyd
THERM_hyd	hydrostatic part of THERM	kg/m3	THERM_hyd
PRES	pressure	Pa	PRES
PT	potential temperature	K	PT
RTOT	Total gas constant	J/kg/K	RTOT
CVTOT	Total heat capacity	J/kg/K	CVTOT
CPTOT	Total heat capacity	J/kg/K	CPTOT
QDRY	dry air	kg/kg	QDRY
PRES_hyd_REF	hydrostatic reference pressure	Pa	PRES_hyd_REF
DENS_tp	DENS_tp	kg/m3/s	DENS_tp
MOMX_tp	MOMX_tp	kg/m2/s	MOMX_tp
MOMY_tp	MOMY_tp	kg/m2/s	MOMY_tp
MOMZ_tp	MOMZ_tp	kg/m2/s	MOMZ_tp
RHOT_tp	RHOT_tp	kg/m3.K/s	RHOT_tp
RHOH_p	RHOH_p	kg/m3.J/s	RHOH_p
DDENS	deviation of density	kg/m3	DDENS
THERM	THERM	-	THERM
MOMZ	momentum z	kg/m2/s	MOMZ
MOMX	momentum x	kg/m2/s	MOMX
MOMY	momentum y	kg/m2/s	MOMY
{TRACER_NAME}	{TRACER_NAME}; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	TRACER_UNIT	TRACER_NAME
{TRACER_NAME}_tp	tendency of physical process for {TRACER_NAME}; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	TRACER_UNIT	TRACER_NAME

Function/Subroutine Documentation

atm_dyn_dgm_nonhydro3d_common_init()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_init

(

class(meshbase3d), intent(in)

mesh

)

Initialize a common module for atmospheric nonhydrostatic dynamical core.

Definition at line 112 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

    implicit none
    class(MeshBase3D), intent(in) :: mesh
 
    type(ElementBase3D), pointer :: elem
    !--------------------------------------------
 
    elem => mesh%refElem3D
    allocate( intrpmat_vpordm1(elem%Np,elem%Np) )
    call mesh%refElem3D%Generate_ModalTruncationMat( elem%PolyOrder_h, elem%PolyOrder_v-1, & ! (in)
      intrpmat_vpordm1 ) ! (out)
    return

References intrpmat_vpordm1.

Referenced by scale_atm_dyn_dgm_globalnonhydro3d_etot_heve::atm_dyn_dgm_globalnonhydro3d_etot_heve_init(), scale_atm_dyn_dgm_globalnonhydro3d_etot_hevi::atm_dyn_dgm_globalnonhydro3d_etot_hevi_init(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_init(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve::atm_dyn_dgm_globalnonhydro3d_rhot_heve_init(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_hevi::atm_dyn_dgm_globalnonhydro3d_rhot_hevi_init(), scale_atm_dyn_dgm_nonhydro3d_etot_heve::atm_dyn_dgm_nonhydro3d_etot_heve_init(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi::atm_dyn_dgm_nonhydro3d_etot_hevi_init(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve_gpu::atm_dyn_dgm_nonhydro3d_rhot_heve_gpu_init(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve::atm_dyn_dgm_nonhydro3d_rhot_heve_init(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve_splitform::atm_dyn_dgm_nonhydro3d_rhot_heve_splitform_init(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi::atm_dyn_dgm_nonhydro3d_rhot_hevi_init(), and scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform::atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform_init().

atm_dyn_dgm_nonhydro3d_common_final()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_final

Finalize a common module for atmospheric nonhydrostatic dynamical core.

Definition at line 128 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

    implicit none
    !--------------------------------------------
    
    deallocate( intrpmat_vpordm1 )
    return

References intrpmat_vpordm1.

Referenced by scale_atm_dyn_dgm_globalnonhydro3d_etot_heve::atm_dyn_dgm_globalnonhydro3d_etot_heve_final(), scale_atm_dyn_dgm_globalnonhydro3d_etot_hevi::atm_dyn_dgm_globalnonhydro3d_etot_hevi_final(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve::atm_dyn_dgm_globalnonhydro3d_rhot_heve_final(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_final(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_hevi::atm_dyn_dgm_globalnonhydro3d_rhot_hevi_final(), scale_atm_dyn_dgm_nonhydro3d_etot_heve::atm_dyn_dgm_nonhydro3d_etot_heve_final(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi::atm_dyn_dgm_nonhydro3d_etot_hevi_final(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve::atm_dyn_dgm_nonhydro3d_rhot_heve_final(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve_gpu::atm_dyn_dgm_nonhydro3d_rhot_heve_gpu_final(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve_splitform::atm_dyn_dgm_nonhydro3d_rhot_heve_splitform_final(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi::atm_dyn_dgm_nonhydro3d_rhot_hevi_final(), and scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform::atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform_final().

atm_dyn_dgm_nonhydro3d_common_get_varinfo()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_get_varinfo	(	type(variableinfo), dimension(prgvar_num), intent(out)	prgvar_info,
		type(variableinfo), dimension(auxvar_num), intent(out)	auxvar_info,
		type(variableinfo), dimension(phytend_num), intent(out), optional	phytend_info )

Get variable information for atmospheric nonhydrostatic dynamical core.

Definition at line 137 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    implicit none
    type(VariableInfo), intent(out) :: prgvar_info(PRGVAR_NUM)
    type(VariableInfo), intent(out) :: auxvar_info(AUXVAR_NUM)
    type(VariableInfo), intent(out), optional :: phytend_info(PHYTEND_NUM)
 
    type(VariableInfo) :: PRGVAR_VARINFO(PRGVAR_NUM)
    DATA prgvar_varinfo / &
      variableinfo( prgvar_ddens_id, 'DDENS', 'deviation of density',         &
                    'kg/m3',  3, 'XYZ',  'air_density'                   ),   &
      variableinfo( prgvar_therm_id, 'THERM', 'THERM',                        &
                    '-', 3, 'XYZ',  ''                                     ), &
      variableinfo( prgvar_momz_id , 'MOMZ', 'momentum z',                    &
                    'kg/m2/s', 3, 'XYZ', 'northward_mass_flux_of_air'      ), &
      variableinfo( prgvar_momx_id , 'MOMX', 'momentum x',                    &
                    'kg/m2/s', 3, 'XYZ', 'upward_mass_flux_of_air'         ), &
      variableinfo( prgvar_momy_id , 'MOMY', 'momentum y',                    &
                    'kg/m2/s', 3, 'XYZ', 'eastward_mass_flux_of_air'       )  /
 
    type(VariableInfo) :: AUXVAR_VARINFO(AUXVAR_NUM)
    DATA auxvar_varinfo / &
      variableinfo( auxvar_preshydro_id, 'PRES_hyd', 'hydrostatic part of pressure',             &
                      'Pa', 3, 'XYZ', ''                                                      ), &
      variableinfo( auxvar_denshydro_id, 'DENS_hyd', 'hydrostatic part of density',              &
                    'kg/m3', 3, 'XYZ', ''                                                     ), &
      variableinfo( auxvar_thermhydro_id, 'THERM_hyd', 'hydrostatic part of THERM',              &
                    'kg/m3', 3, 'XYZ', ''                                                     ), &
      variableinfo( auxvar_pres_id     ,     'PRES', 'pressure',                                 &
                      'Pa', 3, 'XYZ', 'air_pressure'                                          ), &
      variableinfo( auxvar_pt_id       ,       'PT', 'potential temperature',                    &
                        'K', 3, 'XYZ', 'potential_temperature'                                ), &
      variableinfo( auxvar_rtot_id     ,     'RTOT', 'Total gas constant',                       &
                        'J/kg/K', 3, 'XYZ', ''                                                ), &
      variableinfo( auxvar_cvtot_id    ,    'CVTOT', 'Total heat capacity',                      &
                        'J/kg/K', 3, 'XYZ', ''                                                ), &
      variableinfo( auxvar_cptot_id    ,    'CPTOT', 'Total heat capacity',                      &
                        'J/kg/K', 3, 'XYZ', ''                                                ), &
      variableinfo( auxvar_qdry_id     ,     'QDRY', 'dry air',                                  &
                        'kg/kg', 3, 'XYZ', ''                                                 ), &
      variableinfo( auxvar_preshydro_ref_id, 'PRES_hyd_REF', 'hydrostatic reference pressure',   &
                      'Pa', 3, 'XYZ', ''                                                        )/
 
    type(VariableInfo) :: PHYTEND_VARINFO(PHYTEND_NUM)
    DATA phytend_varinfo / &
      variableinfo( phytend_dens_id, 'DENS_tp', 'DENS_tp',                        &
                    'kg/m3/s',  3, 'XYZ',  'tendency of physical process for DENS' ),   &
      variableinfo( phytend_momx_id, 'MOMX_tp', 'MOMX_tp',                        &
                    'kg/m2/s',  3, 'XYZ',  'tendency of physical process for MOMX' ),   &
      variableinfo( phytend_momy_id, 'MOMY_tp', 'MOMY_tp',                        &
                    'kg/m2/s',  3, 'XYZ',  'tendency of physical process for MOMY' ),   &
      variableinfo( phytend_momz_id, 'MOMZ_tp', 'MOMZ_tp',                        &
                    'kg/m2/s',  3, 'XYZ',  'tendency of physical process for MOMZ' ),   &
      variableinfo( phytend_rhot_id, 'RHOT_tp', 'RHOT_tp',                        &
                    'kg/m3.K/s',  3, 'XYZ',  'tendency of physical process for RHOT' ), &
      variableinfo( phytend_rhoh_id,  'RHOH_p',  'RHOH_p',                        &
                    'kg/m3.J/s',  3, 'XYZ',  'heating of physical process for THERM' )   /
 
    !----------------------------------------------------------
 
    prgvar_info(:) = prgvar_varinfo
    auxvar_info(:) = auxvar_varinfo
    if ( present(phytend_info) ) phytend_info(:) = phytend_varinfo
 
    return

References auxvar_cptot_id, auxvar_cvtot_id, auxvar_denshydro_id, auxvar_pres_id, auxvar_preshydro_id, auxvar_preshydro_ref_id, auxvar_pt_id, auxvar_qdry_id, auxvar_rtot_id, auxvar_thermhydro_id, phytend_dens_id, phytend_momx_id, phytend_momy_id, phytend_momz_id, phytend_rhoh_id, phytend_rhot_id, prgvar_ddens_id, prgvar_momx_id, prgvar_momy_id, prgvar_momz_id, and prgvar_therm_id.

Referenced by atm_dyn_dgm_nonhydro3d_common_setup_variables().

atm_dyn_dgm_nonhydro3d_common_setup_variables()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_setup_variables	(	type(meshfield3d), dimension(prgvar_num), intent(inout)	prgvars,
		type(meshfield3d), dimension(0:qa), intent(inout)	qtrcvars,
		type(meshfield3d), dimension(auxvar_num), intent(inout)	auxvars,
		type(meshfield3d), dimension(phytend_num_tot), intent(inout)	phytends,
		type(modelvarmanager), intent(inout)	prgvar_manager,
		type(modelvarmanager), intent(inout)	qtrcvar_manager,
		type(modelvarmanager), intent(inout)	auxvar_manager,
		type(modelvarmanager), intent(inout)	phytend_manager,
		logical, intent(in)	reg_file_hist,
		logical, intent(in)	do_setup_phytend,
		integer, intent(in)	phytend_num_tot,
		class(meshbase3d), intent(in)	mesh3d,
		type(variableinfo), dimension(prgvar_num), intent(out)	prgvar_varinfo )

Setup variable managers for atmospheric nonhydrostatic dynamical core.

Parameters

[in,out]	prgvars	Array of objects to manage prognostic variables
[in,out]	qtrcvars	Array of objects to manage tracer variables
[in,out]	auxvars	Array of objects to manage auxiliary variables
[in,out]	phytends	Array of objects to manage physics tendency variables
[in,out]	prgvar_manager	Object to manage prognostic variables
[in,out]	qtrcvar_manager	Object to manage tracer variables
[in,out]	auxvar_manager	Object to manage auxiliary variables
[in,out]	phytend_manager	Object to manage physics tendency variables
[in]	reg_file_hist	Flag whether variables are registered for history file output or not
[in]	do_setup_phytend	Flag whether to setup physics tendency variables or not
[in]	mesh3d	Object to manage 3D mesh
[out]	prgvar_varinfo	Variable information with prognostic variables

Definition at line 206 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    use scale_atmos_hydrometeor, only: &
      atmos_hydrometeor_dry
    use scale_tracer, only: &
      qa, tracer_name, tracer_desc, tracer_unit    
    use scale_mesh_base3d, only: meshbase3d
    use scale_meshfield_base, only: meshfield3d
    implicit none
    integer, intent(in) :: PHYTEND_NUM_TOT
    type(MeshField3D), intent(inout) :: prgvars(PRGVAR_NUM)       !< Array of objects to manage prognostic variables
    type(MeshField3D), intent(inout) :: qtrcvars(0:QA)            !< Array of objects to manage tracer variables
    type(MeshField3D), intent(inout) :: auxvars(AUXVAR_NUM)       !< Array of objects to manage auxiliary variables
    type(MeshField3D), intent(inout) :: phytends(PHYTEND_NUM_TOT) !< Array of objects to manage physics tendency variables
    type(ModelVarManager), intent(inout) :: prgvar_manager        !< Object to manage prognostic variables
    type(ModelVarManager), intent(inout) :: qtrcvar_manager       !< Object to manage tracer variables
    type(ModelVarManager), intent(inout) :: auxvar_manager        !< Object to manage auxiliary variables
    type(ModelVarManager), intent(inout) :: phytend_manager       !< Object to manage physics tendency variables
    logical, intent(in) :: reg_file_hist                          !< Flag whether variables are registered for history file output or not
    logical, intent(in) :: do_setup_phytend                       !< Flag whether to setup physics tendency variables or not
    class(MeshBase3D), intent(in) :: mesh3D                       !< Object to manage 3D mesh
    type(VariableInfo), intent(out) :: PRGVAR_VARINFO(PRGVAR_NUM) !< Variable information with prognostic variables
 
    type(VariableInfo) :: AUXVAR_VARINFO(AUXVAR_NUM)
    type(VariableInfo) :: PHYTEND_VARINFO(PHYTEND_NUM)
 
    integer :: iv
    integer :: iq
 
    type(VariableInfo) :: qtrc_dry_vinfo_tmp
    type(VariableInfo) :: qtrc_dry_tp_vinfo_tmp
    type(VariableInfo) :: qtrc_vinfo_tmp
    type(VariableInfo) :: qtrc_tp_vinfo_tmp
    !----------------------------------------------------------
 
    call atm_dyn_dgm_nonhydro3d_common_get_varinfo( prgvar_varinfo, auxvar_varinfo, phytend_varinfo ) ! (out)
 
    !- Initialize prognostic variables
 
    do iv = 1, prgvar_num
      call prgvar_manager%Regist(  &
        prgvar_varinfo(iv), mesh3d,                           & ! (in) 
        prgvars(iv),                                          & ! (inout)
        reg_file_hist,  monitor_flag=.true., fill_zero=.true. ) ! (out)
      !$acc update device( prgvars(iv) )
    end do
 
    !- Initialize tracer variables
 
    if ( qa == 0 ) then
      ! Dummy
      qtrc_dry_vinfo_tmp%ndims    = 3
      qtrc_dry_vinfo_tmp%dim_type = 'XYZ'
      qtrc_dry_vinfo_tmp%STDNAME  = ''
 
      qtrc_dry_vinfo_tmp%keyID = 0
      qtrc_dry_vinfo_tmp%NAME  = "QV"
      qtrc_dry_vinfo_tmp%DESC  = "Ratio of Water Vapor mass to total mass (Specific humidity)"
      qtrc_dry_vinfo_tmp%UNIT  = "kg/kg"
      call qtrcvar_manager%Regist( &
        qtrc_dry_vinfo_tmp, mesh3d,                      & ! (in) 
        qtrcvars(0),                                     & ! (inout)
        .false., monitor_flag=.false., fill_zero=.true.  ) ! (in)
    else
      qtrc_vinfo_tmp%ndims    = 3
      qtrc_vinfo_tmp%dim_type = 'XYZ'
      qtrc_vinfo_tmp%STDNAME  = ''
  
      do iq = 1, qa
        qtrc_vinfo_tmp%keyID = iq
        qtrc_vinfo_tmp%NAME  = tracer_name(iq)
        qtrc_vinfo_tmp%DESC  = tracer_desc(iq)
        qtrc_vinfo_tmp%UNIT  = tracer_unit(iq)
        call qtrcvar_manager%Regist( &
          qtrc_vinfo_tmp, mesh3d,                               & ! (in) 
          qtrcvars(iq),                                         & ! (inout)
          reg_file_hist, monitor_flag=.true., fill_zero=.true.  ) ! (in)
      end do
    end if
 
    !- Initialize auxiliary variables
 
    do iv = 1, auxvar_num
      call auxvar_manager%Regist( &
        auxvar_varinfo(iv), mesh3d,       & ! (in) 
        auxvars(iv),                      & ! (inout)
        reg_file_hist, fill_zero=.true.   ) ! (in)
    end do
 
    !- Initialize the tendency of physical processes
 
    if ( do_setup_phytend ) then
 
      do iv = 1, phytend_num    
        call phytend_manager%Regist( &
          phytend_varinfo(iv), mesh3d,     & ! (in) 
          phytends(iv),                    & ! (inout)
          reg_file_hist, fill_zero=.true.  ) ! (in)
      end do
 
      if ( qa == 0 ) then
        ! Dummy
        qtrc_dry_tp_vinfo_tmp%ndims    = 3
        qtrc_dry_tp_vinfo_tmp%dim_type = 'XYZ'
        qtrc_dry_tp_vinfo_tmp%STDNAME  = ''
 
        iv = phytend_num + 1
        qtrc_dry_tp_vinfo_tmp%keyID = iv
        qtrc_dry_tp_vinfo_tmp%NAME  = "QV_tp"
        qtrc_dry_tp_vinfo_tmp%DESC  = "tendency of physical process for QV"
        qtrc_dry_tp_vinfo_tmp%UNIT  = "kg/m3/s"
        call phytend_manager%Regist( &
          qtrc_dry_tp_vinfo_tmp, mesh3d,   & ! (in) 
          phytends(iv),                    & ! (inout)
          .false., fill_zero=.true.        ) ! (in)
      else
        qtrc_tp_vinfo_tmp%ndims    = 3
        qtrc_tp_vinfo_tmp%dim_type = 'XYZ'
        qtrc_tp_vinfo_tmp%STDNAME  = ''
 
        do iq = 1, qa
          iv = phytend_num + iq 
          qtrc_tp_vinfo_tmp%keyID = iv
          qtrc_tp_vinfo_tmp%NAME  = trim(tracer_name(iq))//'_tp'
          qtrc_tp_vinfo_tmp%DESC  = 'tendency of physical process for '//trim(tracer_desc(iq))
          qtrc_tp_vinfo_tmp%UNIT  = trim(tracer_unit(iq))//'/s'
 
          call phytend_manager%Regist( &
            qtrc_tp_vinfo_tmp, mesh3d,       & ! (in) 
            phytends(iv),                    & ! (inout)
            reg_file_hist, fill_zero=.true.  ) ! (in)
        end do    
      end if
 
    end if
 
    return

References atm_dyn_dgm_nonhydro3d_common_get_varinfo(), auxvar_num, phytend_num, and prgvar_num.

atm_dyn_dgm_nonhydro3d_common_calc_pressure()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_calc_pressure	(	class(meshfield3d), intent(inout)	pres,
		class(meshfield3d), intent(inout)	dpres,
		class(meshfield3d), intent(in)	ddens,
		class(meshfield3d), intent(in)	momx,
		class(meshfield3d), intent(in)	momy,
		class(meshfield3d), intent(in)	momz,
		class(meshfield3d), intent(in)	therm,
		class(meshfield3d), intent(in)	pres_hyd,
		class(meshfield3d), intent(in)	dens_hyd,
		class(meshfield3d), intent(in)	therm_hyd,
		class(meshfield3d), intent(in)	rtot,
		class(meshfield3d), intent(in)	cvtot,
		class(meshfield3d), intent(in)	cptot,
		class(meshbase3d), intent(in), target	mesh3d,
		logical, intent(in)	entot_conserve_scheme_flag )

Definition at line 350 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    implicit none
    class(MeshField3D), intent(inout) :: PRES
    class(MeshField3D), intent(inout) :: DPRES
    class(MeshField3D), intent(in) :: DDENS
    class(MeshField3D), intent(in) :: MOMX
    class(MeshField3D), intent(in) :: MOMY
    class(MeshField3D), intent(in) :: MOMZ
    class(MeshField3D), intent(in) :: THERM
    class(MeshField3D), intent(in) :: PRES_hyd
    class(MeshField3D), intent(in) :: DENS_hyd
    class(MeshField3D), intent(in) :: THERM_hyd
    class(MeshField3D), intent(in) :: Rtot
    class(MeshField3D), intent(in) :: CVtot
    class(MeshField3D), intent(in) :: CPtot
    class(MeshBase3D), intent(in), target :: mesh3D
    logical, intent(in) :: ENTOT_CONSERVE_SCHEME_FLAG
    
    integer :: n
    class(LocalMesh3D), pointer :: lcmesh3D
    !---------------------------
 
    do n=1, mesh3d%LOCAL_MESH_NUM
      lcmesh3d => mesh3d%lcmesh_list(n)
 
      if ( entot_conserve_scheme_flag ) then
        call atm_dyn_dgm_nonhydro3d_common_entot2pres( pres%local(n)%val, dpres%local(n)%val,                  &
          ddens%local(n)%val, momx%local(n)%val, momy%local(n)%val, momz%local(n)%val, therm%local(n)%val,     &
          pres_hyd%local(n)%val, dens_hyd%local(n)%val, rtot%local(n)%val, cvtot%local(n)%val,                 &
          lcmesh3d, lcmesh3d%refElem3D )
      else
        call atm_dyn_dgm_nonhydro3d_common_drhot2pres( pres%local(n)%val, dpres%local(n)%val,                                           &
          therm%local(n)%val, pres_hyd%local(n)%val, therm_hyd%local(n)%val, rtot%local(n)%val, cvtot%local(n)%val, cptot%local(n)%val, &
          lcmesh3d, lcmesh3d%refElem3D )
      end if
    end do
 
    return

References atm_dyn_dgm_nonhydro3d_common_drhot2pres(), and atm_dyn_dgm_nonhydro3d_common_entot2pres().

atm_dyn_dgm_nonhydro3d_common_calc_therm_phyd()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_calc_therm_phyd	(	class(meshfield3d), intent(inout)	therm_hyd,
		class(meshfield3d), intent(in)	pres_hyd,
		class(meshfield3d), intent(in)	dens_hyd,
		class(meshbase3d), intent(in), target	mesh3d,
		logical, intent(in)	entot_conserve_scheme_flag )

Definition at line 396 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    implicit none
    class(MeshField3D), intent(inout) :: THERM_hyd
    class(MeshField3D), intent(in) :: PRES_hyd
    class(MeshField3D), intent(in) :: DENS_hyd
    class(MeshBase3D), intent(in), target :: mesh3D
    logical, intent(in) :: ENTOT_CONSERVE_SCHEME_FLAG
    
    integer :: n
    class(LocalMesh3D), pointer :: lcmesh3D
    !---------------------------
 
    do n=1, mesh3d%LOCAL_MESH_NUM
      lcmesh3d => mesh3d%lcmesh_list(n)
 
      if ( entot_conserve_scheme_flag ) then
        ! ToDO
      else
        call atm_dyn_dgm_nonhydro3d_common_calc_rhot_hyd( therm_hyd%local(n)%val, &
          pres_hyd%local(n)%val, lcmesh3d, lcmesh3d%refElem3D                     )
      end if
    end do
 
    return

References atm_dyn_dgm_nonhydro3d_common_calc_rhot_hyd().

atm_dyn_dgm_nonhydro3d_common_drhot2pres()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_drhot2pres	(	real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	pres,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	dpres,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	drhot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	pres_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	therm_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	rtot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	cvtot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	cptot,
		class(localmesh3d), intent(in)	lcmesh,
		class(elementbase3d), intent(in)	elem3d )

Calculate pressure from the deviation of density-weighted potential temperature (DRHOT)

Definition at line 428 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    use scale_const, only: &
      rdry => const_rdry,   &
      cpdry => const_cpdry, &
      cvdry => const_cvdry, &
      pres00 => const_pre00
        
    implicit none
    class(LocalMesh3D), intent(in) :: lcmesh
    class(ElementBase3D), intent(in) :: elem3D
    real(RP), intent(out) :: PRES(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(out) :: DPRES(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: DRHOT(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: PRES_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: THERM_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: Rtot(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: CVtot(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: CPtot(elem3D%Np,lcmesh%NeA) 
 
    integer :: ke, p
    real(RP) :: RHOT(elem3D%Np)
    real(RP) :: RHOT_
    real(RP) :: rP0
    !---------------------------------------------------------------
 
    rp0 = 1.0_rp / pres00
    !$omp parallel do private( RHOT )
    !$acc parallel loop gang present( PRES, DPRES, DRHOT, PRES_hyd, THERM_hyd, Rtot, CVtot, CPtot, lcmesh,elem3D )
    do ke=lcmesh%NeS, lcmesh%NeE
!      RHOT(:) = PRES00 / Rdry * ( PRES_hyd(:,ke) / PRES00 )**(CvDry/CpDry) + DRHOT(:,ke)
 
#ifdef _OPENACC
      !$acc loop vector
      do p=1, elem3d%Np
        rhot_ = pres00 / rdry * ( pres_hyd(p,ke) / pres00 )**(cvdry/cpdry) + drhot(p,ke)
 
        pres(p,ke) = pres00 * ( rtot(p,ke) * rp0 * rhot_ )**( cptot(p,ke) / cvtot(p,ke) )
        dpres(p,ke) = pres(p,ke) - pres_hyd(p,ke)
      end do
#else
      rhot(:) = therm_hyd(:,ke) + drhot(:,ke)
 
      pres(:,ke) = pres00 * ( rtot(:,ke) * rp0 * rhot(:) )**( cptot(:,ke) / cvtot(:,ke) )
      dpres(:,ke) = pres(:,ke) - pres_hyd(:,ke)
#endif
    end do
 
    return

Referenced by atm_dyn_dgm_nonhydro3d_common_calc_pressure(), and atm_dyn_dgm_nonhydro3d_common_drhot2entot().

atm_dyn_dgm_nonhydro3d_common_drhot2entot()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_drhot2entot	(	real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	entot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	ddens,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momx,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momy,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momz,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	drhot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	dens_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	pres_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	therm_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	rtot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	cvtot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	cptot,
		class(localmesh3d), intent(in)	lcmesh,
		class(elementbase3d), intent(in)	elem3d )

Definition at line 482 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    use scale_const, only: &
      grav => const_grav
        
    implicit none
    class(LocalMesh3D), intent(in) :: lcmesh
    class(ElementBase3D), intent(in) :: elem3D
    real(RP), intent(out) :: EnTot(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: DDENS(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMX(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMY(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMZ(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: DRHOT(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: DENS_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: PRES_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: THERM_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: Rtot(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: CVtot(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: CPtot(elem3D%Np,lcmesh%NeA) 
 
    integer :: ke, ke2D
 
    real(RP) :: DENS(elem3D%Np)
    real(RP) :: mom_u1(elem3D%Np), mom_u2(elem3D%Np)
 
    real(RP) :: PRES(elem3D%Np,lcmesh%NeA)
    real(RP) :: DPRES(elem3D%Np,lcmesh%NeA)
    !---------------------------------------------------------------
 
    call atm_dyn_dgm_nonhydro3d_common_drhot2pres( pres, dpres, &
      drhot, pres_hyd, therm_hyd, rtot, cvtot, cptot,           &
      lcmesh, elem3d                                            )
    
    !$omp parallel do private( ke2D, DENS, mom_u1, mom_u2 )
    do ke=lcmesh%NeS, lcmesh%NeE
      ke2d = lcmesh%EMap3Dto2D(ke)
 
      dens(:) = dens_hyd(:,ke) + ddens(:,ke)
      mom_u1(:) = lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,1,1) * momx(:,ke) + lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,1) * momy(:,ke)
      mom_u2(:) = lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,1) * momx(:,ke) + lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,2) * momy(:,ke)
 
      entot(:,ke) = pres(:,ke) * cvtot(:,ke) / rtot(:,ke)                                                  &
                  + 0.5_rp * ( momx(:,ke) * mom_u1(:) + momy(:,ke) * mom_u2(:) + momz(:,ke)**2 ) / dens(:) &
                  + grav * dens(:) * lcmesh%zlev(:,ke)
    end do
 
    return

References atm_dyn_dgm_nonhydro3d_common_drhot2pres().

Referenced by mod_dg_prep::dg_prep().

atm_dyn_dgm_nonhydro3d_common_entot2pres()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_entot2pres	(	real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	pres,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	dpres,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	ddens,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momx,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momy,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	momz,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	entot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	pres_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	dens_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	rtot,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	cvtot,
		class(localmesh3d), intent(in)	lcmesh,
		class(elementbase3d), intent(in)	elem3d )

Definition at line 536 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    use scale_const, only: &
      grav => const_grav
    
    implicit none
    class(LocalMesh3D), intent(in) :: lcmesh
    class(ElementBase3D), intent(in) :: elem3D
    real(RP), intent(out) :: PRES(elem3D%Np,lcmesh%NeA)
    real(RP), intent(out) :: DPRES(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: DDENS(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMX(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMY(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: MOMZ(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: EnTot(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: PRES_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: DENS_hyd(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: Rtot(elem3D%Np,lcmesh%NeA)
    real(RP), intent(in) :: CVtot(elem3D%Np,lcmesh%NeA)
 
    integer :: ke, ke2D
 
    real(RP) :: DENS(elem3D%Np)
    real(RP) :: mom_u1(elem3D%Np), mom_u2(elem3D%Np)
    !---------------------------------------------------------------
 
    !$omp parallel do private( ke2D, DENS, mom_u1, mom_u2 )
    do ke=lcmesh%NeS, lcmesh%NeE
      ke2d = lcmesh%EMap3Dto2D(ke)
 
      dens(:) = dens_hyd(:,ke) + ddens(:,ke)
      mom_u1(:) = lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,1,1) * momx(:,ke) + lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,1) * momy(:,ke)
      mom_u2(:) = lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,1) * momx(:,ke) + lcmesh%G_ij(lcmesh%refElem3D%IndexH2Dto3D,ke2d,2,2) * momy(:,ke)
 
      pres(:,ke) = (  entot(:,ke) - grav * dens(:) * lcmesh%zlev(:,ke)                                        &
                  - 0.5_rp * ( momx(:,ke) * mom_u1(:) + momy(:,ke) * mom_u2(:) + momz(:,ke)**2 ) / dens(:) &
                ) * rtot(:,ke) / cvtot(:,ke)
      
      dpres(:,ke) = pres(:,ke) - pres_hyd(:,ke)
    end do
    
    return

Referenced by atm_dyn_dgm_nonhydro3d_common_calc_pressure().

atm_dyn_dgm_nonhydro3d_common_calc_rhot_hyd()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_calc_rhot_hyd	(	real(rp), dimension(elem3d%np,lcmesh%nea), intent(out)	rhot_hyd,
		real(rp), dimension(elem3d%np,lcmesh%nea), intent(in)	pres_hyd,
		class(localmesh3d), intent(in)	lcmesh,
		class(elementbase3d), intent(in)	elem3d )

Definition at line 584 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

 
    use scale_const, only: &
      rdry => const_rdry,   &
      cpdry => const_cpdry, &
      cvdry => const_cvdry, &
      pres00 => const_pre00
        
    implicit none
    class(LocalMesh3D), intent(in) :: lcmesh
    class(ElementBase3D), intent(in) :: elem3D
    real(RP), intent(out) :: RHOT_hyd(elem3D%Np,lcmesh%NeA) 
    real(RP), intent(in) :: PRES_hyd(elem3D%Np,lcmesh%NeA)
 
    integer :: ke, p
    real(RP) :: rP0
    !---------------------------------------------------------------
 
    rp0 = 1.0_rp / pres00
 
    !$omp parallel do
    !$acc parallel loop gang present( RHOT_hyd, PRES_hyd, lcmesh, elem3D )
    do ke=lcmesh%NeS, lcmesh%NeE
#ifdef _OPENACC
      do p=1, elem3d%Np
        rhot_hyd(p,ke) = pres00 / rdry * ( pres_hyd(p,ke) / pres00 )**(cvdry/cpdry) 
      end do
#else
      rhot_hyd(:,ke) = pres00 / rdry * ( pres_hyd(:,ke) / pres00 )**(cvdry/cpdry) 
#endif
    end do
 
    return

Referenced by atm_dyn_dgm_nonhydro3d_common_calc_therm_phyd(), and mod_dg_prep::dg_prep().

atm_dyn_dgm_nonhydro3d_common_calc_phyd_hgrad_lc()

subroutine, public scale_atm_dyn_dgm_nonhydro3d_common::atm_dyn_dgm_nonhydro3d_common_calc_phyd_hgrad_lc	(	real(rp), dimension(elem%np,lmesh%nea), intent(out)	dphyddx,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	dphyddy,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	pres_hyd,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	pres_hyd_ref,
		class(elementoperationbase3d), intent(in)	element3d_operation,
		class(localmesh3d), intent(in)	lmesh,
		class(elementbase3d), intent(in)	elem )

Calculate horizontal graidient of hydrostatic pressure In this calculation, we assume that PRES_hyd_ref is continuous at element boundaries.

Definition at line 624 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

    implicit none
    class(LocalMesh3D), intent(in) :: lmesh
    class(ElementBase3D), intent(in) :: elem
    real(RP), intent(out)  :: DPhydDx(elem%Np,lmesh%NeA)
    real(RP), intent(out)  :: DPhydDy(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: PRES_hyd(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: PRES_hyd_ref(elem%Np,lmesh%NeA)
    class(ElementOperationBase3D), intent(in) :: element3D_operation
 
    integer :: ke, ke2D, p  
 
    real(RP) :: Flux(elem%Np,3), Fz(elem%Np), DFlux(elem%Np,4,2)
    real(RP) :: del_flux_hyd(elem%NfpTot,2,lmesh%Ne)
    real(RP) :: GsqrtV, RGsqrtV(elem%Np)
 
    real(RP) :: E33
    real(RP) :: GradPhyd_x, GradPhyd_y
    !-----------------------------------------
 
    call get_phyd_hgrad_numflux_generalhvc( del_flux_hyd,                       & ! (out)
      pres_hyd, pres_hyd_ref,                                                   & ! (in)
      lmesh%Gsqrt, lmesh%GsqrtH, lmesh%gam, lmesh%GI3(:,:,1), lmesh%GI3(:,:,2), & ! (in)
      lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3),   & ! (in)
      lmesh%vmapM, lmesh%vmapP, elem%IndexH2Dto3D_bnd,                          & ! (in)
      lmesh, elem, lmesh%lcmesh2D, lmesh%lcmesh2D%refElem2D                     ) ! (in)
 
    !$omp parallel do private( ke, ke2D, p, &
    !$omp Flux, Fz, DFlux, GsqrtV, RGsqrtV, E33, &
    !$omp GradPhyd_x, GradPhyd_y )      
    do ke = lmesh%NeS, lmesh%NeE
      ke2d = lmesh%EMap3Dto2D(ke)
 
      do p=1, elem%Np
        gsqrtv  = lmesh%Gsqrt(p,ke) / ( lmesh%gam(p,ke)**2 * lmesh%GsqrtH(elem%IndexH2Dto3D(p),ke2d) )
 
        rgsqrtv(p) = 1.0_rp / gsqrtv
        flux(p,1) = gsqrtv * ( pres_hyd(p,ke) - pres_hyd_ref(p,ke) )
      end do
 
      do p=1, elem%Np        
        flux(p,2) = flux(p,1)
        flux(p,3) = lmesh%GI3(p,ke,1) * flux(p,1)
        fz(p)   = lmesh%GI3(p,ke,2) * flux(p,1)
      end do   
      
      call element3d_operation%Div( flux, del_flux_hyd(:,1,ke), &
        dflux(:,:,1) )
      call element3d_operation%Dz( fz, dflux(:,3,2) )
      call element3d_operation%Lift( del_flux_hyd(:,2,ke), dflux(:,4,2) )
      
      do p=1, elem%Np
        e33 = lmesh%Escale(p,ke,3,3)
 
        gradphyd_x = lmesh%Escale(p,ke,1,1) * dflux(p,1,1) &
                   + e33                    * dflux(p,3,1) &
                   + dflux(p,4,1)
 
        gradphyd_y = lmesh%Escale(p,ke,2,2) * dflux(p,2,1) &
                   + e33                    * dflux(p,3,2) &
                   + dflux(p,4,2)
 
        dphyddx(p,ke) = gradphyd_x * rgsqrtv(p)
        dphyddy(p,ke) = gradphyd_y * rgsqrtv(p)
      end do      
    end do
 
    return

Variable Documentation

prgvar_ddens_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_ddens_id = 1

Definition at line 68 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_DDENS_ID   = 1

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_atmos_component::atmos_setup(), mod_atmos_dyn::atmosdyn_setup(), mod_atmos_phy_tb::atmosphytb_setup(), and mod_atmos_phy_preproc::physpreprocnone_init().

prgvar_therm_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_therm_id = 2

Definition at line 69 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_THERM_ID   = 2 ! Variable associated with energy equation (DRHOT or ETOT)

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_phy_preproc::physpreprocnone_init().

prgvar_drhot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_drhot_id = 2

Definition at line 70 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_DRHOT_ID   = 2

prgvar_etot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_etot_id = 2

Definition at line 71 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_ETOT_ID    = 2

prgvar_momz_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_momz_id = 3

Definition at line 72 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_MOMZ_ID    = 3

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_atmos_component::atmos_setup(), and mod_atmos_phy_preproc::physpreprocnone_init().

prgvar_momx_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_momx_id = 4

Definition at line 73 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_MOMX_ID    = 4

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_atmos_component::atmos_setup(), and mod_atmos_phy_preproc::physpreprocnone_init().

prgvar_momy_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_momy_id = 5

Definition at line 74 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_MOMY_ID    = 5

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_atmos_component::atmos_setup(), and mod_atmos_phy_preproc::physpreprocnone_init().

prgvar_scalar_num

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_scalar_num = 3

Definition at line 75 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_SCALAR_NUM = 3

prgvar_hvec_num

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_hvec_num = 1

Definition at line 76 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_HVEC_NUM   = 1

prgvar_num

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::prgvar_num = 5

Definition at line 77 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PRGVAR_NUM        = 5

Referenced by atm_dyn_dgm_nonhydro3d_common_setup_variables(), mod_atmos_dyn::atmosdyn_setup(), mod_atmos_phy_preproc::physpreprocmodalbase_operate::physpreprocmodalbase_operate(), and mod_atmos_phy_preproc::physpreprocnone_init().

phytend_dens_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_dens_id = 1

Definition at line 79 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_DENS_ID     = 1

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_momx_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_momx_id = 2

Definition at line 80 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_MOMX_ID     = 2

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_momy_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_momy_id = 3

Definition at line 81 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_MOMY_ID     = 3

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_momz_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_momz_id = 4

Definition at line 82 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_MOMZ_ID     = 4

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_rhot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_rhot_id = 5

Definition at line 83 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_RHOT_ID     = 5

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_rhoh_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_rhoh_id = 6

Definition at line 84 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_RHOH_ID     = 6

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

phytend_num

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::phytend_num = 6

Definition at line 85 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: PHYTEND_NUM         = 6

Referenced by atm_dyn_dgm_nonhydro3d_common_setup_variables(), and mod_atmos_component::atmos_setup().

auxvar_preshydro_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_preshydro_id = 1

Definition at line 88 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_PRESHYDRO_ID     = 1

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_experiment::exp_setinitcond_lc::exp_setinitcond_lc(), and mod_atmos_phy_preproc::physpreprocnone_init().

auxvar_denshydro_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_denshydro_id = 2

Definition at line 89 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_DENSHYDRO_ID     = 2

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), mod_atmos_component::atmos_setup(), mod_experiment::exp_setinitcond_lc::exp_setinitcond_lc(), and mod_atmos_phy_preproc::physpreprocnone_init().

auxvar_thermhydro_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_thermhydro_id = 3

Definition at line 90 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_THERMHYDRO_ID    = 3

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_pres_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_pres_id = 4

Definition at line 91 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_PRES_ID          = 4

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

auxvar_pt_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_pt_id = 5

Definition at line 92 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_PT_ID            = 5

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_rtot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_rtot_id = 6

Definition at line 93 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_Rtot_ID          = 6

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo(), and mod_atmos_component::atmos_setup().

auxvar_cvtot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_cvtot_id = 7

Definition at line 94 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_CVtot_ID         = 7

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_cptot_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_cptot_id = 8

Definition at line 95 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_CPtot_ID         = 8

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_qdry_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_qdry_id = 9

Definition at line 96 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_Qdry_ID          = 9

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_preshydro_ref_id

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_preshydro_ref_id = 10

Definition at line 97 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_PRESHYDRO_REF_ID = 10

Referenced by atm_dyn_dgm_nonhydro3d_common_get_varinfo().

auxvar_num

integer, parameter, public scale_atm_dyn_dgm_nonhydro3d_common::auxvar_num = 10

Definition at line 98 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

  integer, public, parameter :: AUXVAR_NUM              = 10

Referenced by atm_dyn_dgm_nonhydro3d_common_setup_variables(), and mod_atmos_phy_preproc::physpreprocmodalbase_operate::physpreprocmodalbase_operate().

intrpmat_vpordm1

real(rp), dimension(:,:), allocatable, public scale_atm_dyn_dgm_nonhydro3d_common::intrpmat_vpordm1

Definition at line 102 of file scale_atm_dyn_dgm_nonhydro3d_common.F90.

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

Referenced by scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve::atm_dyn_dgm_globalnh3d_rhot_heve_cal_tend_shallow_atm_asis(), scale_atm_dyn_dgm_globalnonhydro3d_etot_heve::atm_dyn_dgm_globalnonhydro3d_etot_heve_cal_tend(), scale_atm_dyn_dgm_globalnonhydro3d_etot_hevi::atm_dyn_dgm_globalnonhydro3d_etot_hevi_cal_vi(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve::atm_dyn_dgm_globalnonhydro3d_rhot_heve_cal_tend_deep_atm(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu::atm_dyn_dgm_globalnonhydro3d_rhot_heve_gpu_cal_tend_deep_atm(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_hevi::atm_dyn_dgm_globalnonhydro3d_rhot_hevi_cal_vi(), scale_atm_dyn_dgm_globalnonhydro3d_rhot_hevi::atm_dyn_dgm_globalnonhydro3d_rhot_hevi_cal_vi_asis(), atm_dyn_dgm_nonhydro3d_common_final(), atm_dyn_dgm_nonhydro3d_common_init(), scale_atm_dyn_dgm_nonhydro3d_etot_heve::atm_dyn_dgm_nonhydro3d_etot_heve_cal_tend(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi::atm_dyn_dgm_nonhydro3d_etot_hevi_cal_vi(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi_common::atm_dyn_dgm_nonhydro3d_etot_hevi_common_construct_matbnd(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi_common::atm_dyn_dgm_nonhydro3d_etot_hevi_common_construct_matbnd_uv(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi_common::atm_dyn_dgm_nonhydro3d_etot_hevi_common_eval_ax(), scale_atm_dyn_dgm_nonhydro3d_etot_hevi_common::atm_dyn_dgm_nonhydro3d_etot_hevi_common_eval_ax_uv(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve::atm_dyn_dgm_nonhydro3d_rhot_heve_cal_tend_asis(), scale_atm_dyn_dgm_nonhydro3d_rhot_heve_splitform::atm_dyn_dgm_nonhydro3d_rhot_heve_splitform_cal_tend(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi::atm_dyn_dgm_nonhydro3d_rhot_hevi_cal_vi(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi::atm_dyn_dgm_nonhydro3d_rhot_hevi_cal_vi_asis(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_common::atm_dyn_dgm_nonhydro3d_rhot_hevi_common_construct_matbnd(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_common::atm_dyn_dgm_nonhydro3d_rhot_hevi_common_construct_matbnd_uv(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_common::atm_dyn_dgm_nonhydro3d_rhot_hevi_common_eval_ax(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_common::atm_dyn_dgm_nonhydro3d_rhot_hevi_common_eval_ax_uv(), scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_common_2::atm_dyn_dgm_nonhydro3d_rhot_hevi_common_solve(), and scale_atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform::atm_dyn_dgm_nonhydro3d_rhot_hevi_splitform_cal_vi().