module ATMOSPHERIC Variables More...

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

subroutine, public	atmos_vars_fillhalo (FILL_BND)
	HALO Communication. More...

subroutine, public	atmos_vars_restart_open
	Open restart file for reading atmospheric variables. More...

subroutine, public	atmos_vars_restart_read
	Read restart of atmospheric variables. More...

subroutine, public	atmos_vars_history_setpres
	Set pressure for history output. More...

subroutine, public	atmos_vars_restart_check
	Check and compare between last data and sample data. More...

subroutine, public	atmos_vars_history
	History output set for atmospheric variables. More...

subroutine, public	atmos_vars_total
	Budget monitor for atmosphere. More...

subroutine, public	atmos_vars_calc_diagnostics
	Calc diagnostic variables. More...

recursive subroutine	atmos_vars_get_diagnostic_3d (vname, var)
	get diagnostic variable 3D More...

recursive subroutine	atmos_vars_get_diagnostic_2d (vname, var)
	get diagnostic variable 2D More...

recursive subroutine	atmos_vars_get_diagnostic_1d (vname, var)
	get diagnostic variable 1D More...

subroutine, public	atmos_vars_monitor
	monitor output More...

subroutine, public	atmos_vars_restart_create
	Create atmospheric restart file. More...

subroutine, public	atmos_vars_restart_enddef
	Exit netCDF define mode. More...

subroutine, public	atmos_vars_restart_close
	Close restart file. More...

subroutine, public	atmos_vars_restart_def_var
	Define atmospheric variables in restart file. More...

subroutine, public	atmos_vars_restart_write
	Write restart of atmospheric variables. More...

subroutine	allocate_3d (ary)

subroutine	allocate_2d (ary)

subroutine	allocate_1d (ary)

Variables
logical, public	atmos_restart_output = .false.
	Output restart file? More...

character(len=h_long), public	atmos_restart_in_basename = ''
	Basename of the input file. More...

logical, public	atmos_restart_in_aggregate
	Switch to use aggregate file. More...

logical, public	atmos_restart_in_postfix_timelabel = .false.
	Add timelabel to the basename of input file? More...

character(len=h_long), public	atmos_restart_out_basename = ''
	Basename of the output file. More...

logical, public	atmos_restart_out_aggregate
	Switch to use aggregate file. More...

logical, public	atmos_restart_out_postfix_timelabel = .true.
	Add timelabel to the basename of output file? More...

character(len=h_mid), public	atmos_restart_out_title = 'ATMOS restart'
	Title of the output file. More...

character(len=h_short), public	atmos_restart_out_dtype = 'DEFAULT'
	REAL4 or REAL8. More...

logical, public	atmos_restart_check = .false.
	Check value consistency? More...

character(len=h_long), public	atmos_restart_check_basename = 'restart_check'

real(rp), public	atmos_restart_check_criterion = 1.E-6_RP

real(rp), dimension(:,:,:), allocatable, target, public	dens

real(rp), dimension(:,:,:), allocatable, target, public	momz

real(rp), dimension(:,:,:), allocatable, target, public	momx

real(rp), dimension(:,:,:), allocatable, target, public	momy

real(rp), dimension(:,:,:), allocatable, target, public	rhot

real(rp), dimension(:,:,:,:), allocatable, target, public	qtrc

real(rp), dimension(:,:,:), allocatable, target, public	dens_avw

real(rp), dimension(:,:,:), allocatable, target, public	momz_avw

real(rp), dimension(:,:,:), allocatable, target, public	momx_avw

real(rp), dimension(:,:,:), allocatable, target, public	momy_avw

real(rp), dimension(:,:,:), allocatable, target, public	rhot_avw

real(rp), dimension(:,:,:,:), allocatable, target, public	qtrc_avw

real(rp), dimension(:,:,:), pointer, public	dens_av

real(rp), dimension(:,:,:), pointer, public	momz_av

real(rp), dimension(:,:,:), pointer, public	momx_av

real(rp), dimension(:,:,:), pointer, public	momy_av

real(rp), dimension(:,:,:), pointer, public	rhot_av

real(rp), dimension(:,:,:,:), pointer, public	qtrc_av

real(rp), dimension(:,:,:), allocatable, pointer, target, public	qv

real(rp), dimension(:,:,:), pointer, public	qc

real(rp), dimension(:,:,:), pointer, public	qr

real(rp), dimension(:,:,:), pointer, public	qi

real(rp), dimension(:,:,:), pointer, public	qs

real(rp), dimension(:,:,:), pointer, public	qg

real(rp), dimension(:,:,:), pointer, public	qh

real(rp), dimension(:,:,:,:), allocatable, target, public	qe

real(rp), allocatable, target, public	mass

real(rp), allocatable, target, public	ratio

real(rp), allocatable, target, public	of

real(rp), allocatable, target, public	hydrometors

real(rp), allocatable, target, public	kg

real(rp), dimension(:,:,:), allocatable, public	dens_ref

real(rp), dimension(:,:,:), allocatable, public	pott_ref

real(rp), dimension(:,:,:), allocatable, public	temp_ref

real(rp), dimension(:,:,:), allocatable, public	pres_ref

real(rp), dimension(:,:,:), allocatable, public	qv_ref

real(rp), dimension(:,:,:), allocatable, public	dens_tp

real(rp), dimension(:,:,:), allocatable, public	momz_tp

real(rp), dimension(:,:,:), allocatable, public	rhou_tp

real(rp), dimension(:,:,:), allocatable, public	rhov_tp

real(rp), dimension(:,:,:), allocatable, public	rhot_tp

real(rp), dimension(:,:,:), allocatable, public	rhoh_p

real(rp), dimension(:,:,:,:), allocatable, public	rhoq_tp

real(rp), dimension(:,:,:), allocatable, public	momx_tp

real(rp), dimension(:,:,:), allocatable, public	momy_tp

real(rp), dimension(:,:,:), allocatable, target, public	w

real(rp), allocatable, target, public	velocity

real(rp), allocatable, target, public	m

real(rp), allocatable, target, public	s

real(rp), dimension(:,:,:), allocatable, target, public	u

real(rp), dimension(:,:,:), allocatable, target, public	v

real(rp), dimension(:,:,:), allocatable, target, public	pott

real(rp), allocatable, target, public	potential

real(rp), allocatable, target, public	temperature

real(rp), allocatable, target, public	k

real(rp), dimension(:,:,:), allocatable, target, public	temp

real(rp), dimension(:,:,:), allocatable, target, public	pres

real(rp), allocatable, target, public	pressure

real(rp), allocatable, target, public	pa =J/m3]

real(rp), dimension(:,:,:), allocatable, target, public	exner

real(rp), dimension(t/pt), allocatable, target, public	function

real(rp), dimension(:,:,:), allocatable, target, public	phyd

real(rp), allocatable, target, public	hydrostatic

real(rp), dimension(:,:,:), allocatable, target, public	phydh

real(rp), allocatable, target, public	layer

real(rp), allocatable, target, public	interface

real(rp), dimension(:,:,:), allocatable, target, public	qdry

real(rp), allocatable, target, public	dry

real(rp), allocatable, target, public	air

real(rp), dimension(:,:,:), allocatable, target, public	rtot

real(rp), allocatable, target, public	specific

real(rp), allocatable, target, public	gass

real(rp), allocatable, target, public	constant

real(rp), allocatable, target, public	j

real(rp), dimension(:,:,:), allocatable, target, public	cvtot

real(rp), allocatable, target, public	heat

real(rp), dimension(:,:,:), allocatable, target, public	cptot

Detailed Description

module ATMOSPHERIC Variables

Description: Container for atmospheric variables

Author: Team SCALE

NAMELIST

PARAM_ATMOS_VARS

name	type	default value	comment
ATMOS_RESTART_IN_BASENAME	character(len=H_LONG)	''	Basename of the input file
ATMOS_RESTART_IN_AGGREGATE	logical		Switch to use aggregate file
ATMOS_RESTART_IN_POSTFIX_TIMELABEL	logical	.false.	Add timelabel to the basename of input file?
ATMOS_RESTART_IN_CHECK_COORDINATES	logical	.true.
ATMOS_RESTART_OUTPUT	logical	.false.	Output restart file?
ATMOS_RESTART_OUT_BASENAME	character(len=H_LONG)	''	Basename of the output file
ATMOS_RESTART_OUT_AGGREGATE	logical		Switch to use aggregate file
ATMOS_RESTART_OUT_POSTFIX_TIMELABEL	logical	.true.	Add timelabel to the basename of output file?
ATMOS_RESTART_OUT_TITLE	character(len=H_MID)	'ATMOS restart'	Title of the output file
ATMOS_RESTART_OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8
ATMOS_RESTART_CHECK	logical	.false.	Check value consistency?
ATMOS_RESTART_CHECK_BASENAME	character(len=H_LONG)	'restart_check'
ATMOS_RESTART_CHECK_CRITERION	real(RP)	1.E-6_RP
ATMOS_VARS_CHECKRANGE	logical	.false.
ATMOS_VARS_CHECKCFL_SOFT	real(RP)	1.0_RP	> if Courant number exceeds this value, put the warning message
ATMOS_VARS_CHECKCFL_HARD	real(RP)	2.0_RP	> if Courant number exceeds this value, put the error message and stop

History Output

name	description	unit	variable
"GPH"	"geopotential height"	"m"	"GPH"
W	velocity w	m/s	W
U	velocity u	m/s	U
V	velocity v	m/s	V
PT	potential temp.	K	PT
T	temperature	K	T
PRES	pressure	Pa	PRES
EXNER	Exner function	1	EXNER
PHYD	hydrostatic pressure	Pa	PHYD
QDRY	dry air	kg/kg	QDRY
RTOT	Total gas constant	J/kg/K	RTOT
CVTOT	Total heat capacity	J/kg/K	CVTOT
CPTOT	Total heat capacity	J/kg/K	CPTOT
LHV	latent heat for vaporization	J/kg	LHV
LHS	latent heat for sublimation	J/kg	LHS
LHF	latent heat for fusion	J/kg	LHF
POTV	virtual potential temp.	K	POTV
TEML	liquid water temperature	K	TEML
POTL	liquid water potential temp.	K	POTL
POTE	equivalent potential temp.	K	POTE
QTOT	total water	kg/kg	QTOT
QHYD	total hydrometeors	kg/kg	QHYD
QLIQ	total liquid water	kg/kg	QLIQ
QICE	total ice water	kg/kg	QICE
LWP	liquid water path	g/m2	LWP
IWP	ice water path	g/m2	IWP
PW	precipitable water	g/m2	PW
PREC	surface precipitation flux	kg/m2/s	PREC
RAIN	surface rain flux	kg/m2/s	RAIN
SNOW	surface snow flux	kg/m2/s	SNOW
QSAT	saturation specific humidity	kg/kg	QSAT
RHA	relative humidity(liq+ice)	%	RHA
RH	relative humidity(liq)	%	RH
RHI	relative humidity(ice)	%	RHI
VOR	vertical vorticity	1/s	VOR
DIV	divergence	1/s	DIV
HDIV	horizontal divergence	1/s	HDIV
Uabs	absolute velocity	m/s	Uabs
N2	squared Brunt-Vaisala frequency	1/s2	N2
PBLH	PBL height	m	PBLH
MSE	moist static energy	m2/s2	MSE
TDEW	dew point	K	TDEW
CAPE	convective avail. pot. energy	m2/s2	CAPE
CIN	convection inhibition	m2/s2	CIN
LCL	lifted condensation level	m	LCL
LFC	level of free convection	m	LFC
LNB	level of neutral buoyancy	m	LNB
ENGT	total energy	J/m3	ENGT
ENGP	potential energy	J/m3	ENGP
ENGK	kinetic energy	J/m3	ENGK
ENGI	internal energy	J/m3	ENGI
DENS_MEAN	horiz. mean of density	kg/m3	DENS_MEAN
W_MEAN	horiz. mean of w	m/s	W_MEAN
U_MEAN	horiz. mean of u	m/s	U_MEAN
V_MEAN	horiz. mean of v	m/s	V_MEAN
PT_MEAN	horiz. mean of pot.	K	PT_MEAN
T_MEAN	horiz. mean of t	K	T_MEAN
QV_MEAN	horiz. mean of QV	1	QV_MEAN
QHYD_MEAN	horiz. mean of QHYD	1	QHYD_MEAN
QLIQ_MEAN	horiz. mean of QLIQ	1	QLIQ_MEAN
QICE_MEAN	horiz. mean of QICE	1	QICE_MEAN
DENS_PRIM	horiz. deviation of density	kg/m3	DENS_PRIM
W_PRIM	horiz. deviation of w	m/s	W_PRIM
U_PRIM	horiz. deviation of u	m/s	U_PRIM
V_PRIM	horiz. deviation of v	m/s	V_PRIM
PT_PRIM	horiz. deviation of pot. temp.	K	PT_PRIM
W_PRIM2	variance of w	m2/s2	W_PRIM2
PT_W_PRIM	resolved scale heat flux	W/s	PT_W_PRIM
W_PRIM3	skewness of w	m3/s3	W_PRIM3
TKE_RS	resolved scale TKE	m2/s2	TKE_RS
VELZ	velocity w at the half level	m/s	VELZ
VELX	velocity u at the half level	m/s	VELX
VELY	velocity v at the half level	m/s	VELY
Umet	eastward velocity	m/s	Umet
Vmet	northward velocity	m/s	Vmet
DENS	density	kg/m3	DENS
MOMZ	momentum z	kg/m2/s	MOMZ
MOMX	momentum x	kg/m2/s	MOMX
MOMY	momentum y	kg/m2/s	MOMY
RHOT	rho * theta	kg/m3*K	RHOT
{TRACER_NAME}	{TRACER_NAME}; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	{TRACER_UNIT}	{TRACER_NAME}

Function/Subroutine Documentation

◆ atmos_vars_setup()

subroutine, public mod_atmos_vars::atmos_vars_setup ( )

Setup.

Definition at line 464 of file mod_atmos_vars.F90.

References mod_atmos_dyn_vars::atmos_dyn_vars_setup(), scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_phy_ae_vars::atmos_phy_ae_vars_setup(), mod_atmos_phy_bl_vars::atmos_phy_bl_vars_setup(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_setup(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_setup(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_setup(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_setup(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_setup(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_setup(), atmos_restart_check, atmos_restart_check_basename, atmos_restart_check_criterion, atmos_restart_in_aggregate, atmos_restart_in_basename, atmos_restart_in_postfix_timelabel, atmos_restart_out_aggregate, atmos_restart_out_basename, atmos_restart_out_dtype, atmos_restart_out_postfix_timelabel, atmos_restart_out_title, atmos_restart_output, mod_atmos_admin::atmos_use_average, scale_const::const_undef, cptot, cvtot, dens, dens_av, dens_avw, dens_tp, exner, scale_file_history::file_history_reg(), scale_atmos_hydrometeor::i_hc, scale_atmos_hydrometeor::i_hg, scale_atmos_hydrometeor::i_hh, scale_atmos_hydrometeor::i_hi, scale_atmos_hydrometeor::i_hr, scale_atmos_hydrometeor::i_hs, scale_atmos_hydrometeor::i_qv, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momx_av, momx_avw, momx_tp, momy, momy_av, momy_avw, momy_tp, momz, momz_av, momz_avw, momz_tp, scale_monitor::monitor_reg(), scale_atmos_hydrometeor::n_hyd, phyd, phydh, pott, scale_prc::prc_abort(), pres, scale_tracer::qa, qc, qdry, qe, qg, qh, qi, qr, qs, qtrc, qtrc_av, qtrc_avw, qv, rhoh_p, rhoq_tp, rhot, rhot_av, rhot_avw, rhot_tp, rhou_tp, rhov_tp, rtot, temp, scale_tracer::tracer_desc, scale_tracer::tracer_name, scale_tracer::tracer_unit, u, v, and w.

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

     use scale_const, only: &
        undef => const_undef
     use scale_prc, only: &
        prc_abort
     use scale_file_history, only: &
        file_history_reg
     use scale_monitor, only: &
        monitor_reg
     use scale_atmos_hydrometeor, only: &
        atmos_hydrometeor_dry, &
        n_hyd, &
        i_qv, &
        i_hc, &
        i_hr, &
        i_hi, &
        i_hs, &
        i_hg, &
        i_hh
     use mod_atmos_admin, only: &
        atmos_use_average
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_setup
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_setup
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_setup
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_setup
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_setup
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_setup
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_setup
     use mod_atmos_phy_bl_vars, only: &
        atmos_phy_bl_vars_setup
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_setup
     implicit none
 
     namelist / param_atmos_vars / &
        atmos_restart_in_basename,           &
        atmos_restart_in_aggregate,          &
        atmos_restart_in_postfix_timelabel,  &
        atmos_restart_in_check_coordinates,  &
        atmos_restart_output,                &
        atmos_restart_out_basename,          &
        atmos_restart_out_aggregate,         &
        atmos_restart_out_postfix_timelabel, &
        atmos_restart_out_title,             &
        atmos_restart_out_dtype,             &
        atmos_restart_check,                 &
        atmos_restart_check_basename,        &
        atmos_restart_check_criterion,       &
        atmos_vars_checkrange,               &
        atmos_vars_checkcfl_soft,            &
        atmos_vars_checkcfl_hard
 
     integer :: ierr
     integer :: iv, iq
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("ATMOS_vars_setup",*) 'Setup'
 
     allocate( dens(ka,ia,ja)    )
     allocate( momz(ka,ia,ja)    )
     allocate( momx(ka,ia,ja)    )
     allocate( momy(ka,ia,ja)    )
     allocate( rhot(ka,ia,ja)    )
     allocate( qtrc(ka,ia,ja,max(qa,1)) )
 
     if ( atmos_use_average ) then
        allocate( dens_avw(ka,ia,ja)    )
        allocate( momz_avw(ka,ia,ja)    )
        allocate( momx_avw(ka,ia,ja)    )
        allocate( momy_avw(ka,ia,ja)    )
        allocate( rhot_avw(ka,ia,ja)    )
        allocate( qtrc_avw(ka,ia,ja,max(qa,1)) )
 
        dens_av => dens_avw
        momz_av => momz_avw
        momx_av => momx_avw
        momy_av => momy_avw
        rhot_av => rhot_avw
        qtrc_av => qtrc_avw
     else
        dens_av => dens
        momz_av => momz
        momx_av => momx
        momy_av => momy
        rhot_av => rhot
        qtrc_av => qtrc
     endif
 
     allocate( dens_tp(ka,ia,ja)    )
     allocate( momz_tp(ka,ia,ja)    )
     allocate( rhou_tp(ka,ia,ja)    )
     allocate( rhov_tp(ka,ia,ja)    )
     allocate( rhot_tp(ka,ia,ja)    )
     allocate( rhoh_p(ka,ia,ja)    )
     allocate( rhoq_tp(ka,ia,ja,max(qa,1)) )
 
     allocate( w(ka,ia,ja) )
     allocate( u(ka,ia,ja) )
     allocate( v(ka,ia,ja) )
     w(:,:,:) = undef
     u(:,:,:) = undef
     v(:,:,:) = undef
 
     allocate( pott(ka,ia,ja) )
     allocate( temp(ka,ia,ja) )
     allocate( pres(ka,ia,ja) )
     allocate( exner(ka,ia,ja) )
     allocate( phyd(ka,ia,ja) )
     allocate( phydh(0:ka,ia,ja) )
     pott(:,:,:) = undef
     temp(:,:,:) = undef
     pres(:,:,:) = undef
     exner(:,:,:) = undef
     phyd(:,:,:) = undef
     phydh(:,:,:) = undef
 
     allocate( qdry(ka,ia,ja) )
     allocate( rtot(ka,ia,ja) )
     allocate( cvtot(ka,ia,ja) )
     allocate( cptot(ka,ia,ja) )
     qdry(:,:,:) = undef
     rtot(:,:,:) = undef
     cvtot(:,:,:) = undef
     cptot(:,:,:) = undef
 
     ! obsolute
     allocate( momx_tp(ka,ia,ja)    )
     allocate( momy_tp(ka,ia,ja)    )
 
 
     momz(1:ks-1,:,:) = 0.0_rp
     momz(ke:ka,:,:) = 0.0_rp
 
     allocate( work3d(ka,ia,ja) )
     allocate( work2d(   ia,ja) )
     allocate( work1d(ka      ) )
 
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_atmos_vars,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("ATMOS_vars_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("ATMOS_vars_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_VARS. Check!'
        call prc_abort
     endif
     log_nml(param_atmos_vars)
 
     log_newline
     log_info("ATMOS_vars_setup",*) 'List of prognostic variables (ATMOS) '
     log_info_cont('(1x,A,A24,A,A48,A,A12,A)') &
                '      |', 'VARNAME                 ','|', &
                'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
     do iv = 1, pv_nmax
        log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                   'NO.',iv,'|',pv_info(iv)%NAME,'|', pv_info(iv)%DESC,'[', pv_info(iv)%UNIT,']'
     enddo
     do iq = 1, qa
        log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                   'NO.',5+iq,'|',tracer_name(iq),'|', tracer_desc(iq),'[', tracer_unit(iq),']'
     enddo
 
     log_newline
     if ( atmos_restart_in_basename /= '' ) then
        log_info("ATMOS_vars_setup",*) 'Restart input?  : YES, file = ', trim(atmos_restart_in_basename)
        log_info("ATMOS_vars_setup",*) 'Add timelabel?  : ', atmos_restart_in_postfix_timelabel
     else
        log_info("ATMOS_vars_setup",*) 'Restart input?  : NO'
     endif
     if (       atmos_restart_output             &
          .AND. atmos_restart_out_basename /= '' ) then
        log_info("ATMOS_vars_setup",*) 'Restart output? : YES, file = ', trim(atmos_restart_out_basename)
        log_info("ATMOS_vars_setup",*) 'Add timelabel?  : ', atmos_restart_out_postfix_timelabel
     else
        log_info("ATMOS_vars_setup",*) 'Restart output? : NO'
        atmos_restart_output = .false.
     endif
 
     if ( atmos_restart_check_basename == '' ) then
        atmos_restart_check = .false.
     endif
 
     if ( atmos_vars_checkcfl_hard > 0.0_rp ) then
        atmos_vars_checkcfl_soft = min( atmos_vars_checkcfl_soft, atmos_vars_checkcfl_hard )
     endif
 
     log_newline
     log_info("ATMOS_vars_setup",*) 'Check restart consistency?          : ', atmos_restart_check
     log_info("ATMOS_vars_setup",*) 'Check value range of variables?     : ', atmos_vars_checkrange
     if ( atmos_vars_checkcfl_soft > 0.0_rp ) then
        log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to warn : ', atmos_vars_checkcfl_soft
     else
        log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to warn : disabled'
     endif
     if ( atmos_vars_checkcfl_hard > 0.0_rp ) then
        log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to stop : ', atmos_vars_checkcfl_hard
     else
        log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to stop : disabled'
     endif
 
     call atmos_dyn_vars_setup
     call atmos_phy_mp_vars_setup
     call atmos_phy_ae_vars_setup
     call atmos_phy_ch_vars_setup
     call atmos_phy_rd_vars_setup
     call atmos_phy_sf_vars_setup
     call atmos_phy_tb_vars_setup
     call atmos_phy_bl_vars_setup
     call atmos_phy_cp_vars_setup
 
 
     ! water content
     if ( atmos_hydrometeor_dry ) then
        allocate( zero(ka,ia,ja) )
 !OCL XFILL
        zero(:,:,:) = 0.0_rp
 
        qv => zero
        qc => zero
        qr => zero
        qi => zero
        qs => zero
        qg => zero
        qh => zero
 
        moist = .false.
     else
        allocate( qe(ka,ia,ja,n_hyd) )
 !OCL XFILL
        qe(:,:,:,:) = undef
 
        qv => qtrc_av(:,:,:,i_qv)
        qc => qe(:,:,:,i_hc)
        qr => qe(:,:,:,i_hr)
        qi => qe(:,:,:,i_hi)
        qs => qe(:,:,:,i_hs)
        qg => qe(:,:,:,i_hg)
        qh => qe(:,:,:,i_hh)
 
        moist = .true.
     end if
 
 
     dv_calculated(dv_nmax) = .false.
 
     !-----< history output setup >-----
     allocate( qp_hist_id( max(qa,1) ) )
     allocate( qp_monit_id( max(qa,1) ) )
     pv_hist_id(:) = -1
     pv_monit_id(:) = -1
     qp_hist_id(:) = -1
     qp_monit_id(:) = -1
     dv_hist_id(:) = -1
     dv_monit_id(:) = -1
 
 
     do iv = 1, pv_nmax
        call file_history_reg( pv_info(iv)%NAME, pv_info(iv)%DESC, pv_info(iv)%UNIT, pv_hist_id(iv), dim_type=pv_info(iv)%dim_type, standard_name=pv_info(iv)%STDNAME )
     end do
 
     do iq = 1, qa
        call file_history_reg( tracer_name(iq), tracer_desc(iq), tracer_unit(iq), qp_hist_id(iq), dim_type='ZXY' )
     enddo
 
     do iv = 1, dv_nmax
        call file_history_reg( dv_info(iv)%NAME, dv_info(iv)%DESC, dv_info(iv)%UNIT, dv_hist_id(iv), dim_type=dv_info(iv)%dim_type, standard_name=dv_info(iv)%STDNAME )
     end do
 
     call file_history_reg( "GPH", "geopotential height", "m", hist_id_gph, dim_type='ZXY', standard_name="geopotential_height" )
 
 
     !-----< monitor output setup >-----
     do iv = 1, pv_nmax
        call monitor_reg( pv_info(iv)%NAME, pv_info(iv)%DESC, trim(pv_info(iv)%UNIT)//"*m3", & ! (in)
                          pv_monit_id(iv),                                                   & ! (out)
                          dim_type=pv_info(iv)%dim_type, isflux=.false.                      ) ! (in)
     end do
     do iq = 1, qa
        call monitor_reg( tracer_name(iq), tracer_desc(iq), tracer_unit(iq)//"*kg", & ! (in)
                          qp_monit_id(iq),                                          & ! (out)
                          dim_type='ZXY', isflux=.false.                            ) ! (in)
     enddo
 
     call monitor_reg( 'QDRY',         'dry air mass',           'kg', & ! (in)
                       dv_monit_id(im_qdry),                           & ! (out)
                       dim_type='ZXY', isflux=.false.                  ) ! (in)
     call monitor_reg( 'QTOT',         'water mass',             'kg', & ! (in)
                       dv_monit_id(im_qtot),                           & ! (out)
                       dim_type='ZXY', isflux=.false.                  ) ! (in)
     call monitor_reg( 'EVAP',         'evaporation at the surface', 'kg', & ! (in)
                       dv_monit_id(im_evap),                           & ! (out)
                       dim_type='XY', isflux=.true.                    ) ! (in)
     call monitor_reg( 'PRCP',         'precipitation',          'kg', & ! (in)
                       dv_monit_id(im_prec),                           & ! (out)
                       dim_type='XY', isflux=.true.                    ) ! (in)
 
     call monitor_reg( 'ENGT',         'total     energy',       'J', & ! (in)
                       dv_monit_id(im_engt),                          & ! (out)
                       dim_type='ZXY', isflux=.false.                 ) ! (in)
     call monitor_reg( 'ENGP',         'potential energy',       'J', & ! (in)
                       dv_monit_id(im_engp),                          & ! (out)
                       dim_type='ZXY', isflux=.false.                 ) ! (in)
     call monitor_reg( 'ENGK',         'kinetic   energy',       'J', & ! (in)
                       dv_monit_id(im_engk),                          & ! (out)
                       dim_type='ZXY', isflux=.false.                 ) ! (in)
     call monitor_reg( 'ENGI',         'internal  energy',       'J', & ! (in)
                       dv_monit_id(im_engi),                          & ! (out)
                       dim_type='ZXY', isflux=.false.                 ) ! (in)
 
     call monitor_reg( 'ENGFLXT',      'total energy flux convergence', 'J', & ! (in)
                       dv_monit_id(im_engflxt),                       & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_SH',    'SFC specific heat flux', 'J', & ! (in)
                       dv_monit_id(im_engsfc_sh),                     & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_LH',    'SFC latent   heat flux', 'J', & ! (in)
                       dv_monit_id(im_engsfc_lh),                     & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_RD',    'SFC net radiation flux', 'J', & ! (in)
                       dv_monit_id(im_engsfc_rd),                     & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGTOA_RD',    'TOA net radiation flux', 'J', & ! (in)
                       dv_monit_id(im_engtoa_rd),                     & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
 
     call monitor_reg( 'ENGSFC_LW_up', 'SFC LW upward   flux',   'J', & ! (in)
                       dv_monit_id(im_engsfc_lw_up),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_LW_dn', 'SFC LW downward flux',   'J', & ! (in)
                       dv_monit_id(im_engsfc_lw_dn),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_SW_up', 'SFC SW upward   flux',   'J', & ! (in)
                       dv_monit_id(im_engsfc_sw_up),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGSFC_SW_dn', 'SFC SW downward flux',   'J', & ! (in)
                       dv_monit_id(im_engsfc_sw_dn),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
 
     call monitor_reg( 'ENGTOA_LW_up', 'TOA LW upward   flux',   'J', & ! (in)
                       dv_monit_id(im_engtoa_lw_up),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGTOA_LW_dn', 'TOA LW downward flux',   'J', & ! (in)
                       dv_monit_id(im_engtoa_lw_dn),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGTOA_SW_up', 'TOA SW upward   flux',   'J', & ! (in)
                       dv_monit_id(im_engtoa_sw_up),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
     call monitor_reg( 'ENGTOA_SW_dn', 'TOA SW downward flux',   'J', & ! (in)
                       dv_monit_id(im_engtoa_sw_dn),                  & ! (out)
                       dim_type='XY', isflux=.true.                   ) ! (in)
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_fillhalo()

subroutine, public mod_atmos_vars::atmos_vars_fillhalo ( logical, intent(in), optional FILL_BND )

HALO Communication.

Definition at line 830 of file mod_atmos_vars.F90.

References dens, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momy, momz, scale_tracer::qa, qtrc, and rhot.

Referenced by atmos_vars_restart_read(), and atmos_vars_restart_write().

     use scale_comm_cartesc, only: &
        comm_vars8, &
        comm_wait
     implicit none
 
     logical, intent(in), optional :: fill_bnd
 
     logical :: fill_bnd_
     integer :: i, j, iq
     !---------------------------------------------------------------------------
 
     fill_bnd_ = .false.
     if ( present(fill_bnd) ) fill_bnd_ = fill_bnd
 
     !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
     do j  = jsb, jeb
     do i  = isb, ieb
        dens(   1:ks-1,i,j) = dens(ks,i,j)
        momz(   1:ks-2,i,j) = momz(ks-1,i,j)
        momx(   1:ks-1,i,j) = momx(ks,i,j)
        momy(   1:ks-1,i,j) = momy(ks,i,j)
        rhot(   1:ks-1,i,j) = rhot(ks,i,j)
        dens(ke+1:ka,  i,j) = dens(ke,i,j)
        momz(ke+1:ka,  i,j) = momz(ke,i,j)
        momx(ke+1:ka,  i,j) = momx(ke,i,j)
        momy(ke+1:ka,  i,j) = momy(ke,i,j)
        rhot(ke+1:ka,  i,j) = rhot(ke,i,j)
     enddo
     enddo
 
     !$omp parallel do private(i,j,iq) OMP_SCHEDULE_ collapse(3)
     do iq = 1, qa
     do j  = jsb, jeb
     do i  = isb, ieb
        qtrc(   1:ks-1,i,j,iq) = qtrc(ks,i,j,iq)
        qtrc(ke+1:ka,  i,j,iq) = qtrc(ke,i,j,iq)
     enddo
     enddo
     enddo
 
     call comm_vars8( dens(:,:,:), 1 )
     call comm_vars8( momz(:,:,:), 2 )
     call comm_vars8( momx(:,:,:), 3 )
     call comm_vars8( momy(:,:,:), 4 )
     call comm_vars8( rhot(:,:,:), 5 )
     call comm_wait ( dens(:,:,:), 1, fill_bnd_ )
     call comm_wait ( momz(:,:,:), 2, fill_bnd_ )
     call comm_wait ( momx(:,:,:), 3, fill_bnd_ )
     call comm_wait ( momy(:,:,:), 4, fill_bnd_ )
     call comm_wait ( rhot(:,:,:), 5, fill_bnd_ )
 
     do iq = 1, qa
        call comm_vars8( qtrc(:,:,:,iq), iq )
     enddo
     do iq = 1, qa
        call comm_wait ( qtrc(:,:,:,iq), iq, fill_bnd_ )
     enddo
 
     return

Here is the caller graph for this function:

◆ atmos_vars_restart_open()

subroutine, public mod_atmos_vars::atmos_vars_restart_open ( )

Open restart file for reading atmospheric variables.

Definition at line 894 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_read().

     use scale_prc, only: &
        prc_abort
     use scale_const, only: &
        grav  => const_grav
     use scale_time, only: &
        time_gettimelabel
     use scale_file_cartesc, only: &
        file_cartesc_open, &
        file_cartesc_check_coordinates
     use mod_atmos_admin, only: &
        atmos_use_average, &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_open
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_open
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_open
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_open
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_open
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_open
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_open
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_open
     use mod_cpl_admin, only: &
        cpl_sw
     implicit none
 
     character(len=19)     :: timelabel
     character(len=H_LONG) :: basename
     !---------------------------------------------------------------------------
 
     log_newline
     log_info("ATMOS_vars_restart_open",*) 'Open restart file (ATMOS) '
 
     if ( atmos_restart_in_basename /= '' ) then
 
        if ( atmos_restart_in_postfix_timelabel ) then
           call time_gettimelabel( timelabel )
           basename = trim(atmos_restart_in_basename)//'_'//trim(timelabel)
        else
           basename = trim(atmos_restart_in_basename)
        endif
 
        log_info("ATMOS_vars_restart_open",*) 'basename: ', trim(basename)
 
        call file_cartesc_open( basename, restart_fid, aggregate=atmos_restart_in_aggregate )
 
        if ( atmos_restart_in_check_coordinates ) then
           call file_cartesc_check_coordinates( restart_fid, atmos=.true. )
        end if
 
     else
        log_error("ATMOS_vars_restart_open",*) 'restart file for atmosphere is not specified. STOP!'
        call prc_abort
     endif
 
     if ( atmos_use_average ) then
        dens_av(:,:,:)   = dens(:,:,:)
        momz_av(:,:,:)   = momz(:,:,:)
        momx_av(:,:,:)   = momx(:,:,:)
        momy_av(:,:,:)   = momy(:,:,:)
        rhot_av(:,:,:)   = rhot(:,:,:)
        qtrc_av(:,:,:,:) = qtrc(:,:,:,:)
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_open
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_open
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_open
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_open
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_open
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_open
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_open
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_open
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_read()

subroutine, public mod_atmos_vars::atmos_vars_restart_read ( )

Read restart of atmospheric variables.

Definition at line 986 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_read().

     use scale_prc, only: &
        prc_abort
     use scale_file, only: &
        file_get_aggregate
     use scale_file_cartesc, only: &
        file_cartesc_read, &
        file_cartesc_flush
     use mod_atmos_admin, only: &
        atmos_use_average, &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_read
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_read
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_read
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_read
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_read
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_read
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_read
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_read
     use mod_cpl_admin, only: &
        cpl_sw
     implicit none
 
     integer  :: i, j, iq
     !---------------------------------------------------------------------------
 
     if ( restart_fid /= -1 ) then
        log_newline
        log_info("ATMOS_vars_restart_read",*) 'Read from restart file (ATMOS) '
 
        call file_cartesc_read( restart_fid, pv_info(i_dens)%NAME, 'ZXY', & ! [IN]
                                dens(:,:,:)                               ) ! [OUT]
        call file_cartesc_read( restart_fid, pv_info(i_momz)%NAME, 'ZHXY', & ! [IN]
                                momz(:,:,:)                                ) ! [OUT]
        call file_cartesc_read( restart_fid, pv_info(i_momx)%NAME, 'ZXHY', & ! [IN]
                                momx(:,:,:)                                ) ! [OUT]
        call file_cartesc_read( restart_fid, pv_info(i_momy)%NAME, 'ZXYH', & ! [IN]
                                momy(:,:,:)                                ) ! [OUT]
        call file_cartesc_read( restart_fid, pv_info(i_rhot)%NAME, 'ZXY', & ! [IN]
                                rhot(:,:,:)                               ) ! [OUT]
 
        do iq = 1, qa
           call file_cartesc_read( restart_fid, tracer_name(iq), 'ZXY', & ! [IN]
                                   qtrc(:,:,:,iq)                       ) ! [OUT]
        enddo
 
        if ( file_get_aggregate(restart_fid) ) then
           call file_cartesc_flush( restart_fid ) ! X/Y halos have been read from file
 
           ! fill k halos
           do j  = 1, ja
           do i  = 1, ia
              dens(   1:ks-1,i,j) = dens(ks,i,j)
              momz(   1:ks-2,i,j) = momz(ks-1,i,j)
              momx(   1:ks-1,i,j) = momx(ks,i,j)
              momy(   1:ks-1,i,j) = momy(ks,i,j)
              rhot(   1:ks-1,i,j) = rhot(ks,i,j)
              dens(ke+1:ka,  i,j) = dens(ke,i,j)
              momz(ke+1:ka,  i,j) = momz(ke,i,j)
              momx(ke+1:ka,  i,j) = momx(ke,i,j)
              momy(ke+1:ka,  i,j) = momy(ke,i,j)
              rhot(ke+1:ka,  i,j) = rhot(ke,i,j)
           enddo
           enddo
        else
           call atmos_vars_fillhalo
        end if
 
        call atmos_vars_total
     else
        log_error("ATMOS_vars_restart_read",*) 'invalid restart file ID for atmosphere. STOP!'
        call prc_abort
     endif
 
     if ( atmos_use_average ) then
        dens_av(:,:,:)   = dens(:,:,:)
        momz_av(:,:,:)   = momz(:,:,:)
        momx_av(:,:,:)   = momx(:,:,:)
        momy_av(:,:,:)   = momy(:,:,:)
        rhot_av(:,:,:)   = rhot(:,:,:)
        qtrc_av(:,:,:,:) = qtrc(:,:,:,:)
     endif
 
     if ( atmos_sw_dyn )    call atmos_dyn_vars_restart_read
     if ( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_read
     if ( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_read
     if ( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_read
     if ( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_read
     if ( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_read
     if ( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_read
     if ( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_read
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_history_setpres()

subroutine, public mod_atmos_vars::atmos_vars_history_setpres ( )

Set pressure for history output.

Definition at line 1098 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_update(), and mod_rm_driver::restart_read().

     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz, &
        real_z1 => atmos_grid_cartesc_real_z1
     use scale_topography, only: &
        topo_zsfc
     use scale_atmos_bottom, only: &
        bottom_estimate => atmos_bottom_estimate
     use scale_file_history_cartesc, only: &
        file_history_cartesc_set_pres
     implicit none
 
     real(RP) :: sfc_dens(ia,ja)
     real(RP) :: sfc_pres(ia,ja)
     !---------------------------------------------------------------------------
 
     call bottom_estimate( ka, ks, ke, ia, isb, ieb, ja, jsb, jeb, &
                           dens_av(:,:,:), pres(:,:,:),                  & ! [IN]
                           real_cz(:,:,:), topo_zsfc(:,:), real_z1(:,:), & ! [IN]
                           sfc_dens(:,:), sfc_pres(:,:)                  ) ! [OUT]
 
     call file_history_cartesc_set_pres( phyd(:,:,:), & ! [IN]
                                         phydh(:,:,:), & ! [IN]
                                         sfc_pres(:,:)    ) ! [IN]
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_check()

subroutine, public mod_atmos_vars::atmos_vars_restart_check ( )

Check and compare between last data and sample data.

Definition at line 1128 of file mod_atmos_vars.F90.

References atmos_restart_check_basename, atmos_restart_check_criterion, dens, scale_file_cartesc::file_cartesc_close(), scale_file_cartesc::file_cartesc_flush(), scale_file_cartesc::file_cartesc_open(), scale_file::file_get_aggregate(), scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momy, momz, scale_prc::prc_myrank, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, qtrc, rhot, and scale_tracer::tracer_name.

Referenced by mod_rm_driver::rm_driver().

     use scale_prc, only: &
        prc_myrank
     use scale_file, only: &
        file_get_aggregate
     use scale_file_cartesc, only: &
        file_cartesc_open, &
        file_cartesc_read, &
        file_cartesc_flush, &
        file_cartesc_close
     implicit none
 
     real(RP) :: dens_check(ka,ia,ja)    ! Density    [kg/m3]
     real(RP) :: momz_check(ka,ia,ja)    ! momentum z [kg/s/m2]
     real(RP) :: momx_check(ka,ia,ja)    ! momentum x [kg/s/m2]
     real(RP) :: momy_check(ka,ia,ja)    ! momentum y [kg/s/m2]
     real(RP) :: rhot_check(ka,ia,ja)    ! DENS * POTT [K*kg/m3]
     real(RP) :: qtrc_check(ka,ia,ja,qa) ! tracer mixing ratio [kg/kg]
 
     character(len=H_LONG) :: basename
 
     logical :: datacheck
     integer :: k, i, j, iq
     integer :: fid
     !---------------------------------------------------------------------------
 
     call prof_rapstart('Debug')
 
     log_info("ATMOS_vars_restart_check",*) 'Compare last Data with ', trim(atmos_restart_check_basename), 'on PE=', prc_myrank
     log_info("ATMOS_vars_restart_check",*) 'criterion = ', atmos_restart_check_criterion
     datacheck = .true.
 
     basename = atmos_restart_check_basename
 
     call file_cartesc_open( basename, fid )
 
     call file_cartesc_read( fid, 'DENS', 'ZXY' , dens_check(:,:,:) )
     call file_cartesc_read( fid, 'MOMZ', 'ZHXY', momz_check(:,:,:) )
     call file_cartesc_read( fid, 'MOMX', 'ZXHY', momx_check(:,:,:) )
     call file_cartesc_read( fid, 'MOMY', 'ZXYH', momy_check(:,:,:) )
     call file_cartesc_read( fid, 'RHOT', 'ZXY' , rhot_check(:,:,:) )
     do iq = 1, qa
        call file_cartesc_read( fid, tracer_name(iq), 'ZXY', qtrc_check(:,:,:,iq) )
     end do
     if ( file_get_aggregate(fid) ) call file_cartesc_flush( fid )
 
     call file_cartesc_close( fid ) ! [IN]
 
     do k = ks, ke
     do j = js, je
     do i = is, ie
        if ( abs( dens(k,i,j)-dens_check(k,i,j) ) > atmos_restart_check_criterion ) then
           log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', dens(k,i,j)-dens_check(k,i,j)
           log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'DENS'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do k = ks-1, ke
     do j = js, je
     do i = is, ie
        if ( abs( momz(k,i,j)-momz_check(k,i,j) ) > atmos_restart_check_criterion ) then
           log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', momz(k,i,j)-momz_check(k,i,j)
           log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMZ'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do k = ks, ke
     do j = js, je
     do i = is, ie
        if ( abs( momx(k,i,j)-momx_check(k,i,j) ) > atmos_restart_check_criterion ) then
           log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', momx(k,i,j)-momx_check(k,i,j)
           log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMX'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do k = ks, ke
     do j = js, je
     do i = is, ie
        if ( abs( momy(k,i,j)-momy_check(k,i,j) ) > atmos_restart_check_criterion ) then
           log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', momy(k,i,j)-momy_check(k,i,j)
           log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMY'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do k = ks, ke
     do j = js, je
     do i = is, ie
        if ( abs( rhot(k,i,j)-rhot_check(k,i,j) ) > atmos_restart_check_criterion ) then
           log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', rhot(k,i,j)-rhot_check(k,i,j)
           log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'RHOT'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do iq = 1, qa
        do k = ks, ke
        do j = js, je
        do i = is, ie
           if ( abs( qtrc(k,i,j,iq)-qtrc_check(k,i,j,iq) ) > atmos_restart_check_criterion ) then
              log_error("ATMOS_vars_restart_check",*) 'there is the difference  : ', qtrc(k,i,j,iq)-qtrc_check(k,i,j,iq)
              log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, tracer_name(iq)
              datacheck = .false.
           endif
        enddo
        enddo
        enddo
     enddo
 
     if (datacheck) then
        log_info("ATMOS_vars_restart_check",*) 'Data Check Clear.'
     else
        log_info("ATMOS_vars_restart_check",*) 'Data Check Failed. See std. output.'
        log_error("ATMOS_vars_restart_check",*) 'Data Check Failed.'
     endif
 
     call prof_rapend('Debug')
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_history()

subroutine, public mod_atmos_vars::atmos_vars_history ( )

History output set for atmospheric variables.

Definition at line 1264 of file mod_atmos_vars.F90.

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, mod_atmos_phy_ae_vars::atmos_phy_ae_vars_history(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_history(), cptot, cvtot, dens, dens_av, exner, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, momx, momy, momz, phyd, pott, pres, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, qdry, qtrc, qtrc_av, rhot, rtot, temp, u, v, and w.

Referenced by mod_atmos_driver::atmos_driver_update(), and mod_rm_driver::restart_read().

     use scale_file_history, only: &
        file_history_query, &
        file_history_put
     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_history
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_history
     implicit none
 
     logical :: do_put
     integer :: iq, iv
     !---------------------------------------------------------------------------
 
     call prof_rapstart('ATM_History', 1)
 
     ! value check for prognostic variables
     if ( atmos_vars_checkrange ) then
        call valcheck( dens(:,:,:),    0.0_rp,    2.0_rp, pv_info(i_dens)%NAME, __file__, __line__ )
        call valcheck( momz(:,:,:), -200.0_rp,  200.0_rp, pv_info(i_momz)%NAME, __file__, __line__ )
        call valcheck( momx(:,:,:), -200.0_rp,  200.0_rp, pv_info(i_momx)%NAME, __file__, __line__ )
        call valcheck( momy(:,:,:), -200.0_rp,  200.0_rp, pv_info(i_momy)%NAME, __file__, __line__ )
        call valcheck( rhot(:,:,:),    0.0_rp, 1000.0_rp, pv_info(i_rhot)%NAME, __file__, __line__ )
     endif
 
     ! history output of prognostic variables
     call file_history_put  ( pv_hist_id(i_dens), dens(:,:,:) )
     call file_history_put  ( pv_hist_id(i_momz), momz(:,:,:) )
     call file_history_put  ( pv_hist_id(i_momx), momx(:,:,:) )
     call file_history_put  ( pv_hist_id(i_momy), momy(:,:,:) )
     call file_history_put  ( pv_hist_id(i_rhot), rhot(:,:,:) )
     do iq = 1, qa
        call file_history_put  ( qp_hist_id(iq), qtrc(:,:,:,iq) )
     enddo
 
 
     ! history output of diagnostic variables
     call file_history_put  ( dv_hist_id(i_w    ), w(:,:,:)     )
     call file_history_put  ( dv_hist_id(i_u    ), u(:,:,:)     )
     call file_history_put  ( dv_hist_id(i_v    ), v(:,:,:)     )
     call file_history_put  ( dv_hist_id(i_pott ), pott(:,:,:)  )
     call file_history_put  ( dv_hist_id(i_temp ), temp(:,:,:)  )
     call file_history_put  ( dv_hist_id(i_pres ), pres(:,:,:)  )
 
     call file_history_put  ( dv_hist_id(i_exner), exner(:,:,:) )
     call file_history_put  ( dv_hist_id(i_phyd ), phyd(:,:,:)  )
 
     call file_history_put  ( dv_hist_id(i_qdry ), qdry(:,:,:)  )
     call file_history_put  ( dv_hist_id(i_rtot ), rtot(:,:,:)  )
     call file_history_put  ( dv_hist_id(i_cvtot), cvtot(:,:,:) )
     call file_history_put  ( dv_hist_id(i_cptot), cptot(:,:,:) )
 
     do iv = i_cptot+1, dv_nmax
        if ( dv_hist_id(iv) > 0 ) then
           call file_history_query( dv_hist_id(iv), do_put )
 
           if ( do_put ) then
              select case( dv_info(iv)%ndims )
              case( 3 )
                 call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work3d(:,:,:) )
                 call file_history_put( dv_hist_id(iv), work3d(:,:,:) )
              case( 2 )
                 call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work2d(:,:) )
                 call file_history_put( dv_hist_id(iv), work2d(:,:) )
              case( 1 )
                 call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work1d(:) )
                 call file_history_put( dv_hist_id(iv), work1d(:) )
              end select
           endif
        endif
     enddo
 
     call file_history_put( hist_id_gph, real_cz(:,:,:) )
 
 
     if ( moist ) &
          call atmos_phy_mp_vars_history( dens_av(:,:,:), temp(:,:,:), qtrc_av(:,:,:,:) )
 !    if ( .false. ) then
 !       call ATMOS_vars_get_diagnostic( "RH", WORK3D(:,:,:) )
 !       call ATMOS_PHY_AE_vars_history( QTRC_av(:,:,:,:), WORK3D(:,:,:) )
 !    end if
 
     call prof_rapend  ('ATM_History', 1)
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_total()

subroutine, public mod_atmos_vars::atmos_vars_total ( )

Budget monitor for atmosphere.

Definition at line 1355 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), atmos_vars_restart_read(), and atmos_vars_restart_write().

     use scale_const, only: &
        grav  => const_grav,  &
        cvdry => const_cvdry
     use scale_statistics, only: &
        statistics_checktotal, &
        statistics_total
     use scale_atmos_grid_cartesc_real, only: &
        atmos_grid_cartesc_real_vol,       &
        atmos_grid_cartesc_real_totvol,    &
        atmos_grid_cartesc_real_volwxy,    &
        atmos_grid_cartesc_real_totvolwxy, &
        atmos_grid_cartesc_real_volzuy,    &
        atmos_grid_cartesc_real_totvolzuy, &
        atmos_grid_cartesc_real_volzxv,    &
        atmos_grid_cartesc_real_totvolzxv
     implicit none
 
     real(RP) :: rhoq(ka,ia,ja)
     integer  :: iq
     !---------------------------------------------------------------------------
 
     if ( statistics_checktotal ) then
 
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               dens(:,:,:), pv_info(i_dens)%NAME,    & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol        ) ! (in)
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                              momz(:,:,:), pv_info(i_momz)%NAME,     & ! (in)
                              atmos_grid_cartesc_real_volwxy(:,:,:), & ! (in)
                              atmos_grid_cartesc_real_totvolwxy      ) ! (in)
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               momx(:,:,:), pv_info(i_momx)%NAME,     & ! (in)
                               atmos_grid_cartesc_real_volzuy(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvolzuy      ) ! (in)
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               momy(:,:,:), pv_info(i_momy)%NAME,     & ! (in)
                               atmos_grid_cartesc_real_volzxv(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvolzxv      ) ! (in)
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               rhot(:,:,:), pv_info(i_rhot)%NAME,    & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol        ) ! (in)
 
        do iq = 1, qa
           rhoq(:,:,:) = dens(:,:,:) * qtrc(:,:,:,iq)
 
           call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                  rhoq(:,:,:), tracer_name(iq),       & ! (in)
                                  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                                  atmos_grid_cartesc_real_totvol      ) ! (in)
        enddo
 
        call atmos_vars_calc_diagnostics
 
 
        rhoq(ks:ke,is:ie,js:je) = dens(ks:ke,is:ie,js:je) * qdry(ks:ke,is:ie,js:je)
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               rhoq(:,:,:), 'QDRY',                & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
 
        rhoq(ks:ke,is:ie,js:je) = dens(ks:ke,is:ie,js:je) * ( 1.0_rp - qdry(ks:ke,is:ie,js:je) ) ! Qtotal
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               rhoq(:,:,:), 'QTOT',                & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
 
 
        call atmos_vars_get_diagnostic( 'ENGT', work3d(:,:,:) )
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               work3d(:,:,:), 'ENGT',              & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
        call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               work3d(:,:,:), 'ENGP',              & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
        call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               work3d(:,:,:), 'ENGK',              & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
        call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
        call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                               work3d(:,:,:), 'ENGI',              & ! (in)
                               atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
                               atmos_grid_cartesc_real_totvol      ) ! (in)
 
     endif
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_calc_diagnostics()

subroutine, public mod_atmos_vars::atmos_vars_calc_diagnostics ( )

Calc diagnostic variables.

Definition at line 1453 of file mod_atmos_vars.F90.

References scale_atmos_diagnostic_cartesc::atmos_diagnostic_cartesc_get_vel(), scale_atmos_diagnostic::atmos_diagnostic_get_phyd(), scale_atmos_diagnostic::atmos_diagnostic_get_therm_rhot(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, mod_atmos_phy_ae_vars::atmos_phy_ae_vars_reset_diagnostics(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_get_diagnostic(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_reset_diagnostics(), cptot, cvtot, dens_av, exner, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx_av, momy_av, momz_av, scale_atmos_hydrometeor::n_hyd, phyd, phydh, pott, pres, scale_tracer::qa, qdry, qe, qtrc_av, rhot_av, rtot, temp, scale_tracer::tracer_cp, scale_tracer::tracer_cv, scale_tracer::tracer_mass, scale_tracer::tracer_r, u, v, and w.

Referenced by mod_atmos_driver::atmos_driver_update(), atmos_vars_total(), and mod_rm_driver::restart_read().

     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz, &
        real_fz => atmos_grid_cartesc_real_fz
     use scale_atmos_thermodyn, only: &
        atmos_thermodyn_specific_heat
     use scale_atmos_diagnostic, only: &
        atmos_diagnostic_get_therm_rhot, &
        atmos_diagnostic_get_phyd
     use scale_atmos_diagnostic_cartesc, only: &
        atmos_diagnostic_cartesc_get_vel
     use scale_comm_cartesc, only: &
        comm_vars8, &
        comm_wait
     use scale_atmos_hydrometeor, only: &
        n_hyd
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_get_diagnostic, &
        atmos_phy_mp_vars_reset_diagnostics
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_reset_diagnostics
     implicit none
 
     integer :: iq
 
     call atmos_thermodyn_specific_heat( &
          ka, ks, ke, ia, 1, ia, ja, 1, ja, qa, &
          qtrc_av(:,:,:,:),                                        & ! (in)
          tracer_mass(:), tracer_r(:), tracer_cv(:), tracer_cp(:), & ! (in)
          qdry(:,:,:), rtot(:,:,:), cvtot(:,:,:), cptot(:,:,:)     ) ! (out)
 
     call atmos_diagnostic_cartesc_get_vel( &
          ka, ks, ke, ia, 1, ia, ja, 1, ja, &
          dens_av(:,:,:), momz_av(:,:,:), momx_av(:,:,:), momy_av(:,:,:), & ! (in)
          w(:,:,:), u(:,:,:), v(:,:,:)                                    ) ! (out)
 
     call atmos_diagnostic_get_therm_rhot( &
          ka, ks, ke, ia, 1, ia, ja, 1, ja, &
          dens_av(:,:,:), rhot_av(:,:,:),                     & ! (in)
          rtot(:,:,:), cvtot(:,:,:), cptot(:,:,:),            & ! (in)
          pott(:,:,:), temp(:,:,:), pres(:,:,:), exner(:,:,:) ) ! (out)
 
     call atmos_diagnostic_get_phyd( &
          ka, ks, ke, ia, 1, ia, ja, 1, ja, &
          dens_av(:,:,:), pres(:,:,:),    & ! (in)
          real_cz(:,:,:), real_fz(:,:,:), & ! (in)
          phyd(:,:,:), phydh(:,:,:)       ) ! (out)
 
     call atmos_phy_mp_vars_reset_diagnostics
     call atmos_phy_ae_vars_reset_diagnostics
 
     if ( moist ) then
        call atmos_phy_mp_vars_get_diagnostic( &
             dens_av(:,:,:), temp(:,:,:), qtrc_av(:,:,:,:), & ! (in)
             qe=qe(:,:,:,:)                                 ) ! (out)
        do iq = 1, n_hyd
           call comm_vars8(qe(:,:,:,iq), iq)
        end do
        do iq = 1, n_hyd
           call comm_wait (qe(:,:,:,iq), iq)
        end do
     end if
 
     ! reset diagnostic variables
     dv_calculated(:) = .false.
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_get_diagnostic_3d()

recursive subroutine mod_atmos_vars::atmos_vars_get_diagnostic_3d	(	character(len=*), intent(in)	vname,
		real(rp), dimension(:,:,:), intent(out)	var
	)

get diagnostic variable 3D

Definition at line 1526 of file mod_atmos_vars.F90.

     use scale_const, only: &
        grav  => const_grav,  &
        rvap  => const_rvap,  &
        cpdry => const_cpdry, &
        cvdry => const_cvdry
     use scale_prc, only: &
        prc_abort
     use scale_atmos_grid_cartesc, only: &
        rcdx => atmos_grid_cartesc_rcdx, &
        rcdy => atmos_grid_cartesc_rcdy
     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz, &
        real_fz => atmos_grid_cartesc_real_fz
     use scale_atmos_grid_cartesc_metric, only: &
        rotc => atmos_grid_cartesc_metric_rotc
     use scale_comm_cartesc, only: &
        comm_vars8, &
        comm_wait
     use scale_atmos_hydrometeor, only: &
        lhvc => lhv, &
        lhfc => lhf, &
        atmos_hydrometeor_lhv, &
        atmos_hydrometeor_lhf, &
        atmos_hydrometeor_lhs
     use scale_atmos_saturation, only: &
        atmos_saturation_dens2qsat_all, &
        atmos_saturation_psat_all, &
        atmos_saturation_psat_liq, &
        atmos_saturation_psat_ice, &
        atmos_saturation_tdew_liq, &
        atmos_saturation_pote
     use scale_atmos_diagnostic, only: &
        atmos_diagnostic_get_potv, &
        atmos_diagnostic_get_teml, &
        atmos_diagnostic_get_n2
     implicit none
     character(len=*), intent(in)  :: vname
     real(RP),         intent(out) :: var(:,:,:)
 
     real(RP) :: uh  (ka,ia,ja)
     real(RP) :: vh  (ka,ia,ja)
 
     real(RP) :: work(ka,ia,ja)
 
     integer :: k, i, j, iq
 
     select case ( vname )
     case ( 'W' )
        var(:,:,:) = w(:,:,:)
 
     case ( 'U' )
        var(:,:,:) = u(:,:,:)
 
     case ( 'V' )
        var(:,:,:) = v(:,:,:)
 
     case ( 'PT' )
        var(:,:,:) = pott(:,:,:)
 
     case ( 'T' )
        var(:,:,:) = temp(:,:,:)
 
     case ( 'EXNER' )
        var(:,:,:) = exner(:,:,:)
 
     case ( 'PHYD' )
        var(:,:,:) = phyd(:,:,:)
 
     case ( 'QDRY' )
        var(:,:,:) = qdry(:,:,:)
 
     case ( 'RTOT' )
        var(:,:,:) = rtot(:,:,:)
 
     case ( 'CVTOT' )
        var(:,:,:) = cvtot(:,:,:)
 
     case ( 'CPTOT' )
        var(:,:,:) = cptot(:,:,:)
 
     case ( 'LHV' )
        if ( .not. dv_calculated(i_lhv) ) then
           call allocate_3d( lhv )
           call atmos_hydrometeor_lhv( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), & ! (in)
                lhv(:,:,:)   ) ! (out)
           dv_calculated(i_lhv) = .true.
        end if
        var(ks:ke,:,:) = lhv(ks:ke,:,:)
 
     case ( 'LHS' )
        if ( .not. dv_calculated(i_lhs) ) then
           call allocate_3d( lhs )
           call atmos_hydrometeor_lhs( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), & ! (in)
                lhs(:,:,:)   ) ! (out)
           dv_calculated(i_lhs) = .true.
        end if
        var(ks:ke,:,:) = lhs(ks:ke,:,:)
 
     case ( 'LHF' )
        if ( .not. dv_calculated(i_lhf) ) then
           call allocate_3d( lhf )
           call atmos_hydrometeor_lhf( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), & ! (in)
                lhf(:,:,:)   ) ! (out)
           dv_calculated(i_lhf) = .true.
        end if
        var(ks:ke,:,:) = lhf(ks:ke,:,:)
 
     case ( 'POTV' )
        if ( .not. dv_calculated(i_potv) ) then
           call allocate_3d( potv )
           call atmos_diagnostic_get_potv( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                pott(:,:,:), rtot(:,:,:), & ! (in)
                potv(:,:,:)               ) ! (out)
           dv_calculated(i_potv) = .true.
        end if
        var(ks:ke,:,:) = potv(ks:ke,:,:)
 
     case ( 'TEML' )
        if ( .not. dv_calculated(i_teml) ) then
           call allocate_3d( teml )
           call atmos_vars_get_diagnostic( 'LHV', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'LHS', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
           call atmos_diagnostic_get_teml( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), lhv(:,:,:), lhs(:,:,:), & ! (in)
                qc(:,:,:), qi(:,:,:), cptot(:,:,:),  & ! (in)
                teml(:,:,:)                          ) ! (out)
           dv_calculated(i_teml) = .true.
        end if
        var(ks:ke,:,:) = teml(ks:ke,:,:)
 
     case ( 'POTL' )
        if ( .not. dv_calculated(i_potl) ) then
           call allocate_3d( potl )
           call atmos_vars_get_diagnostic( 'TEML', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k) &
           !$omp shared(POTL,TEML,EXNER) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              potl(k,i,j) = teml(k,i,j) / exner(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_potl) = .true.
        end if
        var(ks:ke,:,:) = potl(ks:ke,:,:)
 
     case ( 'POTE' )
        if ( .not. dv_calculated(i_pote) ) then
           call allocate_3d( pote )
           call atmos_saturation_pote( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                dens(:,:,:), pott(:,:,:), temp(:,:,:), qv(:,:,:), & ! [IN]
                pote(:,:,:)                                       ) ! [OUT]
        end if
        var(ks:ke,:,:) = pote(ks:ke,:,:)
     case ( 'QTOT' )
        if ( .not. dv_calculated(i_qtot) ) then
           call allocate_3d( qtot )
           if ( moist ) then
              call atmos_vars_get_diagnostic( 'QHYD', work3d(:,:,:) )
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(QTOT,QV,QHYD) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 qtot(k,i,j) = qv(k,i,j) + qhyd(k,i,j)
              enddo
              enddo
              enddo
           else
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(QTOT) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 qtot(k,i,j) = 0.0_rp
              enddo
              enddo
              enddo
           end if
           dv_calculated(i_qtot) = .true.
        end if
        var(ks:ke,:,:) = qtot(ks:ke,:,:)
 
     case ( 'QHYD' )
        if ( .not. dv_calculated(i_qhyd) ) then
           call allocate_3d( qhyd )
           if ( moist ) then
              call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
              call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(QHYD,QLIQ,QICE) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 qhyd(k,i,j) = qliq(k,i,j) + qice(k,i,j)
              enddo
              enddo
              enddo
           else
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(QHYD) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 qhyd(k,i,j) = 0.0_rp
              enddo
              enddo
              enddo
           end if
           dv_calculated(i_qhyd) = .true.
        end if
        var(ks:ke,:,:) = qhyd(ks:ke,:,:)
 
     case ( 'QLIQ' )
        if ( .not. dv_calculated(i_qliq) ) then
           call allocate_3d( qliq )
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k,iq) &
           !$omp shared(QLIQ,QC,QR) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
 !OCL XFILL
              qliq(k,i,j) = qc(k,i,j) + qr(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_qliq) = .true.
        end if
        var(ks:ke,:,:) = qliq(ks:ke,:,:)
 
     case ( 'QICE' )
        if ( .not. dv_calculated(i_qice) ) then
           call allocate_3d( qice )
 !OCL XFILL
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k,iq) &
           !$omp shared(QICE,QI,QS,QG,QH) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              qice(k,i,j) = qi(k,i,j) + qs(k,i,j) + qg(k,i,j) + qh(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_qice) = .true.
        end if
        var(ks:ke,:,:) = qice(ks:ke,:,:)
 
     case ( 'QSAT' )
        if ( .not. dv_calculated(i_qsat) ) then
           call allocate_3d( qsat )
           call atmos_saturation_dens2qsat_all( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), dens_av(:,:,:), & ! (in)
                qsat(:,:,:)                  ) ! (out)
           dv_calculated(i_qsat) = .true.
        end if
        var(ks:ke,:,:) = qsat(ks:ke,:,:)
 
     case ( 'RHA' )
        if ( .not. dv_calculated(i_rha) ) then
           call allocate_3d( rha )
           if ( moist ) then
              call atmos_saturation_psat_all( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                temp(:,:,:), & ! (in)
                work(:,:,:)  ) ! (out)
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHA,DENS_av,QV,WORK,TEMP) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rha(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
                            / work(k,i,j) * rvap * temp(k,i,j) &
                            * 100.0_rp
              enddo
              enddo
              enddo
           else
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHA) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rha(k,i,j) = 0.0_rp
              enddo
              enddo
              enddo
           end if
           dv_calculated(i_rha) = .true.
        end if
        var(ks:ke,:,:) = rha(ks:ke,:,:)
 
     case ( 'RHL', 'RH' )
        if ( .not. dv_calculated(i_rhl) ) then
           call allocate_3d( rhl )
           if ( moist ) then
              call atmos_saturation_psat_liq( &
                   ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                   temp(:,:,:), & ! (in)
                   work(:,:,:)  ) ! (out)
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHL,DENS_av,QV,WORK,TEMP) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rhl(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
                            / work(k,i,j) * rvap * temp(k,i,j) &
                            * 100.0_rp
              enddo
              enddo
              enddo
           else
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHL) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rhl(k,i,j) = 0.0_rp
              enddo
              enddo
              enddo
           end if
           dv_calculated(i_rhl) = .true.
        end if
        var(ks:ke,:,:) = rhl(ks:ke,:,:)
 
     case ( 'RHI' )
        if ( .not. dv_calculated(i_rhi) ) then
           call allocate_3d( rhi )
           if ( moist ) then
              call atmos_saturation_psat_ice( &
                   ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                   temp(:,:,:), & ! (int)
                   work(:,:,:)  ) ! (out)
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHI,DENS_av,QV,WORK,TEMP) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rhi(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
                            / work(k,i,j) * rvap * temp(k,i,j) &
                            * 100.0_rp
              enddo
              enddo
              enddo
           else
 !OCL XFILL
              !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
              !$omp private(i,j,k) &
              !$omp shared(RHI) &
              !$omp shared(KS,KE,IA,JA)
              do j = 1, ja
              do i = 1, ia
              do k = ks, ke
                 rhi(k,i,j) = 0.0_rp
              enddo
              enddo
              enddo
           end if
           dv_calculated(i_rhi) = .true.
        end if
        var(ks:ke,:,:) = rhi(ks:ke,:,:)
 
     case ( 'VOR' )
        if ( .not. dv_calculated(i_vor) ) then
           call allocate_3d( vor )
           !!!  to move to grid !!!
           ! at x, v, layer
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
           do j = 1, ja-1
           do i = 2, ia
           do k = ks, ke
              uh(k,i,j) = 0.5_rp * ( momx_av(k,i  ,j) + momx_av(k,i  ,j+1) &
                                   + momx_av(k,i-1,j) + momx_av(k,i-1,j+1) ) &
                        / ( dens_av(k,i,j) + dens_av(k,i,j+1) )
           enddo
           enddo
           enddo
           ! at u, y, layer
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
           do j = 2, ja
           do i = 1, ia-1
           do k = ks, ke
              vh(k,i,j) = 0.5_rp * ( momy_av(k,i,j  ) + momy_av(k,i+1,j  ) &
                                   + momy_av(k,i,j-1) + momy_av(k,i+1,j-1) ) &
                        / ( dens_av(k,i,j) + dens_av(k,i+1,j) )
           enddo
           enddo
           enddo
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
           do j = 2, ja-1
           do i = 2, ia-1
           do k = ks, ke
              vor(k,i,j) = ( vh(k,i,j  ) - vh(k,i-1,j  ) ) * rcdx(i) &
                         - ( uh(k,i  ,j) - uh(k,i  ,j-1) ) * rcdy(j)
           enddo
           enddo
           enddo
           !$omp parallel do private(j,k) OMP_SCHEDULE_
           do j = 1, ja
           do k = ks, ke
              vor(k,1 ,j) = vor(k,2   ,j)
              vor(k,ia,j) = vor(k,ia-1,j)
           enddo
           enddo
           !$omp parallel do private(i,k) OMP_SCHEDULE_
           do i = 1, ia
           do k = ks, ke
              vor(k,i,1 ) = vor(k,i,2   )
              vor(k,i,ja) = vor(k,i,ja-1)
           enddo
           enddo
           dv_calculated(i_vor) = .true.
        end if
        var(ks:ke,:,:) = vor(ks:ke,:,:)
 
     case ( 'DIV' )
        if ( .not. dv_calculated(i_div) ) then
           call allocate_3d( div )
           call atmos_vars_get_diagnostic( 'HDIV', work3d(:,:,:) )
           !!!! to move to grid !!!!
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              div(k,i,j) = ( momz_av(k,i,j) - momz_av(k-1,i  ,j  ) ) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) &
                         + hdiv(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_div) = .true.
        end if
        var(ks:ke,:,:) = div(ks:ke,:,:)
 
     case ( 'HDIV' )
        if ( .not. dv_calculated(i_hdiv) ) then
           call allocate_3d( hdiv )
           !!!! to move to grid !!!!
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
           do j = 2, ja
           do i = 2, ia
           do k = ks, ke
              hdiv(k,i,j) = ( momx_av(k,i,j) - momx_av(k  ,i-1,j  ) ) * rcdx(i) &
                          + ( momy_av(k,i,j) - momy_av(k  ,i  ,j-1) ) * rcdy(j)
           enddo
           enddo
           enddo
           !$omp parallel do private(i,k) OMP_SCHEDULE_
           do i = 1, ia
           do k = ks, ke
              hdiv(k,i,1) = hdiv(k,i,2)
           enddo
           enddo
           !$omp parallel do private(j,k) OMP_SCHEDULE_
           do j = 1, ja
           do k = ks, ke
              hdiv(k,1,j) = hdiv(k,2,j)
           enddo
           enddo
           dv_calculated(i_hdiv) = .true.
        end if
        var(ks:ke,:,:) = hdiv(ks:ke,:,:)
 
     case ( 'Uabs' )
        if ( .not. dv_calculated(i_uabs) ) then
           call allocate_3d( uabs )
 !OCL XFILL
           !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              uabs(k,i,j) = sqrt( w(k,i,j)**2 + u(k,i,j)**2 + v(k,i,j)**2 )
           enddo
           enddo
           enddo
           dv_calculated(i_uabs) = .true.
        end if
        var(ks:ke,:,:) = uabs(ks:ke,:,:)
 
     case ( 'N2' )
        if ( .not. dv_calculated(i_n2) ) then
           call allocate_3d( n2 )
           call atmos_diagnostic_get_n2( &
                ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                pott(:,:,:), rtot(:,:,:), & !(in)
                real_cz(:,:,:),           & !(in)
                n2(:,:,:)                 ) ! (out)
           dv_calculated(i_n2) = .true.
        end if
        var(ks:ke,:,:) = n2(ks:ke,:,:)
 
     case ( 'MSE' )
        if ( .not. dv_calculated(i_mse) ) then
           call allocate_3d( mse )
           call atmos_vars_get_diagnostic( 'LHV', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              mse(k,i,j) = cptot(k,i,j) * temp(k,i,j)                    &
                         + grav * ( real_cz(k,i,j) - real_fz(ks-1,i,j) ) &
                         + lhv(k,i,j) * qv(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_mse) = .true.
        end if
        var(ks:ke,:,:) = mse(ks:ke,:,:)
 
     case ( 'TDEW' )
        if ( .not. dv_calculated(i_tdew) ) then
           call allocate_3d( tdew )
           call atmos_saturation_tdew_liq( ka, ks, ke, ia, 1, ia, ja, 1, ja, &
                                           dens(:,:,:), temp(:,:,:), qv(:,:,:), & ! [IN]
                                           tdew(:,:,:)                          ) ! [OUT]
           dv_calculated(i_tdew) = .true.
        end if
        var(ks:ke,:,:) = tdew(ks:ke,:,:)
 
     case ( 'ENGP' )
        if ( .not. dv_calculated(i_engp) ) then
           call allocate_3d( engp )
           !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              engp(k,i,j) = dens_av(k,i,j) * grav * real_cz(k,i,j)
           end do
           end do
           end do
           dv_calculated(i_engp) = .true.
        end if
        var(ks:ke,:,:) = engp(ks:ke,:,:)
 
     case ( 'ENGK' )
        if ( .not. dv_calculated(i_engk) ) then
           call allocate_3d( engk )
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              engk(k,i,j) = 0.5_rp * dens_av(k,i,j) &
                          * ( w(k,i,j)**2 + u(k,i,j)**2 + v(k,i,j)**2 )
           end do
           end do
           end do
              dv_calculated(i_engk) = .true.
        end if
        var(ks:ke,:,:) = engk(ks:ke,:,:)
 
     case ( 'ENGI' )
        if ( .not. dv_calculated(i_engi) ) then
           call allocate_3d( engi )
           if ( moist ) then
              call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
           end if
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              engi(k,i,j) = dens_av(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
              do iq = 1, qa
                 engi(k,i,j) = engi(k,i,j) &
                             + dens_av(k,i,j) * qtrc_av(k,i,j,iq) * temp(k,i,j) * tracer_cv(iq)
              enddo
              if ( moist ) then
                 engi(k,i,j) = engi(k,i,j) &
                      + dens_av(k,i,j) * ( qv(k,i,j) * lhvc & ! Latent Heat [vapor->liquid]
                                         - qice(k,i,j) * lhfc ) ! Latent Heat [ice->liquid]
              end if
           end do
           end do
           end do
           dv_calculated(i_engi) = .true.
        end if
        var(ks:ke,:,:) = engi(ks:ke,:,:)
 
     case ( 'ENGT' )
        if ( .not. dv_calculated(i_engt) ) then
           call allocate_3d( engt )
           call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              engt(k,i,j) = engp(k,i,j) + engk(k,i,j) + engi(k,i,j)
           enddo
           enddo
           enddo
           dv_calculated(i_engt) = .true.
        end if
        var(ks:ke,:,:) = engt(ks:ke,:,:)
 
     case ( 'DENS_PRIM' )
        if ( .not. dv_calculated(i_dens_prim) ) then
           call allocate_3d( dens_prim )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              dens_prim(k,i,j) = dens_av(k,i,j) - dens_mean(k)
           enddo
           enddo
           enddo
           dv_calculated(i_dens_prim) = .true.
        end if
        var(ks:ke,:,:) = dens_prim(ks:ke,:,:)
 
     case ( 'W_PRIM' )
        if ( .not. dv_calculated(i_w_prim) ) then
           call allocate_3d( w_prim )
           call atmos_vars_get_diagnostic( 'W_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              w_prim(k,i,j) = w(k,i,j) - w_mean(k)
           enddo
           enddo
           enddo
           dv_calculated(i_w_prim) = .true.
        end if
        var(ks:ke,:,:) = w_prim(ks:ke,:,:)
 
     case ( 'U_PRIM' )
        if ( .not. dv_calculated(i_u_prim) ) then
           call allocate_3d( u_prim )
           call atmos_vars_get_diagnostic( 'U_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              u_prim(k,i,j) = u(k,i,j) - u_mean(k)
           enddo
           enddo
           enddo
           dv_calculated(i_u_prim) = .true.
        end if
        var(ks:ke,:,:) = u_prim(ks:ke,:,:)
 
     case ( 'V_PRIM' )
        if ( .not. dv_calculated(i_v_prim) ) then
           call allocate_3d( v_prim )
           call atmos_vars_get_diagnostic( 'V_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              v_prim(k,i,j) = v(k,i,j) - v_mean(k)
           enddo
           enddo
           enddo
           dv_calculated(i_v_prim) = .true.
        end if
        var(ks:ke,:,:) = v_prim(ks:ke,:,:)
 
     case ( 'PT_PRIM' )
        if ( .not. dv_calculated(i_pt_prim) ) then
           call allocate_3d( pt_prim )
           call atmos_vars_get_diagnostic( 'PT_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              pt_prim(k,i,j) = pott(k,i,j) - pt_mean(k)
           enddo
           enddo
           enddo
           dv_calculated(i_pt_prim) = .true.
        end if
        var(ks:ke,:,:) = pt_prim(ks:ke,:,:)
 
     case ( 'W_PRIM2' )
        if ( .not. dv_calculated(i_w_prim2) ) then
           call allocate_3d( w_prim2 )
           call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              w_prim2(k,i,j) = w_prim(k,i,j)**2
           enddo
           enddo
           enddo
           dv_calculated(i_w_prim2) = .true.
        end if
        var(ks:ke,:,:) = w_prim2(ks:ke,:,:)
 
     case ( 'PT_W_PRIM' )
        if ( .not. dv_calculated(i_pt_w_prim) ) then
           call allocate_3d( pt_w_prim )
           call atmos_vars_get_diagnostic( 'W_PRIM',  work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'PT_PRIM', work3d(:,:,:) )
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              pt_w_prim(k,i,j) = w_prim(k,i,j) * pt_prim(k,i,j) * dens_av(k,i,j) * cpdry
           enddo
           enddo
           enddo
           dv_calculated(i_pt_w_prim) = .true.
        end if
        var(ks:ke,:,:) = pt_w_prim(ks:ke,:,:)
 
     case ( 'W_PRIM3' )
        if ( .not. dv_calculated(i_w_prim3) ) then
           call allocate_3d( w_prim3 )
           call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              w_prim3(k,i,j) = w_prim(k,i,j)**3
           enddo
           enddo
           enddo
           dv_calculated(i_w_prim3) = .true.
        end if
        var(ks:ke,:,:) = w_prim3(ks:ke,:,:)
 
     case ( 'TKE_RS' )
        if ( .not. dv_calculated(i_tke_rs) ) then
           call allocate_3d( tke_rs )
           call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'U_PRIM', work3d(:,:,:) )
           call atmos_vars_get_diagnostic( 'V_PRIM', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              tke_rs(k,i,j) = 0.5_rp * ( w_prim(k,i,j)**2 + u_prim(k,i,j)**2 + v_prim(k,i,j)**2 )
           enddo
           enddo
           enddo
           dv_calculated(i_tke_rs) = .true.
        end if
        var(ks:ke,:,:) = tke_rs(ks:ke,:,:)
 
     case ( 'VELZ' )
        if ( .not. dv_calculated(i_velz) ) then
           call allocate_3d( velz )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
              velz(ks-1,i,j) = 0.0_rp
              do k = ks, ke-1
                 velz(k,i,j) = momz(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k+1,i,j) )
              end do
              velz(ke,i,j) = 0.0_rp
           enddo
           enddo
           dv_calculated(i_velz) = .true.
        end if
        var(ks-1:ke,:,:) = velz(ks-1:ke,:,:)
 
     case ( 'VELX' )
        if ( .not. dv_calculated(i_velx) ) then
           call allocate_3d( velx )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia-1
           do k = ks, ke
              velx(k,i,j) = momx(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k,i+1,j) )
           enddo
           enddo
           enddo
 !OCL XFILL
           !$omp parallel do private(j,k) OMP_SCHEDULE_
           do j = 1, ja
           do k = ks, ke
              velx(k,ia,j) = momx(k,ia,j) / dens(k,ia,j)
           enddo
           enddo
           call comm_vars8( velx(:,:,:), 1 )
           call comm_wait ( velx(:,:,:), 1, .false. )
           dv_calculated(i_velx) = .true.
        end if
        var(ks:ke,:,:) = velx(ks:ke,:,:)
 
     case ( 'VELY' )
        if ( .not. dv_calculated(i_vely) ) then
           call allocate_3d( vely )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja-1
           do i = 1, ia
           do k = ks, ke
              vely(k,i,j) = momy(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k,i,j+1) )
           enddo
           enddo
           enddo
 !OCL XFILL
           !$omp parallel do private(i,k) OMP_SCHEDULE_
           do i = 1, ia
           do k = ks, ke
              vely(k,i,ja) = momy(k,i,ja) / dens(k,i,ja)
           enddo
           enddo
           call comm_vars8( vely(:,:,:), 1 )
           call comm_wait ( vely(:,:,:), 1, .false. )
           dv_calculated(i_vely) = .true.
        end if
        var(ks:ke,:,:) = vely(ks:ke,:,:)
 
     case ( 'Umet' )
        if ( .not. dv_calculated(i_umet) ) then
           call allocate_3d( umet )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              umet(k,i,j) = u(k,i,j) * rotc(i,j,1) - v(k,i,j) * rotc(i,j,2)
           end do
           end do
           end do
           dv_calculated(i_umet) = .true.
        end if
        var(ks:ke,:,:) = umet(ks:ke,:,:)
 
     case ( 'Vmet' )
        if ( .not. dv_calculated(i_vmet) ) then
           call allocate_3d( vmet )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = 1, ja
           do i = 1, ia
           do k = ks, ke
              vmet(k,i,j) = u(k,i,j) * rotc(i,j,2) + v(k,i,j) * rotc(i,j,1)
           end do
           end do
           end do
           dv_calculated(i_vmet) = .true.
        end if
        var(ks:ke,:,:) = vmet(ks:ke,:,:)
 
     case default
        log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_3D: ', trim(vname)
        call prc_abort
     end select
 
 
     return

Here is the call graph for this function:

◆ atmos_vars_get_diagnostic_2d()

recursive subroutine mod_atmos_vars::atmos_vars_get_diagnostic_2d	(	character(len=*), intent(in)	vname,
		real(rp), dimension(:,:), intent(out)	var
	)

get diagnostic variable 2D

Definition at line 2442 of file mod_atmos_vars.F90.

     use scale_prc, only: &
        prc_abort
     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz, &
        real_fz => atmos_grid_cartesc_real_fz
     use scale_atmos_adiabat, only: &
        atmos_adiabat_cape
     use mod_atmos_phy_mp_vars, only: &
        sflx_rain_mp => atmos_phy_mp_sflx_rain, &
        sflx_snow_mp => atmos_phy_mp_sflx_snow
     use mod_atmos_phy_cp_vars, only: &
        sflx_rain_cp => atmos_phy_cp_sflx_rain
     implicit none
 
     character(len=*), intent(in)  :: vname
     real(RP),         intent(out) :: var(:,:)
 
     real(RP) :: fact
     integer  :: k, i, j
     !---------------------------------------------------------------------------
 
     select case ( vname )
     case ( 'LWP' )
        if ( .not. dv_calculated(i_lwp) ) then
           call allocate_2d( lwp )
           call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k) &
           !$omp shared(LWP,QLIQ,DENS_av,REAL_FZ) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
              lwp(i,j) = 0.0_rp
              do k  = ks, ke
                 lwp(i,j) = lwp(i,j) &
                          + qliq(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
              enddo
           enddo
           enddo
           dv_calculated(i_lwp) = .true.
        end if
        var(:,:) = lwp(:,:)
 
     case ( 'IWP' )
        if ( .not. dv_calculated(i_iwp) ) then
           call allocate_2d( iwp )
           call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k) &
           !$omp shared(IWP,QICE,DENS_av,REAL_FZ) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
              iwp(i,j) = 0.0_rp
              do k  = ks, ke
                 iwp(i,j) = iwp(i,j) &
                          + qice(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
              enddo
           enddo
           enddo
           dv_calculated(i_iwp) = .true.
        end if
        var(:,:) = iwp(:,:)
 
     case ( 'PW' )
        if ( .not. dv_calculated(i_pw) ) then
           call allocate_2d( pw )
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(i,j,k) &
           !$omp shared(PW,QV,DENS_av,REAL_FZ) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
              pw(i,j) = 0.0_rp
              do k  = ks, ke
                 pw(i,j) = pw(i,j) &
                         + qv(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
              enddo
           enddo
           enddo
           dv_calculated(i_pw) = .true.
        end if
        var(:,:) = pw(:,:)
 
     case ( 'PBLH' )
        if ( .not. dv_calculated(i_pblh) ) then
           call allocate_2d( pblh )
           call atmos_vars_get_diagnostic( 'POTV', work3d(:,:,:) )
           !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
           !$omp private(k,i,j) &
           !$omp private(fact) &
           !$omp shared(PBLH,POTV,REAL_CZ,REAL_FZ) &
           !$omp shared(KS,KE,IA,JA)
           do j = 1, ja
           do i = 1, ia
              pblh(i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j)
              do k = ks+1, ke
                 if ( potv(k,i,j) > potv(ks,i,j) ) then
                    fact = ( potv(ks,i,j) - potv(k-1,i,j) ) &
                         / ( potv(k,i,j)  - potv(k-1,i,j) )
                    pblh(i,j) = real_cz(k-1,i,j) - real_fz(ks-1,i,j) &
                              + fact * ( real_cz(k,i,j) - real_cz(k-1,i,j) )
 
                    exit
                 endif
              enddo
           enddo
           enddo
           dv_calculated(i_pblh) = .true.
        end if
        var(:,:) = pblh(:,:)
 
     case ( 'CAPE', 'CIN', 'LCL', 'LFC', 'LNB' )
        if ( .not. dv_calculated(i_cape) ) then
           call allocate_2d( cape )
           call allocate_2d( cin )
           call allocate_2d( lcl )
           call allocate_2d( lfc )
           call allocate_2d( lnb )
           call atmos_adiabat_cape( &
                ka, ks, ke, ia, is, ie, ja, js, je, &
                ks,                                               & ! (in)
                dens_av(:,:,:), temp(:,:,:), pres(:,:,:),         & ! (in)
                qv(:,:,:), qc(:,:,:), qdry(:,:,:),                & ! (in)
                rtot(:,:,:), cptot(:,:,:),                        & ! (in)
                real_cz(:,:,:), real_fz(:,:,:),                   & ! (in)
                cape(:,:), cin(:,:), lcl(:,:), lfc(:,:), lnb(:,:) ) ! (out)
           dv_calculated(i_cape) = .true.
        end if
        select case ( vname )
        case ( 'CAPE' )
           var(:,:) = cape(:,:)
        case ( 'CIN' )
           var(:,:) = cin(:,:)
        case ( 'LCL' )
           var(:,:) = lcl(:,:)
        case ( 'LFC' )
           var(:,:) = lfc(:,:)
        case ( 'LNB' )
           var(:,:) = lnb(:,:)
        end select
 
     case ( 'PREC', 'RAIN', 'SNOW' )
        if ( .not. dv_calculated(i_prec) ) then
           call allocate_2d( prec )
           call allocate_2d( rain )
           call allocate_2d( snow )
           !$omp parallel do private(i,j) OMP_SCHEDULE_
           do j = 1, ja
           do i = 1, ia
              rain(i,j) = sflx_rain_mp(i,j) + sflx_rain_cp(i,j)
              snow(i,j) = sflx_snow_mp(i,j)
              prec(i,j) = rain(i,j) + snow(i,j)
           enddo
           enddo
           dv_calculated(i_prec) = .true.
        end if
        select case (vname)
        case ( 'RAIN' )
           var(:,:) = rain(:,:)
        case ( 'SNOW' )
           var(:,:) = snow(:,:)
        case ( 'PREC' )
           var(:,:) = prec(:,:)
        end select
 
     case default
        log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_2D: ', trim(vname)
        call prc_abort
     end select
 
 
     return

Here is the call graph for this function:

◆ atmos_vars_get_diagnostic_1d()

recursive subroutine mod_atmos_vars::atmos_vars_get_diagnostic_1d	(	character(len=*), intent(in)	vname,
		real(rp), dimension(:), intent(out)	var
	)

get diagnostic variable 1D

Definition at line 2622 of file mod_atmos_vars.F90.

     use scale_const, only: &
        cpdry => const_cpdry
     use scale_prc, only: &
        prc_abort
     use scale_statistics, only: &
        statistics_horizontal_mean
     use scale_atmos_grid_cartesc_real, only: &
        area => atmos_grid_cartesc_real_area
     implicit none
 
     character(len=*), intent(in)  :: vname
     real(RP),         intent(out) :: var(:)
 
     real(RP) :: work(ka,ia,ja)
     integer  :: k, i, j
     !---------------------------------------------------------------------------
 
     select case ( vname )
     case ( 'DENS_MEAN' )
        if ( .not. dv_calculated(i_dens_mean) ) then
           call allocate_1d( dens_mean )
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            dens(:,:,:), area(:,:), dens_mean(:) )
           dv_calculated(i_dens_mean) = .true.
        end if
        var(:) = dens_mean(:)
 
     case ( 'W_MEAN' )
        if ( .not. dv_calculated(i_w_mean) ) then
           call allocate_1d( w_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = w(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), w_mean(:) )
           do k = ks, ke
              w_mean(k) = w_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_w_mean) = .true.
        end if
        var(:) = w_mean(:)
 
     case ( 'U_MEAN' )
        if ( .not. dv_calculated(i_u_mean) ) then
           call allocate_1d( u_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = u(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), u_mean(:) )
           do k = ks, ke
              u_mean(k) = u_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_u_mean) = .true.
        end if
        var(:) = u_mean(:)
 
     case ( 'V_MEAN' )
        if ( .not. dv_calculated(i_v_mean) ) then
           call allocate_1d( v_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = v(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), v_mean(:) )
           do k = ks, ke
              v_mean(k) = v_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_v_mean) = .true.
        end if
        var(:) = v_mean(:)
 
     case ( 'PT_MEAN' )
        if ( .not. dv_calculated(i_pt_mean) ) then
           call allocate_1d( pt_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            rhot(:,:,:), area(:,:), pt_mean(:) )
           do k = ks, ke
              pt_mean(k) = pt_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_pt_mean) = .true.
        end if
        var(:) = pt_mean(:)
 
     case ( 'T_MEAN' )
        if ( .not. dv_calculated(i_t_mean) ) then
           call allocate_1d( t_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = temp(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), t_mean(:) )
           do k = ks, ke
              t_mean(k) = t_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_t_mean) = .true.
        end if
        var(:) = t_mean(:)
 
     case ( 'QV_MEAN' )
        if ( .not. dv_calculated(i_qv_mean) ) then
           call allocate_1d( qv_mean )
           if ( moist ) then
              call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
 !OCL XFILL
              !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
              do j = jsb, jeb
              do i = isb, ieb
              do k = ks, ke
                 work(k,i,j) = qv(k,i,j) * dens_av(k,i,j)
              enddo
              enddo
              enddo
              call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                               work(:,:,:), area(:,:), qv_mean(:) )
              do k = ks, ke
                 qv_mean(k) = qv_mean(k) / dens_mean(k)
              enddo
           else
              !$omp parallel do private(k) OMP_SCHEDULE_
              do k = ks, ke
                 qv_mean(k) = 0.0_rp
              enddo
           end if
           dv_calculated(i_qv_mean) = .true.
        end if
        var(:) = qv_mean(:)
 
     case ( 'QHYD_MEAN' )
        if ( .not. dv_calculated(i_qhyd_mean) ) then
           call allocate_1d( qhyd_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
           call atmos_vars_get_diagnostic( 'QHYD', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = qhyd(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), qhyd_mean(:) )
           do k = ks, ke
              qhyd_mean(k) = qhyd_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_qhyd_mean) = .true.
        end if
        var(:) = qhyd_mean(:)
 
     case ( 'QLIQ_MEAN' )
        if ( .not. dv_calculated(i_qliq_mean) ) then
           call allocate_1d( qliq_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
           call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = qliq(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), qliq_mean(:) )
           do k = ks, ke
              qliq_mean(k) = qliq_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_qliq_mean) = .true.
        end if
        var(:) = qliq_mean(:)
 
     case ( 'QICE_MEAN' )
        if ( .not. dv_calculated(i_qice_mean) ) then
           call allocate_1d( qice_mean )
           call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
           call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
 !OCL XFILL
           !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
           do j = jsb, jeb
           do i = isb, ieb
           do k = ks, ke
              work(k,i,j) = qice(k,i,j) * dens_av(k,i,j)
           enddo
           enddo
           enddo
           call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
                                            work(:,:,:), area(:,:), qice_mean(:) )
           do k = ks, ke
              qice_mean(k) = qice_mean(k) / dens_mean(k)
           enddo
           dv_calculated(i_qice_mean) = .true.
        end if
        var(:) = qice_mean(:)
 
     case default
        log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_1D: ', trim(vname)
        call prc_abort
     end select
 
 
     return

Here is the call graph for this function:

◆ atmos_vars_monitor()

subroutine, public mod_atmos_vars::atmos_vars_monitor ( )

monitor output

Definition at line 2860 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_update(), and mod_rm_driver::restart_read().

     use scale_prc, only: &
        prc_myrank, &
        prc_abort
     use scale_const, only: &
        grav  => const_grav,  &
        cvdry => const_cvdry
     use scale_atmos_grid_cartesc, only: &
        rfdx => atmos_grid_cartesc_rfdx, &
        rfdy => atmos_grid_cartesc_rfdy
     use scale_atmos_grid_cartesc_real, only: &
        real_cz => atmos_grid_cartesc_real_cz
     use scale_atmos_grid_cartesc_metric, only: &
        mapf => atmos_grid_cartesc_metric_mapf
     use scale_statistics, only: &
        statistics_checktotal, &
        statistics_total,            &
        statistics_detail
     use scale_monitor, only: &
        monitor_put
     use scale_time, only: &
        time_dtsec_atmos_dyn
     use mod_atmos_admin, only: &
        atmos_dyn_type
     use scale_atmos_hydrometeor, only: &
        i_qv
     use mod_atmos_phy_cp_vars, only: &
        sflx_rain_cp => atmos_phy_cp_sflx_rain
     use mod_atmos_phy_mp_vars, only: &
        sflx_rain_mp => atmos_phy_mp_sflx_rain, &
        sflx_snow_mp => atmos_phy_mp_sflx_snow
     use mod_atmos_phy_rd_vars, only: &
        sflx_lw_up   => atmos_phy_rd_sflx_lw_up,   &
        sflx_lw_dn   => atmos_phy_rd_sflx_lw_dn,   &
        sflx_sw_up   => atmos_phy_rd_sflx_sw_up,   &
        sflx_sw_dn   => atmos_phy_rd_sflx_sw_dn,   &
        toaflx_lw_up => atmos_phy_rd_toaflx_lw_up, &
        toaflx_lw_dn => atmos_phy_rd_toaflx_lw_dn, &
        toaflx_sw_up => atmos_phy_rd_toaflx_sw_up, &
        toaflx_sw_dn => atmos_phy_rd_toaflx_sw_dn
     use mod_atmos_phy_sf_vars, only: &
        sflx_sh   => atmos_phy_sf_sflx_sh, &
        sflx_lh   => atmos_phy_sf_sflx_lh, &
        sflx_qtrc => atmos_phy_sf_sflx_qtrc
     implicit none
 
     real(RP) :: rhoq(ka,ia,ja)
 
     real(RP) :: engflxt    (ia,ja) ! total flux             [J/m2/s]
     real(RP) :: sflx_rd_net(ia,ja) ! net SFC radiation flux [J/m2/s]
     real(RP) :: tflx_rd_net(ia,ja) ! net TOA radiation flux [J/m2/s]
 
     real(RP)               :: work (ka,ia,ja,3)
     character(len=H_SHORT) :: wname(3)
     real(RP)               :: cflmax
 
     integer  :: k, i, j, iq
     !---------------------------------------------------------------------------
 
     call monitor_put( pv_monit_id(i_dens), dens(:,:,:) )
     call monitor_put( pv_monit_id(i_momz), momz(:,:,:) )
     call monitor_put( pv_monit_id(i_momx), momx(:,:,:) )
     call monitor_put( pv_monit_id(i_momy), momy(:,:,:) )
     call monitor_put( pv_monit_id(i_rhot), rhot(:,:,:) )
 
     !##### Mass Budget #####
 
     do iq = 1, qa
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j = js, je
        do i = is, ie
        do k = ks, ke
           rhoq(k,i,j) = dens_av(k,i,j) * qtrc_av(k,i,j,iq)
        enddo
        enddo
        enddo
 
        call monitor_put( qp_monit_id(iq), rhoq(:,:,:) )
     enddo
 
     ! total dry airmass
     if ( dv_monit_id(im_qdry) > 0 ) then
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j = js, je
        do i = is, ie
        do k = ks, ke
           rhoq(k,i,j) = dens(k,i,j) * qdry(k,i,j)
        enddo
        enddo
        enddo
        call monitor_put( dv_monit_id(im_qdry), rhoq(:,:,:) )
     end if
 
     ! total vapor,liquid,solid tracers
     if ( dv_monit_id(im_qtot) > 0 ) then
        call atmos_vars_get_diagnostic( 'QTOT', work3d(:,:,:) )
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j = js, je
        do i = is, ie
        do k = ks, ke
           rhoq(k,i,j) = dens(k,i,j) * qtot(k,i,j)
        enddo
        enddo
        enddo
        call monitor_put( dv_monit_id(im_qtot), rhoq(:,:,:) )
     end if
 
     ! total evapolation
     if ( moist ) call monitor_put( dv_monit_id(im_evap), sflx_qtrc(:,:,i_qv) )
 
     ! total precipitation
     if ( dv_monit_id(im_prec) > 0 ) then
        call atmos_vars_get_diagnostic( 'PREC', work2d(:,:) )
        call monitor_put( dv_monit_id(im_prec), work2d(:,:) )
     end if
 
 
     !##### Energy Budget #####
 
     if ( dv_monit_id(im_engt) > 0 ) then
        call atmos_vars_get_diagnostic( 'ENGT', work3d(:,:,:) )
        call monitor_put( dv_monit_id(im_engt), work3d(:,:,:) )
     end if
     if ( dv_monit_id(im_engp) > 0 ) then
        call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
        call monitor_put( dv_monit_id(im_engp), work3d(:,:,:) )
     end if
     if ( dv_monit_id(im_engk) > 0 ) then
        call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
        call monitor_put( dv_monit_id(im_engk), work3d(:,:,:) )
     end if
     if ( dv_monit_id(im_engi) > 0 ) then
        call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
        call monitor_put( dv_monit_id(im_engi), work3d(:,:,:) )
     end if
 
 
     ! radiation flux
 !OCL XFILL
     !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
     do j = js, je
     do i = is, ie
        sflx_rd_net(i,j) = ( sflx_lw_up(i,j) - sflx_lw_dn(i,j) ) &
                         + ( sflx_sw_up(i,j) - sflx_sw_dn(i,j) )
 
        tflx_rd_net(i,j) = ( toaflx_lw_up(i,j) - toaflx_lw_dn(i,j) ) &
                         + ( toaflx_sw_up(i,j) - toaflx_sw_dn(i,j) )
 
        engflxt(i,j) = sflx_sh(i,j) + sflx_lh(i,j) &
                         + sflx_rd_net(i,j) - tflx_rd_net(i,j)
     enddo
     enddo
 
     call monitor_put( dv_monit_id(im_engflxt),      engflxt(:,:) )
 
     call monitor_put( dv_monit_id(im_engsfc_sh),    sflx_sh(:,:) )
     call monitor_put( dv_monit_id(im_engsfc_lh),    sflx_lh(:,:) )
     call monitor_put( dv_monit_id(im_engsfc_rd),    sflx_rd_net(:,:) )
     call monitor_put( dv_monit_id(im_engtoa_rd),    tflx_rd_net(:,:) )
 
     call monitor_put( dv_monit_id(im_engsfc_lw_up), sflx_lw_up(:,:) )
     call monitor_put( dv_monit_id(im_engsfc_lw_dn), sflx_lw_dn(:,:) )
     call monitor_put( dv_monit_id(im_engsfc_sw_up), sflx_sw_up(:,:) )
     call monitor_put( dv_monit_id(im_engsfc_sw_dn), sflx_sw_dn(:,:) )
 
     call monitor_put( dv_monit_id(im_engtoa_lw_up), toaflx_lw_up(:,:) )
     call monitor_put( dv_monit_id(im_engtoa_lw_dn), toaflx_lw_dn(:,:) )
     call monitor_put( dv_monit_id(im_engtoa_sw_up), toaflx_sw_up(:,:) )
     call monitor_put( dv_monit_id(im_engtoa_sw_dn), toaflx_sw_dn(:,:) )
 
 
 
     if ( atmos_vars_checkrange ) then
 !OCL XFILL
        work(:,:,:,1) = w(:,:,:)
 !OCL XFILL
        work(:,:,:,2) = u(:,:,:)
 !OCL XFILL
        work(:,:,:,3) = v(:,:,:)
 
        wname(1) = "W"
        wname(2) = "U"
        wname(3) = "V"
 
        call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
                                wname(:), work(:,:,:,:)                )
     endif
 
     if (       ( atmos_dyn_type /= 'OFF' .AND. atmos_dyn_type /= 'NONE' )                   &
          .AND. ( atmos_vars_checkcfl_soft > 0.0_rp .OR. atmos_vars_checkcfl_hard > 0.0_rp ) ) then
 !OCL XFILL
        work(:,:,:,:) = 0.0_rp
 
        do j = js, je
        do i = is, ie
        do k = ks, ke
           work(k,i,j,1) = 0.5_rp * abs(momz_av(k,i,j)) / ( dens_av(k+1,i,j) + dens_av(k,i,j) ) &
                         * time_dtsec_atmos_dyn / ( real_cz(k+1,i,j) - real_cz(k,i,j) )
           work(k,i,j,2) = 0.5_rp * abs(momx_av(k,i,j)) / ( dens_av(k,i+1,j) + dens_av(k,i,j) ) &
                         * time_dtsec_atmos_dyn * rfdx(i) * mapf(i,j,1,i_uy)
           work(k,i,j,3) = 0.5_rp * abs(momy_av(k,i,j)) / ( dens_av(k,i,j+1) + dens_av(k,i,j) ) &
                         * time_dtsec_atmos_dyn * rfdy(j) * mapf(i,j,2,i_xv)
        enddo
        enddo
        enddo
 
        cflmax = maxval( work(:,:,:,:) )
 
        if ( atmos_vars_checkcfl_hard > 0.0_rp .AND. cflmax > atmos_vars_checkcfl_hard ) then
           log_info("ATMOS_vars_monitor",*) "Courant number =", cflmax, " exceeded the hard limit =", atmos_vars_checkcfl_hard
           log_error("ATMOS_vars_monitor",*)                     "Courant number =", cflmax, " exceeded the hard limit =", atmos_vars_checkcfl_hard
           log_error_cont(*)                     "Rank =", prc_myrank
           log_error_cont(*)                     "Please set ATMOS_VARS_CHECKCFL_HARD in the namelist PARAM_ATMOS_VARS when you want to change the limit."
 
           wname(1) = "Courant num. Z"
           wname(2) = "Courant num. X"
           wname(3) = "Courant num. Y"
           call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
                                   wname(:), work(:,:,:,:),               &
                                   local=.true.                           )
 
           call prc_abort
        endif
 
        if ( atmos_vars_checkcfl_soft > 0.0_rp .AND. cflmax > atmos_vars_checkcfl_soft ) then
           log_info("ATMOS_vars_monitor",*) "Courant number =", cflmax, " exceeded the soft limit =", atmos_vars_checkcfl_soft
           log_error("ATMOS_vars_monitor",*)                     "Courant number =", cflmax, " exceeded the soft limit =", atmos_vars_checkcfl_soft
           log_error_cont(*)                     "Rank =", prc_myrank
 
           wname(1) = "Courant num. Z"
           wname(2) = "Courant num. X"
           wname(3) = "Courant num. Y"
           call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
                                   wname(:), work(:,:,:,:),               &
                                   local=.true.                           )
        endif
     endif
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_create()

subroutine, public mod_atmos_vars::atmos_vars_restart_create ( )

Create atmospheric restart file.

Definition at line 3106 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_time, only: &
        time_gettimelabel
     use scale_file_cartesc, only: &
        file_cartesc_create
     use mod_atmos_admin, only: &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_cpl_admin, only: &
        cpl_sw
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_create
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_create
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_create
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_create
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_create
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_create
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_create
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_create
 #ifdef SDM
     use scale_atmos_phy_mp_sdm, only: &
        sd_rest_flg_out, &
        atmos_phy_mp_sdm_restart_create
     use scale_time, only: &
        nowdaysec => time_nowdaysec
 #endif
     implicit none
 
     character(len=19)     :: timelabel
     character(len=H_LONG) :: basename
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        log_info("ATMOS_vars_restart_create",*) 'Output random number for SDM '
        call atmos_phy_mp_sdm_restart_create(nowdaysec)
     endif
 #endif
 
     if ( atmos_restart_out_basename /= '' ) then
 
        log_newline
        log_info("ATMOS_vars_restart_create",*) 'Create restart file (ATMOS) '
 
        if ( atmos_restart_out_postfix_timelabel ) then
           call time_gettimelabel( timelabel )
           basename = trim(atmos_restart_out_basename)//'_'//trim(timelabel)
        else
           basename = trim(atmos_restart_out_basename)
        endif
 
        log_info("ATMOS_vars_restart_create",*) 'basename: ', trim(basename)
 
        call file_cartesc_create( &
             basename, atmos_restart_out_title, atmos_restart_out_dtype, & ! [IN]
             restart_fid,                                                & ! [OUT]
             aggregate=atmos_restart_out_aggregate                       ) ! [IN]
 
        allocate( pv_id(pv_nmax+qa) )
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_create
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_create
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_create
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_create
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_create
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_create
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_create
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_create
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_enddef()

subroutine, public mod_atmos_vars::atmos_vars_restart_enddef ( )

Exit netCDF define mode.

Definition at line 3194 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_file_cartesc, only: &
        file_cartesc_enddef
     use mod_atmos_admin, only: &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_cpl_admin, only: &
        cpl_sw
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_enddef
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_enddef
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_enddef
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_enddef
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_enddef
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_enddef
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_enddef
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_enddef
 #ifdef SDM
     use scale_atmos_phy_mp_sdm, only: &
        sd_rest_flg_out, &
        atmos_phy_mp_sdm_restart_enddef
 #endif
     implicit none
 
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        call atmos_phy_mp_sdm_restart_enddef
     endif
 #endif
 
     if ( restart_fid /= -1 ) then
        call file_cartesc_enddef( restart_fid ) ! [IN]
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_enddef
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_enddef
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_enddef
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_enddef
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_enddef
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_enddef
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_enddef
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_enddef
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_close()

subroutine, public mod_atmos_vars::atmos_vars_restart_close ( )

Close restart file.

Definition at line 3257 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_read(), and mod_admin_restart::admin_restart_write().

     use scale_file_cartesc, only: &
        file_cartesc_close
     use mod_atmos_admin, only: &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_cpl_admin, only: &
        cpl_sw
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_close
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_close
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_close
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_close
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_close
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_close
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_close
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_close
 #ifdef SDM
     use scale_atmos_phy_mp_sdm, only: &
        sd_rest_flg_out, &
        atmos_phy_mp_sdm_restart_close
 #endif
     implicit none
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        call atmos_phy_mp_sdm_restart_close
     endif
 #endif
 
     if ( restart_fid /= -1 ) then
        log_newline
        log_info("ATMOS_vars_restart_close",*) 'Close restart file (ATMOS) '
 
        call file_cartesc_close( restart_fid ) ! [IN]
 
        restart_fid = -1
 
        if ( allocated(pv_id) ) deallocate( pv_id )
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_close
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_close
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_close
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_close
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_close
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_close
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_close
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_close
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_def_var()

subroutine, public mod_atmos_vars::atmos_vars_restart_def_var ( )

Define atmospheric variables in restart file.

Definition at line 3326 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_file_cartesc, only: &
        file_cartesc_def_var
     use mod_atmos_admin, only: &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_cpl_admin, only: &
        cpl_sw
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_def_var
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_def_var
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_def_var
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_def_var
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_def_var
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_def_var
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_def_var
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_def_var
 #ifdef SDM
     use scale_atmos_phy_mp_sdm, only: &
        sd_rest_flg_out, &
        atmos_phy_mp_sdm_restart_def_var
 #endif
     implicit none
 
     integer iq
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        call atmos_phy_mp_sdm_restart_def_var
     endif
 #endif
 
     if ( restart_fid /= -1 ) then
 
        call file_cartesc_def_var( restart_fid, pv_info(i_dens)%NAME, pv_info(i_dens)%DESC, pv_info(i_dens)%UNIT, 'ZXY',  atmos_restart_out_dtype, &
                                   pv_id(i_dens), &
                                   standard_name=pv_info(i_dens)%STDNAME )
        call file_cartesc_def_var( restart_fid, pv_info(i_momz)%NAME, pv_info(i_momz)%DESC, pv_info(i_momz)%UNIT, 'ZHXY', atmos_restart_out_dtype, &
                                   pv_id(i_momz), &
                                   standard_name=pv_info(i_momz)%STDNAME )
        call file_cartesc_def_var( restart_fid, pv_info(i_momx)%NAME, pv_info(i_momx)%DESC, pv_info(i_momx)%UNIT, 'ZXHY', atmos_restart_out_dtype, &
                                   pv_id(i_momx), &
                                   standard_name=pv_info(i_momx)%STDNAME )
        call file_cartesc_def_var( restart_fid, pv_info(i_momy)%NAME, pv_info(i_momy)%DESC, pv_info(i_momy)%UNIT, 'ZXYH', atmos_restart_out_dtype, &
                                   pv_id(i_momy), &
                                   standard_name=pv_info(i_momy)%STDNAME )
        call file_cartesc_def_var( restart_fid, pv_info(i_rhot)%NAME, pv_info(i_rhot)%DESC, pv_info(i_rhot)%UNIT, 'ZXY',  atmos_restart_out_dtype, &
                                   pv_id(i_rhot), &
                                   standard_name=pv_info(i_rhot)%STDNAME )
        do iq = 1, qa
           call file_cartesc_def_var( restart_fid, tracer_name(iq), tracer_desc(iq), tracer_unit(iq), 'ZXY',  atmos_restart_out_dtype, &
                                      pv_id(pv_nmax+iq) )
        enddo
 
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_def_var
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_def_var
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_def_var
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_def_var
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_def_var
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_def_var
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_def_var
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_def_var
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_restart_write()

subroutine, public mod_atmos_vars::atmos_vars_restart_write ( )

Write restart of atmospheric variables.

Definition at line 3410 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_file_cartesc, only: &
        file_cartesc_write_var
     use mod_atmos_admin, only: &
        atmos_sw_dyn,      &
        atmos_sw_phy_mp,   &
        atmos_sw_phy_ae,   &
        atmos_sw_phy_ch,   &
        atmos_sw_phy_rd,   &
        atmos_sw_phy_sf,   &
        atmos_sw_phy_tb,   &
        atmos_sw_phy_cp
     use mod_cpl_admin, only: &
        cpl_sw
     use mod_atmos_dyn_vars, only: &
        atmos_dyn_vars_restart_write
     use mod_atmos_phy_mp_vars, only: &
        atmos_phy_mp_vars_restart_write
     use mod_atmos_phy_ae_vars, only: &
        atmos_phy_ae_vars_restart_write
     use mod_atmos_phy_ch_vars, only: &
        atmos_phy_ch_vars_restart_write
     use mod_atmos_phy_rd_vars, only: &
        atmos_phy_rd_vars_restart_write
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_vars_restart_write
     use mod_atmos_phy_tb_vars, only: &
        atmos_phy_tb_vars_restart_write
     use mod_atmos_phy_cp_vars, only: &
        atmos_phy_cp_vars_restart_write
 #ifdef SDM
     use scale_atmos_phy_mp_sdm, only: &
        sd_rest_flg_out, &
        atmos_phy_mp_sdm_restart_write
 #endif
     implicit none
 
     integer iq
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        call atmos_phy_mp_sdm_restart_write
     endif
 #endif
 
     if ( restart_fid /= -1 ) then
 
        call atmos_vars_fillhalo
 
        call atmos_vars_total
 
        call file_cartesc_write_var( restart_fid, pv_id(i_dens), dens(:,:,:), pv_info(i_dens)%NAME, 'ZXY'  ) ! [IN]
        call file_cartesc_write_var( restart_fid, pv_id(i_momz), momz(:,:,:), pv_info(i_momz)%NAME, 'ZHXY' ) ! [IN]
        call file_cartesc_write_var( restart_fid, pv_id(i_momx), momx(:,:,:), pv_info(i_momx)%NAME, 'ZXHY' ) ! [IN]
        call file_cartesc_write_var( restart_fid, pv_id(i_momy), momy(:,:,:), pv_info(i_momy)%NAME, 'ZXYH' ) ! [IN]
        call file_cartesc_write_var( restart_fid, pv_id(i_rhot), rhot(:,:,:), pv_info(i_rhot)%NAME, 'ZXY'  ) ! [IN]
 
        do iq = 1, qa
           call file_cartesc_write_var( restart_fid, pv_id(pv_nmax+iq), qtrc(:,:,:,iq), tracer_name(iq), 'ZXY' ) ! [IN]
        enddo
 
     endif
 
     if( atmos_sw_dyn )    call atmos_dyn_vars_restart_write
     if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_write
     if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_write
     if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_write
     if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_write
     if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_write
     if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_write
     if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_write
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ allocate_3d()

subroutine mod_atmos_vars::allocate_3d ( real(rp), dimension(:,:,:), intent(inout), allocatable ary )

Definition at line 3488 of file mod_atmos_vars.F90.

References scale_const::const_undef, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, and scale_atmos_grid_cartesc_index::ka.

Referenced by atmos_vars_get_diagnostic_3d().

     use scale_const, only: &
        undef => const_undef
     real(RP), intent(inout), allocatable :: ary(:,:,:)
 
     if ( .not. allocated(ary) ) then
        allocate( ary(ka,ia,ja) )
        ary(:,:,:) = undef
     end if
 
     return

Here is the caller graph for this function:

◆ allocate_2d()

subroutine mod_atmos_vars::allocate_2d ( real(rp), dimension(:,:), intent(inout), allocatable ary )

Definition at line 3501 of file mod_atmos_vars.F90.

References scale_const::const_undef, scale_atmos_grid_cartesc_index::ia, and scale_atmos_grid_cartesc_index::ja.

Referenced by atmos_vars_get_diagnostic_2d().

     use scale_const, only: &
        undef => const_undef
     real(RP), intent(inout), allocatable :: ary(:,:)
 
     if ( .not. allocated(ary) ) then
        allocate( ary(ia,ja) )
        ary(:,:) = undef
     end if
 
     return

Here is the caller graph for this function:

◆ allocate_1d()

subroutine mod_atmos_vars::allocate_1d ( real(rp), dimension(:), intent(inout), allocatable ary )

Definition at line 3514 of file mod_atmos_vars.F90.

References scale_const::const_undef, and scale_atmos_grid_cartesc_index::ka.

Referenced by atmos_vars_get_diagnostic_1d().

     use scale_const, only: &
        undef => const_undef
     real(RP), intent(inout), allocatable :: ary(:)
 
     if ( .not. allocated(ary) ) then
        allocate( ary(ka) )
        ary(:) = undef
     end if
 
     return

Here is the caller graph for this function:

Variable Documentation

◆ atmos_restart_output

logical, public mod_atmos_vars::atmos_restart_output = .false.

Output restart file?

Definition at line 59 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), mod_admin_restart::admin_restart_write(), and atmos_vars_setup().

59 logical, public :: atmos_restart_output = .false.

◆ atmos_restart_in_basename

character(len=h_long), public mod_atmos_vars::atmos_restart_in_basename = ''

Basename of the input file.

Definition at line 61 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), mod_admin_time::admin_time_setup(), atmos_vars_restart_open(), and atmos_vars_setup().

61 character(len=H_LONG), public :: atmos_restart_in_basename = ''

◆ atmos_restart_in_aggregate

logical, public mod_atmos_vars::atmos_restart_in_aggregate

Switch to use aggregate file.

Definition at line 62 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_open(), and atmos_vars_setup().

62 logical, public :: atmos_restart_in_aggregate

◆ atmos_restart_in_postfix_timelabel

logical, public mod_atmos_vars::atmos_restart_in_postfix_timelabel = .false.

Add timelabel to the basename of input file?

Definition at line 63 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_open(), and atmos_vars_setup().

63 logical, public :: atmos_restart_in_postfix_timelabel = .false.

◆ atmos_restart_out_basename

character(len=h_long), public mod_atmos_vars::atmos_restart_out_basename = ''

Basename of the output file.

Definition at line 64 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_create(), and atmos_vars_setup().

64 character(len=H_LONG), public :: atmos_restart_out_basename = ''

◆ atmos_restart_out_aggregate

logical, public mod_atmos_vars::atmos_restart_out_aggregate

Switch to use aggregate file.

Definition at line 65 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_create(), and atmos_vars_setup().

65 logical, public :: atmos_restart_out_aggregate

◆ atmos_restart_out_postfix_timelabel

logical, public mod_atmos_vars::atmos_restart_out_postfix_timelabel = .true.

Add timelabel to the basename of output file?

Definition at line 66 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_create(), and atmos_vars_setup().

66 logical, public :: atmos_restart_out_postfix_timelabel = .true.

◆ atmos_restart_out_title

character(len=h_mid), public mod_atmos_vars::atmos_restart_out_title = 'ATMOS restart'

Title of the output file.

Definition at line 67 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_create(), and atmos_vars_setup().

67 character(len=H_MID), public :: atmos_restart_out_title = 'ATMOS restart'

◆ atmos_restart_out_dtype

character(len=h_short), public mod_atmos_vars::atmos_restart_out_dtype = 'DEFAULT'

REAL4 or REAL8.

Definition at line 68 of file mod_atmos_vars.F90.

Referenced by mod_admin_restart::admin_restart_setup(), atmos_vars_restart_create(), atmos_vars_restart_def_var(), and atmos_vars_setup().

68 character(len=H_SHORT), public :: atmos_restart_out_dtype = 'DEFAULT'

◆ atmos_restart_check

logical, public mod_atmos_vars::atmos_restart_check = .false.

Check value consistency?

Definition at line 70 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup(), and mod_rm_driver::rm_driver().

70 logical, public :: atmos_restart_check = .false.

◆ atmos_restart_check_basename

character(len=h_long), public mod_atmos_vars::atmos_restart_check_basename = 'restart_check'

Definition at line 71 of file mod_atmos_vars.F90.

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

71 character(len=H_LONG), public :: atmos_restart_check_basename = 'restart_check'

◆ atmos_restart_check_criterion

real(rp), public mod_atmos_vars::atmos_restart_check_criterion = 1.E-6_RP

Definition at line 72 of file mod_atmos_vars.F90.

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

72 real(RP), public :: atmos_restart_check_criterion = 1.e-6_rp

◆ dens

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::dens

Definition at line 75 of file mod_atmos_vars.F90.

75 real(RP), public, target, allocatable :: dens(:,:,:) ! Density [kg/m3]

mod_atmos_vars::dens

real(rp), dimension(:,:,:), allocatable, target, public dens

Definition: mod_atmos_vars.F90:75

◆ momz

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momz

Definition at line 76 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

76 real(RP), public, target, allocatable :: momz(:,:,:) ! momentum z [kg/m2/s]

mod_atmos_vars::momz

real(rp), dimension(:,:,:), allocatable, target, public momz

Definition: mod_atmos_vars.F90:76

◆ momx

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momx

Definition at line 77 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

77 real(RP), public, target, allocatable :: momx(:,:,:) ! momentum x [kg/m2/s]

mod_atmos_vars::momx

real(rp), dimension(:,:,:), allocatable, target, public momx

Definition: mod_atmos_vars.F90:77

◆ momy

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momy

Definition at line 78 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

78 real(RP), public, target, allocatable :: momy(:,:,:) ! momentum y [kg/m2/s]

mod_atmos_vars::momy

real(rp), dimension(:,:,:), allocatable, target, public momy

Definition: mod_atmos_vars.F90:78

◆ rhot

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::rhot

Definition at line 79 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_1d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

79 real(RP), public, target, allocatable :: rhot(:,:,:) ! DENS * POTT [K*kg/m3]

mod_atmos_vars::rhot

real(rp), dimension(:,:,:), allocatable, target, public rhot

Definition: mod_atmos_vars.F90:79

◆ qtrc

real(rp), dimension(:,:,:,:), allocatable, target, public mod_atmos_vars::qtrc

Definition at line 80 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::mkinit(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), mod_mkinit::rect_setup(), and mod_mkinit::tke_setup().

80 real(RP), public, target, allocatable :: qtrc(:,:,:,:) ! ratio of mass of tracer to total mass[kg/kg]

mod_atmos_vars::qtrc

real(rp), dimension(:,:,:,:), allocatable, target, public qtrc

Definition: mod_atmos_vars.F90:80

◆ dens_avw

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::dens_avw

Definition at line 82 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

82 real(RP), public, target, allocatable :: dens_avw(:,:,:)

◆ momz_avw

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momz_avw

Definition at line 83 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

83 real(RP), public, target, allocatable :: momz_avw(:,:,:)

◆ momx_avw

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momx_avw

Definition at line 84 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

84 real(RP), public, target, allocatable :: momx_avw(:,:,:)

◆ momy_avw

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::momy_avw

Definition at line 85 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

85 real(RP), public, target, allocatable :: momy_avw(:,:,:)

◆ rhot_avw

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::rhot_avw

Definition at line 86 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

86 real(RP), public, target, allocatable :: rhot_avw(:,:,:)

◆ qtrc_avw

real(rp), dimension(:,:,:,:), allocatable, target, public mod_atmos_vars::qtrc_avw

Definition at line 87 of file mod_atmos_vars.F90.

Referenced by atmos_vars_setup().

87 real(RP), public, target, allocatable :: qtrc_avw(:,:,:,:)

◆ dens_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::dens_av

Definition at line 89 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_history_setpres(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

89 real(RP), public, pointer :: dens_av(:,:,:)

mod_atmos_vars::dens_av

real(rp), dimension(:,:,:), pointer, public dens_av

Definition: mod_atmos_vars.F90:89

◆ momz_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::momz_av

Definition at line 90 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

90 real(RP), public, pointer :: momz_av(:,:,:)

mod_atmos_vars::momz_av

real(rp), dimension(:,:,:), pointer, public momz_av

Definition: mod_atmos_vars.F90:90

◆ momx_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::momx_av

Definition at line 91 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

91 real(RP), public, pointer :: momx_av(:,:,:)

mod_atmos_vars::momx_av

real(rp), dimension(:,:,:), pointer, public momx_av

Definition: mod_atmos_vars.F90:91

◆ momy_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::momy_av

Definition at line 92 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

92 real(RP), public, pointer :: momy_av(:,:,:)

mod_atmos_vars::momy_av

real(rp), dimension(:,:,:), pointer, public momy_av

Definition: mod_atmos_vars.F90:92

◆ rhot_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::rhot_av

Definition at line 93 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

93 real(RP), public, pointer :: rhot_av(:,:,:)

mod_atmos_vars::rhot_av

real(rp), dimension(:,:,:), pointer, public rhot_av

Definition: mod_atmos_vars.F90:93

◆ qtrc_av

real(rp), dimension(:,:,:,:), pointer, public mod_atmos_vars::qtrc_av

Definition at line 94 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), and atmos_vars_setup().

94 real(RP), public, pointer :: qtrc_av(:,:,:,:)

mod_atmos_vars::qtrc_av

real(rp), dimension(:,:,:,:), pointer, public qtrc_av

Definition: mod_atmos_vars.F90:94

◆ qv

real(rp), dimension(:,:,:), allocatable, pointer, target, public mod_atmos_vars::qv

Definition at line 96 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_driver::atmos_driver_update(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_setup(), mod_atmos_phy_sf_driver::history_output(), and mod_rm_driver::restart_read().

96 real(RP), public, pointer :: qv(:,:,:)

◆ qc

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qc

Definition at line 97 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

97 real(RP), public, pointer :: qc(:,:,:)

◆ qr

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qr

Definition at line 98 of file mod_atmos_vars.F90.

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

98 real(RP), public, pointer :: qr(:,:,:)

◆ qi

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qi

Definition at line 99 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

99 real(RP), public, pointer :: qi(:,:,:)

◆ qs

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qs

Definition at line 100 of file mod_atmos_vars.F90.

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

100 real(RP), public, pointer :: qs(:,:,:)

◆ qg

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qg

Definition at line 101 of file mod_atmos_vars.F90.

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

101 real(RP), public, pointer :: qg(:,:,:)

◆ qh

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::qh

Definition at line 102 of file mod_atmos_vars.F90.

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

102 real(RP), public, pointer :: qh(:,:,:)

◆ qe

real(rp), dimension(:,:,:,:), allocatable, target, public mod_atmos_vars::qe

Definition at line 104 of file mod_atmos_vars.F90.

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), atmos_vars_calc_diagnostics(), and atmos_vars_setup().

104 real(RP), public, target, allocatable :: qe(:,:,:,:)

◆ mass

real(rp), allocatable, target, public mod_atmos_vars::mass

Definition at line 104 of file mod_atmos_vars.F90.

◆ ratio

real(rp), allocatable, target, public mod_atmos_vars::ratio

Definition at line 104 of file mod_atmos_vars.F90.

◆ of

real(rp), allocatable, target, public mod_atmos_vars::of

Definition at line 104 of file mod_atmos_vars.F90.

◆ hydrometors

real(rp), allocatable, target, public mod_atmos_vars::hydrometors

Definition at line 104 of file mod_atmos_vars.F90.

◆ kg

real(rp), allocatable, target, public mod_atmos_vars::kg

Definition at line 104 of file mod_atmos_vars.F90.

◆ dens_ref

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::dens_ref

Definition at line 107 of file mod_atmos_vars.F90.

107 real(RP), public, allocatable :: dens_ref(:,:,:)

◆ pott_ref

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::pott_ref

Definition at line 108 of file mod_atmos_vars.F90.

108 real(RP), public, allocatable :: pott_ref(:,:,:)

◆ temp_ref

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::temp_ref

Definition at line 109 of file mod_atmos_vars.F90.

109 real(RP), public, allocatable :: temp_ref(:,:,:)

◆ pres_ref

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::pres_ref

Definition at line 110 of file mod_atmos_vars.F90.

110 real(RP), public, allocatable :: pres_ref(:,:,:)

◆ qv_ref

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::qv_ref

Definition at line 111 of file mod_atmos_vars.F90.

111 real(RP), public, allocatable :: qv_ref(:,:,:)

◆ dens_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::dens_tp

Definition at line 114 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), and atmos_vars_setup().

114 real(RP), public, allocatable :: dens_tp(:,:,:)

◆ momz_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::momz_tp

Definition at line 115 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), and atmos_vars_setup().

115 real(RP), public, allocatable :: momz_tp(:,:,:)

mod_atmos_vars::momz_tp

real(rp), dimension(:,:,:), allocatable, public momz_tp

Definition: mod_atmos_vars.F90:115

◆ rhou_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::rhou_tp

Definition at line 116 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), and atmos_vars_setup().

116 real(RP), public, allocatable :: rhou_tp(:,:,:)

◆ rhov_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::rhov_tp

Definition at line 117 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), and atmos_vars_setup().

117 real(RP), public, allocatable :: rhov_tp(:,:,:)

◆ rhot_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::rhot_tp

Definition at line 118 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), and atmos_vars_setup().

118 real(RP), public, allocatable :: rhot_tp(:,:,:)

mod_atmos_vars::rhot_tp

real(rp), dimension(:,:,:), allocatable, public rhot_tp

Definition: mod_atmos_vars.F90:118

◆ rhoh_p

real(rp), dimension (:,:,:), allocatable, public mod_atmos_vars::rhoh_p

Definition at line 119 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), and atmos_vars_setup().

119 real(RP), public, allocatable :: rhoh_p (:,:,:)

◆ rhoq_tp

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_vars::rhoq_tp

Definition at line 120 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), and atmos_vars_setup().

120 real(RP), public, allocatable :: rhoq_tp(:,:,:,:)

mod_atmos_vars::rhoq_tp

real(rp), dimension(:,:,:,:), allocatable, public rhoq_tp

Definition: mod_atmos_vars.F90:120

◆ momx_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::momx_tp

Definition at line 123 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), and atmos_vars_setup().

123 real(RP), public, allocatable :: momx_tp(:,:,:)

mod_atmos_vars::momx_tp

real(rp), dimension(:,:,:), allocatable, public momx_tp

Definition: mod_atmos_vars.F90:123

◆ momy_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::momy_tp

Definition at line 124 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), and atmos_vars_setup().

124 real(RP), public, allocatable :: momy_tp(:,:,:)

mod_atmos_vars::momy_tp

real(rp), dimension(:,:,:), allocatable, public momy_tp

Definition: mod_atmos_vars.F90:124

◆ w

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::w

Definition at line 128 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), and atmos_vars_setup().

128 real(RP), public, allocatable, target :: w (:,:,:)

◆ velocity

real(rp), allocatable, target, public mod_atmos_vars::velocity

Definition at line 128 of file mod_atmos_vars.F90.

◆ m

real(rp), allocatable, target, public mod_atmos_vars::m

Definition at line 128 of file mod_atmos_vars.F90.

◆ s

real(rp), allocatable, target, public mod_atmos_vars::s

Definition at line 128 of file mod_atmos_vars.F90.

◆ u

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::u

Definition at line 129 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), and atmos_vars_setup().

129 real(RP), public, allocatable, target :: u (:,:,:)

◆ v

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::v

Definition at line 130 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), and atmos_vars_setup().

130 real(RP), public, allocatable, target :: v (:,:,:)

◆ pott

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::pott

Definition at line 132 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_update(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_setup(), and mod_rm_driver::restart_read().

132 real(RP), public, allocatable, target :: pott (:,:,:)

mod_atmos_vars::pott

real(rp), dimension(:,:,:), allocatable, target, public pott

Definition: mod_atmos_vars.F90:132

◆ potential

real(rp), allocatable, target, public mod_atmos_vars::potential

Definition at line 132 of file mod_atmos_vars.F90.

◆ temperature

real(rp), allocatable, target, public mod_atmos_vars::temperature

Definition at line 132 of file mod_atmos_vars.F90.

◆ k

real(rp), allocatable, target, public mod_atmos_vars::k

Definition at line 132 of file mod_atmos_vars.F90.

◆ temp

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::temp

Definition at line 133 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_update(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_setup(), mod_atmos_phy_sf_driver::history_output(), and mod_rm_driver::restart_read().

133 real(RP), public, allocatable, target :: temp (:,:,:)

◆ pres

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::pres

Definition at line 134 of file mod_atmos_vars.F90.

Referenced by mod_atmos_driver::atmos_driver_update(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_2d(), atmos_vars_history(), atmos_vars_history_setpres(), atmos_vars_setup(), mod_atmos_phy_sf_driver::history_output(), and mod_rm_driver::restart_read().

134 real(RP), public, allocatable, target :: pres (:,:,:)

◆ pressure

real(rp), allocatable, target, public mod_atmos_vars::pressure

Definition at line 134 of file mod_atmos_vars.F90.

◆ pa

real(rp), allocatable, target, public mod_atmos_vars::pa =J/m3]

Definition at line 134 of file mod_atmos_vars.F90.

◆ exner

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::exner

Definition at line 135 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), and atmos_vars_setup().

135 real(RP), public, allocatable, target :: exner(:,:,:)

◆ function

real(rp), dimension (t/pt), allocatable, target, public mod_atmos_vars::function

Definition at line 135 of file mod_atmos_vars.F90.

◆ phyd

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::phyd

Definition at line 136 of file mod_atmos_vars.F90.

Referenced by atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_history_setpres(), and atmos_vars_setup().

136 real(RP), public, allocatable, target :: phyd (:,:,:)

◆ hydrostatic

real(rp), allocatable, target, public mod_atmos_vars::hydrostatic

Definition at line 136 of file mod_atmos_vars.F90.

◆ phydh

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::phydh

Definition at line 137 of file mod_atmos_vars.F90.

Referenced by atmos_vars_calc_diagnostics(), atmos_vars_history_setpres(), and atmos_vars_setup().

137 real(RP), public, allocatable, target :: phydh(:,:,:)

◆ layer

real(rp), allocatable, target, public mod_atmos_vars::layer

Definition at line 137 of file mod_atmos_vars.F90.

◆ interface

real(rp), allocatable, target, public mod_atmos_vars::interface

Definition at line 137 of file mod_atmos_vars.F90.

◆ qdry

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::qdry

Definition at line 139 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_setup(), and atmos_vars_total().

139 real(RP), public, allocatable, target :: qdry (:,:,:)

◆ dry

real(rp), allocatable, target, public mod_atmos_vars::dry

Definition at line 139 of file mod_atmos_vars.F90.

◆ air

real(rp), allocatable, target, public mod_atmos_vars::air

Definition at line 139 of file mod_atmos_vars.F90.

◆ rtot

real(rp), dimension (:,:,:), allocatable, target, public mod_atmos_vars::rtot

Definition at line 140 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), and atmos_vars_setup().

140 real(RP), public, allocatable, target :: rtot (:,:,:)

◆ specific

real(rp), allocatable, target, public mod_atmos_vars::specific

Definition at line 140 of file mod_atmos_vars.F90.

◆ gass

real(rp), allocatable, target, public mod_atmos_vars::gass

Definition at line 140 of file mod_atmos_vars.F90.

◆ constant

real(rp), allocatable, target, public mod_atmos_vars::constant

Definition at line 140 of file mod_atmos_vars.F90.

◆ j

real(rp), allocatable, target, public mod_atmos_vars::j

Definition at line 140 of file mod_atmos_vars.F90.

◆ cvtot

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::cvtot

Definition at line 141 of file mod_atmos_vars.F90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), and atmos_vars_setup().

141 real(RP), public, allocatable, target :: cvtot(:,:,:)

◆ heat

real(rp), allocatable, target, public mod_atmos_vars::heat

Definition at line 141 of file mod_atmos_vars.F90.

◆ cptot

real(rp), dimension(:,:,:), allocatable, target, public mod_atmos_vars::cptot

Definition at line 142 of file mod_atmos_vars.F90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), atmos_vars_calc_diagnostics(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), and atmos_vars_setup().

142 real(RP), public, allocatable, target :: cptot(:,:,:)

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ atmos_vars_setup()

◆ atmos_vars_fillhalo()

◆ atmos_vars_restart_open()

◆ atmos_vars_restart_read()

◆ atmos_vars_history_setpres()

◆ atmos_vars_restart_check()

◆ atmos_vars_history()

◆ atmos_vars_total()

◆ atmos_vars_calc_diagnostics()

◆ atmos_vars_get_diagnostic_3d()

◆ atmos_vars_get_diagnostic_2d()

◆ atmos_vars_get_diagnostic_1d()

◆ atmos_vars_monitor()

◆ atmos_vars_restart_create()

◆ atmos_vars_restart_enddef()

◆ atmos_vars_restart_close()

◆ atmos_vars_restart_def_var()

◆ atmos_vars_restart_write()

◆ allocate_3d()

◆ allocate_2d()

◆ allocate_1d()

Variable Documentation

◆ atmos_restart_output

◆ atmos_restart_in_basename

◆ atmos_restart_in_aggregate

◆ atmos_restart_in_postfix_timelabel

◆ atmos_restart_out_basename

◆ atmos_restart_out_aggregate

◆ atmos_restart_out_postfix_timelabel

◆ atmos_restart_out_title

◆ atmos_restart_out_dtype

◆ atmos_restart_check

◆ atmos_restart_check_basename

◆ atmos_restart_check_criterion

◆ dens

◆ momz

◆ momx

◆ momy

◆ rhot

◆ qtrc

◆ dens_avw

◆ momz_avw

◆ momx_avw

◆ momy_avw

◆ rhot_avw

◆ qtrc_avw

◆ dens_av

◆ momz_av

◆ momx_av

◆ momy_av

◆ rhot_av

◆ qtrc_av

◆ qv

◆ qc

◆ qr

◆ qi

◆ qs

◆ qg

◆ qh

◆ qe

◆ mass

◆ ratio

◆ of

◆ hydrometors

◆ kg

◆ dens_ref

◆ pott_ref

◆ temp_ref

◆ pres_ref

◆ qv_ref

◆ dens_tp

◆ momz_tp

◆ rhou_tp

◆ rhov_tp

◆ rhot_tp

◆ rhoh_p