mod_atmos_vars Module Reference

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_read
	Read restart of atmospheric variables. More...
subroutine, public	atmos_vars_restart_write
	Write restart of atmospheric variables. 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_var
	Write restart of atmospheric variables. More...

Variables
logical, public	atmos_restart_output = .false.
	output restart file? More...
logical, public	atmos_restart_check = .false.
	check value consistency? More...
character(len=h_long), public	atmos_restart_in_basename = ''
	basename of the restart file More...
character(len=h_long), public	atmos_restart_out_basename = ''
	basename of the output file More...
character(len=h_mid), public	atmos_restart_out_title = 'ATMOS restart'
	title of the output file More...
character(len=h_mid), public	atmos_restart_out_dtype = 'DEFAULT'
	REAL4 or REAL8. More...
logical, public	atmos_restart_in_allowmissingq = .false.
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	temp
real(rp), dimension(:,:,:), allocatable, public	pres
real(rp), dimension(:,:,:), allocatable, public	w
real(rp), dimension(:,:,:), allocatable, public	u
real(rp), dimension(:,:,:), allocatable, public	v
real(rp), dimension(:,:,:), allocatable, public	pott

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 restart file
  ATMOS_RESTART_IN_ALLOWMISSINGQ	logical	.false.
  ATMOS_RESTART_OUTPUT	logical	.false.	output restart file?
  ATMOS_RESTART_OUT_BASENAME	character(len=H_LONG)	''	basename of the output file
  ATMOS_RESTART_OUT_TITLE	character(len=H_MID)	'ATMOS restart'	title of the output file
  ATMOS_RESTART_OUT_DTYPE	character(len=H_MID)	'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
AQ_NAME(iq)	AQ_DESC(iq)	AQ_UNIT(iq)	QTRC
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
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
RH	relative humidity(liq)	%	RHL
RHA	relative humidity(liq+ice)	%	RHA
RHI	relative humidity(ice)	%	RHI
RTOT	Total gas constant	J/kg/K	RTOT
T	temperature	K	TEMP
TKE_RS	resolved scale TKE	m2/s2	TKE_RS
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 247 of file mod_atmos_vars.f90.

References scale_tracer::aq_desc, scale_tracer::aq_name, scale_tracer::aq_unit, 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_allowmissingq, atmos_restart_in_basename, atmos_restart_out_basename, atmos_restart_out_dtype, 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_l, scale_stdio::io_lnml, 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(), pott, scale_process::prc_mpistop(), pres, scale_tracer::qa, qtrc, qtrc_av, qtrc_avw, rhoq_tp, rhot, rhot_av, rhot_avw, rhot_tp, temp, 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_allowmissingq, &
       atmos_restart_output,           &
       atmos_restart_out_basename,     &
       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

    logical :: zinterp ! dummy
    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,qa) )

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

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

    allocate( temp(ka,ia,ja) )
    allocate( pres(ka,ia,ja) )
    allocate( w(ka,ia,ja) )
    allocate( u(ka,ia,ja) )
    allocate( v(ka,ia,ja) )
    allocate( pott(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_lnml ) write(io_fid_log,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,A15,A,A32,3(A))') &
               '***       |','VARNAME        ','|', 'DESCRIPTION                     ','[', 'UNIT            ',']'
    do iv = 1, vmax
       if( io_l ) write(io_fid_log,'(1x,A,i3,A,A15,A,A32,3(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,A15,A,A32,3(A))') &
                  '*** NO.',5+iq,'|',aq_name(iq),'|', aq_desc(iq),'[', aq_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?  : ', trim(atmos_restart_in_basename)
    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? : ', trim(atmos_restart_out_basename)
    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(qa))

    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), zinterp, var_name(i_dens), var_desc(i_dens), var_unit(i_dens), ndim=3 )
    call hist_reg( var_hist_id(i_momz), zinterp, var_name(i_momz), var_desc(i_momz), var_unit(i_momz), ndim=3, zdim='half' )
    call hist_reg( var_hist_id(i_momx), zinterp, var_name(i_momx), var_desc(i_momx), var_unit(i_momx), ndim=3, xdim='half' )
    call hist_reg( var_hist_id(i_momy), zinterp, var_name(i_momy), var_desc(i_momy), var_unit(i_momy), ndim=3, ydim='half' )
    call hist_reg( var_hist_id(i_rhot), zinterp, 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), zinterp, aq_name(iq), aq_desc(iq), aq_unit(iq), ndim=3 )
    enddo

    call hist_reg( ad_hist_id(i_w   ) , zinterp, 'W',     'velocity w',             'm/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_u   ) , zinterp, 'U',     'velocity u',             'm/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_v   ) , zinterp, 'V',     'velocity v',             'm/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_pott) , zinterp, 'PT',    'potential temp.',        'K',      ndim=3 )

    call hist_reg( ad_hist_id(i_qdry) , zinterp, 'QDRY',  'dry air',                'kg/kg',  ndim=3 )
    call hist_reg( ad_hist_id(i_qtot) , zinterp, 'QTOT',  'total water',            'kg/kg',  ndim=3 )
    call hist_reg( ad_hist_id(i_qhyd) , zinterp, 'QHYD',  'total hydrometeors',     'kg/kg',  ndim=3 )
    call hist_reg( ad_hist_id(i_qliq) , zinterp, 'QLIQ',  'total liquid water',     'kg/kg',  ndim=3 )
    call hist_reg( ad_hist_id(i_qice) , zinterp, 'QICE',  'total ice water',        'kg/kg',  ndim=3 )

    call hist_reg( ad_hist_id(i_lwp)  , zinterp, 'LWP',   'liquid water path',      'g/m2',   ndim=2 )
    call hist_reg( ad_hist_id(i_iwp)  , zinterp, 'IWP',   'ice water path',         'g/m2',   ndim=2 )
    call hist_reg( ad_hist_id(i_pw )  , zinterp, 'PW',    'precipitable water',     'g/m2',   ndim=2 )

    call hist_reg( ad_hist_id(i_rtot) , zinterp, 'RTOT',  'Total gas constant',     'J/kg/K', ndim=3 )
    call hist_reg( ad_hist_id(i_cptot), zinterp, 'CPTOT', 'Total heat capacity',    'J/kg/K', ndim=3 )
    call hist_reg( ad_hist_id(i_pres) , zinterp, 'PRES',  'pressure',               'Pa',     ndim=3 )
    call hist_reg( ad_hist_id(i_temp) , zinterp, 'T',     'temperature',            'K',      ndim=3 )

    call hist_reg( ad_hist_id(i_potl) , zinterp, 'LWPT',  'liq. potential temp.',   'K',      ndim=3 )
    call hist_reg( ad_hist_id(i_rha)  , zinterp, 'RHA',   'relative humidity(liq+ice)', '%',  ndim=3 )
    call hist_reg( ad_hist_id(i_rhl)  , zinterp, 'RH',    'relative humidity(liq)', '%',      ndim=3 )
    call hist_reg( ad_hist_id(i_rhi)  , zinterp, 'RHI',   'relative humidity(ice)', '%',      ndim=3 )

    call hist_reg( ad_hist_id(i_vor)  , zinterp, 'VOR',   'vertical vorticity',     '1/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_div)  , zinterp, 'DIV',   'divergence',             '1/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_hdiv) , zinterp, 'HDIV',  'horizontal divergence',  '1/s',    ndim=3 )
    call hist_reg( ad_hist_id(i_uabs) , zinterp, 'Uabs',  'absolute velocity',      'm/s',    ndim=3 )

    call hist_reg( ad_hist_id(i_cape) , zinterp, 'CAPE',  'convection avail. pot. energy', 'm2/s2', ndim=2 )
    call hist_reg( ad_hist_id(i_cin)  , zinterp, 'CIN',   'convection inhibition',         'm2/s2', ndim=2 )
    call hist_reg( ad_hist_id(i_lcl)  , zinterp, 'LCL',   'lifted condensation level',     'm',     ndim=2 )
    call hist_reg( ad_hist_id(i_lfc)  , zinterp, 'LFC',   'level of free convection',      'm',     ndim=2 )
    call hist_reg( ad_hist_id(i_lnb)  , zinterp, 'LNB',   'level of neutral buoyancy',     'm',     ndim=2 )

    call hist_reg( ad_hist_id(i_pblh) , zinterp, 'PBLH',  'PBL height',             'm',      ndim=2 )
    call hist_reg( ad_hist_id(i_mse)  , zinterp, 'MSE',   'moist static energy',    'm2/s2',  ndim=3 )

    call hist_reg( ad_hist_id(i_dens_mean), zinterp, 'DENS_MEAN', 'horiz. mean of density',    'kg/m3', ndim=1 )
    call hist_reg( ad_hist_id(i_w_mean),    zinterp, 'W_MEAN',    'horiz. mean of w',           'm/s',  ndim=1 )
    call hist_reg( ad_hist_id(i_u_mean),    zinterp, 'U_MEAN',    'horiz. mean of u',           'm/s',  ndim=1 )
    call hist_reg( ad_hist_id(i_v_mean),    zinterp, 'V_MEAN',    'horiz. mean of v',           'm/s',  ndim=1 )
    call hist_reg( ad_hist_id(i_pott_mean), zinterp, 'PT_MEAN',   'horiz. mean of pot.',        'K',    ndim=1 )
    call hist_reg( ad_hist_id(i_t_mean),    zinterp, 'T_MEAN',    'horiz. mean of t',           'K',    ndim=1 )
    call hist_reg( ad_hist_id(i_qv_mean),   zinterp, 'QV_MEAN',   'horiz. mean of QV',          '1',    ndim=1 )
    call hist_reg( ad_hist_id(i_qhyd_mean), zinterp, 'QHYD_MEAN', 'horiz. mean of QHYD',        '1',    ndim=1 )
    call hist_reg( ad_hist_id(i_qliq_mean), zinterp, 'QLIQ_MEAN', 'horiz. mean of QLIQ',        '1',    ndim=1 )
    call hist_reg( ad_hist_id(i_qice_mean), zinterp, 'QICE_MEAN', 'horiz. mean of QICE',        '1',    ndim=1 )

    call hist_reg( ad_hist_id(i_dens_prim), zinterp, 'DENS_PRIM', 'horiz. deviation of density',    'kg/m3', ndim=3 )
    call hist_reg( ad_hist_id(i_w_prim   ), zinterp, 'W_PRIM',    'horiz. deviation of w',          'm/s',   ndim=3 )
    call hist_reg( ad_hist_id(i_u_prim   ), zinterp, 'U_PRIM',    'horiz. deviation of u',          'm/s',   ndim=3 )
    call hist_reg( ad_hist_id(i_v_prim   ), zinterp, 'V_PRIM',    'horiz. deviation of v',          'm/s',   ndim=3 )
    call hist_reg( ad_hist_id(i_pott_prim), zinterp, 'PT_PRIM',   'horiz. deviation of pot. temp.', 'K',     ndim=3 )
    call hist_reg( ad_hist_id(i_w_prim2  ), zinterp, 'W_PRIM2',   'variance of w',                  'm2/s2', ndim=3 )
    call hist_reg( ad_hist_id(i_pt_w_prim), zinterp, 'PT_W_PRIM', 'resolved scale heat flux',       'W/s',   ndim=3 )
    call hist_reg( ad_hist_id(i_w_prim3  ), zinterp, 'W_PRIM3',   'skewness of w',                  'm3/s3', ndim=3 )
    call hist_reg( ad_hist_id(i_tke_rs   ), zinterp, 'TKE_RS',    'resolved scale TKE',             'm2/s2', ndim=3 )

    call hist_reg( ad_hist_id(i_engt) , zinterp, 'ENGT',  'total energy',           'J/m3',   ndim=3 )
    call hist_reg( ad_hist_id(i_engp) , zinterp, 'ENGP',  'potential energy',       'J/m3',   ndim=3 )
    call hist_reg( ad_hist_id(i_engk) , zinterp, 'ENGK',  'kinetic energy',         'J/m3',   ndim=3 )
    call hist_reg( ad_hist_id(i_engi) , zinterp, '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', 'total     energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engsfc_lh), 'ENGSFC_LH', 'potential energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engsfc_rd), 'ENGSFC_RD', 'kinetic   energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engtoa_rd), 'ENGTOA_RD', 'internal  energy', 'J',  ndim=2, isflux=.true. )

    call monit_reg( ad_monit_id(i_engsfc_lw_up), 'ENGSFC_LW_up', 'total     energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engsfc_lw_dn), 'ENGSFC_LW_dn', 'potential energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engsfc_sw_up), 'ENGSFC_SW_up', 'kinetic   energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engsfc_sw_dn), 'ENGSFC_SW_dn', 'internal  energy', 'J',  ndim=2, isflux=.true. )

    call monit_reg( ad_monit_id(i_engtoa_lw_up), 'ENGTOA_LW_up', 'total     energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engtoa_lw_dn), 'ENGTOA_LW_dn', 'potential energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engtoa_sw_up), 'ENGTOA_SW_up', 'kinetic   energy', 'J',  ndim=2, isflux=.true. )
    call monit_reg( ad_monit_id(i_engtoa_sw_dn), 'ENGTOA_SW_dn', 'internal  energy', 'J',  ndim=2, isflux=.true. )

    if ( ad_hist_id(i_w) > 0 ) then
       ad_prep_sw(i_w) = 1
    endif
    if ( ad_hist_id(i_u) > 0 ) then
       ad_prep_sw(i_u) = 1
    endif
    if ( ad_hist_id(i_v) > 0 ) then
       ad_prep_sw(i_v) = 1
    endif
    if ( ad_hist_id(i_pott) > 0 ) then
       ad_prep_sw(i_pott) = 1
    endif

    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_pres) > 0 ) then
       ad_prep_sw(i_qdry)  = 1
       ad_prep_sw(i_rtot)  = 1
       ad_prep_sw(i_cptot) = 1
       ad_prep_sw(i_pres)  = 1
    endif
    if ( ad_hist_id(i_temp) > 0 ) then
       ad_prep_sw(i_qdry)  = 1
       ad_prep_sw(i_rtot)  = 1
       ad_prep_sw(i_cptot) = 1
       ad_prep_sw(i_pres)  = 1
       ad_prep_sw(i_temp)  = 1
    endif

    if ( ad_hist_id(i_potl) > 0 ) then
       ad_prep_sw(i_pott)  = 1
       ad_prep_sw(i_qdry)  = 1
       ad_prep_sw(i_rtot)  = 1
       ad_prep_sw(i_cptot) = 1
       ad_prep_sw(i_pres)  = 1
       ad_prep_sw(i_temp)  = 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_pres)  = 1
       ad_prep_sw(i_temp)  = 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_u)    = 1
       ad_prep_sw(i_v)    = 1
       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_pott) = 1
       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_temp)  = 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)      = 1
       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)      = 1
       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)      = 1
       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)      = 1
       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)       = 1
       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)         = 1
       ad_prep_sw(i_w_prim)    = 1
       ad_prep_sw(i_pott)      = 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)       = 1
       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)      = 1
       ad_prep_sw(i_u)      = 1
       ad_prep_sw(i_v)      = 1
       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_w) = 1
       ad_prep_sw(i_u) = 1
       ad_prep_sw(i_v) = 1
       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_pres)  = 1
       ad_prep_sw(i_temp)  = 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_w)  = 1
       ad_prep_sw(i_u)  = 1
       ad_prep_sw(i_v)  = 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_pres)  = 1
       ad_prep_sw(i_temp)  = 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 801 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(), atmos_vars_restart_write(), and atmos_vars_restart_write_var().

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

subroutine, public mod_atmos_vars::atmos_vars_restart_read

(

)

Read restart of atmospheric variables.

Definition at line 865 of file mod_atmos_vars.f90.

References scale_tracer::aq_name, mod_atmos_dyn_vars::atmos_dyn_vars_restart_read(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_read(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_read(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_read(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_read(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_read(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_read(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_read(), atmos_restart_in_basename, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, mod_atmos_admin::atmos_use_average, atmos_vars_fillhalo(), atmos_vars_total(), scale_const::const_grav, dens, dens_av, scale_stdio::io_fid_log, scale_stdio::io_l, momx, momx_av, momy, momy_av, momz, momz_av, scale_process::prc_mpistop(), scale_tracer::qa, qtrc, qtrc_av, rhot, and rhot_av.

Referenced by mod_rm_driver::resume_state().

    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       grav  => const_grav
    use scale_fileio, only: &
       fileio_read
    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  :: iq
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** Input restart file (ATMOS) ***'

    if ( atmos_restart_in_basename /= '' ) then
       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(atmos_restart_in_basename)

       call fileio_read( dens(:,:,:),                                         & ! [OUT]
                         atmos_restart_in_basename, var_name(1), 'ZXY', step=1 ) ! [IN]
       call fileio_read( momz(:,:,:),                                         & ! [OUT]
                         atmos_restart_in_basename, var_name(2), 'ZXY', step=1 ) ! [IN]
       call fileio_read( momx(:,:,:),                                         & ! [OUT]
                         atmos_restart_in_basename, var_name(3), 'ZXY', step=1 ) ! [IN]
       call fileio_read( momy(:,:,:),                                         & ! [OUT]
                         atmos_restart_in_basename, var_name(4), 'ZXY', step=1 ) ! [IN]
       call fileio_read( rhot(:,:,:),                                         & ! [OUT]
                         atmos_restart_in_basename, var_name(5), 'ZXY', step=1 ) ! [IN]

       do iq = 1, qa
          call fileio_read( qtrc(:,:,:,iq),                                       & ! [OUT]
                            atmos_restart_in_basename, aq_name(iq), 'ZXY', step=1 ) ! [IN]
       enddo

       call atmos_vars_fillhalo

       call atmos_vars_total
    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_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_restart_write()

subroutine, public mod_atmos_vars::atmos_vars_restart_write

(

)

Write restart of atmospheric variables.

Definition at line 959 of file mod_atmos_vars.f90.

References scale_tracer::aq_desc, scale_tracer::aq_name, scale_tracer::aq_unit, mod_atmos_dyn_vars::atmos_dyn_vars_restart_write(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_write(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_write(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_write(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_write(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_write(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_write(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_write(), atmos_restart_out_basename, atmos_restart_out_dtype, atmos_restart_out_title, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, atmos_vars_fillhalo(), atmos_vars_total(), dens, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, scale_stdio::io_fid_log, scale_stdio::io_l, momx, momy, momz, scale_tracer::qa, qtrc, rhot, scale_time::time_gettimelabel(), and scale_time::time_nowsec.

Referenced by mod_mkinit::mkinit(), and mod_rm_driver::scalerm().

    use scale_time, only: &
       time_gettimelabel
    use scale_fileio, only: &
       fileio_write
    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_out
    use scale_time, only: &
       nowsec => time_nowsec
#endif
    implicit none

    character(len=20)     :: timelabel
    character(len=H_LONG) :: basename

    integer :: iq
    !---------------------------------------------------------------------------

#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_out(nowsec)
    endif
#endif

    if ( atmos_restart_out_basename /= '' ) then

       call time_gettimelabel( timelabel )
       write(basename,'(A,A,A)') trim(atmos_restart_out_basename), '_', trim(timelabel)

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Output restart file (ATMOS) ***'
       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call atmos_vars_fillhalo

       call atmos_vars_total

       call fileio_write( dens(:,:,:), basename,                                        atmos_restart_out_title, & ! [IN]
                          var_name(i_dens), var_desc(i_dens), var_unit(i_dens), 'ZXY',  atmos_restart_out_dtype  ) ! [IN]
       call fileio_write( momz(:,:,:), basename,                                        atmos_restart_out_title, & ! [IN]
                          var_name(i_momz), var_desc(i_momz), var_unit(i_momz), 'ZHXY', atmos_restart_out_dtype  ) ! [IN]
       call fileio_write( momx(:,:,:), basename,                                        atmos_restart_out_title, & ! [IN]
                          var_name(i_momx), var_desc(i_momx), var_unit(i_momx), 'ZXHY', atmos_restart_out_dtype  ) ! [IN]
       call fileio_write( momy(:,:,:), basename,                                        atmos_restart_out_title, & ! [IN]
                          var_name(i_momy), var_desc(i_momy), var_unit(i_momy), 'ZXYH', atmos_restart_out_dtype  ) ! [IN]
       call fileio_write( rhot(:,:,:), basename,                                        atmos_restart_out_title, & ! [IN]
                          var_name(i_rhot), var_desc(i_rhot), var_unit(i_rhot), 'ZXY',  atmos_restart_out_dtype  ) ! [IN]

       do iq = 1, qa
          call fileio_write( qtrc(:,:,:,iq), basename,                         atmos_restart_out_title, & ! [IN]
                             aq_name(iq), aq_desc(iq), aq_unit(iq), 'ZXY',  atmos_restart_out_dtype     ) ! [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:

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

References scale_tracer::aq_name, atmos_restart_check_basename, atmos_restart_check_criterion, dens, scale_grid_index::ie, scale_grid_index::imax, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::jmax, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::kmax, scale_grid_index::ks, momx, momy, momz, scale_process::prc_myrank, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, qtrc, and rhot.

Referenced by mod_rm_driver::scalerm().

    use scale_process, only: &
       prc_myrank
    use gtool_file, only: &
       fileread
    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]

    real(RP) :: restart_atmos(kmax,imax,jmax)

    character(len=H_LONG) :: basename

    logical :: datacheck
    integer :: k, i, j, iq
    !---------------------------------------------------------------------------

    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 fileread( restart_atmos(:,:,:), basename, 'DENS', 1, prc_myrank )
    dens_check(ks:ke,is:ie,js:je) = restart_atmos(1:kmax,1:imax,1:jmax)
    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

    call fileread( restart_atmos(:,:,:), basename, 'MOMZ', 1, prc_myrank )
    momz_check(ks:ke,is:ie,js:je) = restart_atmos(1:kmax,1:imax,1:jmax)
    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

    call fileread( restart_atmos(:,:,:), basename, 'MOMX', 1, prc_myrank )
    momx_check(ks:ke,is:ie,js:je) = restart_atmos(1:kmax,1:imax,1:jmax)
    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

    call fileread( restart_atmos(:,:,:), basename, 'MOMY', 1, prc_myrank )
    momy_check(ks:ke,is:ie,js:je) = restart_atmos(1:kmax,1:imax,1:jmax)
    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

    call fileread( restart_atmos(:,:,:), basename, 'RHOT', 1, prc_myrank )
    rhot_check(ks:ke,is:ie,js:je) = restart_atmos(1:kmax,1:imax,1:jmax)
    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
       call fileread( restart_atmos(:,:,:), basename, aq_name(iq), 1, prc_myrank )
       qtrc_check(ks:ke,is:ie,js:je,iq) = restart_atmos(1:kmax,1:imax,1:jmax)
       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, aq_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 1187 of file mod_atmos_vars.f90.

References scale_atmos_thermodyn::aq_cp, scale_atmos_thermodyn::aq_cv, scale_tracer::aq_desc, scale_tracer::aq_name, scale_tracer::aq_unit, scale_atmos_adiabat::atmos_adiabat_cape(), 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, scale_comm::comm_horizontal_mean(), scale_const::const_cpdry, scale_const::const_cvdry, scale_const::const_grav, scale_const::const_lhf, scale_const::const_lhv, scale_const::const_pre00, scale_const::const_rdry, scale_const::const_rvap, dens, scale_grid::grid_rcdx, scale_grid::grid_rcdy, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_tracer::i_qc, scale_tracer::i_qv, scale_index::i_rhot, scale_grid_index::ia, scale_grid_index::ie, scale_grid_index::ieb, scale_grid_index::is, scale_grid_index::isb, scale_grid_index::ja, scale_grid_index::je, scale_grid_index::jeb, scale_grid_index::js, scale_grid_index::jsb, scale_grid_index::ke, scale_grid_index::ks, momx, momy, momz, pott, pres, scale_tracer::qa, scale_tracer::qie, scale_tracer::qis, scale_tracer::qqe, scale_tracer::qqs, qtrc, scale_tracer::qwe, scale_tracer::qws, scale_grid_real::real_cz, scale_grid_real::real_fz, rhot, temp, u, v, and w.

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, &
       lhv   => const_lhv,   &
       lhf   => const_lhf,   &
       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,      &
       thermodyn_templhv => atmos_thermodyn_templhv, &
       cpw => aq_cp,                                 &
       cvw => aq_cv
    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) :: 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) :: lhvap (ka,ia,ja) ! latent heat for vaporization [m2/s2]

    real(RP) :: prec  (ia,ja)    ! surface precipitation 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), aq_name(iq), aq_desc(iq), aq_unit(iq) )
    enddo

    ! prepare and history output of diagnostic variables

!    if ( AD_PREP_sw(I_W) > 0 ) then
!       !$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(k-1,i,j)+MOMZ(k,i,j) ) / DENS(k,i,j)
!       enddo
!       enddo
!       enddo
!    endif
!
!    if ( AD_PREP_sw(I_U) > 0 ) then
!       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
!       do j = JS, JE
!       do i = IS, IE
!       do k = KS, KE
!          U(k,i,j) = 0.5_RP * ( MOMX(k,i-1,j)+MOMX(k,i,j) ) / DENS(k,i,j)
!       enddo
!       enddo
!       enddo
!    endif
!
!    if ( AD_PREP_sw(I_V) > 0 ) then
!       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
!       do j = JS, JE
!       do i = IS, IE
!       do k = KS, KE
!          V(k,i,j) = 0.5_RP * ( MOMY(k,i,j-1)+MOMY(k,i,j) ) / DENS(k,i,j)
!       enddo
!       enddo
!       enddo
!    endif
!
!    if ( AD_PREP_sw(I_POTT) > 0 ) then
!       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
!       do j = JS, JE
!       do i = IS, IE
!       do k = KS, KE
!          POTT(k,i,j) = RHOT(k,i,j) / DENS(k,i,j)
!       enddo
!       enddo
!       enddo
!    endif

    if ( ad_prep_sw(i_qdry) > 0 ) then
       call thermodyn_qd( qdry(:,:,:),  & ! [OUT]
                          qtrc(:,:,:,:) ) ! [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 = qws, qwe
         qhyd(:,:,:) = qhyd(:,:,:) + qtrc(:,:,:,iq)
       enddo
       do iq = qis, qie
         qhyd(:,:,:) = qhyd(:,:,:) + qtrc(:,:,:,iq)
       enddo
    endif

    if ( ad_prep_sw(i_qliq) > 0 ) then
!OCL XFILL
       qliq(:,:,:) = 0.0_rp

       do iq = qws, qwe
         qliq(:,:,:) = qliq(:,:,:) + qtrc(:,:,:,iq)
       enddo
    endif

    if ( ad_prep_sw(i_qice) > 0 ) then
!OCL XFILL
       qice(:,:,:) = 0.0_rp

       do iq = qis, qie
         qice(:,:,:) = qice(:,:,:) + qtrc(:,:,:,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(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(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(k,i,j,i_qv) * dens(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
          rtot(k,i,j) = rdry * qdry(k,i,j) + rvap * qtrc(k,i,j,i_qv)
       enddo
       enddo
       enddo
    endif

    if ( ad_prep_sw(i_cptot) > 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
          cptot(k,i,j) = cpdry * qdry(k,i,j)
       enddo
       enddo
       enddo

       do iq = qqs, qqe
       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
       do j  = jsb, jeb
       do i  = isb, ieb
       do k  = ks, ke
          cptot(k,i,j) = cptot(k,i,j) + qtrc(k,i,j,iq) * cpw(iq)
       enddo
       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) / ( cptot(k,i,j) - rtot(k,i,j) )
       enddo
       enddo
       enddo
    endif

!    if ( AD_PREP_sw(I_PRES) > 0 ) then
!       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
!       do j = JS, JE
!       do i = IS, IE
!       do k = KS, KE
!          PRES(k,i,j) = P00 * ( RHOT(k,i,j) * RTOT(k,i,j) / P00 )**CPovCV(k,i,j)
!       enddo
!       enddo
!       enddo
!    endif
!
!    if ( AD_PREP_sw(I_TEMP) > 0 ) then
!       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
!       do j = JS, JE
!       do i = IS, IE
!       do k = KS, KE
!          TEMP(k,i,j) = PRES(k,i,j) / ( DENS(k,i,j) * RTOT(k,i,j) )
!       enddo
!       enddo
!       enddo
!    endif

    if ( ad_prep_sw(i_potl) > 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
          potl(k,i,j) = pott(k,i,j) &
                      - lhv / 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(:,:,:)  ) ! [IN]
    end if

    if ( ad_hist_id(i_rha) > 0 ) then
       call saturation_psat_all( 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
          rha(k,i,j) = dens(k,i,j) * qtrc(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(k,i,j) * qtrc(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(k,i,j) * qtrc(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(k,i,j)+momx(k,i,j+1)+momx(k,i-1,j)+momx(k,i-1,j+1) ) &
                             / ( dens(k,i,j)+dens(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(k,i,j)+momy(k,i+1,j)+momy(k,i,j-1)+momy(k,i+1,j-1) ) &
                             / ( dens(k,i,j)+dens(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(k,i,j) - momx(k  ,i-1,j  ) ) * rcdx(i) &
                      + ( momy(k,i,j) - momy(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(k,i,j) - momz(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(:,:,:) )
    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(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 = ks, ke
          w_prim(k,i,j) = w(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( w_mean(:), w_prim(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          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 = ks, ke
          u_prim(k,i,j) = u(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( u_mean(:), u_prim(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          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 = ks, ke
          v_prim(k,i,j) = v(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( v_mean(:), v_prim(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          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 = ks, ke
          pott_prim(k,i,j) = temp(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( t_mean(:), pott_prim(:,:,:) )
       do k = ks, ke
          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(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          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 = ks, ke
          rhoq(k,i,j) = qtrc(k,i,j,i_qv) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( qv_mean(:), rhoq(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          rhoq(k,i,j) = qhyd(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( qhyd_mean(:), rhoq(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          rhoq(k,i,j) = qliq(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( qliq_mean(:), rhoq(:,:,:) )
       do k = ks, ke
          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 = ks, ke
          rhoq(k,i,j) = qice(k,i,j) * dens(k,i,j)
       enddo
       enddo
       enddo
       call comm_horizontal_mean( qice_mean(:), rhoq(:,:,:) )
       do k = ks, ke
          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(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(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(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(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
       enddo
       enddo
       enddo

       do iq = qqs, qqe
       !$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(k,i,j) * qtrc(k,i,j,iq) * temp(k,i,j) * cvw(iq)
       enddo
       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
          engi(k,i,j) = engi(k,i,j) &
                      + dens(k,i,j) * qtrc(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]
       enddo
       enddo
       enddo

       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(k,i,j) * qtrc(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( 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(:,:,:),   & ! [IN]
                          temp(:,:,:),   & ! [IN]
                          pres(:,:,:),   & ! [IN]
                          qtrc(:,:,:,:), & ! [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
          pottv(k,i,j) = pott(k,i,j) * ( 1.0_rp + 0.61_rp * qtrc(k,i,j,i_qv) - 1.61 * qtrc(k,i,j,i_qc) )
       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 thermodyn_templhv( lhvap(:,:,:), & ! [OUT]
                               temp(:,:,:)  ) ! [IN]

       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) ) &
                     + lhvap(k,i,j) * qtrc(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_cp(i,j)                     &
                 + sflx_rain_mp(i,j) + sflx_snow_mp(i,j)
    enddo
    enddo
    call hist_in( prec(:,:), 'PREC', 'surface precipitation 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 2134 of file mod_atmos_vars.f90.

References scale_atmos_thermodyn::aq_cv, scale_tracer::aq_name, scale_const::const_cvdry, scale_const::const_grav, scale_const::const_lhf, scale_const::const_lhv, dens, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_tracer::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, momx, momy, momz, scale_tracer::qa, scale_tracer::qie, scale_tracer::qis, scale_tracer::qqe, scale_tracer::qqs, qtrc, scale_grid_real::real_cz, rhot, and scale_rm_statistics::statistics_checktotal.

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

    use scale_const, only: &
       grav  => const_grav,  &
       cvdry => const_cvdry, &
       lhv   => const_lhv,   &
       lhf   => const_lhf
    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, &
       cvw => aq_cv
    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(:,:,:) * qtrc(:,:,:,iq)

          call stat_total( total, rhoq(:,:,:), aq_name(iq) )
       enddo

       call thermodyn_qd( qdry(:,:,:),  & ! [OUT]
                          qtrc(:,:,:,:) ) ! [IN]

       call thermodyn_temp_pres( temp(:,:,:),  & ! [OUT]
                                 pres(:,:,:),  & ! [OUT]
                                 dens(:,:,:),  & ! [IN]
                                 rhot(:,:,:),  & ! [IN]
                                 qtrc(:,:,:,:) ) ! [IN]

       rhoq(ks:ke,is:ie,js:je) = dens(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(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(k-1,i,j)+momz(k,i,j) ) / dens(k,i,j)
          u(k,i,j) = 0.5_rp * ( momx(k,i-1,j)+momx(k,i,j) ) / dens(k,i,j)
          v(k,i,j) = 0.5_rp * ( momy(k,i,j-1)+momy(k,i,j) ) / dens(k,i,j)

          engp(k,i,j) = dens(k,i,j) * grav * real_cz(k,i,j)

          engk(k,i,j) = 0.5_rp * dens(k,i,j) * ( w(k,i,j)**2 &
                                               + u(k,i,j)**2 &
                                               + v(k,i,j)**2 )

          engi(k,i,j) = dens(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
          do iq = qqs, qqe
             engi(k,i,j) = engi(k,i,j) &
                         + dens(k,i,j) * qtrc(k,i,j,iq) * temp(k,i,j) * cvw(iq)
          enddo

          engi(k,i,j) = engi(k,i,j) + dens(k,i,j) * qtrc(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]

          do iq = qis, qie
             engi(k,i,j) = engi(k,i,j) - dens(k,i,j) * qtrc(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 2244 of file mod_atmos_vars.f90.

References dens, scale_grid_index::ia, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::ja, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, momx, momy, momz, pott, pres, qtrc, rhot, 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_comm, only: &
       comm_vars8, &
       comm_wait
    use scale_atmos_thermodyn, only: &
       thermodyn_temp_pres => atmos_thermodyn_temp_pres
    implicit none

    integer  :: k, i, j
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Calc diagnostics'

    call thermodyn_temp_pres( temp(:,:,:),  & ! [OUT]
                              pres(:,:,:),  & ! [OUT]
                              dens(:,:,:),  & ! [IN]
                              rhot(:,:,:),  & ! [IN]
                              qtrc(:,:,:,:) ) ! [IN]

!OCL XFILL
    do j = 1, ja
    do i = 1, ia
    do k = ks+1, ke-1
       w(k,i,j) = 0.5_rp * ( momz(k-1,i,j)+momz(k,i,j) ) / dens(k,i,j)
    enddo
    enddo
    enddo
!OCL XFILL
    do j = 1, ja
    do i = 1, ia
       w(ks,i,j) = 0.5_rp * (               momz(ks,i,j) ) / dens(ks,i,j)
    enddo
    enddo
!OCL XFILL
    do j = 1, ja
    do i = 1, ia
       w(ke,i,j) = 0.5_rp * ( momz(ke-1,i,j)             ) / dens(ke,i,j)
    enddo
    enddo

!OCL XFILL
    do j = 1, ja
    do i = 2, ia
    do k = ks, ke
       u(k,i,j) = 0.5_rp * ( momx(k,i-1,j)+momx(k,i,j) ) / dens(k,i,j)
    enddo
    enddo
    enddo
!OCL XFILL
    do j = 1, ja
    do k = ks, ke
       u(k,1,j) = momx(k,1,j) / dens(k,1,j)
    enddo
    enddo

!OCL XFILL
    do j = 2, ja
    do i = 1, ia
    do k = ks, ke
       v(k,i,j) = 0.5_rp * ( momy(k,i,j-1)+momy(k,i,j) ) / dens(k,i,j)
    enddo
    enddo
    enddo
!OCL XFILL
    do i = 1, ia
    do k = ks, ke
       v(k,i,1) = momy(k,i,1) / dens(k,i,1)
    enddo
    enddo

    !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
    do j  = 1, ja
    do i  = 1, ia
       w(   1:ks-1,i,j) = w(ks,i,j)
       u(   1:ks-1,i,j) = u(ks,i,j)
       v(   1:ks-1,i,j) = v(ks,i,j)
       w(ke+1:ka,  i,j) = w(ke,i,j)
       u(ke+1:ka,  i,j) = u(ke,i,j)
       v(ke+1:ka,  i,j) = v(ke,i,j)
    enddo
    enddo

    call comm_vars8( u(:,:,:), 1 )
    call comm_vars8( v(:,:,:), 2 )
    call comm_wait ( u(:,:,:), 1, .false. )
    call comm_wait ( v(:,:,:), 2, .false. )

!OCL XFILL
    do j = 1, ja
    do i = 1, ia
    do k = ks, ke
       pott(k,i,j) = rhot(k,i,j) / dens(k,i,j)
    enddo
    enddo
    enddo

    return

Here is the caller graph for this function:

atmos_vars_monitor()

subroutine, public mod_atmos_vars::atmos_vars_monitor

(

)

monitor output

Definition at line 2345 of file mod_atmos_vars.f90.

References scale_atmos_thermodyn::aq_cv, scale_tracer::aq_desc, scale_tracer::aq_name, scale_tracer::aq_unit, 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, scale_const::const_lhf, scale_const::const_lhv, dens, 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_tracer::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, momx, momy, momz, scale_process::prc_myrank, pres, scale_tracer::qa, scale_tracer::qie, scale_tracer::qis, scale_tracer::qqe, scale_tracer::qqs, qtrc, scale_grid_real::real_cz, scale_grid_real::real_fz, rhot, scale_rm_statistics::stat_detail(), scale_rm_statistics::statistics_checktotal, temp, scale_time::time_dtsec_atmos_dyn, 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, &
       lhv   => const_lhv,   &
       lhf   => const_lhf
    use scale_grid, only: &
       rfdx => grid_rfdx, &
       rfdy => grid_rfdy
    use scale_grid_real, only: &
       real_cz, &
       real_fz
    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, &
       cvw => aq_cv
    use mod_atmos_phy_mp_vars, only: &
       sflx_rain => atmos_phy_mp_sflx_rain, &
       sflx_snow => 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
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Monitor'

    call monit_in( dens(:,:,:), var_name(i_dens), var_desc(i_dens), var_unit(i_dens), ndim=3, isflux=.false. )
    call monit_in( momz(:,:,:), var_name(i_momz), var_desc(i_momz), var_unit(i_momz), ndim=3, isflux=.false. )
    call monit_in( momx(:,:,:), var_name(i_momx), var_desc(i_momx), var_unit(i_momx), ndim=3, isflux=.false. )
    call monit_in( momy(:,:,:), var_name(i_momy), var_desc(i_momy), var_unit(i_momy), ndim=3, isflux=.false. )
    call monit_in( rhot(:,:,:), 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(k,i,j) * qtrc(k,i,j,iq)
       enddo
       enddo
       enddo

       call monit_in( rhoq(:,:,:), aq_name(iq), aq_desc(iq), aq_unit(iq), ndim=3, isflux=.false. )
    enddo

    ! total dry airmass

    call thermodyn_qd( qdry(:,:,:),  & ! [OUT]
                       qtrc(:,:,:,:) ) ! [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(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(k,i,j) * ( 1.0_rp - qdry(k,i,j) )
    enddo
    enddo
    enddo
    call monit_put( ad_monit_id(i_qtot), rhoq(:,:,:) )

    ! total evapolation
    call monit_put( ad_monit_id(i_evap), sflx_qtrc(:,:,i_qv) )

    ! total precipitation
!OCL XFILL
    do j = js, je
    do i = is, ie
       prcp(i,j) = sflx_rain(i,j) + sflx_snow(i,j)
    enddo
    enddo
    call monit_put( ad_monit_id(i_prcp), prcp(:,:) )

    !##### Energy Budget #####

    call thermodyn_temp_pres( temp(:,:,:),  & ! [OUT]
                              pres(:,:,:),  & ! [OUT]
                              dens(:,:,:),  & ! [IN]
                              rhot(:,:,:),  & ! [IN]
                              qtrc(:,:,:,:) ) ! [IN]

    !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       engp(k,i,j) = dens(k,i,j) * grav * real_cz(k,i,j)

       engk(k,i,j) = 0.5_rp * dens(k,i,j) * ( w(k,i,j)**2 &
                                            + u(k,i,j)**2 &
                                            + v(k,i,j)**2 )

       engi(k,i,j) = dens(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
       do iq = qqs, qqe
          engi(k,i,j) = engi(k,i,j) &
                      + dens(k,i,j) * qtrc(k,i,j,iq) * temp(k,i,j) * cvw(iq)
       enddo

       engi(k,i,j) = engi(k,i,j) + dens(k,i,j) * qtrc(k,i,j,i_qv) * lhv ! Latent Heat [vapor->liquid]

       do iq = qis, qie
          engi(k,i,j) = engi(k,i,j) - dens(k,i,j) * qtrc(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_engt), engt(:,:,:) )
    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_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(k,i,j)) / ( dens(k+1,i,j) + dens(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(k,i,j)) / ( dens(k,i+1,j) + dens(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(k,i,j)) / ( dens(k,i,j+1) + dens(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 2617 of file mod_atmos_vars.f90.

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_create(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_create(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_create(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_create(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_create(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_create(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_create(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_create(), atmos_restart_out_basename, atmos_restart_out_dtype, atmos_restart_out_title, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, scale_fileio::fileio_create(), scale_stdio::io_fid_log, scale_stdio::io_l, scale_tracer::qa, scale_time::time_gettimelabel(), and scale_time::time_nowsec.

Referenced by mod_admin_restart::admin_restart().

    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=20)     :: timelabel
    character(len=H_LONG) :: basename

    integer :: iq
    !---------------------------------------------------------------------------

#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

       call time_gettimelabel( timelabel )
       write(basename,'(A,A,A)') trim(atmos_restart_out_basename), '_', trim(timelabel)

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Output restart file (ATMOS) ***'
       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call fileio_create(restart_fid, basename, atmos_restart_out_title, atmos_restart_out_dtype)

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

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_enddef(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_enddef(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_enddef(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_enddef(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_enddef(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_enddef(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_enddef(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_enddef(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, and scale_fileio::fileio_enddef().

Referenced by mod_admin_restart::admin_restart().

    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 .NE. -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 2758 of file mod_atmos_vars.f90.

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_close(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_close(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_close(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_close(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_close(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_close(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_close(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_close(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, and scale_fileio::fileio_close().

Referenced by mod_admin_restart::admin_restart().

    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 .NE. -1 ) then
       call fileio_close( restart_fid ) ! [IN]
       restart_fid = -1

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

References scale_tracer::aq_desc, scale_tracer::aq_name, scale_tracer::aq_unit, mod_atmos_dyn_vars::atmos_dyn_vars_restart_def_var(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_def_var(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_def_var(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_def_var(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_def_var(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_def_var(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_def_var(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_def_var(), atmos_restart_out_dtype, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, scale_fileio::fileio_def_var(), scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, and scale_tracer::qa.

Referenced by mod_admin_restart::admin_restart().

    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 .NE. -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), aq_name(iq), aq_desc(iq), & ! [IN]
                               aq_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_var()

subroutine, public mod_atmos_vars::atmos_vars_restart_write_var

(

)

Write restart of atmospheric variables.

Definition at line 2900 of file mod_atmos_vars.f90.

References scale_tracer::aq_name, mod_atmos_dyn_vars::atmos_dyn_vars_restart_write_var(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_write_var(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_write_var(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_write_var(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_write_var(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_write_var(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_write_var(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_write_var(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, atmos_vars_fillhalo(), atmos_vars_total(), dens, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, momx, momy, momz, scale_tracer::qa, qtrc, and rhot.

Referenced by mod_admin_restart::admin_restart().

    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_var
    use mod_atmos_phy_mp_vars, only: &
       atmos_phy_mp_vars_restart_write_var
    use mod_atmos_phy_ae_vars, only: &
       atmos_phy_ae_vars_restart_write_var
    use mod_atmos_phy_ch_vars, only: &
       atmos_phy_ch_vars_restart_write_var
    use mod_atmos_phy_rd_vars, only: &
       atmos_phy_rd_vars_restart_write_var
    use mod_atmos_phy_sf_vars, only: &
       atmos_phy_sf_vars_restart_write_var
    use mod_atmos_phy_tb_vars, only: &
       atmos_phy_tb_vars_restart_write_var
    use mod_atmos_phy_cp_vars, only: &
       atmos_phy_cp_vars_restart_write_var
#ifdef _SDM
    use scale_atmos_phy_mp_sdm, only: &
       sd_rest_flg_out, &
       atmos_phy_mp_sdm_restart_write_var
#endif
    implicit none

    integer iq
    !---------------------------------------------------------------------------

#ifdef _SDM
    if( sd_rest_flg_out ) then
       call atmos_phy_mp_sdm_restart_write_var
    endif
#endif

    if ( restart_fid .NE. -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), aq_name(iq), 'ZXY' ) ! [IN]
       enddo

    endif

    if( atmos_sw_dyn )    call atmos_dyn_vars_restart_write_var
    if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_write_var
    if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_write_var
    if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_write_var
    if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_write_var
    if( atmos_sw_phy_sf ) call atmos_phy_sf_vars_restart_write_var
    if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_write_var
    if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_write_var

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

Referenced by mod_admin_restart::admin_restart(), mod_admin_restart::admin_restart_setup(), atmos_vars_setup(), mod_mkinit::mkinit(), and mod_rm_driver::scalerm().

  logical,               public :: atmos_restart_output           = .false.         

atmos_restart_check

logical, public mod_atmos_vars::atmos_restart_check = .false.

check value consistency?

Definition at line 57 of file mod_atmos_vars.f90.

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

  logical,               public :: atmos_restart_check            = .false.         

atmos_restart_in_basename

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

basename of the restart file

Definition at line 59 of file mod_atmos_vars.f90.

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

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

atmos_restart_out_basename

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

basename of the output file

Definition at line 60 of file mod_atmos_vars.f90.

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

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

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

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

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

atmos_restart_out_dtype

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

REAL4 or REAL8.

Definition at line 62 of file mod_atmos_vars.f90.

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

  character(len=H_MID),  public :: atmos_restart_out_dtype        = 'DEFAULT'

atmos_restart_in_allowmissingq

logical, public mod_atmos_vars::atmos_restart_in_allowmissingq = .false.

Definition at line 63 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

  logical,               public :: atmos_restart_in_allowmissingq = .false.

atmos_restart_check_basename

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

Definition at line 65 of file mod_atmos_vars.f90.

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

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

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

  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 69 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_diagnostics(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), 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().

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

momz

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

Definition at line 70 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_diagnostics(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

  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 71 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_diagnostics(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

  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 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_tb_driver::atmos_phy_tb_driver(), atmos_vars_diagnostics(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), and mod_realinput::realinput_surface().

  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 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_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_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), 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().

  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 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_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_diagnostics(), atmos_vars_fillhalo(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_restart_write_var(), atmos_vars_setup(), atmos_vars_total(), mod_mkinit::diag_ds(), mod_mkinit::interporation_fact(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_realinput::realinput_surface(), and mod_mkinit::rect_setup().

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

Referenced by atmos_vars_setup().

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

momz_avw

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

Definition at line 77 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

momx_avw

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

Definition at line 78 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

momy_avw

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

Definition at line 79 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

rhot_avw

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

Definition at line 80 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

qtrc_avw

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

Definition at line 81 of file mod_atmos_vars.f90.

Referenced by atmos_vars_setup().

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

dens_av

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

Definition at line 83 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

momz_av

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

Definition at line 84 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

momx_av

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

Definition at line 85 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

momy_av

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

Definition at line 86 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

rhot_av

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

Definition at line 87 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

qtrc_av

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

Definition at line 88 of file mod_atmos_vars.f90.

Referenced by mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), atmos_vars_restart_read(), and atmos_vars_setup().

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

dens_tp

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

Definition at line 91 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().

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

momz_tp

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

Definition at line 92 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().

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

momx_tp

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

Definition at line 93 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().

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

momy_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::momy_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(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver(), and atmos_vars_setup().

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

rhot_tp

real(rp), dimension(:,:,:), allocatable, public mod_atmos_vars::rhot_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_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().

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

rhoq_tp

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_vars::rhoq_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_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().

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

temp

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

Definition at line 99 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_monitor(), and atmos_vars_setup().

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

pres

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

Definition at line 100 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_monitor(), and atmos_vars_setup().

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

w

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

Definition at line 101 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(), and atmos_vars_setup().

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

u

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

Definition at line 102 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(), and atmos_vars_setup().

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

v

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

Definition at line 103 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(), rdma_put8_(), rdma_put_(), and rdma_setup_().

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

pott

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

Definition at line 104 of file mod_atmos_vars.f90.

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

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