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_diagnostics
	Calc diagnostic variables. 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...

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_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_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, public	dens_tp

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

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

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

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

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

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

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

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

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

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

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

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

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

Detailed Description

module ATMOSPHERIC Variables

Description: Container for atmospheric variables

Author: Team SCALE

History

2011-11-11 (H.Yashiro) [new]
2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
2012-06-13 (S.Nishizawa) [mod] follows the change of mod_hist

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_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_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	real(RP)	0.0_RP

History Output

name	description	unit	variable
CAPE	convection avail. pot. energy	m2/s2	CAPE
CIN	convection inhibition	m2/s2	CIN
CPTOT	Total heat capacity	J/kg/K	CPTOT
DENS_MEAN	horiz. mean of density	kg/m3	DENS_MEAN
DENS_PRIM	horiz. deviation of density	kg/m3	DENS_PRIM
DIV	divergence	1/s	DIV
ENGI	internal energy	J/m3	ENGI
ENGK	kinetic energy	J/m3	ENGK
ENGP	potential energy	J/m3	ENGP
ENGT	total energy	J/m3	ENGT
HDIV	horizontal divergence	1/s	HDIV
IWP	ice water path	g/m2	IWP
LCL	lifted condensation level	m	LCL
LFC	level of free convection	m	LFC
LNB	level of neutral buoyancy	m	LNB
LWP	liquid water path	g/m2	LWP
LWPT	liq. potential temp.	K	POTL
MSE	moist static energy	m2/s2	MSE
N2	squared Brunt-Vaisala frequency	1/s2	N2
PBLH	PBL height	m	PBLH
PREC	surface precipitation rate (total)	kg/m2/s	PREC
PRES	pressure	Pa	PRES
PT	potential temp.	K	POTT
PT_MEAN	horiz. mean of pot.	K	PT_MEAN
PT_PRIM	horiz. deviation of pot. temp.	K	POTT_PRIM
PT_W_PRIM	resolved scale heat flux	W/s	PT_W_PRIM
PW	precipitable water	g/m2	PW
QDRY	dry air	kg/kg	QDRY
QHYD	total hydrometeors	kg/kg	QHYD
QHYD_MEAN	horiz. mean of QHYD	1	QHYD_MEAN
QICE	total ice water	kg/kg	QICE
QICE_MEAN	horiz. mean of QICE	1	QICE_MEAN
QLIQ	total liquid water	kg/kg	QLIQ
QLIQ_MEAN	horiz. mean of QLIQ	1	QLIQ_MEAN
QTOT	total water	kg/kg	QTOT
QV_MEAN	horiz. mean of QV	1	QV_MEAN
RAIN	surface rain rate (total)	kg/m2/s	RAIN
RH	relative humidity(liq)	%	RHL
RHA	relative humidity(liq+ice)	%	RHA
RHI	relative humidity(ice)	%	RHI
RTOT	Total gas constant	J/kg/K	RTOT
SNOW	surface snow rate (total)	kg/m2/s	SNOW
T	temperature	K	TEMP
TKE_RS	resolved scale TKE	m2/s2	TKE_RS
TRACER_NAME(iq)	TRACER_DESC(iq)	TRACER_UNIT(iq)	QTRC
T_MEAN	horiz. mean of t	K	T_MEAN
U	velocity u	m/s	U
U_MEAN	horiz. mean of u	m/s	U_MEAN
U_PRIM	horiz. deviation of u	m/s	U_PRIM
Uabs	absolute velocity	m/s	Uabs
V	velocity v	m/s	V
DENS	density	kg/m3	DENS
MOMX	momentum x	kg/m2/s	MOMX
MOMY	momentum y	kg/m2/s	MOMY
MOMZ	momentum z	kg/m2/s	MOMZ
RHOT	rho * theta	kg/m3*K	RHOT
VOR	vertical vorticity	1/s	VOR
V_MEAN	horiz. mean of v	m/s	V_MEAN
V_PRIM	horiz. deviation of v	m/s	V_PRIM
W	velocity w	m/s	W
W_MEAN	horiz. mean of w	m/s	W_MEAN
W_PRIM	horiz. deviation of w	m/s	W_PRIM
W_PRIM2	variance of w	m2/s2	W_PRIM2
W_PRIM3	skewness of w	m3/s3	W_PRIM3

Function/Subroutine Documentation

◆ atmos_vars_setup()

subroutine, public mod_atmos_vars::atmos_vars_setup ( )

Setup.

Definition at line 254 of file mod_atmos_vars.f90.

References mod_atmos_dyn_vars::atmos_dyn_vars_setup(), mod_atmos_phy_ae_vars::atmos_phy_ae_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_basename, atmos_restart_in_postfix_timelabel, 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, dens, dens_av, dens_avw, dens_tp, scale_history::hist_reg(), scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, scale_grid_index::ia, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, scale_grid_index::ja, scale_grid_index::ka, scale_grid_index::ke, scale_grid_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::monit_reg(), n2, phyd, pott, scale_process::prc_mpistop(), pres, scale_tracer::qa, qtrc, qtrc_av, qtrc_avw, rhoq_tp, rhot, rhot_av, rhot_avw, rhot_tp, temp, scale_tracer::tracer_desc, scale_tracer::tracer_name, scale_tracer::tracer_unit, u, v, and w.

Referenced by mod_rm_driver::scalerm(), and mod_rm_prep::scalerm_prep().

     use scale_process, only: &
        prc_mpistop
     use scale_history, only: &
        hist_reg
     use scale_monitor, only: &
        monit_reg
     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_cp_vars, only: &
        atmos_phy_cp_vars_setup
     implicit none
 
     namelist / param_atmos_vars / &
        atmos_restart_in_basename,           &
        atmos_restart_in_postfix_timelabel,  &
        atmos_restart_in_check_coordinates,  &
        atmos_restart_output,                &
        atmos_restart_out_basename,          &
        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
 
     integer :: ierr
     integer :: iv, iq
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[VARS] / Categ[ATMOS] / Origin[SCALE-RM]'
 
     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( momx_tp(ka,ia,ja)    )
     allocate( momy_tp(ka,ia,ja)    )
     allocate( rhot_tp(ka,ia,ja)    )
     allocate( rhoq_tp(ka,ia,ja,max(qa,1)) )
 
     allocate( pott(ka,ia,ja) )
     allocate( temp(ka,ia,ja) )
     allocate( pres(ka,ia,ja) )
     allocate( phyd(ka,ia,ja) )
     allocate( w(ka,ia,ja) )
     allocate( u(ka,ia,ja) )
     allocate( v(ka,ia,ja) )
     allocate( n2(ka,ia,ja) )
 
     momz(1:ks-1,:,:) = 0.0_rp
     momz(ke:ka,:,:) = 0.0_rp
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_atmos_vars,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_VARS. Check!'
        call prc_mpistop
     endif
     if( io_nml ) write(io_fid_nml,nml=param_atmos_vars)
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** List of prognostic variables (ATMOS) ***'
     if( io_l ) write(io_fid_log,'(1x,A,A24,A,A48,A,A12,A)') &
                '***       |', 'VARNAME                 ','|', &
                'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
     do iv = 1, vmax
        if( io_l ) write(io_fid_log,'(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                   '*** NO.',iv,'|',var_name(iv),'|', var_desc(iv),'[', var_unit(iv),']'
     enddo
     do iq = 1, qa
        if( io_l ) write(io_fid_log,'(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                   '*** NO.',5+iq,'|',tracer_name(iq),'|', tracer_desc(iq),'[', tracer_unit(iq),']'
     enddo
 
     if( io_l ) write(io_fid_log,*)
     if ( atmos_restart_in_basename /= '' ) then
        if( io_l ) write(io_fid_log,*) '*** Restart input?  : YES, file = ', trim(atmos_restart_in_basename)
        if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_restart_in_postfix_timelabel
     else
        if( io_l ) write(io_fid_log,*) '*** Restart input?  : NO'
     endif
     if (       atmos_restart_output             &
          .AND. atmos_restart_out_basename /= '' ) then
        if( io_l ) write(io_fid_log,*) '*** Restart output? : YES, file = ', trim(atmos_restart_out_basename)
        if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_restart_out_postfix_timelabel
     else
        if( io_l ) write(io_fid_log,*) '*** Restart output? : NO'
        atmos_restart_output = .false.
     endif
 
     if ( atmos_restart_check_basename == '' ) then
        atmos_restart_check = .false.
     endif
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** Check restart consistency?      : ', atmos_restart_check
     if( io_l ) write(io_fid_log,*) '*** Check value range of variables? : ', atmos_vars_checkrange
     if ( atmos_vars_checkcfl > 0.0_rp ) then
        if( io_l ) write(io_fid_log,*) '*** Check CFL condition?            : YES'
        if( io_l ) write(io_fid_log,*) '*** Limit of Courant number         : ', atmos_vars_checkcfl
     else
        if( io_l ) write(io_fid_log,*) '*** Check CFL condition?            : NO'
     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_cp_vars_setup
 
 
 
     !##### todo: the part below should be moved to the mod_atmos_diag #####
 
     allocate( aq_hist_id(max(qa,1)))
 
     var_hist_id(:) = -1
     aq_hist_id(:) = -1
     ad_hist_id(:) = -1
     ad_monit_id(:) = -1
     ad_prep_sw(:) = -1
 
     call hist_reg( var_hist_id(i_dens), var_name(i_dens), var_desc(i_dens), var_unit(i_dens), ndim=3 )
     call hist_reg( var_hist_id(i_momz), var_name(i_momz), var_desc(i_momz), var_unit(i_momz), ndim=3, zdim='half' )
     call hist_reg( var_hist_id(i_momx), var_name(i_momx), var_desc(i_momx), var_unit(i_momx), ndim=3, xdim='half' )
     call hist_reg( var_hist_id(i_momy), var_name(i_momy), var_desc(i_momy), var_unit(i_momy), ndim=3, ydim='half' )
     call hist_reg( var_hist_id(i_rhot), var_name(i_rhot), var_desc(i_rhot), var_unit(i_rhot), ndim=3 )
     do iq = 1, qa
        call hist_reg( aq_hist_id(iq), tracer_name(iq), tracer_desc(iq), tracer_unit(iq), ndim=3 )
     enddo
 
     call hist_reg( ad_hist_id(i_w)        , 'W',         'velocity w',                     'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_u)        , 'U',         'velocity u',                     'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_v)        , 'V',         'velocity v',                     'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_pott)     , 'PT',        'potential temp.',                'K',      ndim=3 )
 
     call hist_reg( ad_hist_id(i_qdry)     , 'QDRY',      'dry air',                        'kg/kg',  ndim=3 )
     call hist_reg( ad_hist_id(i_qtot)     , 'QTOT',      'total water',                    'kg/kg',  ndim=3 )
     call hist_reg( ad_hist_id(i_qhyd)     , 'QHYD',      'total hydrometeors',             'kg/kg',  ndim=3 )
     call hist_reg( ad_hist_id(i_qliq)     , 'QLIQ',      'total liquid water',             'kg/kg',  ndim=3 )
     call hist_reg( ad_hist_id(i_qice)     , 'QICE',      'total ice water',                'kg/kg',  ndim=3 )
 
     call hist_reg( ad_hist_id(i_lwp)      , 'LWP',       'liquid water path',              'g/m2',   ndim=2 )
     call hist_reg( ad_hist_id(i_iwp)      , 'IWP',       'ice water path',                 'g/m2',   ndim=2 )
     call hist_reg( ad_hist_id(i_pw )      , 'PW',        'precipitable water',             'g/m2',   ndim=2 )
 
     call hist_reg( ad_hist_id(i_rtot)     , 'RTOT',      'Total gas constant',             'J/kg/K', ndim=3 )
     call hist_reg( ad_hist_id(i_cptot)    , 'CPTOT',     'Total heat capacity',            'J/kg/K', ndim=3 )
     call hist_reg( ad_hist_id(i_pres)     , 'PRES',      'pressure',                       'Pa',     ndim=3 )
     call hist_reg( ad_hist_id(i_temp)     , 'T',         'temperature',                    'K',      ndim=3 )
 
     call hist_reg( ad_hist_id(i_potl)     , 'LWPT',      'liq. potential temp.',           'K',      ndim=3 )
     call hist_reg( ad_hist_id(i_rha)      , 'RHA',       'relative humidity(liq+ice)',     '%',      ndim=3 )
     call hist_reg( ad_hist_id(i_rhl)      , 'RH',        'relative humidity(liq)',         '%',      ndim=3 )
     call hist_reg( ad_hist_id(i_rhi)      , 'RHI',       'relative humidity(ice)',         '%',      ndim=3 )
 
     call hist_reg( ad_hist_id(i_vor)      , 'VOR',       'vertical vorticity',             '1/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_div)      , 'DIV',       'divergence',                     '1/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_hdiv)     , 'HDIV',      'horizontal divergence',          '1/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_uabs)     , 'Uabs',      'absolute velocity',              'm/s',    ndim=3 )
 
     call hist_reg( ad_hist_id(i_n2)       , 'N2',        'squared Brunt-Vaisala frequency','1/s2',   ndim=3 )
 
     call hist_reg( ad_hist_id(i_cape)     , 'CAPE',      'convection avail. pot. energy',  'm2/s2',  ndim=2 )
     call hist_reg( ad_hist_id(i_cin)      , 'CIN',       'convection inhibition',          'm2/s2',  ndim=2 )
     call hist_reg( ad_hist_id(i_lcl)      , 'LCL',       'lifted condensation level',      'm',      ndim=2 )
     call hist_reg( ad_hist_id(i_lfc)      , 'LFC',       'level of free convection',       'm',      ndim=2 )
     call hist_reg( ad_hist_id(i_lnb)      , 'LNB',       'level of neutral buoyancy',      'm',      ndim=2 )
 
     call hist_reg( ad_hist_id(i_pblh)     , 'PBLH',      'PBL height',                     'm',      ndim=2 )
     call hist_reg( ad_hist_id(i_mse)      , 'MSE',       'moist static energy',            'm2/s2',  ndim=3 )
 
     call hist_reg( ad_hist_id(i_dens_mean), 'DENS_MEAN', 'horiz. mean of density',         'kg/m3',  ndim=1 )
     call hist_reg( ad_hist_id(i_w_mean)   , 'W_MEAN',    'horiz. mean of w',               'm/s',    ndim=1 )
     call hist_reg( ad_hist_id(i_u_mean)   , 'U_MEAN',    'horiz. mean of u',               'm/s',    ndim=1 )
     call hist_reg( ad_hist_id(i_v_mean)   , 'V_MEAN',    'horiz. mean of v',               'm/s',    ndim=1 )
     call hist_reg( ad_hist_id(i_pott_mean), 'PT_MEAN',   'horiz. mean of pot.',            'K',      ndim=1 )
     call hist_reg( ad_hist_id(i_t_mean)   , 'T_MEAN',    'horiz. mean of t',               'K',      ndim=1 )
     call hist_reg( ad_hist_id(i_qv_mean)  , 'QV_MEAN',   'horiz. mean of QV',              '1',      ndim=1 )
     call hist_reg( ad_hist_id(i_qhyd_mean), 'QHYD_MEAN', 'horiz. mean of QHYD',            '1',      ndim=1 )
     call hist_reg( ad_hist_id(i_qliq_mean), 'QLIQ_MEAN', 'horiz. mean of QLIQ',            '1',      ndim=1 )
     call hist_reg( ad_hist_id(i_qice_mean), 'QICE_MEAN', 'horiz. mean of QICE',            '1',      ndim=1 )
 
     call hist_reg( ad_hist_id(i_dens_prim), 'DENS_PRIM', 'horiz. deviation of density',    'kg/m3',  ndim=3 )
     call hist_reg( ad_hist_id(i_w_prim   ), 'W_PRIM',    'horiz. deviation of w',          'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_u_prim   ), 'U_PRIM',    'horiz. deviation of u',          'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_v_prim   ), 'V_PRIM',    'horiz. deviation of v',          'm/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_pott_prim), 'PT_PRIM',   'horiz. deviation of pot. temp.', 'K',      ndim=3 )
     call hist_reg( ad_hist_id(i_w_prim2  ), 'W_PRIM2',   'variance of w',                  'm2/s2',  ndim=3 )
     call hist_reg( ad_hist_id(i_pt_w_prim), 'PT_W_PRIM', 'resolved scale heat flux',       'W/s',    ndim=3 )
     call hist_reg( ad_hist_id(i_w_prim3  ), 'W_PRIM3',   'skewness of w',                  'm3/s3',  ndim=3 )
     call hist_reg( ad_hist_id(i_tke_rs   ), 'TKE_RS',    'resolved scale TKE',             'm2/s2',  ndim=3 )
 
     call hist_reg( ad_hist_id(i_engt)     , 'ENGT',      'total energy',                   'J/m3',   ndim=3 )
     call hist_reg( ad_hist_id(i_engp)     , 'ENGP',      'potential energy',               'J/m3',   ndim=3 )
     call hist_reg( ad_hist_id(i_engk)     , 'ENGK',      'kinetic energy',                 'J/m3',   ndim=3 )
     call hist_reg( ad_hist_id(i_engi)     , 'ENGI',      'internal energy',                'J/m3',   ndim=3 )
 
     !-----< monitor output setup >-----
 
     call monit_reg( ad_monit_id(i_qdry)        , 'QDRY',         'dry air mass',           'kg', ndim=3, isflux=.false. )
     call monit_reg( ad_monit_id(i_qtot)        , 'QTOT',         'water mass',             'kg', ndim=3, isflux=.false. )
     call monit_reg( ad_monit_id(i_evap)        , 'EVAP',         'evaporation',            'kg', ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_prcp)        , 'PRCP',         'precipitation',          'kg', ndim=2, isflux=.true.  )
 
     call monit_reg( ad_monit_id(i_engt)        , 'ENGT',         'total     energy',       'J',  ndim=3, isflux=.false. )
     call monit_reg( ad_monit_id(i_engp)        , 'ENGP',         'potential energy',       'J',  ndim=3, isflux=.false. )
     call monit_reg( ad_monit_id(i_engk)        , 'ENGK',         'kinetic   energy',       'J',  ndim=3, isflux=.false. )
     call monit_reg( ad_monit_id(i_engi)        , 'ENGI',         'internal  energy',       'J',  ndim=3, isflux=.false. )
 
     call monit_reg( ad_monit_id(i_engflxt)     , 'ENGFLXT',      'total energy flux',      'J',  ndim=2, isflux=.true.  )
 
     call monit_reg( ad_monit_id(i_engsfc_sh)   , 'ENGSFC_SH',    'SFC specific heat flux', 'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engsfc_lh)   , 'ENGSFC_LH',    'SFC latent   heat flux', 'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engsfc_rd)   , 'ENGSFC_RD',    'SFC net radiation flux', 'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engtoa_rd)   , 'ENGTOA_RD',    'TOA net radiation flux', 'J',  ndim=2, isflux=.true.  )
 
     call monit_reg( ad_monit_id(i_engsfc_lw_up), 'ENGSFC_LW_up', 'SFC LW upward   flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engsfc_lw_dn), 'ENGSFC_LW_dn', 'SFC LW downward flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engsfc_sw_up), 'ENGSFC_SW_up', 'SFC SW upward   flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engsfc_sw_dn), 'ENGSFC_SW_dn', 'SFC SW downward flux',   'J',  ndim=2, isflux=.true.  )
 
     call monit_reg( ad_monit_id(i_engtoa_lw_up), 'ENGTOA_LW_up', 'TOA LW upward   flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engtoa_lw_dn), 'ENGTOA_LW_dn', 'TOA LW downward flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engtoa_sw_up), 'ENGTOA_SW_up', 'TOA SW upward   flux',   'J',  ndim=2, isflux=.true.  )
     call monit_reg( ad_monit_id(i_engtoa_sw_dn), 'ENGTOA_SW_dn', 'TOA SW downward flux',   'J',  ndim=2, isflux=.true.  )
 
     if (      ad_hist_id(i_qdry) > 0 &
          .OR. ad_monit_id(i_qdry) > 0 ) then
        ad_prep_sw(i_qdry) = 1
     endif
     if (      ad_hist_id(i_qtot) > 0 &
          .OR. ad_monit_id(i_qtot) > 0 ) then
        ad_prep_sw(i_qdry) = 1
        ad_prep_sw(i_qtot) = 1
     endif
     if ( ad_hist_id(i_qhyd) > 0 ) then
        ad_prep_sw(i_qhyd) = 1
     endif
     if ( ad_hist_id(i_qliq) > 0 ) then
        ad_prep_sw(i_qliq) = 1
     endif
     if ( ad_hist_id(i_qice) > 0 ) then
        ad_prep_sw(i_qice) = 1
     endif
 
     if ( ad_hist_id(i_lwp)  > 0 ) then
        ad_prep_sw(i_qliq) = 1
        ad_prep_sw(i_lwp)  = 1
     endif
     if ( ad_hist_id(i_iwp)  > 0 ) then
        ad_prep_sw(i_qice) = 1
        ad_prep_sw(i_iwp)  = 1
     endif
     if ( ad_hist_id(i_pw)  > 0 ) then
        ad_prep_sw(i_pw)  = 1
     endif
 
     if ( ad_hist_id(i_rtot) > 0 ) then
        ad_prep_sw(i_qdry) = 1
        ad_prep_sw(i_rtot) = 1
     endif
     if ( ad_hist_id(i_cptot) > 0 ) then
        ad_prep_sw(i_qdry)  = 1
        ad_prep_sw(i_cptot) = 1
     endif
 
     if ( ad_hist_id(i_potl) > 0 ) then
        ad_prep_sw(i_qdry)  = 1
        ad_prep_sw(i_rtot)  = 1
        ad_prep_sw(i_cptot) = 1
        ad_prep_sw(i_potl)  = 1
     endif
     if (      ad_hist_id(i_rha) > 0 &
          .OR. ad_hist_id(i_rhl) > 0 &
          .OR. ad_hist_id(i_rhi) > 0 ) then
        ad_prep_sw(i_qdry)  = 1
        ad_prep_sw(i_rtot)  = 1
        ad_prep_sw(i_cptot) = 1
        ad_prep_sw(i_qsat)  = 1
     endif
 
 
     if ( ad_hist_id(i_vor) > 0 ) then
        ad_prep_sw(i_vor) = 1
     endif
 
     if ( ad_prep_sw(i_div) > 0 ) then
        ad_prep_sw(i_hdiv) = 1
     endif
 
     if ( ad_hist_id(i_uabs) > 0 ) then
        ad_prep_sw(i_uabs) = 1
     endif
 
     if (      ad_hist_id(i_cape) > 0 &
          .OR. ad_hist_id(i_cin)  > 0 &
          .OR. ad_hist_id(i_lcl)  > 0 &
          .OR. ad_hist_id(i_lfc)  > 0 &
          .OR. ad_hist_id(i_lnb)  > 0 ) then
        ad_prep_sw(i_cape) = 1
        ad_prep_sw(i_cin)  = 1
        ad_prep_sw(i_lcl)  = 1
        ad_prep_sw(i_lfc)  = 1
        ad_prep_sw(i_lnb)  = 1
     endif
 
     if ( ad_hist_id(i_pblh) > 0 ) then
        ad_prep_sw(i_pblh) = 1
     endif
 
     if ( ad_hist_id(i_mse) > 0 ) then
        ad_prep_sw(i_cptot) = 1
        ad_prep_sw(i_mse)   = 1
     endif
 
     if ( ad_hist_id(i_dens_prim) > 0 ) then
        ad_prep_sw(i_dens_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
     endif
 
     if ( ad_hist_id(i_w_prim) > 0 ) then
        ad_prep_sw(i_w_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_w_mean) = 1
     endif
 
     if ( ad_hist_id(i_u_prim) > 0 ) then
        ad_prep_sw(i_u_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_u_mean) = 1
     endif
 
     if ( ad_hist_id(i_v_prim) > 0 ) then
        ad_prep_sw(i_v_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_v_mean) = 1
     endif
 
     if ( ad_hist_id(i_pott_prim) > 0 ) then
        ad_prep_sw(i_pott_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_pott_mean) = 1
     endif
 
     if ( ad_hist_id(i_dens_mean) > 0 ) then
        ad_prep_sw(i_dens_mean) = 1
     endif
 
     if ( ad_hist_id(i_w_mean) > 0 ) then
        ad_prep_sw(i_w_mean) = 1
     endif
 
     if ( ad_hist_id(i_u_mean) > 0 ) then
        ad_prep_sw(i_u_mean) = 1
     endif
 
     if ( ad_hist_id(i_v_mean) > 0 ) then
        ad_prep_sw(i_v_mean) = 1
     endif
 
     if ( ad_hist_id(i_pott_mean) > 0 ) then
        ad_prep_sw(i_pott_mean) = 1
     end if
 
     if ( ad_hist_id(i_t_mean) > 0 ) then
        ad_prep_sw(i_t_mean) = 1
     end if
 
     if ( ad_hist_id(i_qv_mean) > 0 ) then
        ad_prep_sw(i_qv_mean) = 1
     end if
 
     if ( ad_hist_id(i_qhyd_mean) > 0 ) then
        ad_prep_sw(i_qhyd) = 1
        ad_prep_sw(i_qhyd_mean) = 1
     end if
 
     if ( ad_hist_id(i_qliq_mean) > 0 ) then
        ad_prep_sw(i_qliq) = 1
        ad_prep_sw(i_qliq_mean) = 1
     end if
 
     if ( ad_hist_id(i_qice_mean) > 0 ) then
        ad_prep_sw(i_qice) = 1
        ad_prep_sw(i_qice_mean) = 1
     end if
 
     if ( ad_hist_id(i_w_prim2) > 0 ) then
        ad_prep_sw(i_w_prim)  = 1
        ad_prep_sw(i_w_prim2) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_w_mean)    = 1
     endif
 
     if ( ad_hist_id(i_pt_w_prim) > 0 ) then
        ad_prep_sw(i_w_prim)    = 1
        ad_prep_sw(i_pott_prim) = 1
        ad_prep_sw(i_pt_w_prim) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_w_mean)    = 1
        ad_prep_sw(i_pott_mean) = 1
     endif
 
     if ( ad_hist_id(i_w_prim3) > 0 ) then
        ad_prep_sw(i_w_prim)  = 1
        ad_prep_sw(i_w_prim3) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_w_mean)  = 1
     endif
 
     if ( ad_hist_id(i_tke_rs) > 0 ) then
        ad_prep_sw(i_w_prim) = 1
        ad_prep_sw(i_u_prim) = 1
        ad_prep_sw(i_v_prim) = 1
        ad_prep_sw(i_tke_rs) = 1
        ad_prep_sw(i_dens_mean) = 1
        ad_prep_sw(i_w_mean) = 1
        ad_prep_sw(i_u_mean) = 1
        ad_prep_sw(i_v_mean) = 1
     endif
 
     if (      ad_hist_id(i_engp) > 0 &
          .OR. ad_monit_id(i_engp) > 0 ) then
        ad_prep_sw(i_engp) = 1
     endif
     if (      ad_hist_id(i_engk) > 0 &
          .OR. ad_monit_id(i_engk) > 0 ) then
        ad_prep_sw(i_engk) = 1
     endif
     if (      ad_hist_id(i_engi) > 0 &
          .OR. ad_monit_id(i_engi) > 0 ) then
        ad_prep_sw(i_qdry)  = 1
        ad_prep_sw(i_rtot)  = 1
        ad_prep_sw(i_cptot) = 1
        ad_prep_sw(i_engi)  = 1
     endif
     if (      ad_hist_id(i_engt) > 0 &
          .OR. ad_monit_id(i_engt) > 0 ) then
        ad_prep_sw(i_engp)  = 1
        ad_prep_sw(i_engk)  = 1
        ad_prep_sw(i_qdry)  = 1
        ad_prep_sw(i_rtot)  = 1
        ad_prep_sw(i_cptot) = 1
        ad_prep_sw(i_engi)  = 1
        ad_prep_sw(i_engt)  = 1
     endif
 
     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 760 of file mod_atmos_vars.f90.

References dens, scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ka, scale_grid_index::ke, scale_grid_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, 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-1,i,j) = momz(ks,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 824 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_read().

     use scale_process, only: &
        prc_mpistop
     use scale_const, only: &
        grav  => const_grav
     use scale_time, only: &
        time_gettimelabel
     use scale_fileio, only: &
        fileio_open, &
        fileio_check_coordinates
     use scale_atmos_thermodyn, only: &
        thermodyn_qd        => atmos_thermodyn_qd,        &
        thermodyn_temp_pres => atmos_thermodyn_temp_pres
     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
     implicit none
 
     character(len=19)     :: timelabel
     character(len=H_LONG) :: basename
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '*** 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
 
        if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)
 
        call fileio_open( restart_fid, basename )
 
        if ( atmos_restart_in_check_coordinates ) then
           call fileio_check_coordinates( restart_fid, atmos=.true. )
        end if
 
     else
        write(*,*) '*** restart file for atmosphere is not specified. STOP!'
        call prc_mpistop
     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 ) 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 917 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_read(), and mod_rm_driver::resume_state().

     use scale_process, only: &
        prc_mpistop
     use scale_fileio, only: &
        fileio_read, &
        fileio_flush
     use scale_atmos_thermodyn, only: &
        thermodyn_qd        => atmos_thermodyn_qd,        &
        thermodyn_temp_pres => atmos_thermodyn_temp_pres
     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
     implicit none
 
     integer  :: i, j, iq
     !---------------------------------------------------------------------------
 
     if ( restart_fid /= -1 ) then
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '*** Read from restart file (ATMOS) ***'
 
        call fileio_read( dens(:,:,:),                            & ! [OUT]
                          restart_fid, var_name(1), 'ZXY', step=1 ) ! [IN]
        call fileio_read( momz(:,:,:),                            & ! [OUT]
                          restart_fid, var_name(2), 'ZXY', step=1 ) ! [IN]
        call fileio_read( momx(:,:,:),                            & ! [OUT]
                          restart_fid, var_name(3), 'ZXY', step=1 ) ! [IN]
        call fileio_read( momy(:,:,:),                            & ! [OUT]
                          restart_fid, var_name(4), 'ZXY', step=1 ) ! [IN]
        call fileio_read( rhot(:,:,:),                            & ! [OUT]
                          restart_fid, var_name(5), 'ZXY', step=1 ) ! [IN]
 
        do iq = 1, qa
           call fileio_read( qtrc(:,:,:,iq),                             & ! [OUT]
                             restart_fid, tracer_name(iq), 'ZXY', step=1 ) ! [IN]
        enddo
 
        if ( io_aggregate ) then
           call fileio_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-1,i,j) = momz(ks,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
        write(*,*) '*** invalid restart file ID for atmosphere. STOP!'
        call prc_mpistop
     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 ) 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 1028 of file mod_atmos_vars.f90.

References scale_atmos_bottom::atmos_bottom_estimate(), dens_av, scale_history::hist_setpres(), phyd, pres, scale_grid_real::real_cz, scale_grid_real::real_z1, and scale_topography::topo_zsfc.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), and mod_rm_driver::resume_state().

     use scale_grid_real, only: &
        real_cz, &
        real_z1
     use scale_topography, only: &
        topo_zsfc
     use scale_atmos_bottom, only: &
        bottom_estimate => atmos_bottom_estimate
     use scale_history, only: &
        hist_setpres
     implicit none
 
     real(RP) :: SFC_DENS(IA,JA)
     real(RP) :: SFC_PRES(IA,JA)
     !---------------------------------------------------------------------------
 
     call bottom_estimate( dens_av(:,:,:), & ! [IN]
                           pres(:,:,:), & ! [IN]
                           real_cz(:,:,:), & ! [IN]
                           topo_zsfc(:,:),   & ! [IN]
                           real_z1(:,:),   & ! [IN]
                           sfc_dens(:,:),   & ! [OUT]
                           sfc_pres(:,:)    ) ! [OUT]
 
     call hist_setpres( phyd(:,:,:),  & ! [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 1060 of file mod_atmos_vars.f90.

References atmos_restart_check_basename, atmos_restart_check_criterion, dens, scale_fileio::fileio_close(), scale_fileio::fileio_flush(), scale_fileio::fileio_open(), scale_grid_index::ie, scale_stdio::io_aggregate, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, momx, momy, momz, scale_process::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::scalerm().

     use scale_process, only: &
        prc_myrank
     use scale_fileio, only: &
        fileio_open, &
        fileio_read, &
        fileio_flush, &
        fileio_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')
 
     write(*,*) 'Compare last Data with ', trim(atmos_restart_check_basename), 'on PE=', prc_myrank
     write(*,*) '*** criterion = ', atmos_restart_check_criterion
     datacheck = .true.
 
     basename = atmos_restart_check_basename
 
     call fileio_open( fid, basename )
 
     call fileio_read( dens_check(:,:,:), fid, 'DENS', 'ZXY', step=1 )
     call fileio_read( momz_check(:,:,:), fid, 'MOMZ', 'ZXY', step=1 )
     call fileio_read( momx_check(:,:,:), fid, 'MOMX', 'ZXY', step=1 )
     call fileio_read( momy_check(:,:,:), fid, 'MOMY', 'ZXY', step=1 )
     call fileio_read( rhot_check(:,:,:), fid, 'RHOT', 'ZXY', step=1 )
     do iq = 1, qa
        call fileio_read( qtrc_check(:,:,:,iq), fid, tracer_name(iq), 'ZXY', step=1 )
     end do
     if ( io_aggregate ) call fileio_flush( fid )
 
     call fileio_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
           write(*,*) 'xxx there is the difference  : ', dens(k,i,j)-dens_check(k,i,j)
           write(*,*) 'xxx at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'DENS'
           datacheck = .false.
        endif
     enddo
     enddo
     enddo
 
     do k = ks, 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
           write(*,*) 'xxx there is the difference  : ', momz(k,i,j)-momz_check(k,i,j)
           write(*,*) 'xxx 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
           write(*,*) 'xxx there is the difference  : ', momx(k,i,j)-momx_check(k,i,j)
           write(*,*) 'xxx 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
           write(*,*) 'xxx there is the difference  : ', momy(k,i,j)-momy_check(k,i,j)
           write(*,*) 'xxx 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
           write(*,*) 'xxx there is the difference  : ', rhot(k,i,j)-rhot_check(k,i,j)
           write(*,*) 'xxx 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
              write(*,*) 'xxx there is the difference  : ', qtrc(k,i,j,iq)-qtrc_check(k,i,j,iq)
              write(*,*) 'xxx 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
        if( io_l ) write(io_fid_log,*) 'Data Check Clear.'
        write(*,*) 'Data Check Clear.'
     else
        if( io_l ) write(io_fid_log,*) 'Data Check Failed. See std. output.'
        write(*,*) '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 1195 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), and mod_atmos_driver::atmos_driver_resume2().

     use scale_const, only: &
        grav  => const_grav,  &
        rdry  => const_rdry,  &
        rvap  => const_rvap,  &
        cpdry => const_cpdry, &
        cvdry => const_cvdry, &
        p00   => const_pre00
     use scale_grid, only: &
        rcdx => grid_rcdx, &
        rcdy => grid_rcdy
     use scale_grid_real, only: &
        real_cz, &
        real_fz
     use scale_comm, only: &
        comm_horizontal_mean
     use scale_history, only: &
        hist_in
     use scale_atmos_thermodyn, only: &
        thermodyn_qd => atmos_thermodyn_qd
     use scale_atmos_hydrometeor, only: &
        hydrometeor_lhv => atmos_hydrometeor_lhv, &
        i_qv, &
        i_qc, &
        qhs,  &
        qhe,  &
        qls,  &
        qle,  &
        qis,  &
        qie,  &
        lhv,  &
        lhf
     use scale_atmos_saturation, only: &
        saturation_psat_all => atmos_saturation_psat_all, &
        saturation_psat_liq => atmos_saturation_psat_liq, &
        saturation_psat_ice => atmos_saturation_psat_ice
     use scale_atmos_adiabat, only: &
        adiabat_cape => atmos_adiabat_cape
     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
     implicit none
 
     real(RP) :: QDRY  (KA,IA,JA) ! dry air            [kg/kg]
     real(RP) :: QTOT  (KA,IA,JA) ! total water        [kg/kg]
     real(RP) :: QHYD  (KA,IA,JA) ! total hydrometeor  [kg/kg]
     real(RP) :: QLIQ  (KA,IA,JA) ! total liquid water [kg/kg]
     real(RP) :: QICE  (KA,IA,JA) ! total ice water    [kg/kg]
     real(RP) :: RHOQ  (KA,IA,JA)
 
     real(RP) :: LWP   (IA,JA)    ! liquid water path  [g/m2]
     real(RP) :: IWP   (IA,JA)    ! ice    water path  [g/m2]
     real(RP) :: PW    (IA,JA)    ! precipitable water [g/m2]
 
     real(RP) :: RTOT  (KA,IA,JA) ! Total gas constant  [J/kg/K]
     real(RP) :: CPTOT (KA,IA,JA) ! Total heat capacity [J/kg/K]
     real(RP) :: CVTOT (KA,IA,JA) ! Total heat capacity [J/kg/K]
     real(RP) :: CPovCV(KA,IA,JA) ! Cp/Cv
 
     real(RP) :: POTL  (KA,IA,JA) ! liquid water potential temperature [K]
     real(RP) :: RHA   (KA,IA,JA) ! relative humidity (liquid+ice)      [%]
     real(RP) :: RHL   (KA,IA,JA) ! relative humidity against to liquid [%]
     real(RP) :: RHI   (KA,IA,JA) ! relative humidity against to ice    [%]
 
     real(RP) :: VOR   (KA,IA,JA) ! vertical vorticity    [1/s]
     real(RP) :: DIV   (KA,IA,JA) ! divergence            [1/s]
     real(RP) :: HDIV  (KA,IA,JA) ! horizontal divergence [1/s]
     real(RP) :: Uabs  (KA,IA,JA) ! absolute velocity     [m/s]
 
     real(RP) :: CAPE  (IA,JA)    ! CAPE [m2/s2]
     real(RP) :: CIN   (IA,JA)    ! CIN [m2/s2]
     real(RP) :: LCL   (IA,JA)    ! LCL height [m]
     real(RP) :: LFC   (IA,JA)    ! LFC height [m]
     real(RP) :: LNB   (IA,JA)    ! LNB height [m]
 
     real(RP) :: PBLH  (IA,JA)    ! PBL height [m]
     real(RP) :: POTTv (KA,IA,JA) ! vertual potential temperature [K]
     real(RP) :: fact
 
     real(RP) :: MSE      (KA,IA,JA) ! MSE        [m2/s2]
     real(RP) :: LHV_local(KA,IA,JA) ! latent heat for vaporization [m2/s2]
 
     real(RP) :: PREC  (IA,JA)    ! surface precipitation rate CP+MP(rain+snow) [kg/m2/s]
     real(RP) :: RAIN  (IA,JA)    ! surface rain rate CP+MP [kg/m2/s]
     real(RP) :: SNOW  (IA,JA)    ! surface snow rate CP+MP [kg/m2/s]
 
     real(RP) :: DENS_PRIM(KA,IA,JA) ! horiz. deviation of density    [kg/m3]
     real(RP) :: W_PRIM   (KA,IA,JA) ! horiz. deviation of w          [m/s]
     real(RP) :: U_PRIM   (KA,IA,JA) ! horiz. deviation of u          [m/s]
     real(RP) :: V_PRIM   (KA,IA,JA) ! horiz. deviation of v          [m/s]
     real(RP) :: POTT_PRIM(KA,IA,JA) ! horiz. deviation of pot. temp. [K]
     real(RP) :: W_PRIM2  (KA,IA,JA) ! variance of w                  [m2/s2]
     real(RP) :: PT_W_PRIM(KA,IA,JA) ! resolved scale heat flux       [W/s]
     real(RP) :: W_PRIM3  (KA,IA,JA) ! skewness of w                  [m3/s3]
     real(RP) :: TKE_RS   (KA,IA,JA) ! resolved scale TKE             [m2/s2]
     real(RP) :: DENS_MEAN(KA)       ! horiz. mean of density         [kg/m3]
     real(RP) :: W_MEAN   (KA)       ! horiz. mean of w               [m/s]
     real(RP) :: U_MEAN   (KA)       ! horiz. mean of u               [m/s]
     real(RP) :: V_MEAN   (KA)       ! horiz. mean of v               [m/s]
     real(RP) :: PT_MEAN  (KA)       ! horiz. mean of pot.            [K]
     real(RP) :: T_MEAN   (KA)       ! horiz. mean of t               [K]
     real(RP) :: QV_MEAN  (KA)       ! horiz. mean of QV
     real(RP) :: QHYD_MEAN(KA)       ! horiz. mean of QHYD
     real(RP) :: QLIQ_MEAN(KA)       ! horiz. mean of QLIQ
     real(RP) :: QICE_MEAN(KA)       ! horiz. mean of QICE
 
     real(RP) :: ENGT  (KA,IA,JA) ! total     energy [J/m3]
     real(RP) :: ENGP  (KA,IA,JA) ! potential energy [J/m3]
     real(RP) :: ENGK  (KA,IA,JA) ! kinetic   energy [J/m3]
     real(RP) :: ENGI  (KA,IA,JA) ! internal  energy [J/m3]
 
     real(RP) :: PSAT  (KA,IA,JA)
     real(RP) :: UH    (KA,IA,JA)
     real(RP) :: VH    (KA,IA,JA)
 
     integer :: k, i, j, iq
     !---------------------------------------------------------------------------
 
     ! value check for prognostic variables
     if ( atmos_vars_checkrange ) then
        call valcheck( dens(:,:,:),    0.0_rp,    2.0_rp, var_name(i_dens), __file__, __line__ )
        call valcheck( momz(:,:,:), -200.0_rp,  200.0_rp, var_name(i_momz), __file__, __line__ )
        call valcheck( momx(:,:,:), -200.0_rp,  200.0_rp, var_name(i_momx), __file__, __line__ )
        call valcheck( momy(:,:,:), -200.0_rp,  200.0_rp, var_name(i_momy), __file__, __line__ )
        call valcheck( rhot(:,:,:),    0.0_rp, 1000.0_rp, var_name(i_rhot), __file__, __line__ )
     endif
 
     ! history output of prognostic variables
     call hist_in( dens(:,:,:), var_name(i_dens), var_desc(i_dens), var_unit(i_dens) )
     call hist_in( momz(:,:,:), var_name(i_momz), var_desc(i_momz), var_unit(i_momz) )
     call hist_in( momx(:,:,:), var_name(i_momx), var_desc(i_momx), var_unit(i_momx) )
     call hist_in( momy(:,:,:), var_name(i_momy), var_desc(i_momy), var_unit(i_momy) )
     call hist_in( rhot(:,:,:), var_name(i_rhot), var_desc(i_rhot), var_unit(i_rhot) )
     do iq = 1, qa
        call hist_in( qtrc(:,:,:,iq), tracer_name(iq), tracer_desc(iq), tracer_unit(iq) )
     enddo
 
     ! prepare and history output of diagnostic variables
 
     if ( ad_prep_sw(i_qdry) > 0 ) then
        call thermodyn_qd( qdry(:,:,:),   & ! [OUT]
                           qtrc_av(:,:,:,:), & ! [IN]
                           tracer_mass(:)    ) ! [IN]
     endif
 
     if ( ad_prep_sw(i_qtot) > 0 ) then
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           qtot(k,i,j) = 1.0_rp - qdry(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_qhyd) > 0 ) then
 !OCL XFILL
        qhyd(:,:,:) = 0.0_rp
        do iq = qhs, qhe
          qhyd(:,:,:) = qhyd(:,:,:) + qtrc_av(:,:,:,iq)
        enddo
     endif
 
     if ( ad_prep_sw(i_qliq) > 0 ) then
 !OCL XFILL
        qliq(:,:,:) = 0.0_rp
        do iq = qls, qle
          qliq(:,:,:) = qliq(:,:,:) + qtrc_av(:,:,:,iq)
        enddo
     endif
 
     if ( ad_prep_sw(i_qice) > 0 ) then
 !OCL XFILL
        qice(:,:,:) = 0.0_rp
        do iq = qis, qie
          qice(:,:,:) = qice(:,:,:) + qtrc_av(:,:,:,iq)
        enddo
     endif
 
     if ( ad_prep_sw(i_lwp) > 0 ) then
        do j  = jsb, jeb
        do i  = isb, ieb
           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
     endif
 
     if ( ad_prep_sw(i_iwp) > 0 ) then
        do j  = jsb, jeb
        do i  = isb, ieb
           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
     endif
 
     if ( ad_prep_sw(i_pw) > 0 ) then
        do j  = jsb, jeb
        do i  = isb, ieb
           pw(i,j) = 0.0_rp
           do k  = ks, ke
              pw(i,j) = pw(i,j) &
                       + qtrc_av(k,i,j,i_qv) * 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
     endif
 
     if ( ad_prep_sw(i_rtot) > 0 ) then
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           calc_r(rtot(k,i,j), qdry(k,i,j), qtrc_av, k, i, j, iq, rdry, tracer_r)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_cptot) > 0 ) then
        !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
        !$omp private(iq) &
        !$omp shared(JSB,JEB,ISB,IEB,KS,KE,CPTOT,CPdry,QDRY,QA,CVTOT,CVdry,QTRC_av,TRACER_CP,TRACER_CV) &
        !$omp shared(I_QV,LHV,QIS,QIE,LHF)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           cptot(k,i,j) = cpdry * qdry(k,i,j)
           cvtot(k,i,j) = cvdry * qdry(k,i,j)
           do iq = 1, qa
              cptot(k,i,j) = cptot(k,i,j) + qtrc_av(k,i,j,iq) * tracer_cp(iq)
              cvtot(k,i,j) = cvtot(k,i,j) + qtrc_av(k,i,j,iq) * tracer_cv(iq)
           end do
        enddo
        enddo
        enddo
 
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           cpovcv(k,i,j) = cptot(k,i,j) / cvtot(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_potl) > 0 ) then
        call hydrometeor_lhv( lhv_local(:,:,:), temp(:,:,:) )
 
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k = ks, ke
           potl(k,i,j) = pott(k,i,j) &
                       - lhv_local(k,i,j) / cpdry * qliq(k,i,j) * pott(k,i,j) / temp(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_qsat) > 0 ) then
 !       call SATURATION_dens2qsat_all( QSAT   (:,:,:), & ! [OUT]
 !                                      TEMP   (:,:,:), & ! [IN]
 !                                      DENS_av(:,:,:)  ) ! [IN]
     end if
 
     if ( ad_hist_id(i_rha) > 0 ) then
        call saturation_psat_all( psat(:,:,:), temp(:,:,:) )
 
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           rha(k,i,j) = dens_av(k,i,j) * qtrc_av(k,i,j,i_qv) &
                      / psat(k,i,j) * rvap * temp(k,i,j) &
                      * 100.0_rp
        enddo
        enddo
        enddo
     endif
 
     if ( ad_hist_id(i_rhl) > 0 ) then
        call saturation_psat_liq( psat(:,:,:), & ! [OUT]
                                  temp(:,:,:)  ) ! [IN]
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           rhl(k,i,j) = dens_av(k,i,j) * qtrc_av(k,i,j,i_qv) &
                      / psat(k,i,j) * rvap * temp(k,i,j) &
                      * 100.0_rp
        enddo
        enddo
        enddo
     endif
 
     if ( ad_hist_id(i_rhi) > 0 ) then
        call saturation_psat_ice( psat(:,:,:), & ! [OUT]
                                  temp(:,:,:)  ) ! [IN]
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           rhi(k,i,j) = dens_av(k,i,j) * qtrc_av(k,i,j,i_qv) &
                      / psat(k,i,j) * rvap * temp(k,i,j) &
                      * 100.0_rp
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_vor) > 0 ) then
        ! 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
     endif
 
     if ( ad_prep_sw(i_hdiv) > 0 ) then
        !$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
     endif
 
     if ( ad_prep_sw(i_div) > 0 ) then
        !$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
     endif
 
     if ( ad_prep_sw(i_uabs) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = ks, ke
           uabs(k,i,j) = sqrt( u(k,i,j) * u(k,i,j) &
                             + v(k,i,j) * v(k,i,j) )
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_dens_mean) > 0 ) then
        call comm_horizontal_mean( dens_mean(:), dens_av(:,:,:) )
     end if
 
     if ( ad_prep_sw(i_dens_prim) > 0 ) then
 !OCL XFILL
        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
     endif
 
     if ( ad_prep_sw(i_w_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           w_prim(k,i,j) = w(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( w_mean(:), w_prim(:,:,:) )
        do k = 1, ka
           w_mean(k) = w_mean(k) / dens_mean(k)
        enddo
     end if
     if ( ad_prep_sw(i_w_prim) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           w_prim(k,i,j) = w(k,i,j) - w_mean(k)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_u_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           u_prim(k,i,j) = u(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( u_mean(:), u_prim(:,:,:) )
        do k = 1, ka
           u_mean(k) = u_mean(k) / dens_mean(k)
        enddo
     end if
     if ( ad_prep_sw(i_u_prim) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           u_prim(k,i,j) = u(k,i,j) - u_mean(k)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_v_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           v_prim(k,i,j) = v(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( v_mean(:), v_prim(:,:,:) )
        do k = 1, ka
           v_mean(k) = v_mean(k) / dens_mean(k)
        enddo
     end if
     if ( ad_prep_sw(i_v_prim) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           v_prim(k,i,j) = v(k,i,j) - v_mean(k)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_t_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           pott_prim(k,i,j) = temp(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( t_mean(:), pott_prim(:,:,:) )
        do k = 1, ka
           t_mean(k) = t_mean(k) / dens_mean(k)
        enddo
     end if
 
     if ( ad_prep_sw(i_pott_mean) > 0 ) then
        call comm_horizontal_mean( pt_mean(:), rhot_av(:,:,:) )
        do k = 1, ka
           pt_mean(k) = pt_mean(k) / dens_mean(k)
        enddo
     end if
     if ( ad_prep_sw(i_pott_prim) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           pott_prim(k,i,j) = pott(k,i,j) - pt_mean(k)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_qv_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           rhoq(k,i,j) = qtrc_av(k,i,j,i_qv) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( qv_mean(:), rhoq(:,:,:) )
        do k = 1, ka
           qv_mean(k) = qv_mean(k) / dens_mean(k)
        enddo
     end if
 
     if ( ad_prep_sw(i_qhyd_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           rhoq(k,i,j) = qhyd(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( qhyd_mean(:), rhoq(:,:,:) )
        do k = 1, ka
           qhyd_mean(k) = qhyd_mean(k) / dens_mean(k)
        enddo
     end if
 
     if ( ad_prep_sw(i_qliq_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           rhoq(k,i,j) = qliq(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( qliq_mean(:), rhoq(:,:,:) )
        do k = 1, ka
           qliq_mean(k) = qliq_mean(k) / dens_mean(k)
        enddo
     end if
 
     if ( ad_prep_sw(i_qice_mean) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = 1, ka
           rhoq(k,i,j) = qice(k,i,j) * dens_av(k,i,j)
        enddo
        enddo
        enddo
        call comm_horizontal_mean( qice_mean(:), rhoq(:,:,:) )
        do k = 1, ka
           qice_mean(k) = qice_mean(k) / dens_mean(k)
        enddo
     end if
 
     if ( ad_prep_sw(i_w_prim2) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = ks, ke
           w_prim2(k,i,j) = w_prim(k,i,j) * w_prim(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_pt_w_prim) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = ks, ke
           pt_w_prim(k,i,j) = w_prim(k,i,j) * pott_prim(k,i,j) * dens_av(k,i,j) * cpdry
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_w_prim3) > 0 ) then
 !OCL XFILL
        do j = 1, ja
        do i = 1, ia
        do k = ks, ke
           w_prim3(k,i,j) = w_prim(k,i,j) * w_prim(k,i,j) * w_prim(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_tke_rs) > 0 ) then
 !OCL XFILL
        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) * w_prim(k,i,j) &
                                    + u_prim(k,i,j) * u_prim(k,i,j) &
                                    + v_prim(k,i,j) * v_prim(k,i,j) )
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_engp) > 0 ) then
        !$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
           engp(k,i,j) = dens_av(k,i,j) * grav * real_cz(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_engk) > 0 ) then
        !$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
           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 )
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_engi) > 0 ) then
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k = ks, ke
           engi(k,i,j) = dens_av(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
        enddo
        enddo
        enddo
 
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
        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
        enddo
        enddo
        enddo
 
        if ( i_qv > 0 ) then
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           engi(k,i,j) = engi(k,i,j) &
                       + dens_av(k,i,j) * qtrc_av(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]
        enddo
        enddo
        enddo
        end if
 
        do iq = qis, qie
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           engi(k,i,j) = engi(k,i,j) &
                       - dens_av(k,i,j) * qtrc_av(k,i,j,iq) * lhf ! Latent Heat [ice->liquid]
        enddo
        enddo
        enddo
        enddo
     endif
 
     if ( ad_prep_sw(i_engt) > 0 ) then
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
        do j  = jsb, jeb
        do i  = isb, ieb
        do k  = ks, ke
           engt(k,i,j) = engp(k,i,j) + engk(k,i,j) + engi(k,i,j)
        enddo
        enddo
        enddo
     endif
 
     call hist_in( w(:,:,:), 'W',     'velocity w',             'm/s'    )
     call hist_in( u(:,:,:), 'U',     'velocity u',             'm/s'    )
     call hist_in( v(:,:,:), 'V',     'velocity v',             'm/s'    )
     call hist_in( pott(:,:,:), 'PT',    'potential temp.',        'K'      )
 
     call hist_in( qdry(:,:,:), 'QDRY',  'dry air',                'kg/kg'  )
     call hist_in( qtot(:,:,:), 'QTOT',  'total water',            'kg/kg'  )
     call hist_in( qhyd(:,:,:), 'QHYD',  'total hydrometeors',     'kg/kg'  )
     call hist_in( qliq(:,:,:), 'QLIQ',  'total liquid water',     'kg/kg'  )
     call hist_in( qice(:,:,:), 'QICE',  'total ice water',        'kg/kg'  )
 
     call hist_in( lwp(:,:),   'LWP',   'liquid water path',      'g/m2'   )
     call hist_in( iwp(:,:),   'IWP',   'ice    water path',      'g/m2'   )
     call hist_in( pw(:,:),   'PW',    'precipitable water',     'g/m2'   )
 
     call hist_in( rtot(:,:,:), 'RTOT',  'Total gas constant',     'J/kg/K' )
     call hist_in( cptot(:,:,:), 'CPTOT', 'Total heat capacity',    'J/kg/K' )
     call hist_in( pres(:,:,:), 'PRES',  'pressure',               'Pa'     )
     call hist_in( temp(:,:,:), 'T',     'temperature',            'K'      )
 
     call hist_in( potl(:,:,:), 'LWPT',  'liq. potential temp.',   'K'      )
     call hist_in( rha(:,:,:), 'RHA',   'relative humidity(liq+ice)','%'   )
     call hist_in( rhl(:,:,:), 'RH' ,   'relative humidity(liq)', '%'      )
     call hist_in( rhi(:,:,:), 'RHI',   'relative humidity(ice)', '%'      )
 
     call hist_in( vor(:,:,:), 'VOR',   'vertical vorticity',     '1/s'    )
     call hist_in( div(:,:,:), 'DIV',   'divergence',             '1/s'    )
     call hist_in( hdiv(:,:,:), 'HDIV',  'horizontal divergence',  '1/s'    )
     call hist_in( uabs(:,:,:), 'Uabs',  'absolute velocity',      'm/s'    )
 
     call hist_in( n2(:,:,:), 'N2',    'squared Brunt-Vaisala frequency','1/s2' )
 
     call hist_in( dens_mean(:), 'DENS_MEAN', 'horiz. mean of density',    'kg/m3' )
     call hist_in( w_mean(:), 'W_MEAN',    'horiz. mean of w',          'm/s' )
     call hist_in( u_mean(:), 'U_MEAN',    'horiz. mean of u',          'm/s' )
     call hist_in( v_mean(:), 'V_MEAN',    'horiz. mean of v',          'm/s' )
     call hist_in( pt_mean(:), 'PT_MEAN',   'horiz. mean of pot.',       'K' )
     call hist_in( t_mean(:), 'T_MEAN',    'horiz. mean of t',          'K' )
     call hist_in( qv_mean(:), 'QV_MEAN',   'horiz. mean of QV',         '1' )
     call hist_in( qhyd_mean(:), 'QHYD_MEAN', 'horiz. mean of QHYD',       '1' )
     call hist_in( qliq_mean(:), 'QLIQ_MEAN', 'horiz. mean of QLIQ',       '1' )
     call hist_in( qice_mean(:), 'QICE_MEAN', 'horiz. mean of QICE',       '1' )
 
     call hist_in( dens_prim(:,:,:), 'DENS_PRIM', 'horiz. deviation of density',    'kg/m3' )
     call hist_in( w_prim(:,:,:), 'W_PRIM',    'horiz. deviation of w',          'm/s'   )
     call hist_in( u_prim(:,:,:), 'U_PRIM',    'horiz. deviation of u',          'm/s'   )
     call hist_in( v_prim(:,:,:), 'V_PRIM',    'horiz. deviation of v',          'm/s'   )
     call hist_in( pott_prim(:,:,:), 'PT_PRIM',   'horiz. deviation of pot. temp.', 'K'     )
     call hist_in( w_prim2(:,:,:), 'W_PRIM2',   'variance of w',                  'm2/s2' )
     call hist_in( pt_w_prim(:,:,:), 'PT_W_PRIM', 'resolved scale heat flux',       'W/s'   )
     call hist_in( w_prim3(:,:,:), 'W_PRIM3',   'skewness of w',                  'm3/s3' )
     call hist_in( tke_rs(:,:,:), 'TKE_RS',    'resolved scale TKE',             'm2/s2' )
 
     call hist_in( engt(:,:,:), 'ENGT',  'total energy',           'J/m3'   )
     call hist_in( engp(:,:,:), 'ENGP',  'potential energy',       'J/m3'   )
     call hist_in( engk(:,:,:), 'ENGK',  'kinetic energy',         'J/m3'   )
     call hist_in( engi(:,:,:), 'ENGI',  'internal energy',        'J/m3'   )
 
     if (      ad_prep_sw(i_cape) > 0 &
          .OR. ad_prep_sw(i_cin)  > 0 &
          .OR. ad_prep_sw(i_lcl)  > 0 &
          .OR. ad_prep_sw(i_lfc)  > 0 &
          .OR. ad_prep_sw(i_lnb)  > 0 ) then
 
        call adiabat_cape( ks,               & ! [IN]
                           dens_av(:,:,:),   & ! [IN]
                           temp(:,:,:),   & ! [IN]
                           pres(:,:,:),   & ! [IN]
                           qtrc_av(:,:,:,:), & ! [IN]
                           real_cz(:,:,:),   & ! [IN]
                           real_fz(:,:,:),   & ! [IN]
                           cape(:,:),     & ! [OUT]
                           cin(:,:),     & ! [OUT]
                           lcl(:,:),     & ! [OUT]
                           lfc(:,:),     & ! [OUT]
                           lnb(:,:)      ) ! [OUT]
 
     endif
 
     call hist_in( cape(:,:), 'CAPE', 'convection avail. pot. energy', 'm2/s2' )
     call hist_in( cin(:,:), 'CIN',  'convection inhibition',         'm2/s2' )
     call hist_in( lcl(:,:), 'LCL',  'lifted condensation level',     'm'     )
     call hist_in( lfc(:,:), 'LFC',  'level of free convection',      'm'     )
     call hist_in( lnb(:,:), 'LNB',  'level of neutral buoyancy',     'm'     )
 
     if ( ad_prep_sw(i_pblh) > 0 ) then
        do j = js, je
        do i = is, ie
        do k = ks, ke
           fact = 1.0_rp
           if ( i_qv > 0 ) fact = fact + 0.61_rp * qtrc_av(k,i,j,i_qv)
           if ( i_qc > 0 ) fact = fact - 1.61_rp * qtrc_av(k,i,j,i_qc)
           pottv(k,i,j) = pott(k,i,j) * fact
        enddo
        enddo
        enddo
 
        do j = js, je
        do i = is, ie
           pblh(i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j)
 
           do k = ks+1, ke
              if ( pottv(k,i,j) > pottv(ks,i,j) ) then
                 fact = ( pottv(ks,i,j) - pottv(k-1,i,j) ) &
                      / ( pottv(k,i,j)  - pottv(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
     endif
     call hist_in( pblh(:,:), 'PBLH', 'PBL height', 'm' )
 
     if ( ad_prep_sw(i_mse) > 0 ) then
        call hydrometeor_lhv( lhv_local(:,:,:), temp(:,:,:) )
 
        do j = js, je
        do i = is, ie
        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_local(k,i,j) * qtrc_av(k,i,j,i_qv)
        enddo
        enddo
        enddo
     endif
     call hist_in( mse(:,:,:), 'MSE', 'moist static energy', 'm2/s2' )
 
     do j = js, je
     do i = is, ie
        prec(i,j) = sflx_rain_mp(i,j) + sflx_snow_mp(i,j) &
                  + sflx_rain_cp(i,j)
        rain(i,j) = sflx_rain_mp(i,j)                     &
                  + sflx_rain_cp(i,j)
        snow(i,j) = sflx_snow_mp(i,j)
     enddo
     enddo
     call hist_in( prec(:,:), 'PREC', 'surface precipitation rate (total)', 'kg/m2/s' )
     call hist_in( rain(:,:), 'RAIN', 'surface rain rate (total)', 'kg/m2/s' )
     call hist_in( snow(:,:), 'SNOW', 'surface snow rate (total)', 'kg/m2/s' )
 
     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 2091 of file mod_atmos_vars.f90.

References scale_const::const_cvdry, scale_const::const_grav, dens, dens_av, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_atmos_hydrometeor::i_qv, scale_index::i_rhot, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, scale_atmos_hydrometeor::lhf, scale_atmos_hydrometeor::lhv, momx, momx_av, momy, momy_av, momz, momz_av, pres, scale_tracer::qa, scale_atmos_hydrometeor::qie, scale_atmos_hydrometeor::qis, qtrc, qtrc_av, scale_grid_real::real_cz, rhot, rhot_av, scale_rm_statistics::statistics_checktotal, temp, scale_tracer::tracer_cv, scale_tracer::tracer_mass, scale_tracer::tracer_name, scale_tracer::tracer_r, u, v, and w.

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_grid_real, only: &
        real_cz
     use scale_rm_statistics, only: &
        statistics_checktotal, &
        stat_total
     use scale_atmos_thermodyn, only: &
        thermodyn_qd        => atmos_thermodyn_qd,        &
        thermodyn_temp_pres => atmos_thermodyn_temp_pres
     use scale_atmos_hydrometeor, only: &
        i_qv, &
        qis,  &
        qie,  &
        lhv,  &
        lhf
     implicit none
 
     real(RP) :: W   (KA,IA,JA) ! velocity w at cell center [m/s]
     real(RP) :: U   (KA,IA,JA) ! velocity u at cell center [m/s]
     real(RP) :: V   (KA,IA,JA) ! velocity v at cell center [m/s]
 
     real(RP) :: QDRY(KA,IA,JA) ! dry air     [kg/kg]
     real(RP) :: PRES(KA,IA,JA) ! pressure    [Pa]
     real(RP) :: TEMP(KA,IA,JA) ! temperature [K]
 
     real(RP) :: ENGT(KA,IA,JA) ! total     energy [J/m3]
     real(RP) :: ENGP(KA,IA,JA) ! potential energy [J/m3]
     real(RP) :: ENGK(KA,IA,JA) ! kinetic   energy [J/m3]
     real(RP) :: ENGI(KA,IA,JA) ! internal  energy [J/m3]
 
     real(RP) :: RHOQ(KA,IA,JA)
 
     real(RP) :: total ! dummy
     integer  :: i, j, k, iq
     !---------------------------------------------------------------------------
 
     if ( statistics_checktotal ) then
 
        call stat_total( total, dens(:,:,:), var_name(i_dens) )
        call stat_total( total, momz(:,:,:), var_name(i_momz) )
        call stat_total( total, momx(:,:,:), var_name(i_momx) )
        call stat_total( total, momy(:,:,:), var_name(i_momy) )
        call stat_total( total, rhot(:,:,:), var_name(i_rhot) )
 
        do iq = 1, qa
           rhoq(:,:,:) = dens_av(:,:,:) * qtrc_av(:,:,:,iq)
 
           call stat_total( total, rhoq(:,:,:), tracer_name(iq) )
        enddo
 
        call thermodyn_qd( qdry(:,:,:),   & ! [OUT]
                           qtrc_av(:,:,:,:), & ! [IN]
                           tracer_mass(:) ) ! [IN]
 
        call thermodyn_temp_pres( temp(:,:,:),   & ! [OUT]
                                  pres(:,:,:),   & ! [OUT]
                                  dens_av(:,:,:),   & ! [IN]
                                  rhot_av(:,:,:),   & ! [IN]
                                  qtrc_av(:,:,:,:), & ! [IN]
                                  tracer_cv(:),  & ! [IN]
                                  tracer_r(:),   & ! [IN]
                                  tracer_mass(:) ) ! [IN]
 
        rhoq(ks:ke,is:ie,js:je) = dens_av(ks:ke,is:ie,js:je) * qdry(ks:ke,is:ie,js:je)
 
        call stat_total( total, rhoq(:,:,:), 'QDRY' )
 
        rhoq(ks:ke,is:ie,js:je) = dens_av(ks:ke,is:ie,js:je) * ( 1.0_rp - qdry(ks:ke,is:ie,js:je) ) ! Qtotal
 
        call stat_total( total, rhoq(:,:,:), 'QTOT' )
 
        !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
        do j = js, je
        do i = is, ie
        do k = ks, ke
           w(k,i,j) = 0.5_rp * ( momz_av(k-1,i,j)+momz_av(k,i,j) ) / dens_av(k,i,j)
           u(k,i,j) = 0.5_rp * ( momx_av(k,i-1,j)+momx_av(k,i,j) ) / dens_av(k,i,j)
           v(k,i,j) = 0.5_rp * ( momy_av(k,i,j-1)+momy_av(k,i,j) ) / dens_av(k,i,j)
 
           engp(k,i,j) = dens_av(k,i,j) * grav * real_cz(k,i,j)
 
           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 )
 
           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 ( i_qv > 0 ) then
              engi(k,i,j) = engi(k,i,j) + dens(k,i,j) * qtrc(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]
           end if
 
           do iq = qis, qie
              engi(k,i,j) = engi(k,i,j) - dens_av(k,i,j) * qtrc_av(k,i,j,iq) * lhf ! Latent Heat [ice->liquid]
           enddo
 
           engt(k,i,j) = engp(k,i,j) + engk(k,i,j) + engi(k,i,j)
        enddo
        enddo
        enddo
 
        call stat_total( total, engp(:,:,:), 'ENGP' )
        call stat_total( total, engk(:,:,:), 'ENGK' )
        call stat_total( total, engi(:,:,:), 'ENGI' )
        call stat_total( total, engt(:,:,:), 'ENGT' )
 
     endif
 
     return

Here is the caller graph for this function:

◆ atmos_vars_diagnostics()

subroutine, public mod_atmos_vars::atmos_vars_diagnostics ( )

Calc diagnostic variables.

Definition at line 2210 of file mod_atmos_vars.f90.

References scale_atmos_diagnostic::atmos_diagnostic_get(), dens_av, momx_av, momy_av, momz_av, n2, phyd, pott, pres, qtrc_av, rhot_av, temp, u, v, and w.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), and mod_rm_driver::resume_state().

     use scale_atmos_diagnostic, only: &
        atmos_diagnostic_get
     implicit none
 
     call atmos_diagnostic_get( pott, temp, pres, phyd, w, u, v, n2,                 & ! (out)
                                dens_av, momz_av, momx_av, momy_av, rhot_av, qtrc_av ) ! (in)
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_vars_monitor()

subroutine, public mod_atmos_vars::atmos_vars_monitor ( )

monitor output

Definition at line 2223 of file mod_atmos_vars.f90.

References mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_up, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_up, mod_atmos_phy_rd_vars::atmos_phy_rd_toaflx_lw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_toaflx_lw_up, mod_atmos_phy_rd_vars::atmos_phy_rd_toaflx_sw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_toaflx_sw_up, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, scale_const::const_cvdry, scale_const::const_grav, dens_av, scale_grid::grid_rfdx, scale_grid::grid_rfdy, scale_gridtrans::gtrans_mapf, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_atmos_hydrometeor::i_qv, scale_index::i_rhot, scale_gridtrans::i_uy, scale_gridtrans::i_xv, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, scale_atmos_hydrometeor::lhf, scale_atmos_hydrometeor::lhv, momx_av, momy_av, momz_av, scale_process::prc_myrank, scale_tracer::qa, scale_atmos_hydrometeor::qie, scale_atmos_hydrometeor::qis, qtrc_av, scale_grid_real::real_cz, rhot_av, scale_rm_statistics::stat_detail(), scale_rm_statistics::statistics_checktotal, temp, scale_time::time_dtsec_atmos_dyn, scale_tracer::tracer_cv, scale_tracer::tracer_desc, scale_tracer::tracer_mass, scale_tracer::tracer_name, scale_tracer::tracer_unit, u, v, and w.

Referenced by mod_atmos_driver::atmos_driver(), and mod_atmos_driver::atmos_driver_resume2().

     use scale_process, only: &
        prc_myrank
     use scale_const, only: &
        grav  => const_grav,  &
        cvdry => const_cvdry
     use scale_grid, only: &
        rfdx => grid_rfdx, &
        rfdy => grid_rfdy
     use scale_grid_real, only: &
        real_cz
     use scale_gridtrans, only: &
        mapf => gtrans_mapf, &
        i_uy, &
        i_xv
     use scale_comm, only: &
        comm_vars8, &
        comm_wait
     use scale_rm_statistics, only: &
        statistics_checktotal, &
        stat_total,            &
        stat_detail
     use scale_monitor, only: &
        monit_put, &
        monit_in
     use scale_time, only: &
        time_dtsec_atmos_dyn
     use scale_atmos_thermodyn, only: &
        thermodyn_qd        => atmos_thermodyn_qd,        &
        thermodyn_temp_pres => atmos_thermodyn_temp_pres
     use scale_atmos_hydrometeor, only: &
        i_qv, &
        qis,  &
        qie,  &
        lhv,  &
        lhf
     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) :: QDRY(KA,IA,JA) ! dry air         [kg/kg]
     real(RP) :: RHOQ(KA,IA,JA) ! DENS * tracer   [kg/m3]
     real(RP) :: PRCP(IA,JA)    ! rain + snow     [kg/m2/s]
 
     real(RP) :: ENGT(KA,IA,JA) ! total     energy [J/m3]
     real(RP) :: ENGP(KA,IA,JA) ! potential energy [J/m3]
     real(RP) :: ENGK(KA,IA,JA) ! kinetic   energy [J/m3]
     real(RP) :: ENGI(KA,IA,JA) ! internal  energy [J/m3]
 
     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) :: TOAFLX_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 monit_in( dens_av(:,:,:), var_name(i_dens), var_desc(i_dens), var_unit(i_dens), ndim=3, isflux=.false. )
     call monit_in( momz_av(:,:,:), var_name(i_momz), var_desc(i_momz), var_unit(i_momz), ndim=3, isflux=.false. )
     call monit_in( momx_av(:,:,:), var_name(i_momx), var_desc(i_momx), var_unit(i_momx), ndim=3, isflux=.false. )
     call monit_in( momy_av(:,:,:), var_name(i_momy), var_desc(i_momy), var_unit(i_momy), ndim=3, isflux=.false. )
     call monit_in( rhot_av(:,:,:), var_name(i_rhot), var_desc(i_rhot), var_unit(i_rhot), ndim=3, isflux=.false. )
 
     !##### 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 monit_in( rhoq(:,:,:), tracer_name(iq), tracer_desc(iq), tracer_unit(iq), ndim=3, isflux=.false. )
     enddo
 
     ! total dry airmass
 
     call thermodyn_qd( qdry(:,:,:),   & ! [OUT]
                        qtrc_av(:,:,:,:), & ! [IN]
                        tracer_mass(:) ) ! [IN]
 
     !$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) * qdry(k,i,j)
     enddo
     enddo
     enddo
     call monit_put( ad_monit_id(i_qdry), rhoq(:,:,:) )
 
     ! total vapor,liquid,solid tracers
     !$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) * ( 1.0_rp - qdry(k,i,j) )
     enddo
     enddo
     enddo
     call monit_put( ad_monit_id(i_qtot), rhoq(:,:,:) )
 
     ! total evapolation
     if ( i_qv > 0 ) then
        call monit_put( ad_monit_id(i_evap), sflx_qtrc(:,:,i_qv) )
     endif
 
     ! total precipitation
 !OCL XFILL
     do j = js, je
     do i = is, ie
        prcp(i,j) = sflx_rain_cp(i,j) &
                  + sflx_rain_mp(i,j) + sflx_snow_mp(i,j)
     enddo
     enddo
     call monit_put( ad_monit_id(i_prcp), prcp(:,:) )
 
     !##### Energy Budget #####
 
 
     !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
     !$omp private(iq) &
     !$omp shared(JS,JE,IS,IE,KS,KE,ENGP,DENS_av,GRAV,REAL_CZ,W,U,V,ENGI,ENGK,QDRY,TEMP,CVdry,QA) &
     !$omp shared(TRACER_CV,QTRC_av,I_QV,LHV,QIS,QIE,LHF)
     do j = js, je
     do i = is, ie
     do k = ks, ke
        engp(k,i,j) = dens_av(k,i,j) * grav * real_cz(k,i,j)
 
        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 )
 
        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 ( i_qv > 0 ) then
           engi(k,i,j) = engi(k,i,j) + dens_av(k,i,j) * qtrc_av(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]
        end if
 
        do iq = qis, qie
           engi(k,i,j) = engi(k,i,j) - dens_av(k,i,j) * qtrc_av(k,i,j,iq) * lhf ! Latent Heat [ice->liquid]
        enddo
     enddo
     enddo
     enddo
 
     !$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
        engt(k,i,j) = engp(k,i,j) + engk(k,i,j) + engi(k,i,j)
     enddo
     enddo
     enddo
 
 !OCL XFILL
     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) )
 
        toaflx_rd_net(i,j) = ( toaflx_lw_up(i,j) - toaflx_lw_dn(i,j) ) &
                           + ( toaflx_sw_up(i,j) - toaflx_sw_dn(i,j) )
     enddo
     enddo
 
     !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
 !OCL XFILL
     do j = js, je
     do i = is, ie
        engflxt(i,j) = sflx_sh(i,j) + sflx_lh(i,j) &
                     + sflx_rd_net(i,j) - toaflx_rd_net(i,j)
     enddo
     enddo
 
     call monit_put( ad_monit_id(i_engp), engp(:,:,:) )
     call monit_put( ad_monit_id(i_engk), engk(:,:,:) )
     call monit_put( ad_monit_id(i_engi), engi(:,:,:) )
     call monit_put( ad_monit_id(i_engt), engt(:,:,:) )
 
     call monit_put( ad_monit_id(i_engflxt), engflxt(:,:) )
 
 
     call monit_put( ad_monit_id(i_engsfc_sh), sflx_sh(:,:) )
     call monit_put( ad_monit_id(i_engsfc_lh), sflx_lh(:,:) )
     call monit_put( ad_monit_id(i_engsfc_rd), sflx_rd_net(:,:) )
     call monit_put( ad_monit_id(i_engtoa_rd), toaflx_rd_net(:,:) )
 
     call monit_put( ad_monit_id(i_engsfc_lw_up), sflx_lw_up(:,:) )
     call monit_put( ad_monit_id(i_engsfc_lw_dn), sflx_lw_dn(:,:) )
     call monit_put( ad_monit_id(i_engsfc_sw_up), sflx_sw_up(:,:) )
     call monit_put( ad_monit_id(i_engsfc_sw_dn), sflx_sw_dn(:,:) )
 
     call monit_put( ad_monit_id(i_engtoa_lw_up), toaflx_lw_up(:,:) )
     call monit_put( ad_monit_id(i_engtoa_lw_dn), toaflx_lw_dn(:,:) )
     call monit_put( ad_monit_id(i_engtoa_sw_up), toaflx_sw_up(:,:) )
     call monit_put( ad_monit_id(i_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 stat_detail( work(:,:,:,:), wname(:) )
     endif
 
     if ( atmos_vars_checkcfl > 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 ( cflmax > atmos_vars_checkcfl ) then
           if( io_l ) write(io_fid_log,*) "*** [ATMOS_vars_monitor] Courant number exceeded the upper limit. : ", cflmax
                      write(*,*)          "*** [ATMOS_vars_monitor] Courant number exceeded the upper limit. : ", cflmax, &
                                          ", rank = ", prc_myrank
 
           wname(1) = "Courant num. Z"
           wname(2) = "Courant num. X"
           wname(3) = "Courant num. Y"
 
           call stat_detail( work(:,:,:,:), wname(:), supress_globalcomm=.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 2502 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_time, only: &
        time_gettimelabel
     use scale_fileio, only: &
        fileio_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_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: &
        nowsec => time_nowsec
 #endif
     implicit none
 
     character(len=19)     :: timelabel
     character(len=H_LONG) :: basename
     !---------------------------------------------------------------------------
 
 #ifdef SDM
     if( sd_rest_flg_out ) then
        if( io_l ) write(io_fid_log,*) '*** Output random number for SDM ***'
        call atmos_phy_mp_sdm_restart_create(nowsec)
     endif
 #endif
 
     if ( atmos_restart_out_basename /= '' ) then
 
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '*** 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
 
        if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)
 
        call fileio_create( restart_fid,                                               & ! [OUT]
                            basename, atmos_restart_out_title, atmos_restart_out_dtype ) ! [IN]
 
        allocate( var_id(vmax+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 ) 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 2586 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_fileio, only: &
        fileio_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_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 fileio_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 ) 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 2647 of file mod_atmos_vars.f90.

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

     use scale_fileio, only: &
        fileio_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_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
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '*** Close restart file (ATMOS) ***'
 
        call fileio_close( restart_fid ) ! [IN]
 
        restart_fid = -1
 
        if ( allocated(var_id) ) deallocate( var_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 ) 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 2714 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_fileio, only: &
        fileio_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_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 fileio_def_var( restart_fid, var_id(i_dens), var_name(i_dens), var_desc(i_dens), & ! [IN]
                             var_unit(i_dens), 'ZXY',  atmos_restart_out_dtype                ) ! [IN]
        call fileio_def_var( restart_fid, var_id(i_momz), var_name(i_momz), var_desc(i_momz), & ! [IN]
                             var_unit(i_momz), 'ZHXY', atmos_restart_out_dtype                ) ! [IN]
        call fileio_def_var( restart_fid, var_id(i_momx), var_name(i_momx), var_desc(i_momx), & ! [IN]
                             var_unit(i_momx), 'ZXHY', atmos_restart_out_dtype                ) ! [IN]
        call fileio_def_var( restart_fid, var_id(i_momy), var_name(i_momy), var_desc(i_momy), & ! [IN]
                             var_unit(i_momy), 'ZXYH', atmos_restart_out_dtype                ) ! [IN]
        call fileio_def_var( restart_fid, var_id(i_rhot), var_name(i_rhot), var_desc(i_rhot), & ! [IN]
                             var_unit(i_rhot), 'ZXY',  atmos_restart_out_dtype                ) ! [IN]
 
        do iq = 1, qa
           call fileio_def_var( restart_fid, var_id(vmax+iq), tracer_name(iq), tracer_desc(iq), & ! [IN]
                                tracer_unit(iq), 'ZXY',  atmos_restart_out_dtype       ) ! [IN]
        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 ) 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 2792 of file mod_atmos_vars.f90.

Referenced by mod_admin_restart::admin_restart_write().

     use scale_fileio, only: &
        fileio_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_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 fileio_write_var( restart_fid, var_id(i_dens), dens(:,:,:), var_name(i_dens), 'ZXY'  ) ! [IN]
        call fileio_write_var( restart_fid, var_id(i_momz), momz(:,:,:), var_name(i_momz), 'ZHXY' ) ! [IN]
        call fileio_write_var( restart_fid, var_id(i_momx), momx(:,:,:), var_name(i_momx), 'ZXHY' ) ! [IN]
        call fileio_write_var( restart_fid, var_id(i_momy), momy(:,:,:), var_name(i_momy), 'ZXYH' ) ! [IN]
        call fileio_write_var( restart_fid, var_id(i_rhot), rhot(:,:,:), var_name(i_rhot), 'ZXY'  ) ! [IN]
 
        do iq = 1, qa
           call fileio_write_var( restart_fid, var_id(vmax+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 ) 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:

Variable Documentation

◆ atmos_restart_output

logical, public mod_atmos_vars::atmos_restart_output = .false.

Output restart file?

Definition at line 58 of file mod_atmos_vars.f90.

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

58 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 60 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().

60 character(len=H_LONG), public :: ATMOS_RESTART_IN_BASENAME = ''

◆ 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 61 of file mod_atmos_vars.f90.

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

61 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 62 of file mod_atmos_vars.f90.

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

62 character(len=H_LONG), public :: ATMOS_RESTART_OUT_BASENAME = ''

◆ 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 63 of file mod_atmos_vars.f90.

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

63 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 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_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 65 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().

65 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 67 of file mod_atmos_vars.f90.

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

67 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 68 of file mod_atmos_vars.f90.

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

68 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 69 of file mod_atmos_vars.f90.

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

69 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 72 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), and mod_mkinit::rect_setup().

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

◆ momz

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

Definition at line 73 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

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

◆ momx

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

Definition at line 74 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

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

◆ momy

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

Definition at line 75 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

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

◆ rhot

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

Definition at line 76 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), and mod_mkinit::rect_setup().

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

◆ qtrc

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

Definition at line 77 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), 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::interporation_fact(), mod_mkinit::mkinit(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), mod_mkinit::rect_setup(), and mod_mkinit::tke_setup().

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

◆ dens_avw

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

Definition at line 79 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

79 real(RP), public, target, allocatable :: DENS_avw(:,:,:)

◆ momz_avw

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

Definition at line 80 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

80 real(RP), public, target, allocatable :: MOMZ_avw(:,:,:)

◆ momx_avw

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

Definition at line 81 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

81 real(RP), public, target, allocatable :: MOMX_avw(:,:,:)

◆ momy_avw

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

Definition at line 82 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

◆ rhot_avw

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

Definition at line 83 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

◆ qtrc_avw

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

Definition at line 84 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

◆ dens_av

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

Definition at line 86 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_history_setpres(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

86 real(RP), public, pointer :: DENS_av(:,:,:)

◆ momz_av

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

Definition at line 87 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

87 real(RP), public, pointer :: MOMZ_av(:,:,:)

◆ momx_av

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

Definition at line 88 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

88 real(RP), public, pointer :: MOMX_av(:,:,:)

◆ momy_av

real(rp), dimension(:,:,:), pointer, public mod_atmos_vars::momy_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(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

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

◆ rhot_av

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

Definition at line 90 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

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

◆ qtrc_av

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

Definition at line 91 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_setup(), and atmos_vars_total().

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

◆ dens_tp

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

Definition at line 94 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), and atmos_vars_setup().

94 real(RP), public, allocatable :: DENS_tp(:,:,:)

◆ momz_tp

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

Definition at line 95 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

95 real(RP), public, allocatable :: MOMZ_tp(:,:,:)

◆ momx_tp

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

Definition at line 96 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

96 real(RP), public, allocatable :: MOMX_tp(:,:,:)

◆ momy_tp

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

Definition at line 97 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

97 real(RP), public, allocatable :: MOMY_tp(:,:,:)

◆ rhot_tp

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

Definition at line 98 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

98 real(RP), public, allocatable :: RHOT_tp(:,:,:)

◆ rhoq_tp

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

Definition at line 99 of file mod_atmos_vars.f90.

Referenced by mod_atmos_driver::atmos_driver(), mod_atmos_driver::atmos_driver_resume1(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

99 real(RP), public, allocatable :: RHOQ_tp(:,:,:,:)

◆ pott

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

Definition at line 102 of file mod_atmos_vars.f90.

Referenced by atmos_vars_diagnostics(), atmos_vars_history(), and atmos_vars_setup().

102 real(RP), public, allocatable :: POTT(:,:,:) ! potential temperature [K]

◆ temp

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

Definition at line 103 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_mp_driver::atmos_phy_mp_driver(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_driver::atmos_surface_set(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_setup(), and atmos_vars_total().

103 real(RP), public, allocatable :: TEMP(:,:,:) ! temperature [K]

◆ pres

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

Definition at line 104 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_driver::atmos_surface_set(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_history_setpres(), atmos_vars_setup(), and atmos_vars_total().

104 real(RP), public, allocatable :: PRES(:,:,:) ! pressure [Pa=J/m3]

◆ phyd

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

Definition at line 105 of file mod_atmos_vars.f90.

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

105 real(RP), public, allocatable :: PHYD(:,:,:) ! hydrostatic pressure [Pa=J/m3]

◆ w

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

Definition at line 106 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_driver::atmos_surface_set(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_setup(), and atmos_vars_total().

106 real(RP), public, allocatable :: W (:,:,:) ! velocity w [m/s]

◆ u

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

Definition at line 107 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_driver::atmos_surface_set(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_setup(), and atmos_vars_total().

107 real(RP), public, allocatable :: U (:,:,:) ! velocity u [m/s]

◆ v

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

Definition at line 108 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_sf_driver::atmos_phy_sf_driver(), mod_atmos_driver::atmos_surface_set(), atmos_vars_diagnostics(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_setup(), and atmos_vars_total().

108 real(RP), public, allocatable :: V (:,:,:) ! velocity v [m/s]

◆ n2

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

Definition at line 109 of file mod_atmos_vars.f90.

Referenced by mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_history(), and atmos_vars_setup().

109 real(RP), public, allocatable :: N2 (:,:,:) ! squared Brunt-Vaisala frequency [/s2]

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

◆ atmos_vars_monitor()

◆ atmos_vars_restart_create()

◆ atmos_vars_restart_enddef()

◆ atmos_vars_restart_close()

◆ atmos_vars_restart_def_var()

◆ atmos_vars_restart_write()

Variable Documentation

◆ atmos_restart_output

◆ atmos_restart_in_basename

◆ atmos_restart_in_postfix_timelabel

◆ atmos_restart_out_basename

◆ 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

◆ dens_tp

◆ momz_tp

◆ momx_tp

◆ momy_tp

◆ rhot_tp

◆ rhoq_tp

◆ pott

◆ temp

◆ pres

◆ phyd

◆ w

◆ u

◆ v

◆ n2