SCALE-RM
Data Types | Functions/Subroutines | Variables
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_open
 Open restart file for reading atmospheric variables. More...
 
subroutine, public atmos_vars_restart_read
 Read restart of atmospheric variables. More...
 
subroutine, public atmos_vars_history_setpres
 Set pressure for history output. More...
 
subroutine, public atmos_vars_restart_check
 Check and compare between last data and sample data. More...
 
subroutine, public atmos_vars_history
 History output set for atmospheric variables. More...
 
subroutine, public atmos_vars_check (force)
 Check variables for atmosphere. More...
 
subroutine, public atmos_vars_calc_diagnostics
 Calc diagnostic variables. More...
 
recursive subroutine atmos_vars_get_diagnostic_3d (vname, var)
 get diagnostic variable 3D More...
 
recursive subroutine atmos_vars_get_diagnostic_2d (vname, var)
 get diagnostic variable 2D More...
 
recursive subroutine atmos_vars_get_diagnostic_1d (vname, var)
 get diagnostic variable 1D More...
 
subroutine, public atmos_vars_monitor
 monitor output More...
 
subroutine, public atmos_vars_restart_create
 Create atmospheric restart file. More...
 
subroutine, public atmos_vars_restart_enddef
 Exit netCDF define mode. More...
 
subroutine, public atmos_vars_restart_close
 Close restart file. More...
 
subroutine, public atmos_vars_restart_def_var
 Define atmospheric variables in restart file. More...
 
subroutine, public atmos_vars_restart_write
 Write restart of atmospheric variables. More...
 
subroutine allocate_3d (ary)
 
subroutine allocate_2d (ary)
 
subroutine allocate_1d (ary)
 

Variables

logical, public atmos_restart_output = .false.
 Output restart file? More...
 
character(len=h_long), public atmos_restart_in_basename = ''
 Basename of the input file. More...
 
logical, public atmos_restart_in_aggregate
 Switch to use aggregate file. More...
 
logical, public atmos_restart_in_postfix_timelabel = .false.
 Add timelabel to the basename of input file? More...
 
character(len=h_long), public atmos_restart_out_basename = ''
 Basename of the output file. More...
 
logical, public atmos_restart_out_aggregate
 Switch to use aggregate file. More...
 
logical, public atmos_restart_out_postfix_timelabel = .true.
 Add timelabel to the basename of output file? More...
 
character(len=h_mid), public atmos_restart_out_title = 'ATMOS restart'
 Title of the output file. More...
 
character(len=h_short), public atmos_restart_out_dtype = 'DEFAULT'
 REAL4 or REAL8. More...
 
logical, public atmos_restart_check = .false.
 Check value consistency? More...
 
character(len=h_long), public atmos_restart_check_basename = 'restart_check'
 
real(rp), public atmos_restart_check_criterion = 1.E-6_RP
 
real(rp), dimension(:,:,:), allocatable, target, public dens
 
real(rp), dimension(:,:,:), allocatable, target, public momz
 
real(rp), dimension(:,:,:), allocatable, target, public momx
 
real(rp), dimension(:,:,:), allocatable, target, public momy
 
real(rp), dimension(:,:,:), allocatable, target, public rhot
 
real(rp), dimension(:,:,:,:), allocatable, target, public qtrc
 
real(rp), dimension(:,:,:), allocatable, target, public dens_avw
 
real(rp), dimension(:,:,:), allocatable, target, public momz_avw
 
real(rp), dimension(:,:,:), allocatable, target, public momx_avw
 
real(rp), dimension(:,:,:), allocatable, target, public momy_avw
 
real(rp), dimension(:,:,:), allocatable, target, public rhot_avw
 
real(rp), dimension(:,:,:,:), allocatable, target, public qtrc_avw
 
real(rp), dimension(:,:,:), pointer, public dens_av
 
real(rp), dimension(:,:,:), pointer, public momz_av
 
real(rp), dimension(:,:,:), pointer, public momx_av
 
real(rp), dimension(:,:,:), pointer, public momy_av
 
real(rp), dimension(:,:,:), pointer, public rhot_av
 
real(rp), dimension(:,:,:,:), pointer, public qtrc_av
 
real(rp), dimension(:,:,:), allocatable, pointer, target, public qv
 
real(rp), dimension(:,:,:), pointer, public qc
 
real(rp), dimension(:,:,:), pointer, public qr
 
real(rp), dimension(:,:,:), pointer, public qi
 
real(rp), dimension(:,:,:), pointer, public qs
 
real(rp), dimension(:,:,:), pointer, public qg
 
real(rp), dimension(:,:,:), pointer, public qh
 
real(rp), dimension(:,:,:,:), allocatable, target, public qe
 
real(rp), allocatable, target, public mass
 
real(rp), allocatable, target, public ratio
 
real(rp), allocatable, target, public of
 
real(rp), allocatable, target, public hydrometors
 
real(rp), allocatable, target, public kg
 
real(rp), dimension(:,:,:), allocatable, public dens_ref
 
real(rp), dimension(:,:,:), allocatable, public pott_ref
 
real(rp), dimension(:,:,:), allocatable, public temp_ref
 
real(rp), dimension(:,:,:), allocatable, public pres_ref
 
real(rp), dimension(:,:,:), allocatable, public qv_ref
 
real(rp), dimension(:,:,:), allocatable, public dens_tp
 
real(rp), dimension(:,:,:), allocatable, public momz_tp
 
real(rp), dimension(:,:,:), allocatable, public rhou_tp
 
real(rp), dimension(:,:,:), allocatable, public rhov_tp
 
real(rp), dimension(:,:,:), allocatable, public rhot_tp
 
real(rp), dimension(:,:,:), allocatable, public rhoh_p
 
real(rp), dimension(:,:,:,:), allocatable, public rhoq_tp
 
real(rp), dimension(:,:,:), allocatable, public momx_tp
 
real(rp), dimension(:,:,:), allocatable, public momy_tp
 
real(rp), dimension(:,:,:), allocatable, target, public w
 
real(rp), allocatable, target, public velocity
 
real(rp), allocatable, target, public m
 
real(rp), allocatable, target, public s
 
real(rp), dimension(:,:,:), allocatable, target, public u
 
real(rp), dimension(:,:,:), allocatable, target, public v
 
real(rp), dimension(:,:,:), allocatable, target, public pott
 
real(rp), allocatable, target, public potential
 
real(rp), allocatable, target, public temperature
 
real(rp), allocatable, target, public k
 
real(rp), dimension(:,:,:), allocatable, target, public temp
 
real(rp), dimension(:,:,:), allocatable, target, public pres
 
real(rp), allocatable, target, public pressure
 
real(rp), allocatable, target, public pa =J/m3]
 
real(rp), dimension(:,:,:), allocatable, target, public exner
 
real(rp), dimension(t/pt), allocatable, target, public function
 
real(rp), dimension(:,:,:), allocatable, target, public phyd
 
real(rp), allocatable, target, public hydrostatic
 
real(rp), dimension(:,:,:), allocatable, target, public phydh
 
real(rp), allocatable, target, public layer
 
real(rp), allocatable, target, public interface
 
real(rp), dimension(:,:,:), allocatable, target, public qdry
 
real(rp), allocatable, target, public dry
 
real(rp), allocatable, target, public air
 
real(rp), dimension(:,:,:), allocatable, target, public rtot
 
real(rp), allocatable, target, public specific
 
real(rp), allocatable, target, public gass
 
real(rp), allocatable, target, public constant
 
real(rp), allocatable, target, public j
 
real(rp), dimension(:,:,:), allocatable, target, public cvtot
 
real(rp), allocatable, target, public heat
 
real(rp), dimension(:,:,:), allocatable, target, public cptot
 
real(rp), dimension(:,:), allocatable, target, public prec
 
real(rp), allocatable, target, public total
 
real(rp), allocatable, target, public precipitation
 
real(rp), allocatable, target, public m2
 
real(rp), dimension(:,:), allocatable, public prec_engi
 
real(rp), allocatable, target, public internal
 
real(rp), allocatable, target, public energy
 

Detailed Description

module ATMOSPHERIC Variables

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

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

Function/Subroutine Documentation

◆ atmos_vars_setup()

subroutine, public mod_atmos_vars::atmos_vars_setup

Setup.

Definition at line 467 of file mod_atmos_vars.F90.

467  use scale_const, only: &
468  undef => const_undef
469  use scale_prc, only: &
470  prc_abort
471  use scale_file_history, only: &
472  file_history_reg
473  use scale_monitor, only: &
474  monitor_reg
475  use scale_atmos_hydrometeor, only: &
476  atmos_hydrometeor_dry, &
477  n_hyd, &
478  i_qv, &
479  i_hc, &
480  i_hr, &
481  i_hi, &
482  i_hs, &
483  i_hg, &
484  i_hh
485  use mod_atmos_admin, only: &
486  atmos_use_average
487  use mod_atmos_dyn_vars, only: &
488  atmos_dyn_vars_setup
489  use mod_atmos_phy_mp_vars, only: &
490  atmos_phy_mp_vars_setup
491  use mod_atmos_phy_ae_vars, only: &
492  atmos_phy_ae_vars_setup
493  use mod_atmos_phy_ch_vars, only: &
494  atmos_phy_ch_vars_setup
495  use mod_atmos_phy_rd_vars, only: &
496  atmos_phy_rd_vars_setup
497  use mod_atmos_phy_sf_vars, only: &
498  atmos_phy_sf_vars_setup
499  use mod_atmos_phy_tb_vars, only: &
500  atmos_phy_tb_vars_setup
501  use mod_atmos_phy_bl_vars, only: &
502  atmos_phy_bl_vars_setup
503  use mod_atmos_phy_cp_vars, only: &
504  atmos_phy_cp_vars_setup
505  use mod_atmos_phy_lt_vars, only: &
506  atmos_phy_lt_vars_setup
507  implicit none
508 
509  namelist / param_atmos_vars / &
510  atmos_restart_in_basename, &
511  atmos_restart_in_aggregate, &
512  atmos_restart_in_postfix_timelabel, &
513  atmos_restart_in_check_coordinates, &
514  atmos_restart_output, &
515  atmos_restart_out_basename, &
516  atmos_restart_out_aggregate, &
517  atmos_restart_out_postfix_timelabel, &
518  atmos_restart_out_title, &
519  atmos_restart_out_dtype, &
520  atmos_restart_check, &
521  atmos_restart_check_basename, &
522  atmos_restart_check_criterion, &
523  atmos_vars_checkrange, &
524  atmos_vars_checkcfl_soft, &
525  atmos_vars_checkcfl_hard
526 
527  integer :: ierr
528  integer :: iv, iq
529  !---------------------------------------------------------------------------
530 
531  log_newline
532  log_info("ATMOS_vars_setup",*) 'Setup'
533 
534  allocate( dens(ka,ia,ja) )
535  allocate( momz(ka,ia,ja) )
536  allocate( momx(ka,ia,ja) )
537  allocate( momy(ka,ia,ja) )
538  allocate( rhot(ka,ia,ja) )
539  allocate( qtrc(ka,ia,ja,max(qa,1)) )
540 
541  if ( atmos_use_average ) then
542  allocate( dens_avw(ka,ia,ja) )
543  allocate( momz_avw(ka,ia,ja) )
544  allocate( momx_avw(ka,ia,ja) )
545  allocate( momy_avw(ka,ia,ja) )
546  allocate( rhot_avw(ka,ia,ja) )
547  allocate( qtrc_avw(ka,ia,ja,max(qa,1)) )
548 
549  dens_av => dens_avw
550  momz_av => momz_avw
551  momx_av => momx_avw
552  momy_av => momy_avw
553  rhot_av => rhot_avw
554  qtrc_av => qtrc_avw
555  else
556  dens_av => dens
557  momz_av => momz
558  momx_av => momx
559  momy_av => momy
560  rhot_av => rhot
561  qtrc_av => qtrc
562  endif
563 
564  allocate( dens_tp(ka,ia,ja) )
565  allocate( momz_tp(ka,ia,ja) )
566  allocate( rhou_tp(ka,ia,ja) )
567  allocate( rhov_tp(ka,ia,ja) )
568  allocate( rhot_tp(ka,ia,ja) )
569  allocate( rhoh_p(ka,ia,ja) )
570  allocate( rhoq_tp(ka,ia,ja,max(qa,1)) )
571 
572  allocate( w(ka,ia,ja) )
573  allocate( u(ka,ia,ja) )
574  allocate( v(ka,ia,ja) )
575  w(:,:,:) = undef
576  u(:,:,:) = undef
577  v(:,:,:) = undef
578 
579  allocate( pott(ka,ia,ja) )
580  allocate( temp(ka,ia,ja) )
581  allocate( pres(ka,ia,ja) )
582  allocate( exner(ka,ia,ja) )
583  allocate( phyd(ka,ia,ja) )
584  allocate( phydh(0:ka,ia,ja) )
585  pott(:,:,:) = undef
586  temp(:,:,:) = undef
587  pres(:,:,:) = undef
588  exner(:,:,:) = undef
589  phyd(:,:,:) = undef
590  phydh(:,:,:) = undef
591 
592  allocate( qdry(ka,ia,ja) )
593  allocate( rtot(ka,ia,ja) )
594  allocate( cvtot(ka,ia,ja) )
595  allocate( cptot(ka,ia,ja) )
596  qdry(:,:,:) = undef
597  rtot(:,:,:) = undef
598  cvtot(:,:,:) = undef
599  cptot(:,:,:) = undef
600 
601  allocate( prec(ia,ja) )
602  allocate( prec_engi(ia,ja) )
603  prec(:,:) = undef
604  prec_engi(:,:) = undef
605 
606  ! obsolute
607  allocate( momx_tp(ka,ia,ja) )
608  allocate( momy_tp(ka,ia,ja) )
609 
610 
611  momz(1:ks-1,:,:) = 0.0_rp
612  momz(ke:ka,:,:) = 0.0_rp
613 
614  allocate( work3d(ka,ia,ja) )
615  allocate( work2d( ia,ja) )
616  allocate( work1d(ka ) )
617 
618 
619  !--- read namelist
620  rewind(io_fid_conf)
621  read(io_fid_conf,nml=param_atmos_vars,iostat=ierr)
622  if( ierr < 0 ) then !--- missing
623  log_info("ATMOS_vars_setup",*) 'Not found namelist. Default used.'
624  elseif( ierr > 0 ) then !--- fatal error
625  log_error("ATMOS_vars_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_VARS. Check!'
626  call prc_abort
627  endif
628  log_nml(param_atmos_vars)
629 
630  log_newline
631  log_info("ATMOS_vars_setup",*) 'List of prognostic variables (ATMOS) '
632  log_info_cont('(1x,A,A24,A,A48,A,A12,A)') &
633  ' |', 'VARNAME ','|', &
634  'DESCRIPTION ', '[', 'UNIT ', ']'
635  do iv = 1, pv_nmax
636  log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
637  'NO.',iv,'|',pv_info(iv)%NAME,'|', pv_info(iv)%DESC,'[', pv_info(iv)%UNIT,']'
638  enddo
639  do iq = 1, qa
640  log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
641  'NO.',5+iq,'|',tracer_name(iq),'|', tracer_desc(iq),'[', tracer_unit(iq),']'
642  enddo
643 
644  log_newline
645  if ( atmos_restart_in_basename /= '' ) then
646  log_info("ATMOS_vars_setup",*) 'Restart input? : YES, file = ', trim(atmos_restart_in_basename)
647  log_info("ATMOS_vars_setup",*) 'Add timelabel? : ', atmos_restart_in_postfix_timelabel
648  else
649  log_info("ATMOS_vars_setup",*) 'Restart input? : NO'
650  endif
651  if ( atmos_restart_output &
652  .AND. atmos_restart_out_basename /= '' ) then
653  log_info("ATMOS_vars_setup",*) 'Restart output? : YES, file = ', trim(atmos_restart_out_basename)
654  log_info("ATMOS_vars_setup",*) 'Add timelabel? : ', atmos_restart_out_postfix_timelabel
655  else
656  log_info("ATMOS_vars_setup",*) 'Restart output? : NO'
657  atmos_restart_output = .false.
658  endif
659 
660  if ( atmos_restart_check_basename == '' ) then
661  atmos_restart_check = .false.
662  endif
663 
664  if ( atmos_vars_checkcfl_hard > 0.0_rp ) then
665  atmos_vars_checkcfl_soft = min( atmos_vars_checkcfl_soft, atmos_vars_checkcfl_hard )
666  endif
667 
668  log_newline
669  log_info("ATMOS_vars_setup",*) 'Check restart consistency? : ', atmos_restart_check
670  log_info("ATMOS_vars_setup",*) 'Check value range of variables? : ', atmos_vars_checkrange
671  if ( atmos_vars_checkcfl_soft > 0.0_rp ) then
672  log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to warn : ', atmos_vars_checkcfl_soft
673  else
674  log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to warn : disabled'
675  endif
676  if ( atmos_vars_checkcfl_hard > 0.0_rp ) then
677  log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to stop : ', atmos_vars_checkcfl_hard
678  else
679  log_info("ATMOS_vars_setup",*) 'Threshold of Courant number to stop : disabled'
680  endif
681 
682  call atmos_dyn_vars_setup
683  call atmos_phy_mp_vars_setup
684  call atmos_phy_ae_vars_setup
685  call atmos_phy_ch_vars_setup
686  call atmos_phy_rd_vars_setup
687  call atmos_phy_sf_vars_setup
688  call atmos_phy_tb_vars_setup
689  call atmos_phy_bl_vars_setup
690  call atmos_phy_cp_vars_setup
691  call atmos_phy_lt_vars_setup
692 
693 
694  ! water content
695  if ( atmos_hydrometeor_dry ) then
696  allocate( zero(ka,ia,ja) )
697 !OCL XFILL
698  zero(:,:,:) = 0.0_rp
699 
700  qv => zero
701  qc => zero
702  qr => zero
703  qi => zero
704  qs => zero
705  qg => zero
706  qh => zero
707 
708  moist = .false.
709  else
710  allocate( qe(ka,ia,ja,n_hyd) )
711 !OCL XFILL
712  qe(:,:,:,:) = undef
713 
714  qv => qtrc_av(:,:,:,i_qv)
715  qc => qe(:,:,:,i_hc)
716  qr => qe(:,:,:,i_hr)
717  qi => qe(:,:,:,i_hi)
718  qs => qe(:,:,:,i_hs)
719  qg => qe(:,:,:,i_hg)
720  qh => qe(:,:,:,i_hh)
721 
722  moist = .true.
723  end if
724 
725 
726  dv_calculated(dv_nmax) = .false.
727 
728  !-----< history output setup >-----
729  allocate( qp_hist_id( max(qa,1) ) )
730  allocate( qp_monit_id( max(qa,1) ) )
731  pv_hist_id(:) = -1
732  pv_monit_id(:) = -1
733  qp_hist_id(:) = -1
734  qp_monit_id(:) = -1
735  dv_hist_id(:) = -1
736  dv_monit_id(:) = -1
737 
738 
739  do iv = 1, pv_nmax
740  call file_history_reg( pv_info(iv)%NAME, pv_info(iv)%DESC, pv_info(iv)%UNIT, pv_hist_id(iv), dim_type=pv_info(iv)%dim_type, standard_name=pv_info(iv)%STDNAME )
741  end do
742 
743  do iq = 1, qa
744  call file_history_reg( tracer_name(iq), tracer_desc(iq), tracer_unit(iq), qp_hist_id(iq), dim_type='ZXY' )
745  enddo
746 
747  do iv = 1, dv_nmax
748  call file_history_reg( dv_info(iv)%NAME, dv_info(iv)%DESC, dv_info(iv)%UNIT, dv_hist_id(iv), dim_type=dv_info(iv)%dim_type, standard_name=dv_info(iv)%STDNAME )
749  end do
750 
751  call file_history_reg( "GPH", "geopotential height", "m", hist_id_gph, dim_type='ZXY', standard_name="geopotential_height" )
752 
753 
754  !-----< monitor output setup >-----
755  do iv = 1, pv_nmax
756  call monitor_reg( pv_info(iv)%NAME, pv_info(iv)%DESC, trim(pv_info(iv)%UNIT)//"*m3", & ! (in)
757  pv_monit_id(iv), & ! (out)
758  dim_type=pv_info(iv)%dim_type, is_tendency=.false. ) ! (in)
759  end do
760  do iq = 1, qa
761  call monitor_reg( tracer_name(iq), tracer_desc(iq), tracer_unit(iq)//"*kg", & ! (in)
762  qp_monit_id(iq), & ! (out)
763  dim_type='ZXY', is_tendency=.false. ) ! (in)
764  enddo
765 
766  call monitor_reg( 'QDRY', 'dry air mass', 'kg', & ! (in)
767  dv_monit_id(im_qdry), & ! (out)
768  dim_type='ZXY', is_tendency=.false. ) ! (in)
769  call monitor_reg( 'QTOT', 'water mass', 'kg', & ! (in)
770  dv_monit_id(im_qtot), & ! (out)
771  dim_type='ZXY', is_tendency=.false. ) ! (in)
772  call monitor_reg( 'EVAP', 'evaporation at the surface', 'kg', & ! (in)
773  dv_monit_id(im_evap), & ! (out)
774  dim_type='XY', is_tendency=.true. ) ! (in)
775  call monitor_reg( 'PREC', 'precipitation', 'kg', & ! (in)
776  dv_monit_id(im_prec), & ! (out)
777  dim_type='XY', is_tendency=.true. ) ! (in)
778 
779  call monitor_reg( 'ENGT', 'total energy', 'J', & ! (in)
780  dv_monit_id(im_engt), & ! (out)
781  dim_type='ZXY', is_tendency=.false. ) ! (in)
782  call monitor_reg( 'ENGP', 'potential energy', 'J', & ! (in)
783  dv_monit_id(im_engp), & ! (out)
784  dim_type='ZXY', is_tendency=.false. ) ! (in)
785  call monitor_reg( 'ENGK', 'kinetic energy', 'J', & ! (in)
786  dv_monit_id(im_engk), & ! (out)
787  dim_type='ZXY', is_tendency=.false. ) ! (in)
788  call monitor_reg( 'ENGI', 'internal energy', 'J', & ! (in)
789  dv_monit_id(im_engi), & ! (out)
790  dim_type='ZXY', is_tendency=.false. ) ! (in)
791 
792  call monitor_reg( 'ENGFLXT', 'total energy flux convergence', 'J', & ! (in)
793  dv_monit_id(im_engflxt), & ! (out)
794  dim_type='XY', is_tendency=.true. ) ! (in)
795  call monitor_reg( 'ENGSFC_SH', 'SFC sensible heat flux', 'J', & ! (in)
796  dv_monit_id(im_engsfc_sh), & ! (out)
797  dim_type='XY', is_tendency=.true. ) ! (in)
798  call monitor_reg( 'ENGSFC_LH', 'SFC latent heat flux', 'J', & ! (in)
799  dv_monit_id(im_engsfc_lh), & ! (out)
800  dim_type='XY', is_tendency=.true. ) ! (in)
801  call monitor_reg( 'ENGSFC_EVAP', 'SFC internal energy flux of the evapolation', 'J', & ! (in)
802  dv_monit_id(im_engsfc_evap), & ! (out)
803  dim_type='XY', is_tendency=.true. ) ! (in)
804  call monitor_reg( 'ENGSFC_PREC', 'SFC internal energy flux of the precipitation', 'J', & ! (in)
805  dv_monit_id(im_engsfc_prec), & ! (out)
806  dim_type='XY', is_tendency=.true. ) ! (in)
807  call monitor_reg( 'ENGSFC_RD', 'SFC net radiation flux', 'J', & ! (in)
808  dv_monit_id(im_engsfc_rd), & ! (out)
809  dim_type='XY', is_tendency=.true. ) ! (in)
810  call monitor_reg( 'ENGTOM_RD', 'TOM net radiation flux', 'J', & ! (in)
811  dv_monit_id(im_engtom_rd), & ! (out)
812  dim_type='XY', is_tendency=.true. ) ! (in)
813 
814  call monitor_reg( 'ENGSFC_LW_up', 'SFC LW upward flux', 'J', & ! (in)
815  dv_monit_id(im_engsfc_lw_up), & ! (out)
816  dim_type='XY', is_tendency=.true. ) ! (in)
817  call monitor_reg( 'ENGSFC_LW_dn', 'SFC LW downward flux', 'J', & ! (in)
818  dv_monit_id(im_engsfc_lw_dn), & ! (out)
819  dim_type='XY', is_tendency=.true. ) ! (in)
820  call monitor_reg( 'ENGSFC_SW_up', 'SFC SW upward flux', 'J', & ! (in)
821  dv_monit_id(im_engsfc_sw_up), & ! (out)
822  dim_type='XY', is_tendency=.true. ) ! (in)
823  call monitor_reg( 'ENGSFC_SW_dn', 'SFC SW downward flux', 'J', & ! (in)
824  dv_monit_id(im_engsfc_sw_dn), & ! (out)
825  dim_type='XY', is_tendency=.true. ) ! (in)
826 
827  call monitor_reg( 'ENGTOM_LW_up', 'TOM LW upward flux', 'J', & ! (in)
828  dv_monit_id(im_engtom_lw_up), & ! (out)
829  dim_type='XY', is_tendency=.true. ) ! (in)
830  call monitor_reg( 'ENGTOM_LW_dn', 'TOM LW downward flux', 'J', & ! (in)
831  dv_monit_id(im_engtom_lw_dn), & ! (out)
832  dim_type='XY', is_tendency=.true. ) ! (in)
833  call monitor_reg( 'ENGTOM_SW_up', 'TOM SW upward flux', 'J', & ! (in)
834  dv_monit_id(im_engtom_sw_up), & ! (out)
835  dim_type='XY', is_tendency=.true. ) ! (in)
836  call monitor_reg( 'ENGTOM_SW_dn', 'TOM SW downward flux', 'J', & ! (in)
837  dv_monit_id(im_engtom_sw_dn), & ! (out)
838  dim_type='XY', is_tendency=.true. ) ! (in)
839 
840  return

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

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

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 847 of file mod_atmos_vars.F90.

847  use scale_comm_cartesc, only: &
848  comm_vars8, &
849  comm_wait
850  implicit none
851 
852  logical, intent(in), optional :: FILL_BND
853 
854  logical :: FILL_BND_
855  integer :: i, j, iq
856  !---------------------------------------------------------------------------
857 
858  fill_bnd_ = .false.
859  if ( present(fill_bnd) ) fill_bnd_ = fill_bnd
860 
861  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
862  do j = jsb, jeb
863  do i = isb, ieb
864  dens( 1:ks-1,i,j) = dens(ks,i,j)
865  momz( 1:ks-2,i,j) = momz(ks-1,i,j)
866  momx( 1:ks-1,i,j) = momx(ks,i,j)
867  momy( 1:ks-1,i,j) = momy(ks,i,j)
868  rhot( 1:ks-1,i,j) = rhot(ks,i,j)
869  dens(ke+1:ka, i,j) = dens(ke,i,j)
870  momz(ke+1:ka, i,j) = momz(ke,i,j)
871  momx(ke+1:ka, i,j) = momx(ke,i,j)
872  momy(ke+1:ka, i,j) = momy(ke,i,j)
873  rhot(ke+1:ka, i,j) = rhot(ke,i,j)
874  enddo
875  enddo
876 
877  !$omp parallel do private(i,j,iq) OMP_SCHEDULE_ collapse(3)
878  do iq = 1, qa
879  do j = jsb, jeb
880  do i = isb, ieb
881  qtrc( 1:ks-1,i,j,iq) = qtrc(ks,i,j,iq)
882  qtrc(ke+1:ka, i,j,iq) = qtrc(ke,i,j,iq)
883  enddo
884  enddo
885  enddo
886 
887  call comm_vars8( dens(:,:,:), 1 )
888  call comm_vars8( momz(:,:,:), 2 )
889  call comm_vars8( momx(:,:,:), 3 )
890  call comm_vars8( momy(:,:,:), 4 )
891  call comm_vars8( rhot(:,:,:), 5 )
892  call comm_wait ( dens(:,:,:), 1, fill_bnd_ )
893  call comm_wait ( momz(:,:,:), 2, fill_bnd_ )
894  call comm_wait ( momx(:,:,:), 3, fill_bnd_ )
895  call comm_wait ( momy(:,:,:), 4, fill_bnd_ )
896  call comm_wait ( rhot(:,:,:), 5, fill_bnd_ )
897 
898  do iq = 1, qa
899  call comm_vars8( qtrc(:,:,:,iq), iq )
900  enddo
901  do iq = 1, qa
902  call comm_wait ( qtrc(:,:,:,iq), iq, fill_bnd_ )
903  enddo
904 
905  return

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

Referenced by mod_atmos_driver::atmos_driver_update(), atmos_vars_restart_read(), and atmos_vars_restart_write().

Here is the caller graph for this function:

◆ atmos_vars_restart_open()

subroutine, public mod_atmos_vars::atmos_vars_restart_open

Open restart file for reading atmospheric variables.

Definition at line 911 of file mod_atmos_vars.F90.

911  use scale_prc, only: &
912  prc_abort
913  use scale_const, only: &
914  grav => const_grav
915  use scale_time, only: &
916  time_gettimelabel
917  use scale_file_cartesc, only: &
918  file_cartesc_open, &
919  file_cartesc_check_coordinates
920  use mod_atmos_admin, only: &
921  atmos_use_average, &
922  atmos_sw_dyn, &
923  atmos_sw_phy_mp, &
924  atmos_sw_phy_ae, &
925  atmos_sw_phy_ch, &
926  atmos_sw_phy_rd, &
927  atmos_sw_phy_sf, &
928  atmos_sw_phy_tb, &
929  atmos_sw_phy_bl, &
930  atmos_sw_phy_cp, &
931  atmos_sw_phy_lt
932  use mod_atmos_dyn_vars, only: &
933  atmos_dyn_vars_restart_open
934  use mod_atmos_phy_mp_vars, only: &
935  atmos_phy_mp_vars_restart_open
936  use mod_atmos_phy_ae_vars, only: &
937  atmos_phy_ae_vars_restart_open
938  use mod_atmos_phy_ch_vars, only: &
939  atmos_phy_ch_vars_restart_open
940  use mod_atmos_phy_rd_vars, only: &
941  atmos_phy_rd_vars_restart_open
942  use mod_atmos_phy_sf_vars, only: &
943  atmos_phy_sf_vars_restart_open
944  use mod_atmos_phy_tb_vars, only: &
945  atmos_phy_tb_vars_restart_open
946  use mod_atmos_phy_bl_vars, only: &
947  atmos_phy_bl_vars_restart_open
948  use mod_atmos_phy_cp_vars, only: &
949  atmos_phy_cp_vars_restart_open
950  use mod_atmos_phy_lt_vars, only: &
951  atmos_phy_lt_vars_restart_open
952  use mod_cpl_admin, only: &
953  cpl_sw
954  implicit none
955 
956  character(len=19) :: timelabel
957  character(len=H_LONG) :: basename
958  !---------------------------------------------------------------------------
959 
960  log_newline
961  log_info("ATMOS_vars_restart_open",*) 'Open restart file (ATMOS) '
962 
963  if ( atmos_restart_in_basename /= '' ) then
964 
965  if ( atmos_restart_in_postfix_timelabel ) then
966  call time_gettimelabel( timelabel )
967  basename = trim(atmos_restart_in_basename)//'_'//trim(timelabel)
968  else
969  basename = trim(atmos_restart_in_basename)
970  endif
971 
972  log_info("ATMOS_vars_restart_open",*) 'basename: ', trim(basename)
973 
974  call file_cartesc_open( basename, restart_fid, aggregate=atmos_restart_in_aggregate )
975 
976  if ( atmos_restart_in_check_coordinates ) then
977  call file_cartesc_check_coordinates( restart_fid, atmos=.true. )
978  end if
979 
980  else
981  log_error("ATMOS_vars_restart_open",*) 'restart file for atmosphere is not specified. STOP!'
982  call prc_abort
983  endif
984 
985  if ( atmos_use_average ) then
986  dens_av(:,:,:) = dens(:,:,:)
987  momz_av(:,:,:) = momz(:,:,:)
988  momx_av(:,:,:) = momx(:,:,:)
989  momy_av(:,:,:) = momy(:,:,:)
990  rhot_av(:,:,:) = rhot(:,:,:)
991  qtrc_av(:,:,:,:) = qtrc(:,:,:,:)
992  endif
993 
994  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_open
995  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_open
996  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_open
997  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_open
998  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_open
999  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_open
1000  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_open
1001  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_open
1002  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_open
1003  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_open
1004 
1005  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_open(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_open(), mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_open(), mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_open(), mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_open(), mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_open(), mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_open(), mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_open(), mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_open(), mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_open(), atmos_restart_in_aggregate, atmos_restart_in_basename, atmos_restart_in_postfix_timelabel, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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, scale_const::const_grav, mod_cpl_admin::cpl_sw, dens, dens_av, scale_file_cartesc::file_cartesc_open(), momx, momx_av, momy, momy_av, momz, momz_av, scale_prc::prc_abort(), qtrc, qtrc_av, rhot, rhot_av, and scale_time::time_gettimelabel().

Referenced by mod_admin_restart::admin_restart_read().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_vars_restart_read()

subroutine, public mod_atmos_vars::atmos_vars_restart_read

Read restart of atmospheric variables.

Definition at line 1011 of file mod_atmos_vars.F90.

1011  use scale_prc, only: &
1012  prc_abort
1013  use scale_file, only: &
1014  file_get_aggregate
1015  use scale_file_cartesc, only: &
1016  file_cartesc_read, &
1017  file_cartesc_flush
1018  use mod_atmos_admin, only: &
1019  atmos_use_average, &
1020  atmos_sw_dyn, &
1021  atmos_sw_phy_mp, &
1022  atmos_sw_phy_ae, &
1023  atmos_sw_phy_ch, &
1024  atmos_sw_phy_rd, &
1025  atmos_sw_phy_sf, &
1026  atmos_sw_phy_tb, &
1027  atmos_sw_phy_bl, &
1028  atmos_sw_phy_cp, &
1029  atmos_sw_phy_lt
1030  use mod_atmos_dyn_vars, only: &
1031  atmos_dyn_vars_restart_read
1032  use mod_atmos_phy_mp_vars, only: &
1033  atmos_phy_mp_vars_restart_read
1034  use mod_atmos_phy_ae_vars, only: &
1035  atmos_phy_ae_vars_restart_read
1036  use mod_atmos_phy_ch_vars, only: &
1037  atmos_phy_ch_vars_restart_read
1038  use mod_atmos_phy_rd_vars, only: &
1039  atmos_phy_rd_vars_restart_read
1040  use mod_atmos_phy_sf_vars, only: &
1041  atmos_phy_sf_vars_restart_read
1042  use mod_atmos_phy_tb_vars, only: &
1043  atmos_phy_tb_vars_restart_read
1044  use mod_atmos_phy_bl_vars, only: &
1045  atmos_phy_bl_vars_restart_read
1046  use mod_atmos_phy_cp_vars, only: &
1047  atmos_phy_cp_vars_restart_read
1048  use mod_atmos_phy_lt_vars, only: &
1049  atmos_phy_lt_vars_restart_read
1050  use mod_cpl_admin, only: &
1051  cpl_sw
1052  implicit none
1053 
1054  integer :: i, j, iq
1055  !---------------------------------------------------------------------------
1056 
1057  if ( restart_fid /= -1 ) then
1058  log_newline
1059  log_info("ATMOS_vars_restart_read",*) 'Read from restart file (ATMOS) '
1060 
1061  call file_cartesc_read( restart_fid, pv_info(i_dens)%NAME, 'ZXY', & ! [IN]
1062  dens(:,:,:) ) ! [OUT]
1063  call file_cartesc_read( restart_fid, pv_info(i_momz)%NAME, 'ZHXY', & ! [IN]
1064  momz(:,:,:) ) ! [OUT]
1065  call file_cartesc_read( restart_fid, pv_info(i_momx)%NAME, 'ZXHY', & ! [IN]
1066  momx(:,:,:) ) ! [OUT]
1067  call file_cartesc_read( restart_fid, pv_info(i_momy)%NAME, 'ZXYH', & ! [IN]
1068  momy(:,:,:) ) ! [OUT]
1069  call file_cartesc_read( restart_fid, pv_info(i_rhot)%NAME, 'ZXY', & ! [IN]
1070  rhot(:,:,:) ) ! [OUT]
1071 
1072  do iq = 1, qa
1073  call file_cartesc_read( restart_fid, tracer_name(iq), 'ZXY', & ! [IN]
1074  qtrc(:,:,:,iq) ) ! [OUT]
1075  enddo
1076 
1077  if ( file_get_aggregate(restart_fid) ) then
1078  call file_cartesc_flush( restart_fid ) ! X/Y halos have been read from file
1079  ! fill k halos
1080  do j = 1, ja
1081  do i = 1, ia
1082  dens( 1:ks-1,i,j) = dens(ks,i,j)
1083  momz( 1:ks-2,i,j) = momz(ks-1,i,j)
1084  momx( 1:ks-1,i,j) = momx(ks,i,j)
1085  momy( 1:ks-1,i,j) = momy(ks,i,j)
1086  rhot( 1:ks-1,i,j) = rhot(ks,i,j)
1087  dens(ke+1:ka, i,j) = dens(ke,i,j)
1088  momz(ke+1:ka, i,j) = momz(ke,i,j)
1089  momx(ke+1:ka, i,j) = momx(ke,i,j)
1090  momy(ke+1:ka, i,j) = momy(ke,i,j)
1091  rhot(ke+1:ka, i,j) = rhot(ke,i,j)
1092  do iq = 1, qa
1093  qtrc( 1:ks-1,i,j,iq) = qtrc(ks,i,j,iq)
1094  qtrc(ke+1:ka ,i,j,iq) = qtrc(ke,i,j,iq)
1095  end do
1096  enddo
1097  enddo
1098  else
1099  call atmos_vars_fillhalo
1100  end if
1101 
1102  call atmos_vars_calc_diagnostics
1103  call atmos_vars_check( force = .true. )
1104  else
1105  log_error("ATMOS_vars_restart_read",*) 'invalid restart file ID for atmosphere. STOP!'
1106  call prc_abort
1107  endif
1108 
1109  if ( atmos_use_average ) then
1110  dens_av(:,:,:) = dens(:,:,:)
1111  momz_av(:,:,:) = momz(:,:,:)
1112  momx_av(:,:,:) = momx(:,:,:)
1113  momy_av(:,:,:) = momy(:,:,:)
1114  rhot_av(:,:,:) = rhot(:,:,:)
1115  qtrc_av(:,:,:,:) = qtrc(:,:,:,:)
1116  endif
1117 
1118  if ( atmos_sw_dyn ) call atmos_dyn_vars_restart_read
1119  if ( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_read
1120  if ( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_read
1121  if ( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_read
1122  if ( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_read
1123  if ( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_read
1124  if ( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_read
1125  if ( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_read
1126  if ( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_read
1127  if ( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_read
1128 
1129  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_read(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_read(), mod_atmos_phy_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_calc_diagnostics(), atmos_vars_check(), atmos_vars_fillhalo(), mod_cpl_admin::cpl_sw, dens, dens_av, scale_file_cartesc::file_cartesc_flush(), scale_file::file_get_aggregate(), scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, scale_atmos_grid_cartesc_index::ia, j, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momx_av, momy, momy_av, momz, momz_av, scale_prc::prc_abort(), scale_tracer::qa, qtrc, qtrc_av, rhot, rhot_av, and scale_tracer::tracer_name.

Referenced by mod_admin_restart::admin_restart_read().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_vars_history_setpres()

subroutine, public mod_atmos_vars::atmos_vars_history_setpres

Set pressure for history output.

Definition at line 1135 of file mod_atmos_vars.F90.

1135  use scale_atmos_grid_cartesc_real, only: &
1136  real_fz => atmos_grid_cartesc_real_fz
1137  use scale_atmos_bottom, only: &
1138  bottom_estimate => atmos_bottom_estimate
1139  use scale_file_history_cartesc, only: &
1140  file_history_cartesc_set_pres
1141  use mod_atmos_phy_sf_vars, only: &
1142  sfc_temp => atmos_phy_sf_sfc_temp
1143  implicit none
1144 
1145  real(RP) :: SFC_DENS(IA,JA)
1146  real(RP) :: SFC_PRES(IA,JA)
1147  !---------------------------------------------------------------------------
1148 
1149  call bottom_estimate( ka, ks, ke, ia, isb, ieb, ja, jsb, jeb, &
1150  dens_av(:,:,:), pres(:,:,:), qv(:,:,:), & ! [IN]
1151  sfc_temp(:,:), & ! [IN]
1152  real_fz(:,:,:), & ! [IN]
1153  sfc_dens(:,:), sfc_pres(:,:) ) ! [OUT]
1154 
1155  call file_history_cartesc_set_pres( phyd(:,:,:), & ! [IN]
1156  phydh(:,:,:), & ! [IN]
1157  sfc_pres(:,:) ) ! [IN]
1158 
1159  return

References scale_atmos_bottom::atmos_bottom_estimate(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, dens_av, scale_file_history_cartesc::file_history_cartesc_set_pres(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, phyd, phydh, pres, and qv.

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

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 1165 of file mod_atmos_vars.F90.

1165  use scale_prc, only: &
1166  prc_myrank
1167  use scale_file, only: &
1168  file_get_aggregate
1169  use scale_file_cartesc, only: &
1170  file_cartesc_open, &
1171  file_cartesc_read, &
1172  file_cartesc_flush, &
1173  file_cartesc_close
1174  implicit none
1175 
1176  real(RP) :: DENS_check(KA,IA,JA) ! Density [kg/m3]
1177  real(RP) :: MOMZ_check(KA,IA,JA) ! momentum z [kg/s/m2]
1178  real(RP) :: MOMX_check(KA,IA,JA) ! momentum x [kg/s/m2]
1179  real(RP) :: MOMY_check(KA,IA,JA) ! momentum y [kg/s/m2]
1180  real(RP) :: RHOT_check(KA,IA,JA) ! DENS * POTT [K*kg/m3]
1181  real(RP) :: QTRC_check(KA,IA,JA,QA) ! tracer mixing ratio [kg/kg]
1182 
1183  character(len=H_LONG) :: basename
1184 
1185  logical :: datacheck
1186  integer :: k, i, j, iq
1187  integer :: fid
1188  !---------------------------------------------------------------------------
1189 
1190  call prof_rapstart('Debug')
1191 
1192  log_info("ATMOS_vars_restart_check",*) 'Compare last Data with ', trim(atmos_restart_check_basename), 'on PE=', prc_myrank
1193  log_info("ATMOS_vars_restart_check",*) 'criterion = ', atmos_restart_check_criterion
1194  datacheck = .true.
1195 
1196  basename = atmos_restart_check_basename
1197 
1198  call file_cartesc_open( basename, fid )
1199 
1200  call file_cartesc_read( fid, 'DENS', 'ZXY' , dens_check(:,:,:) )
1201  call file_cartesc_read( fid, 'MOMZ', 'ZHXY', momz_check(:,:,:) )
1202  call file_cartesc_read( fid, 'MOMX', 'ZXHY', momx_check(:,:,:) )
1203  call file_cartesc_read( fid, 'MOMY', 'ZXYH', momy_check(:,:,:) )
1204  call file_cartesc_read( fid, 'RHOT', 'ZXY' , rhot_check(:,:,:) )
1205  do iq = 1, qa
1206  call file_cartesc_read( fid, tracer_name(iq), 'ZXY', qtrc_check(:,:,:,iq) )
1207  end do
1208  if ( file_get_aggregate(fid) ) call file_cartesc_flush( fid )
1209 
1210  call file_cartesc_close( fid ) ! [IN]
1211 
1212  do k = ks, ke
1213  do j = js, je
1214  do i = is, ie
1215  if ( abs( dens(k,i,j)-dens_check(k,i,j) ) > atmos_restart_check_criterion ) then
1216  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', dens(k,i,j)-dens_check(k,i,j)
1217  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'DENS'
1218  datacheck = .false.
1219  endif
1220  enddo
1221  enddo
1222  enddo
1223 
1224  do k = ks-1, ke
1225  do j = js, je
1226  do i = is, ie
1227  if ( abs( momz(k,i,j)-momz_check(k,i,j) ) > atmos_restart_check_criterion ) then
1228  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', momz(k,i,j)-momz_check(k,i,j)
1229  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMZ'
1230  datacheck = .false.
1231  endif
1232  enddo
1233  enddo
1234  enddo
1235 
1236  do k = ks, ke
1237  do j = js, je
1238  do i = is, ie
1239  if ( abs( momx(k,i,j)-momx_check(k,i,j) ) > atmos_restart_check_criterion ) then
1240  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', momx(k,i,j)-momx_check(k,i,j)
1241  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMX'
1242  datacheck = .false.
1243  endif
1244  enddo
1245  enddo
1246  enddo
1247 
1248  do k = ks, ke
1249  do j = js, je
1250  do i = is, ie
1251  if ( abs( momy(k,i,j)-momy_check(k,i,j) ) > atmos_restart_check_criterion ) then
1252  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', momy(k,i,j)-momy_check(k,i,j)
1253  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'MOMY'
1254  datacheck = .false.
1255  endif
1256  enddo
1257  enddo
1258  enddo
1259 
1260  do k = ks, ke
1261  do j = js, je
1262  do i = is, ie
1263  if ( abs( rhot(k,i,j)-rhot_check(k,i,j) ) > atmos_restart_check_criterion ) then
1264  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', rhot(k,i,j)-rhot_check(k,i,j)
1265  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, 'RHOT'
1266  datacheck = .false.
1267  endif
1268  enddo
1269  enddo
1270  enddo
1271 
1272  do iq = 1, qa
1273  do k = ks, ke
1274  do j = js, je
1275  do i = is, ie
1276  if ( abs( qtrc(k,i,j,iq)-qtrc_check(k,i,j,iq) ) > atmos_restart_check_criterion ) then
1277  log_error("ATMOS_vars_restart_check",*) 'there is the difference : ', qtrc(k,i,j,iq)-qtrc_check(k,i,j,iq)
1278  log_error_cont(*) 'at (PE-id,k,i,j,varname) : ', prc_myrank, k, i, j, tracer_name(iq)
1279  datacheck = .false.
1280  endif
1281  enddo
1282  enddo
1283  enddo
1284  enddo
1285 
1286  if (datacheck) then
1287  log_info("ATMOS_vars_restart_check",*) 'Data Check Clear.'
1288  else
1289  log_info("ATMOS_vars_restart_check",*) 'Data Check Failed. See std. output.'
1290  log_error("ATMOS_vars_restart_check",*) 'Data Check Failed.'
1291  endif
1292 
1293  call prof_rapend('Debug')
1294 
1295  return

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

Referenced by mod_rm_driver::rm_driver().

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 1301 of file mod_atmos_vars.F90.

1301  use scale_file_history, only: &
1302  file_history_query, &
1303  file_history_put
1304  use scale_atmos_grid_cartesc_real, only: &
1305  real_cz => atmos_grid_cartesc_real_cz
1306  use mod_atmos_phy_mp_vars, only: &
1307  atmos_phy_mp_vars_history
1308  use mod_atmos_phy_ae_vars, only: &
1309  atmos_phy_ae_vars_history
1310  implicit none
1311 
1312  logical :: do_put
1313  integer :: iq, iv
1314  !---------------------------------------------------------------------------
1315 
1316  call prof_rapstart('ATM_History', 1)
1317 
1318  ! history output of prognostic variables
1319  call file_history_put ( pv_hist_id(i_dens), dens(:,:,:) )
1320  call file_history_put ( pv_hist_id(i_momz), momz(:,:,:) )
1321  call file_history_put ( pv_hist_id(i_momx), momx(:,:,:) )
1322  call file_history_put ( pv_hist_id(i_momy), momy(:,:,:) )
1323  call file_history_put ( pv_hist_id(i_rhot), rhot(:,:,:) )
1324  do iq = 1, qa
1325  call file_history_put ( qp_hist_id(iq), qtrc(:,:,:,iq) )
1326  enddo
1327 
1328 
1329  ! history output of diagnostic variables
1330  call file_history_put ( dv_hist_id(i_w ), w(:,:,:) )
1331  call file_history_put ( dv_hist_id(i_u ), u(:,:,:) )
1332  call file_history_put ( dv_hist_id(i_v ), v(:,:,:) )
1333  call file_history_put ( dv_hist_id(i_pott ), pott(:,:,:) )
1334  call file_history_put ( dv_hist_id(i_temp ), temp(:,:,:) )
1335  call file_history_put ( dv_hist_id(i_pres ), pres(:,:,:) )
1336 
1337  call file_history_put ( dv_hist_id(i_exner), exner(:,:,:) )
1338  call file_history_put ( dv_hist_id(i_phyd ), phyd(:,:,:) )
1339 
1340  call file_history_put ( dv_hist_id(i_qdry ), qdry(:,:,:) )
1341  call file_history_put ( dv_hist_id(i_rtot ), rtot(:,:,:) )
1342  call file_history_put ( dv_hist_id(i_cvtot), cvtot(:,:,:) )
1343  call file_history_put ( dv_hist_id(i_cptot), cptot(:,:,:) )
1344 
1345  do iv = i_cptot+1, dv_nmax
1346  if ( dv_hist_id(iv) > 0 ) then
1347  call file_history_query( dv_hist_id(iv), do_put )
1348 
1349  if ( do_put ) then
1350  select case( dv_info(iv)%ndims )
1351  case( 3 )
1352  call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work3d(:,:,:) )
1353  call file_history_put( dv_hist_id(iv), work3d(:,:,:) )
1354  case( 2 )
1355  call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work2d(:,:) )
1356  call file_history_put( dv_hist_id(iv), work2d(:,:) )
1357  case( 1 )
1358  call atmos_vars_get_diagnostic( dv_info(iv)%NAME, work1d(:) )
1359  call file_history_put( dv_hist_id(iv), work1d(:) )
1360  end select
1361  endif
1362  endif
1363  enddo
1364 
1365  call file_history_put( hist_id_gph, real_cz(:,:,:) )
1366 
1367 
1368  if ( moist ) &
1369  call atmos_phy_mp_vars_history( dens_av(:,:,:), temp(:,:,:), qtrc_av(:,:,:,:) )
1370 ! if ( .false. ) then
1371 ! call ATMOS_vars_get_diagnostic( "RH", WORK3D(:,:,:) )
1372 ! call ATMOS_PHY_AE_vars_history( QTRC_av(:,:,:,:), WORK3D(:,:,:) )
1373 ! end if
1374 
1375  call prof_rapend ('ATM_History', 1)
1376 
1377  return

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

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

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_vars_check()

subroutine, public mod_atmos_vars::atmos_vars_check ( logical, intent(in), optional  force)

Check variables for atmosphere.

Definition at line 1383 of file mod_atmos_vars.F90.

1383  use scale_prc, only: &
1384  prc_myrank, &
1385  prc_abort
1386  use scale_prc_cartesc, only: &
1387  prc_twod
1388  use scale_const, only: &
1389  grav => const_grav, &
1390  cvdry => const_cvdry
1391  use scale_statistics, only: &
1392  statistics_checktotal, &
1393  statistics_total, &
1394  statistics_detail
1395  use scale_atmos_grid_cartesc_real, only: &
1396  atmos_grid_cartesc_real_vol, &
1397  atmos_grid_cartesc_real_totvol, &
1398  atmos_grid_cartesc_real_volwxy, &
1399  atmos_grid_cartesc_real_totvolwxy, &
1400  atmos_grid_cartesc_real_volzuy, &
1401  atmos_grid_cartesc_real_totvolzuy, &
1402  atmos_grid_cartesc_real_volzxv, &
1403  atmos_grid_cartesc_real_totvolzxv
1404  use mod_atmos_admin, only: &
1405  atmos_dyn_type
1406  use scale_atmos_grid_cartesc, only: &
1407  rfdx => atmos_grid_cartesc_rfdx, &
1408  rfdy => atmos_grid_cartesc_rfdy
1409  use scale_atmos_grid_cartesc_real, only: &
1410  real_cz => atmos_grid_cartesc_real_cz
1411  use scale_atmos_grid_cartesc_metric, only: &
1412  mapf => atmos_grid_cartesc_metric_mapf
1413  use scale_time, only: &
1414  time_dtsec_atmos_dyn
1415  implicit none
1416  logical, intent(in), optional :: force
1417 
1418  real(RP) :: RHOQ(KA,IA,JA)
1419 
1420  real(RP) :: WORK (KA,IA,JA,3)
1421  character(len=H_SHORT) :: WNAME(3)
1422  real(RP) :: CFLMAX
1423 
1424  integer :: k, i, j, iq
1425  logical :: check
1426  !---------------------------------------------------------------------------
1427 
1428  if ( present(force) ) then
1429  check = force
1430  else
1431  check = atmos_vars_checkrange
1432  end if
1433 
1434  ! value check for prognostic variables
1435  if ( check ) then
1436  call valcheck( ka, ks, ke, ia, is, ie, ja, js, je, &
1437  dens(:,:,:), 0.0_rp, 2.0_rp, pv_info(i_dens)%NAME, __file__, __line__ )
1438  call valcheck( ka, ks, ke, ia, is, ie, ja, js, je, &
1439  momz(:,:,:), -200.0_rp, 200.0_rp, pv_info(i_momz)%NAME, __file__, __line__ )
1440  call valcheck( ka, ks, ke, ia, is, ie, ja, js, je, &
1441  momx(:,:,:), -200.0_rp, 200.0_rp, pv_info(i_momx)%NAME, __file__, __line__ )
1442  call valcheck( ka, ks, ke, ia, is, ie, ja, js, je, &
1443  momy(:,:,:), -200.0_rp, 200.0_rp, pv_info(i_momy)%NAME, __file__, __line__ )
1444  call valcheck( ka, ks, ke, ia, is, ie, ja, js, je, &
1445  rhot(:,:,:), 0.0_rp, 1000.0_rp, pv_info(i_rhot)%NAME, __file__, __line__ )
1446 
1447  !$omp parallel workshare
1448 !OCL XFILL
1449  work(:,:,:,1) = w(:,:,:)
1450 !OCL XFILL
1451  work(:,:,:,2) = u(:,:,:)
1452 !OCL XFILL
1453  work(:,:,:,3) = v(:,:,:)
1454  !$omp end parallel workshare
1455 
1456  wname(1) = "W"
1457  wname(2) = "U"
1458  wname(3) = "V"
1459 
1460  call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
1461  wname(:), work(:,:,:,:) )
1462  endif
1463 
1464  if ( present(force) ) then
1465  check = force
1466  else
1467  check = statistics_checktotal
1468  end if
1469 
1470  if ( check ) then
1471 
1472  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1473  dens(:,:,:), pv_info(i_dens)%NAME, & ! (in)
1474  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1475  atmos_grid_cartesc_real_totvol ) ! (in)
1476  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1477  momz(:,:,:), pv_info(i_momz)%NAME, & ! (in)
1478  atmos_grid_cartesc_real_volwxy(:,:,:), & ! (in)
1479  atmos_grid_cartesc_real_totvolwxy ) ! (in)
1480  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1481  momx(:,:,:), pv_info(i_momx)%NAME, & ! (in)
1482  atmos_grid_cartesc_real_volzuy(:,:,:), & ! (in)
1483  atmos_grid_cartesc_real_totvolzuy ) ! (in)
1484  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1485  momy(:,:,:), pv_info(i_momy)%NAME, & ! (in)
1486  atmos_grid_cartesc_real_volzxv(:,:,:), & ! (in)
1487  atmos_grid_cartesc_real_totvolzxv ) ! (in)
1488  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1489  rhot(:,:,:), pv_info(i_rhot)%NAME, & ! (in)
1490  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1491  atmos_grid_cartesc_real_totvol ) ! (in)
1492 
1493  do iq = 1, qa
1494  rhoq(ks:ke,is:ie,js:je) = dens(ks:ke,is:ie,js:je) * qtrc(ks:ke,is:ie,js:je,iq)
1495 
1496  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1497  rhoq(:,:,:), tracer_name(iq), & ! (in)
1498  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1499  atmos_grid_cartesc_real_totvol ) ! (in)
1500  enddo
1501 
1502  call atmos_vars_calc_diagnostics
1503 
1504 
1505  rhoq(ks:ke,is:ie,js:je) = dens(ks:ke,is:ie,js:je) * qdry(ks:ke,is:ie,js:je)
1506  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1507  rhoq(:,:,:), 'QDRY', & ! (in)
1508  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1509  atmos_grid_cartesc_real_totvol ) ! (in)
1510 
1511  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
1512  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1513  rhoq(:,:,:), 'QTOT', & ! (in)
1514  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1515  atmos_grid_cartesc_real_totvol ) ! (in)
1516 
1517 
1518  call atmos_vars_get_diagnostic( 'ENGT', work3d(:,:,:) )
1519  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1520  work3d(:,:,:), 'ENGT', & ! (in)
1521  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1522  atmos_grid_cartesc_real_totvol ) ! (in)
1523  call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
1524  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1525  work3d(:,:,:), 'ENGP', & ! (in)
1526  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1527  atmos_grid_cartesc_real_totvol ) ! (in)
1528  call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
1529  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1530  work3d(:,:,:), 'ENGK', & ! (in)
1531  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1532  atmos_grid_cartesc_real_totvol ) ! (in)
1533  call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
1534  call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
1535  work3d(:,:,:), 'ENGI', & ! (in)
1536  atmos_grid_cartesc_real_vol(:,:,:), & ! (in)
1537  atmos_grid_cartesc_real_totvol ) ! (in)
1538 
1539  end if
1540 
1541 
1542  ! CFL condition check
1543  if ( ( atmos_dyn_type /= 'OFF' .AND. atmos_dyn_type /= 'NONE' ) &
1544  .AND. ( atmos_vars_checkcfl_soft > 0.0_rp .OR. atmos_vars_checkcfl_hard > 0.0_rp ) ) then
1545  !$omp parallel workshare
1546 !OCL XFILL
1547  work(:,:,:,:) = 0.0_rp
1548  !$omp end parallel workshare
1549 
1550  !$omp parallel do
1551  do j = js, je
1552  do i = is, ie
1553  do k = ks, ke
1554  work(k,i,j,1) = 0.5_rp * abs(momz_av(k,i,j)) / ( dens_av(k+1,i,j) + dens_av(k,i,j) ) &
1555  * time_dtsec_atmos_dyn / ( real_cz(k+1,i,j) - real_cz(k,i,j) )
1556  work(k,i,j,3) = 0.5_rp * abs(momy_av(k,i,j)) / ( dens_av(k,i,j+1) + dens_av(k,i,j) ) &
1557  * time_dtsec_atmos_dyn * rfdy(j) * mapf(i,j,2,i_xv)
1558  enddo
1559  enddo
1560  enddo
1561  if ( prc_twod ) then
1562  !$omp parallel do
1563  do j = js, je
1564  do k = ks, ke
1565  work(k,is,j,2) = 0.0_rp
1566  enddo
1567  enddo
1568  else
1569  !$omp parallel do
1570  do j = js, je
1571  do i = is, ie
1572  do k = ks, ke
1573  work(k,i,j,2) = 0.5_rp * abs(momx_av(k,i,j)) / ( dens_av(k,i+1,j) + dens_av(k,i,j) ) &
1574  * time_dtsec_atmos_dyn * rfdx(i) * mapf(i,j,1,i_uy)
1575  enddo
1576  enddo
1577  enddo
1578  end if
1579 
1580  cflmax = maxval( work(:,:,:,:) )
1581 
1582  if ( atmos_vars_checkcfl_hard > 0.0_rp .AND. cflmax > atmos_vars_checkcfl_hard ) then
1583  log_info("ATMOS_vars_check",*) "Courant number =", cflmax, " exceeded the hard limit =", atmos_vars_checkcfl_hard
1584  log_error("ATMOS_vars_check",*) "Courant number =", cflmax, " exceeded the hard limit =", atmos_vars_checkcfl_hard
1585  log_error_cont(*) "Rank =", prc_myrank
1586  log_error_cont(*) "Please set ATMOS_VARS_CHECKCFL_HARD in the namelist PARAM_ATMOS_VARS when you want to change the limit."
1587 
1588  wname(1) = "Courant num. Z"
1589  wname(2) = "Courant num. X"
1590  wname(3) = "Courant num. Y"
1591  call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
1592  wname(:), work(:,:,:,:), &
1593  local=.true. )
1594 
1595  call prc_abort
1596  endif
1597 
1598  if ( atmos_vars_checkcfl_soft > 0.0_rp .AND. cflmax > atmos_vars_checkcfl_soft ) then
1599  log_info("ATMOS_vars_check",*) "Courant number =", cflmax, " exceeded the soft limit =", atmos_vars_checkcfl_soft
1600  log_error("ATMOS_vars_check",*) "Courant number =", cflmax, " exceeded the soft limit =", atmos_vars_checkcfl_soft
1601  log_error_cont(*) "Rank =", prc_myrank
1602 
1603  wname(1) = "Courant num. Z"
1604  wname(2) = "Courant num. X"
1605  wname(3) = "Courant num. Y"
1606  call statistics_detail( ka, ks, ke, ia, is, ie, ja, js, je, 3, &
1607  wname(:), work(:,:,:,:), &
1608  local=.true. )
1609  endif
1610 
1611  endif
1612 
1613 
1614  return

References mod_atmos_admin::atmos_dyn_type, scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_mapf, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolwxy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolzuy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolzxv, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volwxy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volzuy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volzxv, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdy, atmos_vars_calc_diagnostics(), scale_debug::check(), scale_const::const_cvdry, scale_const::const_grav, dens, dens_av, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, scale_atmos_grid_cartesc_index::i_uy, scale_atmos_grid_cartesc_index::i_xv, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, j, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momx_av, momy, momy_av, momz, momz_av, scale_prc::prc_abort(), scale_prc::prc_myrank, scale_prc_cartesc::prc_twod, scale_tracer::qa, qdry, qtrc, rhot, scale_statistics::statistics_checktotal, scale_time::time_dtsec_atmos_dyn, scale_tracer::tracer_name, u, v, and w.

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

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_vars_calc_diagnostics()

subroutine, public mod_atmos_vars::atmos_vars_calc_diagnostics

Calc diagnostic variables.

Definition at line 1620 of file mod_atmos_vars.F90.

1620  use scale_atmos_grid_cartesc_real, only: &
1621  real_cz => atmos_grid_cartesc_real_cz, &
1622  real_fz => atmos_grid_cartesc_real_fz
1623  use scale_atmos_thermodyn, only: &
1624  atmos_thermodyn_specific_heat
1625  use scale_atmos_diagnostic, only: &
1626  atmos_diagnostic_get_therm_rhot, &
1627  atmos_diagnostic_get_phyd
1628  use scale_atmos_diagnostic_cartesc, only: &
1629  atmos_diagnostic_cartesc_get_vel
1630  use scale_comm_cartesc, only: &
1631  comm_vars8, &
1632  comm_wait
1633  use scale_atmos_hydrometeor, only: &
1634  n_hyd
1635  use mod_atmos_phy_mp_vars, only: &
1636  atmos_phy_mp_vars_get_diagnostic, &
1637  atmos_phy_mp_vars_reset_diagnostics
1638  use mod_atmos_phy_ae_vars, only: &
1639  atmos_phy_ae_vars_reset_diagnostics
1640  implicit none
1641 
1642  integer :: iq
1643 
1644  call atmos_thermodyn_specific_heat( &
1645  ka, ks, ke, ia, 1, ia, ja, 1, ja, qa, &
1646  qtrc_av(:,:,:,:), & ! (in)
1647  tracer_mass(:), tracer_r(:), tracer_cv(:), tracer_cp(:), & ! (in)
1648  qdry(:,:,:), rtot(:,:,:), cvtot(:,:,:), cptot(:,:,:) ) ! (out)
1649 
1650  call atmos_diagnostic_cartesc_get_vel( &
1651  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1652  dens_av(:,:,:), momz_av(:,:,:), momx_av(:,:,:), momy_av(:,:,:), & ! (in)
1653  w(:,:,:), u(:,:,:), v(:,:,:) ) ! (out)
1654 
1655  call atmos_diagnostic_get_therm_rhot( &
1656  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1657  dens_av(:,:,:), rhot_av(:,:,:), & ! (in)
1658  rtot(:,:,:), cvtot(:,:,:), cptot(:,:,:), & ! (in)
1659  pott(:,:,:), temp(:,:,:), pres(:,:,:), exner(:,:,:) ) ! (out)
1660 
1661  call atmos_diagnostic_get_phyd( &
1662  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1663  dens_av(:,:,:), pres(:,:,:), & ! (in)
1664  real_cz(:,:,:), real_fz(:,:,:), & ! (in)
1665  phyd(:,:,:), phydh(:,:,:) ) ! (out)
1666 
1667  call atmos_phy_mp_vars_reset_diagnostics
1668  call atmos_phy_ae_vars_reset_diagnostics
1669 
1670  if ( moist ) then
1671  call atmos_phy_mp_vars_get_diagnostic( &
1672  dens_av(:,:,:), temp(:,:,:), qtrc_av(:,:,:,:), & ! (in)
1673  qe=qe(:,:,:,:) ) ! (out)
1674  do iq = 1, n_hyd
1675  call comm_vars8(qe(:,:,:,iq), iq)
1676  end do
1677  do iq = 1, n_hyd
1678  call comm_wait (qe(:,:,:,iq), iq)
1679  end do
1680  end if
1681 
1682  ! reset diagnostic variables
1683  dv_calculated(:) = .false.
1684 
1685  return

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

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

Here is the call graph for this function:
Here is the caller graph for this function:

◆ atmos_vars_get_diagnostic_3d()

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

get diagnostic variable 3D

Definition at line 1693 of file mod_atmos_vars.F90.

1693  use scale_const, only: &
1694  grav => const_grav, &
1695  rvap => const_rvap, &
1696  cpdry => const_cpdry, &
1697  cvdry => const_cvdry
1698  use scale_prc, only: &
1699  prc_abort
1700  use scale_prc_cartesc, only: &
1701  prc_twod
1702  use scale_atmos_grid_cartesc, only: &
1703  rcdx => atmos_grid_cartesc_rcdx, &
1704  rcdy => atmos_grid_cartesc_rcdy
1705  use scale_atmos_grid_cartesc_real, only: &
1706  real_cz => atmos_grid_cartesc_real_cz, &
1707  real_fz => atmos_grid_cartesc_real_fz, &
1708  f2h => atmos_grid_cartesc_real_f2h
1709  use scale_atmos_grid_cartesc_metric, only: &
1710  rotc => atmos_grid_cartesc_metric_rotc
1711  use scale_comm_cartesc, only: &
1712  comm_vars8, &
1713  comm_wait
1714  use scale_atmos_hydrometeor, only: &
1715  lhvc => lhv, &
1716  lhfc => lhf, &
1717  atmos_hydrometeor_lhv, &
1718  atmos_hydrometeor_lhf, &
1719  atmos_hydrometeor_lhs
1720  use scale_atmos_saturation, only: &
1721  atmos_saturation_dens2qsat_all, &
1722  atmos_saturation_psat_all, &
1723  atmos_saturation_psat_liq, &
1724  atmos_saturation_psat_ice, &
1725  atmos_saturation_tdew_liq, &
1726  atmos_saturation_pote
1727  use scale_atmos_diagnostic, only: &
1728  atmos_diagnostic_get_potv, &
1729  atmos_diagnostic_get_teml, &
1730  atmos_diagnostic_get_n2
1731  implicit none
1732  character(len=*), intent(in) :: vname
1733  real(RP), intent(out) :: var(:,:,:)
1734 
1735  real(RP) :: UH (KA,IA,JA)
1736  real(RP) :: VH (KA,IA,JA)
1737 
1738  real(RP) :: WORK(KA,IA,JA)
1739 
1740  integer :: k, i, j, iq
1741 
1742  select case ( vname )
1743  case ( 'W' )
1744  var(:,:,:) = w(:,:,:)
1745 
1746  case ( 'U' )
1747  var(:,:,:) = u(:,:,:)
1748 
1749  case ( 'V' )
1750  var(:,:,:) = v(:,:,:)
1751 
1752  case ( 'PT' )
1753  var(:,:,:) = pott(:,:,:)
1754 
1755  case ( 'T' )
1756  var(:,:,:) = temp(:,:,:)
1757 
1758  case ( 'EXNER' )
1759  var(:,:,:) = exner(:,:,:)
1760 
1761  case ( 'PHYD' )
1762  var(:,:,:) = phyd(:,:,:)
1763 
1764  case ( 'QDRY' )
1765  var(:,:,:) = qdry(:,:,:)
1766 
1767  case ( 'RTOT' )
1768  var(:,:,:) = rtot(:,:,:)
1769 
1770  case ( 'CVTOT' )
1771  var(:,:,:) = cvtot(:,:,:)
1772 
1773  case ( 'CPTOT' )
1774  var(:,:,:) = cptot(:,:,:)
1775 
1776  case ( 'LHV' )
1777  if ( .not. dv_calculated(i_lhv) ) then
1778  call allocate_3d( lhv )
1779  call atmos_hydrometeor_lhv( &
1780  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1781  temp(:,:,:), & ! (in)
1782  lhv(:,:,:) ) ! (out)
1783  dv_calculated(i_lhv) = .true.
1784  end if
1785  var(ks:ke,:,:) = lhv(ks:ke,:,:)
1786 
1787  case ( 'LHS' )
1788  if ( .not. dv_calculated(i_lhs) ) then
1789  call allocate_3d( lhs )
1790  call atmos_hydrometeor_lhs( &
1791  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1792  temp(:,:,:), & ! (in)
1793  lhs(:,:,:) ) ! (out)
1794  dv_calculated(i_lhs) = .true.
1795  end if
1796  var(ks:ke,:,:) = lhs(ks:ke,:,:)
1797 
1798  case ( 'LHF' )
1799  if ( .not. dv_calculated(i_lhf) ) then
1800  call allocate_3d( lhf )
1801  call atmos_hydrometeor_lhf( &
1802  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1803  temp(:,:,:), & ! (in)
1804  lhf(:,:,:) ) ! (out)
1805  dv_calculated(i_lhf) = .true.
1806  end if
1807  var(ks:ke,:,:) = lhf(ks:ke,:,:)
1808 
1809  case ( 'POTV' )
1810  if ( .not. dv_calculated(i_potv) ) then
1811  call allocate_3d( potv )
1812  call atmos_diagnostic_get_potv( &
1813  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1814  pott(:,:,:), rtot(:,:,:), & ! (in)
1815  potv(:,:,:) ) ! (out)
1816  dv_calculated(i_potv) = .true.
1817  end if
1818  var(ks:ke,:,:) = potv(ks:ke,:,:)
1819 
1820  case ( 'TEML' )
1821  if ( .not. dv_calculated(i_teml) ) then
1822  call allocate_3d( teml )
1823  call atmos_vars_get_diagnostic( 'LHV', work3d(:,:,:) )
1824  call atmos_vars_get_diagnostic( 'LHS', work3d(:,:,:) )
1825  call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
1826  call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
1827  call atmos_diagnostic_get_teml( &
1828  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1829  temp(:,:,:), lhv(:,:,:), lhs(:,:,:), & ! (in)
1830  qc(:,:,:), qi(:,:,:), cptot(:,:,:), & ! (in)
1831  teml(:,:,:) ) ! (out)
1832  dv_calculated(i_teml) = .true.
1833  end if
1834  var(ks:ke,:,:) = teml(ks:ke,:,:)
1835 
1836  case ( 'POTL' )
1837  if ( .not. dv_calculated(i_potl) ) then
1838  call allocate_3d( potl )
1839  call atmos_vars_get_diagnostic( 'TEML', work3d(:,:,:) )
1840 !OCL XFILL
1841  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1842  !$omp private(i,j,k) &
1843  !$omp shared(POTL,TEML,EXNER) &
1844  !$omp shared(KS,KE,IA,JA)
1845  do j = 1, ja
1846  do i = 1, ia
1847  do k = ks, ke
1848  potl(k,i,j) = teml(k,i,j) / exner(k,i,j)
1849  enddo
1850  enddo
1851  enddo
1852  dv_calculated(i_potl) = .true.
1853  end if
1854  var(ks:ke,:,:) = potl(ks:ke,:,:)
1855 
1856  case ( 'POTE' )
1857  if ( .not. dv_calculated(i_pote) ) then
1858  call allocate_3d( pote )
1859  call atmos_saturation_pote( &
1860  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1861  dens(:,:,:), pott(:,:,:), temp(:,:,:), qv(:,:,:), & ! [IN]
1862  pote(:,:,:) ) ! [OUT]
1863  end if
1864  var(ks:ke,:,:) = pote(ks:ke,:,:)
1865  case ( 'QTOT' )
1866  if ( .not. dv_calculated(i_qtot) ) then
1867  call allocate_3d( qtot )
1868  if ( moist ) then
1869  call atmos_vars_get_diagnostic( 'QHYD', work3d(:,:,:) )
1870 !OCL XFILL
1871  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1872  !$omp private(i,j,k) &
1873  !$omp shared(QTOT,QV,QHYD) &
1874  !$omp shared(KS,KE,IA,JA)
1875  do j = 1, ja
1876  do i = 1, ia
1877  do k = ks, ke
1878  qtot(k,i,j) = qv(k,i,j) + qhyd(k,i,j)
1879  enddo
1880  enddo
1881  enddo
1882  else
1883 !OCL XFILL
1884  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1885  !$omp private(i,j,k) &
1886  !$omp shared(QTOT) &
1887  !$omp shared(KS,KE,IA,JA)
1888  do j = 1, ja
1889  do i = 1, ia
1890  do k = ks, ke
1891  qtot(k,i,j) = 0.0_rp
1892  enddo
1893  enddo
1894  enddo
1895  end if
1896  dv_calculated(i_qtot) = .true.
1897  end if
1898  var(ks:ke,:,:) = qtot(ks:ke,:,:)
1899 
1900  case ( 'QHYD' )
1901  if ( .not. dv_calculated(i_qhyd) ) then
1902  call allocate_3d( qhyd )
1903  if ( moist ) then
1904  call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
1905  call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
1906 !OCL XFILL
1907  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1908  !$omp private(i,j,k) &
1909  !$omp shared(QHYD,QLIQ,QICE) &
1910  !$omp shared(KS,KE,IA,JA)
1911  do j = 1, ja
1912  do i = 1, ia
1913  do k = ks, ke
1914  qhyd(k,i,j) = qliq(k,i,j) + qice(k,i,j)
1915  enddo
1916  enddo
1917  enddo
1918  else
1919 !OCL XFILL
1920  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1921  !$omp private(i,j,k) &
1922  !$omp shared(QHYD) &
1923  !$omp shared(KS,KE,IA,JA)
1924  do j = 1, ja
1925  do i = 1, ia
1926  do k = ks, ke
1927  qhyd(k,i,j) = 0.0_rp
1928  enddo
1929  enddo
1930  enddo
1931  end if
1932  dv_calculated(i_qhyd) = .true.
1933  end if
1934  var(ks:ke,:,:) = qhyd(ks:ke,:,:)
1935 
1936  case ( 'QLIQ' )
1937  if ( .not. dv_calculated(i_qliq) ) then
1938  call allocate_3d( qliq )
1939 !OCL XFILL
1940  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1941  !$omp private(i,j,k) &
1942  !$omp shared(QLIQ,QC,QR) &
1943  !$omp shared(KS,KE,IA,JA)
1944  do j = 1, ja
1945  do i = 1, ia
1946  do k = ks, ke
1947  qliq(k,i,j) = qc(k,i,j) + qr(k,i,j)
1948  enddo
1949  enddo
1950  enddo
1951  dv_calculated(i_qliq) = .true.
1952  end if
1953  var(ks:ke,:,:) = qliq(ks:ke,:,:)
1954 
1955  case ( 'QICE' )
1956  if ( .not. dv_calculated(i_qice) ) then
1957  call allocate_3d( qice )
1958 !OCL XFILL
1959  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1960  !$omp private(i,j,k) &
1961  !$omp shared(QICE,QI,QS,QG,QH) &
1962  !$omp shared(KS,KE,IA,JA)
1963  do j = 1, ja
1964  do i = 1, ia
1965  do k = ks, ke
1966  qice(k,i,j) = qi(k,i,j) + qs(k,i,j) + qg(k,i,j) + qh(k,i,j)
1967  enddo
1968  enddo
1969  enddo
1970  dv_calculated(i_qice) = .true.
1971  end if
1972  var(ks:ke,:,:) = qice(ks:ke,:,:)
1973 
1974  case ( 'QSAT' )
1975  if ( .not. dv_calculated(i_qsat) ) then
1976  call allocate_3d( qsat )
1977  call atmos_saturation_dens2qsat_all( &
1978  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1979  temp(:,:,:), dens_av(:,:,:), & ! (in)
1980  qsat(:,:,:) ) ! (out)
1981  dv_calculated(i_qsat) = .true.
1982  end if
1983  var(ks:ke,:,:) = qsat(ks:ke,:,:)
1984 
1985  case ( 'RHA' )
1986  if ( .not. dv_calculated(i_rha) ) then
1987  call allocate_3d( rha )
1988  if ( moist ) then
1989  call atmos_saturation_psat_all( &
1990  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
1991  temp(:,:,:), & ! (in)
1992  work(:,:,:) ) ! (out)
1993 !OCL XFILL
1994  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
1995  !$omp private(i,j,k) &
1996  !$omp shared(RHA,DENS_av,QV,WORK,TEMP) &
1997  !$omp shared(KS,KE,IA,JA)
1998  do j = 1, ja
1999  do i = 1, ia
2000  do k = ks, ke
2001  rha(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
2002  / work(k,i,j) * rvap * temp(k,i,j) &
2003  * 100.0_rp
2004  enddo
2005  enddo
2006  enddo
2007  else
2008 !OCL XFILL
2009  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2010  !$omp private(i,j,k) &
2011  !$omp shared(RHA) &
2012  !$omp shared(KS,KE,IA,JA)
2013  do j = 1, ja
2014  do i = 1, ia
2015  do k = ks, ke
2016  rha(k,i,j) = 0.0_rp
2017  enddo
2018  enddo
2019  enddo
2020  end if
2021  dv_calculated(i_rha) = .true.
2022  end if
2023  var(ks:ke,:,:) = rha(ks:ke,:,:)
2024 
2025  case ( 'RHL', 'RH' )
2026  if ( .not. dv_calculated(i_rhl) ) then
2027  call allocate_3d( rhl )
2028  if ( moist ) then
2029  call atmos_saturation_psat_liq( &
2030  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
2031  temp(:,:,:), & ! (in)
2032  work(:,:,:) ) ! (out)
2033 !OCL XFILL
2034  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2035  !$omp private(i,j,k) &
2036  !$omp shared(RHL,DENS_av,QV,WORK,TEMP) &
2037  !$omp shared(KS,KE,IA,JA)
2038  do j = 1, ja
2039  do i = 1, ia
2040  do k = ks, ke
2041  rhl(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
2042  / work(k,i,j) * rvap * temp(k,i,j) &
2043  * 100.0_rp
2044  enddo
2045  enddo
2046  enddo
2047  else
2048 !OCL XFILL
2049  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2050  !$omp private(i,j,k) &
2051  !$omp shared(RHL) &
2052  !$omp shared(KS,KE,IA,JA)
2053  do j = 1, ja
2054  do i = 1, ia
2055  do k = ks, ke
2056  rhl(k,i,j) = 0.0_rp
2057  enddo
2058  enddo
2059  enddo
2060  end if
2061  dv_calculated(i_rhl) = .true.
2062  end if
2063  var(ks:ke,:,:) = rhl(ks:ke,:,:)
2064 
2065  case ( 'RHI' )
2066  if ( .not. dv_calculated(i_rhi) ) then
2067  call allocate_3d( rhi )
2068  if ( moist ) then
2069  call atmos_saturation_psat_ice( &
2070  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
2071  temp(:,:,:), & ! (int)
2072  work(:,:,:) ) ! (out)
2073 !OCL XFILL
2074  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2075  !$omp private(i,j,k) &
2076  !$omp shared(RHI,DENS_av,QV,WORK,TEMP) &
2077  !$omp shared(KS,KE,IA,JA)
2078  do j = 1, ja
2079  do i = 1, ia
2080  do k = ks, ke
2081  rhi(k,i,j) = dens_av(k,i,j) * qv(k,i,j) &
2082  / work(k,i,j) * rvap * temp(k,i,j) &
2083  * 100.0_rp
2084  enddo
2085  enddo
2086  enddo
2087  else
2088 !OCL XFILL
2089  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2090  !$omp private(i,j,k) &
2091  !$omp shared(RHI) &
2092  !$omp shared(KS,KE,IA,JA)
2093  do j = 1, ja
2094  do i = 1, ia
2095  do k = ks, ke
2096  rhi(k,i,j) = 0.0_rp
2097  enddo
2098  enddo
2099  enddo
2100  end if
2101  dv_calculated(i_rhi) = .true.
2102  end if
2103  var(ks:ke,:,:) = rhi(ks:ke,:,:)
2104 
2105  case ( 'VOR' )
2106  if ( .not. dv_calculated(i_vor) ) then
2107  call allocate_3d( vor )
2108  !!! to move to grid !!!
2109  ! at x, v, layer
2110  ! at u, y, layer
2111  if ( prc_twod ) then
2112  !$omp parallel do private(j,k) OMP_SCHEDULE_
2113 !OCL XFILL
2114  do j = 1, ja-1
2115  do k = ks, ke
2116  uh(k,is,j) = ( momx_av(k,is,j) + momx_av(k,is,j+1) ) &
2117  / ( dens_av(k,is,j) + dens_av(k,is,j+1) )
2118  enddo
2119  enddo
2120  !$omp parallel do private(j,k) OMP_SCHEDULE_
2121 !OCL XFILL
2122  do j = 2, ja-1
2123  do k = ks, ke
2124  vor(k,is,j) = - ( uh(k,is,j) - uh(k,is,j-1) ) * rcdy(j)
2125  enddo
2126  enddo
2127  !$omp parallel do private(k) OMP_SCHEDULE_
2128  do k = ks, ke
2129  vor(k,is,1 ) = vor(k,is,2 )
2130  vor(k,is,ja) = vor(k,is,ja-1)
2131  enddo
2132  else
2133  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2134 !OCL XFILL
2135  do j = 1, ja-1
2136  do i = 2, ia
2137  do k = ks, ke
2138  uh(k,i,j) = 0.5_rp * ( momx_av(k,i ,j) + momx_av(k,i ,j+1) &
2139  + momx_av(k,i-1,j) + momx_av(k,i-1,j+1) ) &
2140  / ( dens_av(k,i,j) + dens_av(k,i,j+1) )
2141  enddo
2142  enddo
2143  enddo
2144  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2145 !OCL XFILL
2146  do j = 2, ja
2147  do i = 1, ia-1
2148  do k = ks, ke
2149  vh(k,i,j) = 0.5_rp * ( momy_av(k,i,j ) + momy_av(k,i+1,j ) &
2150  + momy_av(k,i,j-1) + momy_av(k,i+1,j-1) ) &
2151  / ( dens_av(k,i,j) + dens_av(k,i+1,j) )
2152  enddo
2153  enddo
2154  enddo
2155  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2156 !OCL XFILL
2157  do j = 2, ja-1
2158  do i = 2, ia-1
2159  do k = ks, ke
2160  vor(k,i,j) = ( vh(k,i,j ) - vh(k,i-1,j ) ) * rcdx(i) &
2161  - ( uh(k,i ,j) - uh(k,i ,j-1) ) * rcdy(j)
2162  enddo
2163  enddo
2164  enddo
2165  !$omp parallel do private(j,k) OMP_SCHEDULE_
2166  do j = 1, ja
2167  do k = ks, ke
2168  vor(k,1 ,j) = vor(k,2 ,j)
2169  vor(k,ia,j) = vor(k,ia-1,j)
2170  enddo
2171  enddo
2172  !$omp parallel do private(i,k) OMP_SCHEDULE_
2173  do i = 1, ia
2174  do k = ks, ke
2175  vor(k,i,1 ) = vor(k,i,2 )
2176  vor(k,i,ja) = vor(k,i,ja-1)
2177  enddo
2178  enddo
2179  end if
2180  dv_calculated(i_vor) = .true.
2181  end if
2182  call comm_vars8( vor(:,:,:), 1 )
2183  call comm_wait ( vor(:,:,:), 1, .false. )
2184  var(ks:ke,:,:) = vor(ks:ke,:,:)
2185 
2186  case ( 'DIV' )
2187  if ( .not. dv_calculated(i_div) ) then
2188  call allocate_3d( div )
2189  call atmos_vars_get_diagnostic( 'HDIV', work3d(:,:,:) )
2190  !!!! to move to grid !!!!
2191  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2192 !OCL XFILL
2193  do j = 1, ja
2194  do i = 1, ia
2195  do k = ks, ke
2196  div(k,i,j) = ( momz_av(k,i,j) - momz_av(k-1,i ,j ) ) / ( real_fz(k,i,j)-real_fz(k-1,i,j) ) &
2197  + hdiv(k,i,j)
2198  enddo
2199  enddo
2200  enddo
2201  dv_calculated(i_div) = .true.
2202  end if
2203  var(ks:ke,:,:) = div(ks:ke,:,:)
2204 
2205  case ( 'HDIV' )
2206  if ( .not. dv_calculated(i_hdiv) ) then
2207  call allocate_3d( hdiv )
2208  !!!! to move to grid !!!!
2209  if ( prc_twod ) then
2210  !$omp parallel do private(j,k) OMP_SCHEDULE_
2211 !OCL XFILL
2212  do j = 2, ja
2213  do k = ks, ke
2214  hdiv(k,is,j) = ( momy_av(k,is,j) - momy_av(k ,is,j-1) ) * rcdy(j)
2215  enddo
2216  enddo
2217  !$omp parallel do private(k) OMP_SCHEDULE_
2218  do k = ks, ke
2219  hdiv(k,is,1) = hdiv(k,is,2)
2220  enddo
2221  else
2222  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2223 !OCL XFILL
2224  do j = 2, ja
2225  do i = 2, ia
2226  do k = ks, ke
2227  hdiv(k,i,j) = ( momx_av(k,i,j) - momx_av(k ,i-1,j ) ) * rcdx(i) &
2228  + ( momy_av(k,i,j) - momy_av(k ,i ,j-1) ) * rcdy(j)
2229  enddo
2230  enddo
2231  enddo
2232  !$omp parallel do private(i,k) OMP_SCHEDULE_
2233  do i = 1, ia
2234  do k = ks, ke
2235  hdiv(k,i,1) = hdiv(k,i,2)
2236  enddo
2237  enddo
2238  !$omp parallel do private(j,k) OMP_SCHEDULE_
2239  do j = 1, ja
2240  do k = ks, ke
2241  hdiv(k,1,j) = hdiv(k,2,j)
2242  enddo
2243  enddo
2244  end if
2245  call comm_vars8( hdiv(:,:,:), 1 )
2246  call comm_wait ( hdiv(:,:,:), 1, .false. )
2247  dv_calculated(i_hdiv) = .true.
2248  end if
2249  var(ks:ke,:,:) = hdiv(ks:ke,:,:)
2250 
2251  case ( 'Uabs' )
2252  if ( .not. dv_calculated(i_uabs) ) then
2253  call allocate_3d( uabs )
2254 !OCL XFILL
2255  !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
2256  do j = 1, ja
2257  do i = 1, ia
2258  do k = ks, ke
2259  uabs(k,i,j) = sqrt( w(k,i,j)**2 + u(k,i,j)**2 + v(k,i,j)**2 )
2260  enddo
2261  enddo
2262  enddo
2263  dv_calculated(i_uabs) = .true.
2264  end if
2265  var(ks:ke,:,:) = uabs(ks:ke,:,:)
2266 
2267  case ( 'N2' )
2268  if ( .not. dv_calculated(i_n2) ) then
2269  call allocate_3d( n2 )
2270  call atmos_diagnostic_get_n2( &
2271  ka, ks, ke, ia, 1, ia, ja, 1, ja, &
2272  pott(:,:,:), rtot(:,:,:), & !(in)
2273  real_cz(:,:,:), real_fz(:,:,:), & !(in)
2274  f2h(:,:,:,:), & !(in)
2275  n2(:,:,:) ) ! (out)
2276  dv_calculated(i_n2) = .true.
2277  end if
2278  var(ks:ke,:,:) = n2(ks:ke,:,:)
2279 
2280  case ( 'MSE' )
2281  if ( .not. dv_calculated(i_mse) ) then
2282  call allocate_3d( mse )
2283  call atmos_vars_get_diagnostic( 'LHV', work3d(:,:,:) )
2284 !OCL XFILL
2285  !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
2286  do j = 1, ja
2287  do i = 1, ia
2288  do k = ks, ke
2289  mse(k,i,j) = cptot(k,i,j) * temp(k,i,j) &
2290  + grav * ( real_cz(k,i,j) - real_fz(ks-1,i,j) ) &
2291  + lhv(k,i,j) * qv(k,i,j)
2292  enddo
2293  enddo
2294  enddo
2295  dv_calculated(i_mse) = .true.
2296  end if
2297  var(ks:ke,:,:) = mse(ks:ke,:,:)
2298 
2299  case ( 'TDEW' )
2300  if ( .not. dv_calculated(i_tdew) ) then
2301  call allocate_3d( tdew )
2302  call atmos_saturation_tdew_liq( ka, ks, ke, ia, 1, ia, ja, 1, ja, &
2303  dens(:,:,:), temp(:,:,:), qv(:,:,:), & ! [IN]
2304  tdew(:,:,:) ) ! [OUT]
2305  dv_calculated(i_tdew) = .true.
2306  end if
2307  var(ks:ke,:,:) = tdew(ks:ke,:,:)
2308 
2309  case ( 'ENGP' )
2310  if ( .not. dv_calculated(i_engp) ) then
2311  call allocate_3d( engp )
2312  !$omp parallel do private(k,i,j) OMP_SCHEDULE_ collapse(2)
2313  do j = 1, ja
2314  do i = 1, ia
2315  do k = ks, ke
2316  engp(k,i,j) = dens_av(k,i,j) * grav * real_cz(k,i,j)
2317  end do
2318  end do
2319  end do
2320  dv_calculated(i_engp) = .true.
2321  end if
2322  var(ks:ke,:,:) = engp(ks:ke,:,:)
2323 
2324  case ( 'ENGK' )
2325  if ( .not. dv_calculated(i_engk) ) then
2326  call allocate_3d( engk )
2327  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2328  do j = 1, ja
2329  do i = 1, ia
2330  do k = ks, ke
2331  engk(k,i,j) = 0.5_rp * dens_av(k,i,j) &
2332  * ( w(k,i,j)**2 + u(k,i,j)**2 + v(k,i,j)**2 )
2333  end do
2334  end do
2335  end do
2336  dv_calculated(i_engk) = .true.
2337  end if
2338  var(ks:ke,:,:) = engk(ks:ke,:,:)
2339 
2340  case ( 'ENGI' )
2341  if ( .not. dv_calculated(i_engi) ) then
2342  call allocate_3d( engi )
2343  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2344  do j = 1, ja
2345  do i = 1, ia
2346  do k = ks, ke
2347  engi(k,i,j) = dens_av(k,i,j) * qdry(k,i,j) * temp(k,i,j) * cvdry
2348  do iq = 1, qa
2349  engi(k,i,j) = engi(k,i,j) &
2350  + dens_av(k,i,j) * qtrc_av(k,i,j,iq) * ( temp(k,i,j) * tracer_cv(iq) + tracer_engi0(iq) )
2351  enddo
2352  end do
2353  end do
2354  end do
2355  dv_calculated(i_engi) = .true.
2356  end if
2357  var(ks:ke,:,:) = engi(ks:ke,:,:)
2358 
2359  case ( 'ENGT' )
2360  if ( .not. dv_calculated(i_engt) ) then
2361  call allocate_3d( engt )
2362  call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
2363  call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
2364  call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
2365  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2366  do j = 1, ja
2367  do i = 1, ia
2368  do k = ks, ke
2369  engt(k,i,j) = engp(k,i,j) + engk(k,i,j) + engi(k,i,j)
2370  enddo
2371  enddo
2372  enddo
2373  dv_calculated(i_engt) = .true.
2374  end if
2375  var(ks:ke,:,:) = engt(ks:ke,:,:)
2376 
2377  case ( 'DENS_PRIM' )
2378  if ( .not. dv_calculated(i_dens_prim) ) then
2379  call allocate_3d( dens_prim )
2380  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2381 !OCL XFILL
2382  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2383  do j = 1, ja
2384  do i = 1, ia
2385  do k = ks, ke
2386  dens_prim(k,i,j) = dens_av(k,i,j) - dens_mean(k)
2387  enddo
2388  enddo
2389  enddo
2390  dv_calculated(i_dens_prim) = .true.
2391  end if
2392  var(ks:ke,:,:) = dens_prim(ks:ke,:,:)
2393 
2394  case ( 'W_PRIM' )
2395  if ( .not. dv_calculated(i_w_prim) ) then
2396  call allocate_3d( w_prim )
2397  call atmos_vars_get_diagnostic( 'W_MEAN', work1d(:) )
2398 !OCL XFILL
2399  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2400  do j = 1, ja
2401  do i = 1, ia
2402  do k = ks, ke
2403  w_prim(k,i,j) = w(k,i,j) - w_mean(k)
2404  enddo
2405  enddo
2406  enddo
2407  dv_calculated(i_w_prim) = .true.
2408  end if
2409  var(ks:ke,:,:) = w_prim(ks:ke,:,:)
2410 
2411  case ( 'U_PRIM' )
2412  if ( .not. dv_calculated(i_u_prim) ) then
2413  call allocate_3d( u_prim )
2414  call atmos_vars_get_diagnostic( 'U_MEAN', work1d(:) )
2415 !OCL XFILL
2416  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2417  do j = 1, ja
2418  do i = 1, ia
2419  do k = ks, ke
2420  u_prim(k,i,j) = u(k,i,j) - u_mean(k)
2421  enddo
2422  enddo
2423  enddo
2424  dv_calculated(i_u_prim) = .true.
2425  end if
2426  var(ks:ke,:,:) = u_prim(ks:ke,:,:)
2427 
2428  case ( 'V_PRIM' )
2429  if ( .not. dv_calculated(i_v_prim) ) then
2430  call allocate_3d( v_prim )
2431  call atmos_vars_get_diagnostic( 'V_MEAN', work1d(:) )
2432 !OCL XFILL
2433  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2434  do j = 1, ja
2435  do i = 1, ia
2436  do k = ks, ke
2437  v_prim(k,i,j) = v(k,i,j) - v_mean(k)
2438  enddo
2439  enddo
2440  enddo
2441  dv_calculated(i_v_prim) = .true.
2442  end if
2443  var(ks:ke,:,:) = v_prim(ks:ke,:,:)
2444 
2445  case ( 'PT_PRIM' )
2446  if ( .not. dv_calculated(i_pt_prim) ) then
2447  call allocate_3d( pt_prim )
2448  call atmos_vars_get_diagnostic( 'PT_MEAN', work1d(:) )
2449 !OCL XFILL
2450  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2451  do j = 1, ja
2452  do i = 1, ia
2453  do k = ks, ke
2454  pt_prim(k,i,j) = pott(k,i,j) - pt_mean(k)
2455  enddo
2456  enddo
2457  enddo
2458  dv_calculated(i_pt_prim) = .true.
2459  end if
2460  var(ks:ke,:,:) = pt_prim(ks:ke,:,:)
2461 
2462  case ( 'W_PRIM2' )
2463  if ( .not. dv_calculated(i_w_prim2) ) then
2464  call allocate_3d( w_prim2 )
2465  call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
2466 !OCL XFILL
2467  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2468  do j = 1, ja
2469  do i = 1, ia
2470  do k = ks, ke
2471  w_prim2(k,i,j) = w_prim(k,i,j)**2
2472  enddo
2473  enddo
2474  enddo
2475  dv_calculated(i_w_prim2) = .true.
2476  end if
2477  var(ks:ke,:,:) = w_prim2(ks:ke,:,:)
2478 
2479  case ( 'PT_W_PRIM' )
2480  if ( .not. dv_calculated(i_pt_w_prim) ) then
2481  call allocate_3d( pt_w_prim )
2482  call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
2483  call atmos_vars_get_diagnostic( 'PT_PRIM', work3d(:,:,:) )
2484  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2485  do j = 1, ja
2486  do i = 1, ia
2487  do k = ks, ke
2488  pt_w_prim(k,i,j) = w_prim(k,i,j) * pt_prim(k,i,j) * dens_av(k,i,j) * cpdry
2489  enddo
2490  enddo
2491  enddo
2492  dv_calculated(i_pt_w_prim) = .true.
2493  end if
2494  var(ks:ke,:,:) = pt_w_prim(ks:ke,:,:)
2495 
2496  case ( 'W_PRIM3' )
2497  if ( .not. dv_calculated(i_w_prim3) ) then
2498  call allocate_3d( w_prim3 )
2499  call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
2500 !OCL XFILL
2501  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2502  do j = 1, ja
2503  do i = 1, ia
2504  do k = ks, ke
2505  w_prim3(k,i,j) = w_prim(k,i,j)**3
2506  enddo
2507  enddo
2508  enddo
2509  dv_calculated(i_w_prim3) = .true.
2510  end if
2511  var(ks:ke,:,:) = w_prim3(ks:ke,:,:)
2512 
2513  case ( 'TKE_RS' )
2514  if ( .not. dv_calculated(i_tke_rs) ) then
2515  call allocate_3d( tke_rs )
2516  call atmos_vars_get_diagnostic( 'W_PRIM', work3d(:,:,:) )
2517  call atmos_vars_get_diagnostic( 'U_PRIM', work3d(:,:,:) )
2518  call atmos_vars_get_diagnostic( 'V_PRIM', work3d(:,:,:) )
2519 !OCL XFILL
2520  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2521  do j = 1, ja
2522  do i = 1, ia
2523  do k = ks, ke
2524  tke_rs(k,i,j) = 0.5_rp * ( w_prim(k,i,j)**2 + u_prim(k,i,j)**2 + v_prim(k,i,j)**2 )
2525  enddo
2526  enddo
2527  enddo
2528  dv_calculated(i_tke_rs) = .true.
2529  end if
2530  var(ks:ke,:,:) = tke_rs(ks:ke,:,:)
2531 
2532  case ( 'VELZ' )
2533  if ( .not. dv_calculated(i_velz) ) then
2534  call allocate_3d( velz )
2535 !OCL XFILL
2536  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2537  do j = 1, ja
2538  do i = 1, ia
2539  velz(ks-1,i,j) = 0.0_rp
2540  do k = ks, ke-1
2541  velz(k,i,j) = momz(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k+1,i,j) )
2542  end do
2543  velz(ke,i,j) = 0.0_rp
2544  enddo
2545  enddo
2546  dv_calculated(i_velz) = .true.
2547  end if
2548  var(ks-1:ke,:,:) = velz(ks-1:ke,:,:)
2549 
2550  case ( 'VELX' )
2551  if ( .not. dv_calculated(i_velx) ) then
2552  call allocate_3d( velx )
2553  if ( prc_twod ) then
2554 !OCL XFILL
2555  !$omp parallel do private(j,k) OMP_SCHEDULE_
2556  do j = 1, ja
2557  do k = ks, ke
2558  velx(k,is,j) = momx(k,is,j) / dens(k,is,j)
2559  enddo
2560  enddo
2561  else
2562  !OCL XFILL
2563  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2564  do j = 1, ja
2565  do i = 1, ia-1
2566  do k = ks, ke
2567  velx(k,i,j) = momx(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k,i+1,j) )
2568  enddo
2569  enddo
2570  enddo
2571 !OCL XFILL
2572  !$omp parallel do private(j,k) OMP_SCHEDULE_
2573  do j = 1, ja
2574  do k = ks, ke
2575  velx(k,ia,j) = momx(k,ia,j) / dens(k,ia,j)
2576  enddo
2577  enddo
2578  call comm_vars8( velx(:,:,:), 1 )
2579  call comm_wait ( velx(:,:,:), 1, .false. )
2580  end if
2581  dv_calculated(i_velx) = .true.
2582  end if
2583  var(ks:ke,:,:) = velx(ks:ke,:,:)
2584 
2585  case ( 'VELY' )
2586  if ( .not. dv_calculated(i_vely) ) then
2587  call allocate_3d( vely )
2588 !OCL XFILL
2589  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2590  do j = 1, ja-1
2591  do i = 1, ia
2592  do k = ks, ke
2593  vely(k,i,j) = momy(k,i,j) * 2.0_rp / ( dens(k,i,j) + dens(k,i,j+1) )
2594  enddo
2595  enddo
2596  enddo
2597 !OCL XFILL
2598  !$omp parallel do private(i,k) OMP_SCHEDULE_
2599  do i = 1, ia
2600  do k = ks, ke
2601  vely(k,i,ja) = momy(k,i,ja) / dens(k,i,ja)
2602  enddo
2603  enddo
2604  call comm_vars8( vely(:,:,:), 1 )
2605  call comm_wait ( vely(:,:,:), 1, .false. )
2606  dv_calculated(i_vely) = .true.
2607  end if
2608  var(ks:ke,:,:) = vely(ks:ke,:,:)
2609 
2610  case ( 'Umet' )
2611  if ( .not. dv_calculated(i_umet) ) then
2612  call allocate_3d( umet )
2613 !OCL XFILL
2614  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2615  do j = 1, ja
2616  do i = 1, ia
2617  do k = ks, ke
2618  umet(k,i,j) = u(k,i,j) * rotc(i,j,1) - v(k,i,j) * rotc(i,j,2)
2619  end do
2620  end do
2621  end do
2622  dv_calculated(i_umet) = .true.
2623  end if
2624  var(ks:ke,:,:) = umet(ks:ke,:,:)
2625 
2626  case ( 'Vmet' )
2627  if ( .not. dv_calculated(i_vmet) ) then
2628  call allocate_3d( vmet )
2629 !OCL XFILL
2630  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2631  do j = 1, ja
2632  do i = 1, ia
2633  do k = ks, ke
2634  vmet(k,i,j) = u(k,i,j) * rotc(i,j,2) + v(k,i,j) * rotc(i,j,1)
2635  end do
2636  end do
2637  end do
2638  dv_calculated(i_vmet) = .true.
2639  end if
2640  var(ks:ke,:,:) = vmet(ks:ke,:,:)
2641 
2642  case default
2643  log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_3D: ', trim(vname)
2644  call prc_abort
2645  end select
2646 
2647 
2648  return

References allocate_3d(), scale_atmos_diagnostic::atmos_diagnostic_get_n2(), scale_atmos_diagnostic::atmos_diagnostic_get_potv(), scale_atmos_diagnostic::atmos_diagnostic_get_teml(), scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdy, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_f2h, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_const::const_cpdry, scale_const::const_cvdry, scale_const::const_grav, scale_const::const_rvap, cptot, cvtot, dens, dens_av, exner, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::is, j, scale_atmos_grid_cartesc_index::ja, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_atmos_hydrometeor::lhf, scale_atmos_hydrometeor::lhv, momx, momx_av, momy, momy_av, momz, momz_av, phyd, pott, scale_prc::prc_abort(), scale_prc_cartesc::prc_twod, scale_tracer::qa, qc, qdry, qg, qh, qi, qr, qs, qtrc_av, qv, rtot, temp, scale_tracer::tracer_cv, scale_tracer::tracer_engi0, u, v, and w.

Here is the call graph for this function:

◆ atmos_vars_get_diagnostic_2d()

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

get diagnostic variable 2D

Definition at line 2656 of file mod_atmos_vars.F90.

2656  use scale_prc, only: &
2657  prc_abort
2658  use scale_atmos_grid_cartesc_real, only: &
2659  real_cz => atmos_grid_cartesc_real_cz, &
2660  real_fz => atmos_grid_cartesc_real_fz
2661  use scale_atmos_adiabat, only: &
2662  atmos_adiabat_cape
2663  use mod_atmos_phy_mp_vars, only: &
2664  sflx_rain_mp => atmos_phy_mp_sflx_rain, &
2665  sflx_snow_mp => atmos_phy_mp_sflx_snow
2666  use mod_atmos_phy_cp_vars, only: &
2667  sflx_rain_cp => atmos_phy_cp_sflx_rain, &
2668  sflx_snow_cp => atmos_phy_cp_sflx_snow
2669  implicit none
2670 
2671  character(len=*), intent(in) :: vname
2672  real(RP), intent(out) :: var(:,:)
2673 
2674  real(RP) :: fact
2675  integer :: k, i, j
2676  !---------------------------------------------------------------------------
2677 
2678  select case ( vname )
2679  case ( 'LWP' )
2680  if ( .not. dv_calculated(i_lwp) ) then
2681  call allocate_2d( lwp )
2682  call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
2683  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2684  !$omp private(i,j,k) &
2685  !$omp shared(LWP,QLIQ,DENS_av,REAL_FZ) &
2686  !$omp shared(KS,KE,IA,JA)
2687  do j = 1, ja
2688  do i = 1, ia
2689  lwp(i,j) = 0.0_rp
2690  do k = ks, ke
2691  lwp(i,j) = lwp(i,j) &
2692  + qliq(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
2693  enddo
2694  enddo
2695  enddo
2696  dv_calculated(i_lwp) = .true.
2697  end if
2698  var(:,:) = lwp(:,:)
2699 
2700  case ( 'IWP' )
2701  if ( .not. dv_calculated(i_iwp) ) then
2702  call allocate_2d( iwp )
2703  call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
2704  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2705  !$omp private(i,j,k) &
2706  !$omp shared(IWP,QICE,DENS_av,REAL_FZ) &
2707  !$omp shared(KS,KE,IA,JA)
2708  do j = 1, ja
2709  do i = 1, ia
2710  iwp(i,j) = 0.0_rp
2711  do k = ks, ke
2712  iwp(i,j) = iwp(i,j) &
2713  + qice(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
2714  enddo
2715  enddo
2716  enddo
2717  dv_calculated(i_iwp) = .true.
2718  end if
2719  var(:,:) = iwp(:,:)
2720 
2721  case ( 'PW' )
2722  if ( .not. dv_calculated(i_pw) ) then
2723  call allocate_2d( pw )
2724  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2725  !$omp private(i,j,k) &
2726  !$omp shared(PW,QV,DENS_av,REAL_FZ) &
2727  !$omp shared(KS,KE,IA,JA)
2728  do j = 1, ja
2729  do i = 1, ia
2730  pw(i,j) = 0.0_rp
2731  do k = ks, ke
2732  pw(i,j) = pw(i,j) &
2733  + qv(k,i,j) * dens_av(k,i,j) * ( real_fz(k,i,j)-real_fz(k-1,i,j) ) * 1.e3_rp ! [kg/m2->g/m2]
2734  enddo
2735  enddo
2736  enddo
2737  dv_calculated(i_pw) = .true.
2738  end if
2739  var(:,:) = pw(:,:)
2740 
2741  case ( 'PBLH' )
2742  if ( .not. dv_calculated(i_pblh) ) then
2743  call allocate_2d( pblh )
2744  call atmos_vars_get_diagnostic( 'POTV', work3d(:,:,:) )
2745  !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
2746  !$omp private(k,i,j) &
2747  !$omp private(fact) &
2748  !$omp shared(PBLH,POTV,REAL_CZ,REAL_FZ) &
2749  !$omp shared(KS,KE,IA,JA)
2750  do j = 1, ja
2751  do i = 1, ia
2752  pblh(i,j) = real_cz(ks,i,j) - real_fz(ks-1,i,j)
2753  do k = ks+1, ke
2754  if ( potv(k,i,j) > potv(ks,i,j) ) then
2755  fact = ( potv(ks,i,j) - potv(k-1,i,j) ) &
2756  / ( potv(k,i,j) - potv(k-1,i,j) )
2757  pblh(i,j) = real_cz(k-1,i,j) - real_fz(ks-1,i,j) &
2758  + fact * ( real_cz(k,i,j) - real_cz(k-1,i,j) )
2759 
2760  exit
2761  endif
2762  enddo
2763  enddo
2764  enddo
2765  dv_calculated(i_pblh) = .true.
2766  end if
2767  var(:,:) = pblh(:,:)
2768 
2769  case ( 'CAPE', 'CIN', 'LCL', 'LFC', 'LNB' )
2770  if ( .not. dv_calculated(i_cape) ) then
2771  call allocate_2d( cape )
2772  call allocate_2d( cin )
2773  call allocate_2d( lcl )
2774  call allocate_2d( lfc )
2775  call allocate_2d( lnb )
2776  call atmos_adiabat_cape( &
2777  ka, ks, ke, ia, is, ie, ja, js, je, &
2778  ks, & ! (in)
2779  dens_av(:,:,:), temp(:,:,:), pres(:,:,:), & ! (in)
2780  qv(:,:,:), qc(:,:,:), qdry(:,:,:), & ! (in)
2781  rtot(:,:,:), cptot(:,:,:), & ! (in)
2782  real_cz(:,:,:), real_fz(:,:,:), & ! (in)
2783  cape(:,:), cin(:,:), lcl(:,:), lfc(:,:), lnb(:,:) ) ! (out)
2784  dv_calculated(i_cape) = .true.
2785  end if
2786  select case ( vname )
2787  case ( 'CAPE' )
2788  !$omp parallel do private(i,j) OMP_SCHEDULE_
2789  do j = js, je
2790  do i = is, ie
2791  var(i,j) = cape(i,j)
2792  end do
2793  end do
2794  case ( 'CIN' )
2795  !$omp parallel do private(i,j) OMP_SCHEDULE_
2796  do j = js, je
2797  do i = is, ie
2798  var(i,j) = cin(i,j)
2799  end do
2800  end do
2801  case ( 'LCL' )
2802  !$omp parallel do private(i,j) OMP_SCHEDULE_
2803  do j = js, je
2804  do i = is, ie
2805  var(i,j) = lcl(i,j)
2806  end do
2807  end do
2808  case ( 'LFC' )
2809  !$omp parallel do private(i,j) OMP_SCHEDULE_
2810  do j = js, je
2811  do i = is, ie
2812  var(i,j) = lfc(i,j)
2813  end do
2814  end do
2815  case ( 'LNB' )
2816  !$omp parallel do private(i,j) OMP_SCHEDULE_
2817  do j = js, je
2818  do i = is, ie
2819  var(i,j) = lnb(i,j)
2820  end do
2821  end do
2822  end select
2823 
2824  case ( 'PREC' )
2825  !$omp parallel do private(i,j) OMP_SCHEDULE_
2826  do j = js, je
2827  do i = is, ie
2828  var(i,j) = prec(i,j)
2829  end do
2830  end do
2831 
2832  case ( 'RAIN', 'SNOW' )
2833  if ( .not. dv_calculated(i_rain) ) then
2834  call allocate_2d( rain )
2835  call allocate_2d( snow )
2836  !$omp parallel do private(i,j) OMP_SCHEDULE_
2837  do j = js, je
2838  do i = is, ie
2839  rain(i,j) = sflx_rain_mp(i,j) + sflx_rain_cp(i,j)
2840  snow(i,j) = sflx_snow_mp(i,j) + sflx_snow_cp(i,j)
2841  enddo
2842  enddo
2843  dv_calculated(i_rain) = .true.
2844  end if
2845  select case (vname)
2846  case ( 'RAIN' )
2847  !$omp parallel do private(i,j) OMP_SCHEDULE_
2848  do j = js, je
2849  do i = is, ie
2850  var(i,j) = rain(i,j)
2851  end do
2852  end do
2853  case ( 'SNOW' )
2854  !$omp parallel do private(i,j) OMP_SCHEDULE_
2855  do j = js, je
2856  do i = is, ie
2857  var(i,j) = snow(i,j)
2858  end do
2859  end do
2860  end select
2861 
2862  case default
2863  log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_2D: ', trim(vname)
2864  call prc_abort
2865  end select
2866 
2867 
2868  return

References allocate_2d(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_rain, mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_snow, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, cptot, dens_av, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, j, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_prc::prc_abort(), prec, pres, qc, qdry, qv, rtot, and temp.

Here is the call graph for this function:

◆ atmos_vars_get_diagnostic_1d()

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

get diagnostic variable 1D

Definition at line 2876 of file mod_atmos_vars.F90.

2876  use scale_const, only: &
2877  cpdry => const_cpdry
2878  use scale_prc, only: &
2879  prc_abort
2880  use scale_statistics, only: &
2881  statistics_horizontal_mean
2882  use scale_atmos_grid_cartesc_real, only: &
2883  area => atmos_grid_cartesc_real_area
2884  implicit none
2885 
2886  character(len=*), intent(in) :: vname
2887  real(RP), intent(out) :: var(:)
2888 
2889  real(RP) :: WORK(KA,IA,JA)
2890  integer :: k, i, j
2891  !---------------------------------------------------------------------------
2892 
2893  select case ( vname )
2894  case ( 'DENS_MEAN' )
2895  if ( .not. dv_calculated(i_dens_mean) ) then
2896  call allocate_1d( dens_mean )
2897  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2898  dens(:,:,:), area(:,:), dens_mean(:) )
2899  dv_calculated(i_dens_mean) = .true.
2900  end if
2901  var(:) = dens_mean(:)
2902 
2903  case ( 'W_MEAN' )
2904  if ( .not. dv_calculated(i_w_mean) ) then
2905  call allocate_1d( w_mean )
2906  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2907 !OCL XFILL
2908  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2909  do j = jsb, jeb
2910  do i = isb, ieb
2911  do k = ks, ke
2912  work(k,i,j) = w(k,i,j) * dens_av(k,i,j)
2913  enddo
2914  enddo
2915  enddo
2916  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2917  work(:,:,:), area(:,:), w_mean(:) )
2918  do k = ks, ke
2919  w_mean(k) = w_mean(k) / dens_mean(k)
2920  enddo
2921  dv_calculated(i_w_mean) = .true.
2922  end if
2923  var(:) = w_mean(:)
2924 
2925  case ( 'U_MEAN' )
2926  if ( .not. dv_calculated(i_u_mean) ) then
2927  call allocate_1d( u_mean )
2928  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2929 !OCL XFILL
2930  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2931  do j = jsb, jeb
2932  do i = isb, ieb
2933  do k = ks, ke
2934  work(k,i,j) = u(k,i,j) * dens_av(k,i,j)
2935  enddo
2936  enddo
2937  enddo
2938  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2939  work(:,:,:), area(:,:), u_mean(:) )
2940  do k = ks, ke
2941  u_mean(k) = u_mean(k) / dens_mean(k)
2942  enddo
2943  dv_calculated(i_u_mean) = .true.
2944  end if
2945  var(:) = u_mean(:)
2946 
2947  case ( 'V_MEAN' )
2948  if ( .not. dv_calculated(i_v_mean) ) then
2949  call allocate_1d( v_mean )
2950  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2951 !OCL XFILL
2952  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2953  do j = jsb, jeb
2954  do i = isb, ieb
2955  do k = ks, ke
2956  work(k,i,j) = v(k,i,j) * dens_av(k,i,j)
2957  enddo
2958  enddo
2959  enddo
2960  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2961  work(:,:,:), area(:,:), v_mean(:) )
2962  do k = ks, ke
2963  v_mean(k) = v_mean(k) / dens_mean(k)
2964  enddo
2965  dv_calculated(i_v_mean) = .true.
2966  end if
2967  var(:) = v_mean(:)
2968 
2969  case ( 'PT_MEAN' )
2970  if ( .not. dv_calculated(i_pt_mean) ) then
2971  call allocate_1d( pt_mean )
2972  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2973  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2974  rhot(:,:,:), area(:,:), pt_mean(:) )
2975  do k = ks, ke
2976  pt_mean(k) = pt_mean(k) / dens_mean(k)
2977  enddo
2978  dv_calculated(i_pt_mean) = .true.
2979  end if
2980  var(:) = pt_mean(:)
2981 
2982  case ( 'T_MEAN' )
2983  if ( .not. dv_calculated(i_t_mean) ) then
2984  call allocate_1d( t_mean )
2985  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
2986 !OCL XFILL
2987  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
2988  do j = jsb, jeb
2989  do i = isb, ieb
2990  do k = ks, ke
2991  work(k,i,j) = temp(k,i,j) * dens_av(k,i,j)
2992  enddo
2993  enddo
2994  enddo
2995  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
2996  work(:,:,:), area(:,:), t_mean(:) )
2997  do k = ks, ke
2998  t_mean(k) = t_mean(k) / dens_mean(k)
2999  enddo
3000  dv_calculated(i_t_mean) = .true.
3001  end if
3002  var(:) = t_mean(:)
3003 
3004  case ( 'QV_MEAN' )
3005  if ( .not. dv_calculated(i_qv_mean) ) then
3006  call allocate_1d( qv_mean )
3007  if ( moist ) then
3008  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
3009 !OCL XFILL
3010  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3011  do j = jsb, jeb
3012  do i = isb, ieb
3013  do k = ks, ke
3014  work(k,i,j) = qv(k,i,j) * dens_av(k,i,j)
3015  enddo
3016  enddo
3017  enddo
3018  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
3019  work(:,:,:), area(:,:), qv_mean(:) )
3020  do k = ks, ke
3021  qv_mean(k) = qv_mean(k) / dens_mean(k)
3022  enddo
3023  else
3024  !$omp parallel do private(k) OMP_SCHEDULE_
3025  do k = ks, ke
3026  qv_mean(k) = 0.0_rp
3027  enddo
3028  end if
3029  dv_calculated(i_qv_mean) = .true.
3030  end if
3031  var(:) = qv_mean(:)
3032 
3033  case ( 'QHYD_MEAN' )
3034  if ( .not. dv_calculated(i_qhyd_mean) ) then
3035  call allocate_1d( qhyd_mean )
3036  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
3037  call atmos_vars_get_diagnostic( 'QHYD', work3d(:,:,:) )
3038 !OCL XFILL
3039  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3040  do j = jsb, jeb
3041  do i = isb, ieb
3042  do k = ks, ke
3043  work(k,i,j) = qhyd(k,i,j) * dens_av(k,i,j)
3044  enddo
3045  enddo
3046  enddo
3047  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
3048  work(:,:,:), area(:,:), qhyd_mean(:) )
3049  do k = ks, ke
3050  qhyd_mean(k) = qhyd_mean(k) / dens_mean(k)
3051  enddo
3052  dv_calculated(i_qhyd_mean) = .true.
3053  end if
3054  var(:) = qhyd_mean(:)
3055 
3056  case ( 'QLIQ_MEAN' )
3057  if ( .not. dv_calculated(i_qliq_mean) ) then
3058  call allocate_1d( qliq_mean )
3059  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
3060  call atmos_vars_get_diagnostic( 'QLIQ', work3d(:,:,:) )
3061 !OCL XFILL
3062  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3063  do j = jsb, jeb
3064  do i = isb, ieb
3065  do k = ks, ke
3066  work(k,i,j) = qliq(k,i,j) * dens_av(k,i,j)
3067  enddo
3068  enddo
3069  enddo
3070  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
3071  work(:,:,:), area(:,:), qliq_mean(:) )
3072  do k = ks, ke
3073  qliq_mean(k) = qliq_mean(k) / dens_mean(k)
3074  enddo
3075  dv_calculated(i_qliq_mean) = .true.
3076  end if
3077  var(:) = qliq_mean(:)
3078 
3079  case ( 'QICE_MEAN' )
3080  if ( .not. dv_calculated(i_qice_mean) ) then
3081  call allocate_1d( qice_mean )
3082  call atmos_vars_get_diagnostic( 'DENS_MEAN', work1d(:) )
3083  call atmos_vars_get_diagnostic( 'QICE', work3d(:,:,:) )
3084 !OCL XFILL
3085  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3086  do j = jsb, jeb
3087  do i = isb, ieb
3088  do k = ks, ke
3089  work(k,i,j) = qice(k,i,j) * dens_av(k,i,j)
3090  enddo
3091  enddo
3092  enddo
3093  call statistics_horizontal_mean( ka, ks, ke, ia, is, ie, ja, js, je, &
3094  work(:,:,:), area(:,:), qice_mean(:) )
3095  do k = ks, ke
3096  qice_mean(k) = qice_mean(k) / dens_mean(k)
3097  enddo
3098  dv_calculated(i_qice_mean) = .true.
3099  end if
3100  var(:) = qice_mean(:)
3101 
3102  case default
3103  log_error("ATMOS_vars_calc_diagnostics",*) 'name is invalid for ATMOS_vars_get_diagnostic_1D: ', trim(vname)
3104  call prc_abort
3105  end select
3106 
3107 
3108  return

References allocate_1d(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_area, scale_const::const_cpdry, dens, dens_av, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::isb, j, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::jsb, scale_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_prc::prc_abort(), qv, rhot, temp, u, v, and w.

Here is the call graph for this function:

◆ atmos_vars_monitor()

subroutine, public mod_atmos_vars::atmos_vars_monitor

monitor output

Definition at line 3114 of file mod_atmos_vars.F90.

3114  use scale_monitor, only: &
3115  monitor_put
3116  use scale_atmos_hydrometeor, only: &
3117  i_qv
3118  use mod_atmos_phy_rd_vars, only: &
3119  sflx_lw_up => atmos_phy_rd_sflx_lw_up, &
3120  sflx_lw_dn => atmos_phy_rd_sflx_lw_dn, &
3121  sflx_sw_up => atmos_phy_rd_sflx_sw_up, &
3122  sflx_sw_dn => atmos_phy_rd_sflx_sw_dn, &
3123  tomflx_lw_up => atmos_phy_rd_tomflx_lw_up, &
3124  tomflx_lw_dn => atmos_phy_rd_tomflx_lw_dn, &
3125  tomflx_sw_up => atmos_phy_rd_tomflx_sw_up, &
3126  tomflx_sw_dn => atmos_phy_rd_tomflx_sw_dn
3127  use mod_atmos_phy_sf_vars, only: &
3128  sflx_sh => atmos_phy_sf_sflx_sh, &
3129  sflx_lh => atmos_phy_sf_sflx_lh, &
3130  sflx_engi => atmos_phy_sf_sflx_engi, &
3131  sflx_qtrc => atmos_phy_sf_sflx_qtrc
3132  implicit none
3133 
3134  real(RP) :: RHOQ(KA,IA,JA)
3135 
3136  real(RP) :: ENGFLXT (IA,JA) ! total flux [J/m2/s]
3137  real(RP) :: SFLX_RD_net(IA,JA) ! net SFC radiation flux [J/m2/s]
3138  real(RP) :: TFLX_RD_net(IA,JA) ! net TOM radiation flux [J/m2/s]
3139 
3140  integer :: k, i, j, iq
3141  !---------------------------------------------------------------------------
3142 
3143  call monitor_put( pv_monit_id(i_dens), dens(:,:,:) )
3144  call monitor_put( pv_monit_id(i_momz), momz(:,:,:) )
3145  call monitor_put( pv_monit_id(i_momx), momx(:,:,:) )
3146  call monitor_put( pv_monit_id(i_momy), momy(:,:,:) )
3147  call monitor_put( pv_monit_id(i_rhot), rhot(:,:,:) )
3148 
3149  !##### Mass Budget #####
3150 
3151  do iq = 1, qa
3152  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3153 !OCL XFILL
3154  do j = js, je
3155  do i = is, ie
3156  do k = ks, ke
3157  rhoq(k,i,j) = dens_av(k,i,j) * qtrc_av(k,i,j,iq)
3158  enddo
3159  enddo
3160  enddo
3161 
3162  call monitor_put( qp_monit_id(iq), rhoq(:,:,:) )
3163  enddo
3164 
3165  ! total dry airmass
3166  if ( dv_monit_id(im_qdry) > 0 ) then
3167  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3168 !OCL XFILL
3169  do j = js, je
3170  do i = is, ie
3171  do k = ks, ke
3172  rhoq(k,i,j) = dens(k,i,j) * qdry(k,i,j)
3173  enddo
3174  enddo
3175  enddo
3176  call monitor_put( dv_monit_id(im_qdry), rhoq(:,:,:) )
3177  end if
3178 
3179  ! total vapor,liquid,solid tracers
3180  if ( dv_monit_id(im_qtot) > 0 ) then
3181  call atmos_vars_get_diagnostic( 'QTOT', work3d(:,:,:) )
3182  !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
3183 !OCL XFILL
3184  do j = js, je
3185  do i = is, ie
3186  do k = ks, ke
3187  rhoq(k,i,j) = dens(k,i,j) * qtot(k,i,j)
3188  enddo
3189  enddo
3190  enddo
3191  call monitor_put( dv_monit_id(im_qtot), rhoq(:,:,:) )
3192  end if
3193 
3194  ! total evapolation
3195  if ( moist ) call monitor_put( dv_monit_id(im_evap), sflx_qtrc(:,:,i_qv) )
3196 
3197  ! total precipitation
3198  if ( dv_monit_id(im_prec) > 0 ) then
3199  call monitor_put( dv_monit_id(im_prec), prec(:,:) )
3200  end if
3201 
3202 
3203  !##### Energy Budget #####
3204 
3205  if ( dv_monit_id(im_engt) > 0 ) then
3206  call atmos_vars_get_diagnostic( 'ENGT', work3d(:,:,:) )
3207  call monitor_put( dv_monit_id(im_engt), work3d(:,:,:) )
3208  end if
3209  if ( dv_monit_id(im_engp) > 0 ) then
3210  call atmos_vars_get_diagnostic( 'ENGP', work3d(:,:,:) )
3211  call monitor_put( dv_monit_id(im_engp), work3d(:,:,:) )
3212  end if
3213  if ( dv_monit_id(im_engk) > 0 ) then
3214  call atmos_vars_get_diagnostic( 'ENGK', work3d(:,:,:) )
3215  call monitor_put( dv_monit_id(im_engk), work3d(:,:,:) )
3216  end if
3217  if ( dv_monit_id(im_engi) > 0 ) then
3218  call atmos_vars_get_diagnostic( 'ENGI', work3d(:,:,:) )
3219  call monitor_put( dv_monit_id(im_engi), work3d(:,:,:) )
3220  end if
3221 
3222  ! radiation flux
3223 !OCL XFILL
3224  !$omp parallel do private(i,j) OMP_SCHEDULE_ collapse(2)
3225  do j = js, je
3226  do i = is, ie
3227  sflx_rd_net(i,j) = ( sflx_lw_up(i,j) - sflx_lw_dn(i,j) ) &
3228  + ( sflx_sw_up(i,j) - sflx_sw_dn(i,j) )
3229 
3230  tflx_rd_net(i,j) = ( tomflx_lw_up(i,j) - tomflx_lw_dn(i,j) ) &
3231  + ( tomflx_sw_up(i,j) - tomflx_sw_dn(i,j) )
3232 
3233  engflxt(i,j) = sflx_sh(i,j) &
3234  + sflx_engi(i,j) - prec_engi(i,j) &
3235  + sflx_rd_net(i,j) - tflx_rd_net(i,j)
3236  enddo
3237  enddo
3238 
3239  call monitor_put( dv_monit_id(im_engflxt), engflxt(:,:) )
3240 
3241  call monitor_put( dv_monit_id(im_engsfc_sh), sflx_sh(:,:) )
3242  call monitor_put( dv_monit_id(im_engsfc_lh), sflx_lh(:,:) )
3243  call monitor_put( dv_monit_id(im_engsfc_evap), sflx_engi(:,:) )
3244  call monitor_put( dv_monit_id(im_engsfc_prec), prec_engi(:,:) )
3245  call monitor_put( dv_monit_id(im_engsfc_rd), sflx_rd_net(:,:) )
3246  call monitor_put( dv_monit_id(im_engtom_rd), tflx_rd_net(:,:) )
3247 
3248  call monitor_put( dv_monit_id(im_engsfc_lw_up), sflx_lw_up(:,:) )
3249  call monitor_put( dv_monit_id(im_engsfc_lw_dn), sflx_lw_dn(:,:) )
3250  call monitor_put( dv_monit_id(im_engsfc_sw_up), sflx_sw_up(:,:) )
3251  call monitor_put( dv_monit_id(im_engsfc_sw_dn), sflx_sw_dn(:,:) )
3252 
3253  call monitor_put( dv_monit_id(im_engtom_lw_up), tomflx_lw_up(:,:) )
3254  call monitor_put( dv_monit_id(im_engtom_lw_dn), tomflx_lw_dn(:,:) )
3255  call monitor_put( dv_monit_id(im_engtom_sw_up), tomflx_sw_up(:,:) )
3256  call monitor_put( dv_monit_id(im_engtom_sw_dn), tomflx_sw_dn(:,:) )
3257 
3258  return

References 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_tomflx_lw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_lw_up, mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_sw_dn, mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_sw_up, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi, 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, dens, dens_av, scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_atmos_hydrometeor::i_qv, scale_index::i_rhot, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, j, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, momx, momy, momz, prec, prec_engi, scale_tracer::qa, qdry, qtrc_av, and rhot.

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

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 3264 of file mod_atmos_vars.F90.

3264  use scale_time, only: &
3265  time_gettimelabel
3266  use scale_file_cartesc, only: &
3267  file_cartesc_create
3268  use mod_atmos_admin, only: &
3269  atmos_sw_dyn, &
3270  atmos_sw_phy_mp, &
3271  atmos_sw_phy_ae, &
3272  atmos_sw_phy_ch, &
3273  atmos_sw_phy_rd, &
3274  atmos_sw_phy_sf, &
3275  atmos_sw_phy_tb, &
3276  atmos_sw_phy_bl, &
3277  atmos_sw_phy_cp, &
3278  atmos_sw_phy_lt
3279  use mod_cpl_admin, only: &
3280  cpl_sw
3281  use mod_atmos_dyn_vars, only: &
3282  atmos_dyn_vars_restart_create
3283  use mod_atmos_phy_mp_vars, only: &
3284  atmos_phy_mp_vars_restart_create
3285  use mod_atmos_phy_ae_vars, only: &
3286  atmos_phy_ae_vars_restart_create
3287  use mod_atmos_phy_ch_vars, only: &
3288  atmos_phy_ch_vars_restart_create
3289  use mod_atmos_phy_rd_vars, only: &
3290  atmos_phy_rd_vars_restart_create
3291  use mod_atmos_phy_sf_vars, only: &
3292  atmos_phy_sf_vars_restart_create
3293  use mod_atmos_phy_tb_vars, only: &
3294  atmos_phy_tb_vars_restart_create
3295  use mod_atmos_phy_bl_vars, only: &
3296  atmos_phy_bl_vars_restart_create
3297  use mod_atmos_phy_cp_vars, only: &
3298  atmos_phy_cp_vars_restart_create
3299  use mod_atmos_phy_lt_vars, only: &
3300  atmos_phy_lt_vars_restart_create
3301 #ifdef SDM
3302  use scale_atmos_phy_mp_sdm, only: &
3303  sd_rest_flg_out, &
3304  atmos_phy_mp_sdm_restart_create
3305  use scale_time, only: &
3306  nowdaysec => time_nowdaysec
3307 #endif
3308  implicit none
3309 
3310  character(len=19) :: timelabel
3311  character(len=H_LONG) :: basename
3312  !---------------------------------------------------------------------------
3313 
3314 #ifdef SDM
3315  if( sd_rest_flg_out ) then
3316  log_info("ATMOS_vars_restart_create",*) 'Output random number for SDM '
3317  call atmos_phy_mp_sdm_restart_create(nowdaysec)
3318  endif
3319 #endif
3320 
3321  if ( atmos_restart_out_basename /= '' ) then
3322 
3323  log_newline
3324  log_info("ATMOS_vars_restart_create",*) 'Create restart file (ATMOS) '
3325 
3326  if ( atmos_restart_out_postfix_timelabel ) then
3327  call time_gettimelabel( timelabel )
3328  basename = trim(atmos_restart_out_basename)//'_'//trim(timelabel)
3329  else
3330  basename = trim(atmos_restart_out_basename)
3331  endif
3332 
3333  log_info("ATMOS_vars_restart_create",*) 'basename: ', trim(basename)
3334 
3335  call file_cartesc_create( &
3336  basename, atmos_restart_out_title, atmos_restart_out_dtype, & ! [IN]
3337  restart_fid, & ! [OUT]
3338  aggregate=atmos_restart_out_aggregate ) ! [IN]
3339 
3340  allocate( pv_id(pv_nmax+qa) )
3341  endif
3342 
3343  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_create
3344  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_create
3345  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_create
3346  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_create
3347  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_create
3348  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_create
3349  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_create
3350  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_create
3351  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_create
3352  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_create
3353 
3354  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_create(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_create(), mod_atmos_phy_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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_aggregate, atmos_restart_out_basename, atmos_restart_out_dtype, atmos_restart_out_postfix_timelabel, atmos_restart_out_title, mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_cpl_admin::cpl_sw, scale_file_cartesc::file_cartesc_create(), scale_tracer::qa, scale_time::time_gettimelabel(), and scale_time::time_nowdaysec.

Referenced by mod_admin_restart::admin_restart_write().

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 3360 of file mod_atmos_vars.F90.

3360  use scale_file_cartesc, only: &
3361  file_cartesc_enddef
3362  use mod_atmos_admin, only: &
3363  atmos_sw_dyn, &
3364  atmos_sw_phy_mp, &
3365  atmos_sw_phy_ae, &
3366  atmos_sw_phy_ch, &
3367  atmos_sw_phy_rd, &
3368  atmos_sw_phy_sf, &
3369  atmos_sw_phy_tb, &
3370  atmos_sw_phy_bl, &
3371  atmos_sw_phy_cp, &
3372  atmos_sw_phy_lt
3373  use mod_cpl_admin, only: &
3374  cpl_sw
3375  use mod_atmos_dyn_vars, only: &
3376  atmos_dyn_vars_restart_enddef
3377  use mod_atmos_phy_mp_vars, only: &
3378  atmos_phy_mp_vars_restart_enddef
3379  use mod_atmos_phy_ae_vars, only: &
3380  atmos_phy_ae_vars_restart_enddef
3381  use mod_atmos_phy_ch_vars, only: &
3382  atmos_phy_ch_vars_restart_enddef
3383  use mod_atmos_phy_rd_vars, only: &
3384  atmos_phy_rd_vars_restart_enddef
3385  use mod_atmos_phy_sf_vars, only: &
3386  atmos_phy_sf_vars_restart_enddef
3387  use mod_atmos_phy_tb_vars, only: &
3388  atmos_phy_tb_vars_restart_enddef
3389  use mod_atmos_phy_bl_vars, only: &
3390  atmos_phy_bl_vars_restart_enddef
3391  use mod_atmos_phy_cp_vars, only: &
3392  atmos_phy_cp_vars_restart_enddef
3393  use mod_atmos_phy_lt_vars, only: &
3394  atmos_phy_lt_vars_restart_enddef
3395 #ifdef SDM
3396  use scale_atmos_phy_mp_sdm, only: &
3397  sd_rest_flg_out, &
3398  atmos_phy_mp_sdm_restart_enddef
3399 #endif
3400  implicit none
3401 
3402  !---------------------------------------------------------------------------
3403 
3404 #ifdef SDM
3405  if( sd_rest_flg_out ) then
3406  call atmos_phy_mp_sdm_restart_enddef
3407  endif
3408 #endif
3409 
3410  if ( restart_fid /= -1 ) then
3411  call file_cartesc_enddef( restart_fid ) ! [IN]
3412  endif
3413 
3414  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_enddef
3415  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_enddef
3416  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_enddef
3417  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_enddef
3418  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_enddef
3419  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_enddef
3420  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_enddef
3421  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_enddef
3422  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_enddef
3423  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_enddef
3424 
3425  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_enddef(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_enddef(), mod_atmos_phy_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_cpl_admin::cpl_sw, and scale_file_cartesc::file_cartesc_enddef().

Referenced by mod_admin_restart::admin_restart_write().

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 3431 of file mod_atmos_vars.F90.

3431  use scale_file_cartesc, only: &
3432  file_cartesc_close
3433  use mod_atmos_admin, only: &
3434  atmos_sw_dyn, &
3435  atmos_sw_phy_mp, &
3436  atmos_sw_phy_ae, &
3437  atmos_sw_phy_ch, &
3438  atmos_sw_phy_rd, &
3439  atmos_sw_phy_sf, &
3440  atmos_sw_phy_tb, &
3441  atmos_sw_phy_bl, &
3442  atmos_sw_phy_cp, &
3443  atmos_sw_phy_lt
3444  use mod_cpl_admin, only: &
3445  cpl_sw
3446  use mod_atmos_dyn_vars, only: &
3447  atmos_dyn_vars_restart_close
3448  use mod_atmos_phy_mp_vars, only: &
3449  atmos_phy_mp_vars_restart_close
3450  use mod_atmos_phy_ae_vars, only: &
3451  atmos_phy_ae_vars_restart_close
3452  use mod_atmos_phy_ch_vars, only: &
3453  atmos_phy_ch_vars_restart_close
3454  use mod_atmos_phy_rd_vars, only: &
3455  atmos_phy_rd_vars_restart_close
3456  use mod_atmos_phy_sf_vars, only: &
3457  atmos_phy_sf_vars_restart_close
3458  use mod_atmos_phy_tb_vars, only: &
3459  atmos_phy_tb_vars_restart_close
3460  use mod_atmos_phy_bl_vars, only: &
3461  atmos_phy_bl_vars_restart_close
3462  use mod_atmos_phy_cp_vars, only: &
3463  atmos_phy_cp_vars_restart_close
3464  use mod_atmos_phy_lt_vars, only: &
3465  atmos_phy_lt_vars_restart_close
3466 #ifdef SDM
3467  use scale_atmos_phy_mp_sdm, only: &
3468  sd_rest_flg_out, &
3469  atmos_phy_mp_sdm_restart_close
3470 #endif
3471  implicit none
3472  !---------------------------------------------------------------------------
3473 
3474 #ifdef SDM
3475  if( sd_rest_flg_out ) then
3476  call atmos_phy_mp_sdm_restart_close
3477  endif
3478 #endif
3479 
3480  if ( restart_fid /= -1 ) then
3481  log_newline
3482  log_info("ATMOS_vars_restart_close",*) 'Close restart file (ATMOS) '
3483 
3484  call file_cartesc_close( restart_fid ) ! [IN]
3485 
3486  restart_fid = -1
3487 
3488  if ( allocated(pv_id) ) deallocate( pv_id )
3489  endif
3490 
3491  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_close
3492  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_close
3493  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_close
3494  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_close
3495  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_close
3496  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_close
3497  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_close
3498  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_close
3499  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_close
3500  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_close
3501 
3502  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_close(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_close(), mod_atmos_phy_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_cpl_admin::cpl_sw, and scale_file_cartesc::file_cartesc_close().

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

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 3508 of file mod_atmos_vars.F90.

3508  use scale_file_cartesc, only: &
3509  file_cartesc_def_var
3510  use mod_atmos_admin, only: &
3511  atmos_sw_dyn, &
3512  atmos_sw_phy_mp, &
3513  atmos_sw_phy_ae, &
3514  atmos_sw_phy_ch, &
3515  atmos_sw_phy_rd, &
3516  atmos_sw_phy_sf, &
3517  atmos_sw_phy_tb, &
3518  atmos_sw_phy_bl, &
3519  atmos_sw_phy_cp, &
3520  atmos_sw_phy_lt
3521  use mod_cpl_admin, only: &
3522  cpl_sw
3523  use mod_atmos_dyn_vars, only: &
3524  atmos_dyn_vars_restart_def_var
3525  use mod_atmos_phy_mp_vars, only: &
3526  atmos_phy_mp_vars_restart_def_var
3527  use mod_atmos_phy_ae_vars, only: &
3528  atmos_phy_ae_vars_restart_def_var
3529  use mod_atmos_phy_ch_vars, only: &
3530  atmos_phy_ch_vars_restart_def_var
3531  use mod_atmos_phy_rd_vars, only: &
3532  atmos_phy_rd_vars_restart_def_var
3533  use mod_atmos_phy_sf_vars, only: &
3534  atmos_phy_sf_vars_restart_def_var
3535  use mod_atmos_phy_tb_vars, only: &
3536  atmos_phy_tb_vars_restart_def_var
3537  use mod_atmos_phy_bl_vars, only: &
3538  atmos_phy_bl_vars_restart_def_var
3539  use mod_atmos_phy_cp_vars, only: &
3540  atmos_phy_cp_vars_restart_def_var
3541  use mod_atmos_phy_lt_vars, only: &
3542  atmos_phy_lt_vars_restart_def_var
3543 #ifdef SDM
3544  use scale_atmos_phy_mp_sdm, only: &
3545  sd_rest_flg_out, &
3546  atmos_phy_mp_sdm_restart_def_var
3547 #endif
3548  implicit none
3549 
3550  integer iq
3551  !---------------------------------------------------------------------------
3552 
3553 #ifdef SDM
3554  if( sd_rest_flg_out ) then
3555  call atmos_phy_mp_sdm_restart_def_var
3556  endif
3557 #endif
3558 
3559  if ( restart_fid /= -1 ) then
3560 
3561  call file_cartesc_def_var( restart_fid, pv_info(i_dens)%NAME, pv_info(i_dens)%DESC, pv_info(i_dens)%UNIT, 'ZXY', atmos_restart_out_dtype, &
3562  pv_id(i_dens), &
3563  standard_name=pv_info(i_dens)%STDNAME )
3564  call file_cartesc_def_var( restart_fid, pv_info(i_momz)%NAME, pv_info(i_momz)%DESC, pv_info(i_momz)%UNIT, 'ZHXY', atmos_restart_out_dtype, &
3565  pv_id(i_momz), &
3566  standard_name=pv_info(i_momz)%STDNAME )
3567  call file_cartesc_def_var( restart_fid, pv_info(i_momx)%NAME, pv_info(i_momx)%DESC, pv_info(i_momx)%UNIT, 'ZXHY', atmos_restart_out_dtype, &
3568  pv_id(i_momx), &
3569  standard_name=pv_info(i_momx)%STDNAME )
3570  call file_cartesc_def_var( restart_fid, pv_info(i_momy)%NAME, pv_info(i_momy)%DESC, pv_info(i_momy)%UNIT, 'ZXYH', atmos_restart_out_dtype, &
3571  pv_id(i_momy), &
3572  standard_name=pv_info(i_momy)%STDNAME )
3573  call file_cartesc_def_var( restart_fid, pv_info(i_rhot)%NAME, pv_info(i_rhot)%DESC, pv_info(i_rhot)%UNIT, 'ZXY', atmos_restart_out_dtype, &
3574  pv_id(i_rhot), &
3575  standard_name=pv_info(i_rhot)%STDNAME )
3576  do iq = 1, qa
3577  call file_cartesc_def_var( restart_fid, tracer_name(iq), tracer_desc(iq), tracer_unit(iq), 'ZXY', atmos_restart_out_dtype, &
3578  pv_id(pv_nmax+iq) )
3579  enddo
3580 
3581  endif
3582 
3583  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_def_var
3584  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_def_var
3585  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_def_var
3586  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_def_var
3587  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_def_var
3588  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_def_var
3589  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_def_var
3590  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_def_var
3591  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_def_var
3592  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_def_var
3593 
3594  return

References 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_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_cpl_admin::cpl_sw, scale_file_cartesc::file_cartesc_def_var(), scale_index::i_dens, scale_index::i_momx, scale_index::i_momy, scale_index::i_momz, scale_index::i_rhot, scale_tracer::qa, scale_tracer::tracer_desc, scale_tracer::tracer_name, and scale_tracer::tracer_unit.

Referenced by mod_admin_restart::admin_restart_write().

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 3600 of file mod_atmos_vars.F90.

3600  use scale_file_cartesc, only: &
3601  file_cartesc_write_var
3602  use mod_atmos_admin, only: &
3603  atmos_sw_dyn, &
3604  atmos_sw_phy_mp, &
3605  atmos_sw_phy_ae, &
3606  atmos_sw_phy_ch, &
3607  atmos_sw_phy_rd, &
3608  atmos_sw_phy_sf, &
3609  atmos_sw_phy_tb, &
3610  atmos_sw_phy_bl, &
3611  atmos_sw_phy_cp, &
3612  atmos_sw_phy_lt
3613  use mod_cpl_admin, only: &
3614  cpl_sw
3615  use mod_atmos_dyn_vars, only: &
3616  atmos_dyn_vars_restart_write
3617  use mod_atmos_phy_mp_vars, only: &
3618  atmos_phy_mp_vars_restart_write
3619  use mod_atmos_phy_ae_vars, only: &
3620  atmos_phy_ae_vars_restart_write
3621  use mod_atmos_phy_ch_vars, only: &
3622  atmos_phy_ch_vars_restart_write
3623  use mod_atmos_phy_rd_vars, only: &
3624  atmos_phy_rd_vars_restart_write
3625  use mod_atmos_phy_sf_vars, only: &
3626  atmos_phy_sf_vars_restart_write
3627  use mod_atmos_phy_tb_vars, only: &
3628  atmos_phy_tb_vars_restart_write
3629  use mod_atmos_phy_bl_vars, only: &
3630  atmos_phy_bl_vars_restart_write
3631  use mod_atmos_phy_cp_vars, only: &
3632  atmos_phy_cp_vars_restart_write
3633  use mod_atmos_phy_lt_vars, only: &
3634  atmos_phy_lt_vars_restart_write
3635 #ifdef SDM
3636  use scale_atmos_phy_mp_sdm, only: &
3637  sd_rest_flg_out, &
3638  atmos_phy_mp_sdm_restart_write
3639 #endif
3640  implicit none
3641 
3642  integer iq
3643  !---------------------------------------------------------------------------
3644 
3645 #ifdef SDM
3646  if( sd_rest_flg_out ) then
3647  call atmos_phy_mp_sdm_restart_write
3648  endif
3649 #endif
3650 
3651  if ( restart_fid /= -1 ) then
3652 
3653  call atmos_vars_fillhalo
3654 
3655  call atmos_vars_check( force = .true. )
3656 
3657  call file_cartesc_write_var( restart_fid, pv_id(i_dens), dens(:,:,:), pv_info(i_dens)%NAME, 'ZXY' ) ! [IN]
3658  call file_cartesc_write_var( restart_fid, pv_id(i_momz), momz(:,:,:), pv_info(i_momz)%NAME, 'ZHXY' ) ! [IN]
3659  call file_cartesc_write_var( restart_fid, pv_id(i_momx), momx(:,:,:), pv_info(i_momx)%NAME, 'ZXHY' ) ! [IN]
3660  call file_cartesc_write_var( restart_fid, pv_id(i_momy), momy(:,:,:), pv_info(i_momy)%NAME, 'ZXYH' ) ! [IN]
3661  call file_cartesc_write_var( restart_fid, pv_id(i_rhot), rhot(:,:,:), pv_info(i_rhot)%NAME, 'ZXY' ) ! [IN]
3662 
3663  do iq = 1, qa
3664  call file_cartesc_write_var( restart_fid, pv_id(pv_nmax+iq), qtrc(:,:,:,iq), tracer_name(iq), 'ZXY' ) ! [IN]
3665  enddo
3666 
3667  endif
3668 
3669  if( atmos_sw_dyn ) call atmos_dyn_vars_restart_write
3670  if( atmos_sw_phy_mp ) call atmos_phy_mp_vars_restart_write
3671  if( atmos_sw_phy_ae ) call atmos_phy_ae_vars_restart_write
3672  if( atmos_sw_phy_ch ) call atmos_phy_ch_vars_restart_write
3673  if( atmos_sw_phy_rd ) call atmos_phy_rd_vars_restart_write
3674  if( atmos_sw_phy_sf .and. (.not. cpl_sw) ) call atmos_phy_sf_vars_restart_write
3675  if( atmos_sw_phy_tb ) call atmos_phy_tb_vars_restart_write
3676  if( atmos_sw_phy_bl ) call atmos_phy_bl_vars_restart_write
3677  if( atmos_sw_phy_cp ) call atmos_phy_cp_vars_restart_write
3678  if( atmos_sw_phy_lt ) call atmos_phy_lt_vars_restart_write
3679 
3680  return

References mod_atmos_dyn_vars::atmos_dyn_vars_restart_write(), mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_write(), mod_atmos_phy_bl_vars::atmos_phy_bl_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_lt_vars::atmos_phy_lt_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(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_lt, 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_check(), atmos_vars_fillhalo(), mod_cpl_admin::cpl_sw, 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, rhot, and scale_tracer::tracer_name.

Referenced by mod_admin_restart::admin_restart_write().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ allocate_3d()

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

Definition at line 3686 of file mod_atmos_vars.F90.

3686  use scale_const, only: &
3687  undef => const_undef
3688  real(RP), intent(inout), allocatable :: ary(:,:,:)
3689 
3690  if ( .not. allocated(ary) ) then
3691  allocate( ary(ka,ia,ja) )
3692  ary(:,:,:) = undef
3693  end if
3694 
3695  return

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

Referenced by atmos_vars_get_diagnostic_3d().

Here is the caller graph for this function:

◆ allocate_2d()

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

Definition at line 3699 of file mod_atmos_vars.F90.

3699  use scale_const, only: &
3700  undef => const_undef
3701  real(RP), intent(inout), allocatable :: ary(:,:)
3702 
3703  if ( .not. allocated(ary) ) then
3704  allocate( ary(ia,ja) )
3705  ary(:,:) = undef
3706  end if
3707 
3708  return

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

Referenced by atmos_vars_get_diagnostic_2d().

Here is the caller graph for this function:

◆ allocate_1d()

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

Definition at line 3712 of file mod_atmos_vars.F90.

3712  use scale_const, only: &
3713  undef => const_undef
3714  real(RP), intent(inout), allocatable :: ary(:)
3715 
3716  if ( .not. allocated(ary) ) then
3717  allocate( ary(ka) )
3718  ary(:) = undef
3719  end if
3720 
3721  return

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

Referenced by atmos_vars_get_diagnostic_1d().

Here is the caller graph for this function:

Variable Documentation

◆ atmos_restart_output

logical, public mod_atmos_vars::atmos_restart_output = .false.

Output restart file?

Definition at line 59 of file mod_atmos_vars.F90.

59  logical, public :: ATMOS_RESTART_OUTPUT = .false.

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

◆ atmos_restart_in_basename

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

Basename of the input file.

Definition at line 61 of file mod_atmos_vars.F90.

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

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

◆ atmos_restart_in_aggregate

logical, public mod_atmos_vars::atmos_restart_in_aggregate

Switch to use aggregate file.

Definition at line 62 of file mod_atmos_vars.F90.

62  logical, public :: ATMOS_RESTART_IN_AGGREGATE

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

◆ atmos_restart_in_postfix_timelabel

logical, public mod_atmos_vars::atmos_restart_in_postfix_timelabel = .false.

Add timelabel to the basename of input file?

Definition at line 63 of file mod_atmos_vars.F90.

63  logical, public :: ATMOS_RESTART_IN_POSTFIX_TIMELABEL = .false.

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

◆ atmos_restart_out_basename

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

Basename of the output file.

Definition at line 64 of file mod_atmos_vars.F90.

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

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

◆ atmos_restart_out_aggregate

logical, public mod_atmos_vars::atmos_restart_out_aggregate

Switch to use aggregate file.

Definition at line 65 of file mod_atmos_vars.F90.

65  logical, public :: ATMOS_RESTART_OUT_AGGREGATE

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

◆ atmos_restart_out_postfix_timelabel

logical, public mod_atmos_vars::atmos_restart_out_postfix_timelabel = .true.

Add timelabel to the basename of output file?

Definition at line 66 of file mod_atmos_vars.F90.

66  logical, public :: ATMOS_RESTART_OUT_POSTFIX_TIMELABEL = .true.

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

◆ atmos_restart_out_title

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

Title of the output file.

Definition at line 67 of file mod_atmos_vars.F90.

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

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

◆ atmos_restart_out_dtype

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

REAL4 or REAL8.

Definition at line 68 of file mod_atmos_vars.F90.

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

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

◆ atmos_restart_check

logical, public mod_atmos_vars::atmos_restart_check = .false.

Check value consistency?

Definition at line 70 of file mod_atmos_vars.F90.

70  logical, public :: ATMOS_RESTART_CHECK = .false.

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

◆ atmos_restart_check_basename

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

Definition at line 71 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

◆ atmos_restart_check_criterion

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

Definition at line 72 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_restart_check(), and atmos_vars_setup().

◆ dens

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

Definition at line 75 of file mod_atmos_vars.F90.

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

Referenced by mod_atmos_driver::atmos_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_adjustment(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), mod_atmos_driver::atmos_surface_set(), atmos_vars_check(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_3d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), mod_atmos_phy_lt_driver::history(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_mkinit::rect_setup(), mod_rm_driver::restart_read(), and mod_realinput::urban_input().

◆ momz

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

Definition at line 76 of file mod_atmos_vars.F90.

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

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

◆ momx

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

Definition at line 77 of file mod_atmos_vars.F90.

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

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

◆ momy

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

Definition at line 78 of file mod_atmos_vars.F90.

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

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

◆ rhot

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

Definition at line 79 of file mod_atmos_vars.F90.

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

Referenced by mod_atmos_bnd_driver::atmos_boundary_driver_send(), mod_atmos_bnd_driver::atmos_boundary_driver_set(), mod_atmos_dyn_driver::atmos_dyn_driver(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_adjustment(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_vars_check(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_1d(), atmos_vars_history(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_open(), atmos_vars_restart_read(), atmos_vars_restart_write(), atmos_vars_setup(), mod_mkinit::read_sounding(), mod_realinput::realinput_atmos(), mod_atmos_bnd_driver::set_boundary(), and mod_realinput::urban_input().

◆ qtrc

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

Definition at line 80 of file mod_atmos_vars.F90.

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

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

◆ dens_avw

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

Definition at line 82 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ momz_avw

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

Definition at line 83 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ momx_avw

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

Definition at line 84 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ momy_avw

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

Definition at line 85 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ rhot_avw

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

Definition at line 86 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ qtrc_avw

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

Definition at line 87 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_setup().

◆ dens_av

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

Definition at line 89 of file mod_atmos_vars.F90.

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

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

◆ momz_av

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

Definition at line 90 of file mod_atmos_vars.F90.

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

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

◆ momx_av

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

Definition at line 91 of file mod_atmos_vars.F90.

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

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

◆ momy_av

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

Definition at line 92 of file mod_atmos_vars.F90.

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

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

◆ rhot_av

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

Definition at line 93 of file mod_atmos_vars.F90.

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

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

◆ qtrc_av

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

Definition at line 94 of file mod_atmos_vars.F90.

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

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

◆ qv

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

Definition at line 96 of file mod_atmos_vars.F90.

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

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

◆ qc

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

Definition at line 97 of file mod_atmos_vars.F90.

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

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

◆ qr

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

Definition at line 98 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

◆ qi

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

Definition at line 99 of file mod_atmos_vars.F90.

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

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

◆ qs

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

Definition at line 100 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

◆ qg

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

Definition at line 101 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

◆ qh

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

Definition at line 102 of file mod_atmos_vars.F90.

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

Referenced by atmos_vars_get_diagnostic_3d(), and atmos_vars_setup().

◆ qe

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

Definition at line 104 of file mod_atmos_vars.F90.

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

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

◆ mass

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

Definition at line 104 of file mod_atmos_vars.F90.

◆ ratio

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

Definition at line 104 of file mod_atmos_vars.F90.

◆ of

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

Definition at line 104 of file mod_atmos_vars.F90.

◆ hydrometors

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

Definition at line 104 of file mod_atmos_vars.F90.

◆ kg

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

Definition at line 104 of file mod_atmos_vars.F90.

◆ dens_ref

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

Definition at line 107 of file mod_atmos_vars.F90.

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

◆ pott_ref

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

Definition at line 108 of file mod_atmos_vars.F90.

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

◆ temp_ref

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

Definition at line 109 of file mod_atmos_vars.F90.

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

◆ pres_ref

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

Definition at line 110 of file mod_atmos_vars.F90.

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

◆ qv_ref

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

Definition at line 111 of file mod_atmos_vars.F90.

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

◆ dens_tp

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

Definition at line 114 of file mod_atmos_vars.F90.

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

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

◆ momz_tp

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

Definition at line 115 of file mod_atmos_vars.F90.

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

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

◆ rhou_tp

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

Definition at line 116 of file mod_atmos_vars.F90.

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

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

◆ rhov_tp

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

Definition at line 117 of file mod_atmos_vars.F90.

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

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

◆ rhot_tp

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

Definition at line 118 of file mod_atmos_vars.F90.

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

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

◆ rhoh_p

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

Definition at line 119 of file mod_atmos_vars.F90.

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

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

◆ rhoq_tp

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

Definition at line 120 of file mod_atmos_vars.F90.

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

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

◆ momx_tp

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

Definition at line 123 of file mod_atmos_vars.F90.

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

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

◆ momy_tp

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

Definition at line 124 of file mod_atmos_vars.F90.

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

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

◆ w

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

Definition at line 128 of file mod_atmos_vars.F90.

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

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

◆ velocity

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

Definition at line 128 of file mod_atmos_vars.F90.

◆ m

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

Definition at line 128 of file mod_atmos_vars.F90.

◆ s

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

Definition at line 128 of file mod_atmos_vars.F90.

◆ u

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

Definition at line 129 of file mod_atmos_vars.F90.

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

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

◆ v

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

Definition at line 130 of file mod_atmos_vars.F90.

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

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

◆ pott

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

Definition at line 132 of file mod_atmos_vars.F90.

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

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

◆ potential

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

Definition at line 132 of file mod_atmos_vars.F90.

◆ temperature

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

Definition at line 132 of file mod_atmos_vars.F90.

◆ k

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

Definition at line 132 of file mod_atmos_vars.F90.

◆ temp

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

Definition at line 133 of file mod_atmos_vars.F90.

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

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

◆ pres

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

Definition at line 134 of file mod_atmos_vars.F90.

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

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

◆ pressure

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

Definition at line 134 of file mod_atmos_vars.F90.

◆ pa

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

Definition at line 134 of file mod_atmos_vars.F90.

◆ exner

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

Definition at line 135 of file mod_atmos_vars.F90.

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

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

◆ function

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

Definition at line 135 of file mod_atmos_vars.F90.

◆ phyd

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

Definition at line 136 of file mod_atmos_vars.F90.

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

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

◆ hydrostatic

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

Definition at line 136 of file mod_atmos_vars.F90.

◆ phydh

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

Definition at line 137 of file mod_atmos_vars.F90.

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

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

◆ layer

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

Definition at line 137 of file mod_atmos_vars.F90.

◆ interface

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

Definition at line 137 of file mod_atmos_vars.F90.

◆ qdry

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

Definition at line 139 of file mod_atmos_vars.F90.

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

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

◆ dry

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

Definition at line 139 of file mod_atmos_vars.F90.

◆ air

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

Definition at line 139 of file mod_atmos_vars.F90.

◆ rtot

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

Definition at line 140 of file mod_atmos_vars.F90.

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

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

◆ specific

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

Definition at line 140 of file mod_atmos_vars.F90.

◆ gass

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

Definition at line 140 of file mod_atmos_vars.F90.

◆ constant

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

Definition at line 140 of file mod_atmos_vars.F90.

◆ j

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

Definition at line 140 of file mod_atmos_vars.F90.

Referenced by atmos_vars_check(), atmos_vars_fillhalo(), atmos_vars_get_diagnostic_1d(), atmos_vars_get_diagnostic_2d(), atmos_vars_get_diagnostic_3d(), atmos_vars_monitor(), atmos_vars_restart_check(), atmos_vars_restart_read(), and file_add_variable_c().

◆ cvtot

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

Definition at line 141 of file mod_atmos_vars.F90.

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

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

◆ heat

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

Definition at line 141 of file mod_atmos_vars.F90.

◆ cptot

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

Definition at line 142 of file mod_atmos_vars.F90.

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

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

◆ prec

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

Definition at line 144 of file mod_atmos_vars.F90.

144  real(RP), public, allocatable, target :: PREC (:,:)

Referenced by mod_atmos_driver::atmos_surface_set(), atmos_vars_get_diagnostic_2d(), atmos_vars_monitor(), and atmos_vars_setup().

◆ total

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

Definition at line 144 of file mod_atmos_vars.F90.

◆ precipitation

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

Definition at line 144 of file mod_atmos_vars.F90.

◆ m2

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

Definition at line 144 of file mod_atmos_vars.F90.

◆ prec_engi

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

Definition at line 145 of file mod_atmos_vars.F90.

145  real(RP), public, allocatable :: PREC_ENGI(:,:)

Referenced by mod_atmos_driver::atmos_surface_set(), atmos_vars_monitor(), and atmos_vars_setup().

◆ internal

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

Definition at line 145 of file mod_atmos_vars.F90.

◆ energy

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

Definition at line 145 of file mod_atmos_vars.F90.

scale_const::const_grav
real(rp), public const_grav
standard acceleration of gravity [m/s2]
Definition: scale_const.F90:46
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_write
subroutine, public atmos_phy_cp_vars_restart_write
Write restart.
Definition: mod_atmos_phy_cp_vars.F90:543
mod_atmos_vars::momz_av
real(rp), dimension(:,:,:), pointer, public momz_av
Definition: mod_atmos_vars.F90:90
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_close
subroutine, public atmos_phy_tb_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_tb_vars.F90:311
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_create
subroutine, public atmos_phy_ch_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_ch_vars.F90:265
scale_statistics
module Statistics
Definition: scale_statistics.F90:11
mod_atmos_phy_mp_vars
module Atmosphere / Physics Cloud Microphysics
Definition: mod_atmos_phy_mp_vars.F90:12
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_create
subroutine, public atmos_phy_rd_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_rd_vars.F90:339
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_close
subroutine, public atmos_phy_ae_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_ae_vars.F90:372
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_close
subroutine, public atmos_phy_ch_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_ch_vars.F90:316
scale_atmos_diagnostic::atmos_diagnostic_get_teml
subroutine, public atmos_diagnostic_get_teml(KA, KS, KE, IA, IS, IE, JA, JS, JE, TEMP, LHV, LHS, QC, QI, CPtot, TEML)
ATMOS_DIAGNOSTIC_get_teml liqued water temperature.
Definition: scale_atmos_diagnostic.F90:320
scale_time::time_nowdaysec
real(dp), public time_nowdaysec
second of current time [sec]
Definition: scale_time.F90:70
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_open
subroutine, public atmos_phy_ch_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_ch_vars.F90:187
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_write
subroutine, public atmos_phy_lt_vars_restart_write
Write restart.
Definition: mod_atmos_phy_lt_vars.F90:368
scale_prc::prc_abort
subroutine, public prc_abort
Abort Process.
Definition: scale_prc.F90:342
mod_atmos_dyn_vars
module Atmosphere / Dynamics
Definition: mod_atmos_dyn_vars.F90:12
scale_atmos_hydrometeor::i_hr
integer, parameter, public i_hr
liquid water rain
Definition: scale_atmos_hydrometeor.F90:82
mod_atmos_admin::atmos_sw_phy_tb
logical, public atmos_sw_phy_tb
Definition: mod_atmos_admin.F90:57
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_open
subroutine, public atmos_phy_sf_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_sf_vars.F90:392
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_enddef
subroutine, public atmos_phy_ch_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_ch_vars.F90:302
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_get_diagnostic
subroutine, public atmos_phy_mp_vars_get_diagnostic(DENS, TEMP, QTRC, CLDFRAC, Re, Qe, Ne)
Definition: mod_atmos_phy_mp_vars.F90:573
mod_atmos_vars::momx_av
real(rp), dimension(:,:,:), pointer, public momx_av
Definition: mod_atmos_vars.F90:91
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolzxv
real(rp), public atmos_grid_cartesc_real_totvolzxv
total volume (zxv, local) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:90
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_cz
geopotential height [m] (zxy)
Definition: scale_atmos_grid_cartesC_real.F90:38
mod_atmos_vars::rhoq_tp
real(rp), dimension(:,:,:,:), allocatable, public rhoq_tp
Definition: mod_atmos_vars.F90:120
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_enddef
subroutine, public atmos_phy_lt_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_lt_vars.F90:318
mod_atmos_vars::pott
real(rp), dimension(:,:,:), allocatable, target, public pott
Definition: mod_atmos_vars.F90:132
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_close
subroutine, public atmos_phy_rd_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_rd_vars.F90:390
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_enddef
subroutine, public atmos_phy_tb_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_tb_vars.F90:297
scale_atmos_hydrometeor::i_hs
integer, parameter, public i_hs
snow
Definition: scale_atmos_hydrometeor.F90:84
scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdx
real(rp), dimension(:), allocatable, public atmos_grid_cartesc_rfdx
reciprocal of face-dx
Definition: scale_atmos_grid_cartesC.F90:67
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volwxy
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_volwxy
control volume (wxy) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:84
scale_file_cartesc::file_cartesc_enddef
subroutine, public file_cartesc_enddef(fid)
Exit netCDF file define mode.
Definition: scale_file_cartesC.F90:943
mod_atmos_admin::atmos_use_average
logical, public atmos_use_average
Definition: mod_atmos_admin.F90:49
scale_file_cartesc::file_cartesc_def_var
subroutine, public file_cartesc_def_var(fid, varname, desc, unit, dim_type, datatype, vid, standard_name, timeintv, nsteps, cell_measures)
Define a variable to file.
Definition: scale_file_cartesC.F90:3307
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_reset_diagnostics
subroutine, public atmos_phy_ae_vars_reset_diagnostics
Definition: mod_atmos_phy_ae_vars.F90:532
mod_atmos_vars::qtrc_av
real(rp), dimension(:,:,:,:), pointer, public qtrc_av
Definition: mod_atmos_vars.F90:94
scale_atmos_diagnostic::atmos_diagnostic_get_potv
subroutine, public atmos_diagnostic_get_potv(KA, KS, KE, IA, IS, IE, JA, JS, JE, POTT, Rtot, POTV)
ATMOS_DIAGNOSTIC_get_potv virtual potential temperature.
Definition: scale_atmos_diagnostic.F90:277
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_def_var
subroutine, public atmos_phy_bl_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_bl_vars.F90:336
mod_atmos_admin
module ATMOS admin
Definition: mod_atmos_admin.F90:11
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_def_var
subroutine, public atmos_phy_sf_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_sf_vars.F90:549
mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_lw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_tomflx_lw_dn
Definition: mod_atmos_phy_rd_vars.F90:67
scale_file_history_cartesc::file_history_cartesc_set_pres
subroutine, public file_history_cartesc_set_pres(PRES, PRESH, SFC_PRES)
set hydrostatic pressure for pressure coordinate
Definition: scale_file_history_cartesC.F90:260
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_def_var
subroutine, public atmos_phy_ch_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_ch_vars.F90:336
mod_atmos_phy_lt_vars
module Atmosphere / Physics Chemistry
Definition: mod_atmos_phy_lt_vars.F90:12
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_create
subroutine, public atmos_phy_ae_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_ae_vars.F90:321
scale_const::const_rvap
real(rp), parameter, public const_rvap
specific gas constant (water vapor) [J/kg/K]
Definition: scale_const.F90:63
scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdx
real(rp), dimension(:), allocatable, public atmos_grid_cartesc_rcdx
reciprocal of center-dx
Definition: scale_atmos_grid_cartesC.F90:65
mod_atmos_admin::atmos_sw_dyn
logical, public atmos_sw_dyn
Definition: mod_atmos_admin.F90:51
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_open
subroutine, public atmos_phy_cp_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_cp_vars.F90:323
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_create
subroutine, public atmos_phy_cp_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_cp_vars.F90:432
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_read
subroutine, public atmos_phy_tb_vars_restart_read
Read restart.
Definition: mod_atmos_phy_tb_vars.F90:218
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolwxy
real(rp), public atmos_grid_cartesc_real_totvolwxy
total volume (wxy, local) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:88
scale_atmos_grid_cartesc_metric
module Atmosphere Grid CartesianC metirc
Definition: scale_atmos_grid_cartesC_metric.F90:12
mod_atmos_admin::atmos_sw_phy_mp
logical, public atmos_sw_phy_mp
Definition: mod_atmos_admin.F90:52
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_open
subroutine, public atmos_phy_rd_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_rd_vars.F90:265
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_setup
subroutine, public atmos_phy_lt_vars_setup
Setup.
Definition: mod_atmos_phy_lt_vars.F90:91
mod_atmos_vars::rhot_av
real(rp), dimension(:,:,:), pointer, public rhot_av
Definition: mod_atmos_vars.F90:93
scale_atmos_hydrometeor
module atmosphere / hydrometeor
Definition: scale_atmos_hydrometeor.F90:12
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_setup
subroutine, public atmos_phy_ch_vars_setup
Setup.
Definition: mod_atmos_phy_ch_vars.F90:90
scale_atmos_diagnostic
module atmosphere / diagnostic
Definition: scale_atmos_diagnostic.F90:12
mod_atmos_phy_rd_vars
module Atmosphere / Physics Radiation
Definition: mod_atmos_phy_rd_vars.F90:12
scale_prc::prc_myrank
integer, public prc_myrank
process num in local communicator
Definition: scale_prc.F90:90
scale_atmos_diagnostic::atmos_diagnostic_get_phyd
subroutine, public atmos_diagnostic_get_phyd(KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, PRES, CZ, FZ, PHYD, PHYDH)
ATMOS_DIAGNOSTIC_get_phyd hydrostatic pressure.
Definition: scale_atmos_diagnostic.F90:157
mod_atmos_dyn_vars::atmos_dyn_vars_restart_read
subroutine, public atmos_dyn_vars_restart_read
Read restart.
Definition: mod_atmos_dyn_vars.F90:240
mod_atmos_phy_sf_vars
module ATMOSPHERIC Surface Variables
Definition: mod_atmos_phy_sf_vars.F90:12
mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain
real(rp), dimension(:,:), allocatable, public atmos_phy_mp_sflx_rain
Definition: mod_atmos_phy_mp_vars.F90:73
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_enddef
subroutine, public atmos_phy_cp_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_cp_vars.F90:469
scale_atmos_hydrometeor::i_hh
integer, parameter, public i_hh
hail
Definition: scale_atmos_hydrometeor.F90:86
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_create
subroutine, public atmos_phy_lt_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_lt_vars.F90:281
scale_atmos_hydrometeor::atmos_hydrometeor_dry
logical, public atmos_hydrometeor_dry
Definition: scale_atmos_hydrometeor.F90:97
mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_sw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_tomflx_sw_up
Definition: mod_atmos_phy_rd_vars.F90:68
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_def_var
subroutine, public atmos_phy_mp_vars_restart_def_var
Define variables in restart file.
Definition: mod_atmos_phy_mp_vars.F90:428
mod_atmos_vars::rhot
real(rp), dimension(:,:,:), allocatable, target, public rhot
Definition: mod_atmos_vars.F90:79
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_create
subroutine, public atmos_phy_tb_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_tb_vars.F90:261
scale_atmos_grid_cartesc_real
module Atmosphere GRID CartesC Real(real space)
Definition: scale_atmos_grid_cartesC_real.F90:11
mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp
real(rp), dimension(:,:), allocatable, public atmos_phy_sf_sfc_temp
Definition: mod_atmos_phy_sf_vars.F90:64
mod_atmos_vars::qtrc
real(rp), dimension(:,:,:,:), allocatable, target, public qtrc
Definition: mod_atmos_vars.F90:80
scale_file_history
module file_history
Definition: scale_file_history.F90:15
scale_file
module file
Definition: scale_file.F90:15
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_write
subroutine, public atmos_phy_tb_vars_restart_write
Write restart.
Definition: mod_atmos_phy_tb_vars.F90:351
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volzuy
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_volzuy
control volume (zuy) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:85
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_create
subroutine, public atmos_phy_mp_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_mp_vars.F90:357
scale_prc
module PROCESS
Definition: scale_prc.F90:11
mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_snow
real(rp), dimension(:,:), allocatable, public atmos_phy_cp_sflx_snow
Definition: mod_atmos_phy_cp_vars.F90:71
mod_atmos_admin::atmos_sw_phy_ae
logical, public atmos_sw_phy_ae
Definition: mod_atmos_admin.F90:53
mod_atmos_dyn_vars::atmos_dyn_vars_restart_open
subroutine, public atmos_dyn_vars_restart_open
Open restart file for read.
Definition: mod_atmos_dyn_vars.F90:203
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_history
subroutine, public atmos_phy_ae_vars_history(QTRC, RH)
Definition: mod_atmos_phy_ae_vars.F90:428
scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_mapf
real(rp), dimension(:,:,:,:), allocatable, public atmos_grid_cartesc_metric_mapf
map factor
Definition: scale_atmos_grid_cartesC_metric.F90:34
scale_atmos_diagnostic_cartesc::atmos_diagnostic_cartesc_get_vel
subroutine, public atmos_diagnostic_cartesc_get_vel(KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, MOMZ, MOMX, MOMY, W, U, V)
ATMOS_DIAGNOSTIC_CARTESC_get_vel W, U, V.
Definition: scale_atmos_diagnostic_cartesC.F90:58
scale_atmos_hydrometeor::i_hi
integer, parameter, public i_hi
ice water cloud
Definition: scale_atmos_hydrometeor.F90:83
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_setup
subroutine, public atmos_phy_tb_vars_setup
Setup.
Definition: mod_atmos_phy_tb_vars.F90:83
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_lw_dn
Definition: mod_atmos_phy_rd_vars.F90:61
scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdy
real(rp), dimension(:), allocatable, public atmos_grid_cartesc_rcdy
reciprocal of center-dy
Definition: scale_atmos_grid_cartesC.F90:66
mod_atmos_phy_bl_vars
module atmosphere / physics / PBL
Definition: mod_atmos_phy_bl_vars.F90:12
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_sw_dn
Definition: mod_atmos_phy_rd_vars.F90:63
mod_atmos_vars::dens
real(rp), dimension(:,:,:), allocatable, target, public dens
Definition: mod_atmos_vars.F90:75
mod_atmos_dyn_vars::atmos_dyn_vars_restart_create
subroutine, public atmos_dyn_vars_restart_create
Create restart file.
Definition: mod_atmos_dyn_vars.F90:292
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_reset_diagnostics
subroutine, public atmos_phy_mp_vars_reset_diagnostics
Definition: mod_atmos_phy_mp_vars.F90:768
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_write
subroutine, public atmos_phy_sf_vars_restart_write
Write variables to restart file.
Definition: mod_atmos_phy_sf_vars.F90:574
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_close
subroutine, public atmos_phy_sf_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_sf_vars.F90:529
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_volzxv
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_volzxv
control volume (zxv) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:86
scale_const
module CONSTANT
Definition: scale_const.F90:11
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_vol
control volume (zxy) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:83
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_enddef
subroutine, public atmos_phy_sf_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_sf_vars.F90:515
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_sw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_sw_up
Definition: mod_atmos_phy_rd_vars.F90:62
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_read
subroutine, public atmos_phy_lt_vars_restart_read
Read restart.
Definition: mod_atmos_phy_lt_vars.F90:228
mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_lw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_sflx_lw_up
Definition: mod_atmos_phy_rd_vars.F90:60
mod_atmos_vars::momz
real(rp), dimension(:,:,:), allocatable, target, public momz
Definition: mod_atmos_vars.F90:76
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_setup
subroutine, public atmos_phy_sf_vars_setup
Setup.
Definition: mod_atmos_phy_sf_vars.F90:189
scale_atmos_hydrometeor::lhv
real(rp), public lhv
latent heat of vaporization for use [J/kg]
Definition: scale_atmos_hydrometeor.F90:126
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_def_var
subroutine, public atmos_phy_rd_vars_restart_def_var
Define variables in restart file.
Definition: mod_atmos_phy_rd_vars.F90:410
mod_atmos_admin::atmos_sw_phy_lt
logical, public atmos_sw_phy_lt
Definition: mod_atmos_admin.F90:60
scale_atmos_adiabat
module atmosphere / adiabat
Definition: scale_atmos_adiabat.F90:11
mod_atmos_phy_tb_vars
module Atmosphere / Physics Turbulence
Definition: mod_atmos_phy_tb_vars.F90:12
scale_const::const_cvdry
real(rp), public const_cvdry
specific heat (dry air,constant volume) [J/kg/K]
Definition: scale_const.F90:57
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_close
subroutine, public atmos_phy_cp_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_cp_vars.F90:483
mod_atmos_vars::momy_av
real(rp), dimension(:,:,:), pointer, public momy_av
Definition: mod_atmos_vars.F90:92
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_write
subroutine, public atmos_phy_rd_vars_restart_write
Write variables to restart file.
Definition: mod_atmos_phy_rd_vars.F90:434
scale_file_cartesc::file_cartesc_close
subroutine, public file_cartesc_close(fid)
Close a netCDF file.
Definition: scale_file_cartesC.F90:1023
scale_prc_cartesc
module process / cartesC
Definition: scale_prc_cartesC.F90:11
scale_const::const_cpdry
real(rp), public const_cpdry
specific heat (dry air,constant pressure) [J/kg/K]
Definition: scale_const.F90:56
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_read
subroutine, public atmos_phy_ae_vars_restart_read
Read restart.
Definition: mod_atmos_phy_ae_vars.F90:278
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_def_var
subroutine, public atmos_phy_lt_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_lt_vars.F90:352
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_create
subroutine, public atmos_phy_bl_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_bl_vars.F90:266
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_setup
subroutine, public atmos_phy_cp_vars_setup
Setup.
Definition: mod_atmos_phy_cp_vars.F90:127
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_area
real(rp), dimension(:,:), allocatable, public atmos_grid_cartesc_real_area
horizontal area ( xy, normal z) [m2]
Definition: scale_atmos_grid_cartesC_real.F90:65
mod_atmos_phy_ch_vars
module Atmosphere / Physics Chemistry
Definition: mod_atmos_phy_ch_vars.F90:12
mod_atmos_vars::momz_tp
real(rp), dimension(:,:,:), allocatable, public momz_tp
Definition: mod_atmos_vars.F90:115
mod_atmos_vars::momx
real(rp), dimension(:,:,:), allocatable, target, public momx
Definition: mod_atmos_vars.F90:77
mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_sw_dn
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_tomflx_sw_dn
Definition: mod_atmos_phy_rd_vars.F90:69
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_create
subroutine, public atmos_phy_sf_vars_restart_create
Create restart file.
Definition: mod_atmos_phy_sf_vars.F90:479
mod_atmos_admin::atmos_sw_phy_cp
logical, public atmos_sw_phy_cp
Definition: mod_atmos_admin.F90:59
mod_atmos_admin::atmos_sw_phy_ch
logical, public atmos_sw_phy_ch
Definition: mod_atmos_admin.F90:54
scale_monitor::monitor_reg
subroutine, public monitor_reg(name, desc, unit, itemid, ndims, dim_type, is_tendency)
Search existing item, or matching check between requested and registered item.
Definition: scale_monitor.F90:241
scale_time::time_dtsec_atmos_dyn
real(dp), public time_dtsec_atmos_dyn
time interval of dynamics [sec]
Definition: scale_time.F90:35
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvolzuy
real(rp), public atmos_grid_cartesc_real_totvolzuy
total volume (zuy, local) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:89
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_write
subroutine, public atmos_phy_mp_vars_restart_write
Write restart.
Definition: mod_atmos_phy_mp_vars.F90:450
mod_atmos_vars::momy
real(rp), dimension(:,:,:), allocatable, target, public momy
Definition: mod_atmos_vars.F90:78
scale_time
module TIME
Definition: scale_time.F90:11
mod_atmos_dyn_vars::atmos_dyn_vars_restart_write
subroutine, public atmos_dyn_vars_restart_write
Write variables to restart file.
Definition: mod_atmos_dyn_vars.F90:387
mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc
real(rp), dimension(:,:,:), allocatable, target, public atmos_phy_sf_sflx_qtrc
Definition: mod_atmos_phy_sf_vars.F90:84
scale_atmos_hydrometeor::i_hc
integer, parameter, public i_hc
liquid water cloud
Definition: scale_atmos_hydrometeor.F90:81
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_open
subroutine, public atmos_phy_bl_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_bl_vars.F90:196
scale_atmos_hydrometeor::i_qv
integer, public i_qv
Definition: scale_atmos_hydrometeor.F90:77
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_def_var
subroutine, public atmos_phy_tb_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_tb_vars.F90:331
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_setup
subroutine, public atmos_phy_rd_vars_setup
Setup.
Definition: mod_atmos_phy_rd_vars.F90:121
scale_atmos_diagnostic_cartesc
module atmosphere / diagnostic / CartesianC
Definition: scale_atmos_diagnostic_cartesC.F90:12
scale_atmos_diagnostic::atmos_diagnostic_get_therm_rhot
subroutine, public atmos_diagnostic_get_therm_rhot(KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, RHOT, Rtot, CVtot, CPtot, POTT, TEMP, PRES, EXNER)
ATMOS_DIAGNOSTIC_get_therm_rhot potential temperature, temperature, pressure.
Definition: scale_atmos_diagnostic.F90:59
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_write
subroutine, public atmos_phy_ae_vars_restart_write
Write restart.
Definition: mod_atmos_phy_ae_vars.F90:408
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_def_var
subroutine, public atmos_phy_ae_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_ae_vars.F90:392
mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi
real(rp), dimension(:,:), allocatable, public atmos_phy_sf_sflx_engi
Definition: mod_atmos_phy_sf_vars.F90:85
mod_atmos_admin::atmos_sw_phy_rd
logical, public atmos_sw_phy_rd
Definition: mod_atmos_admin.F90:55
mod_atmos_vars::dens_av
real(rp), dimension(:,:,:), pointer, public dens_av
Definition: mod_atmos_vars.F90:89
mod_atmos_dyn_vars::atmos_dyn_vars_restart_enddef
subroutine, public atmos_dyn_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_dyn_vars.F90:331
mod_atmos_phy_rd_vars::atmos_phy_rd_tomflx_lw_up
real(rp), dimension(:,:), allocatable, public atmos_phy_rd_tomflx_lw_up
Definition: mod_atmos_phy_rd_vars.F90:66
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_read
subroutine, public atmos_phy_rd_vars_restart_read
Read restart.
Definition: mod_atmos_phy_rd_vars.F90:300
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_f2h
real(rp), dimension(:,:,:,:), allocatable, public atmos_grid_cartesc_real_f2h
coefficient for interpolation from full to half levels
Definition: scale_atmos_grid_cartesC_real.F90:46
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_write
subroutine, public atmos_phy_ch_vars_restart_write
Write restart.
Definition: mod_atmos_phy_ch_vars.F90:352
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_setup
subroutine, public atmos_phy_mp_vars_setup
Setup.
Definition: mod_atmos_phy_mp_vars.F90:126
mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh
real(rp), dimension(:,:), allocatable, public atmos_phy_sf_sflx_lh
Definition: mod_atmos_phy_sf_vars.F90:80
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_read
subroutine, public atmos_phy_mp_vars_restart_read
Read restart.
Definition: mod_atmos_phy_mp_vars.F90:323
mod_atmos_dyn_vars::atmos_dyn_vars_restart_close
subroutine, public atmos_dyn_vars_restart_close
Close restart file.
Definition: mod_atmos_dyn_vars.F90:345
scale_file_history_cartesc
module file / history_cartesC
Definition: scale_file_history_cartesC.F90:12
scale_file_cartesc::file_cartesc_create
subroutine, public file_cartesc_create(basename, title, datatype, fid, date, subsec, haszcoord, append, aggregate, single)
Create/open a netCDF file.
Definition: scale_file_cartesC.F90:780
mod_atmos_dyn_vars::atmos_dyn_vars_restart_def_var
subroutine, public atmos_dyn_vars_restart_def_var
Define variables in restart file.
Definition: mod_atmos_dyn_vars.F90:365
mod_atmos_admin::atmos_dyn_type
character(len=h_short), public atmos_dyn_type
Definition: mod_atmos_admin.F90:35
scale_file_cartesc::file_cartesc_open
subroutine, public file_cartesc_open(basename, fid, aggregate)
open a netCDF file for read
Definition: scale_file_cartesC.F90:746
scale_statistics::statistics_checktotal
logical, public statistics_checktotal
calc&report variable totals to logfile?
Definition: scale_statistics.F90:64
scale_file_cartesc::file_cartesc_flush
subroutine, public file_cartesc_flush(fid)
Flush all pending requests to a netCDF file (PnetCDF only)
Definition: scale_file_cartesC.F90:997
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_setup
subroutine, public atmos_phy_ae_vars_setup
Setup.
Definition: mod_atmos_phy_ae_vars.F90:109
scale_time::time_gettimelabel
subroutine, public time_gettimelabel(timelabel)
generate time label
Definition: scale_time.F90:91
scale_comm_cartesc
module COMMUNICATION
Definition: scale_comm_cartesC.F90:11
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_read
subroutine, public atmos_phy_bl_vars_restart_read
Read restart.
Definition: mod_atmos_phy_bl_vars.F90:231
scale_atmos_bottom
module atmosphere / bottom boundary extrapolation
Definition: scale_atmos_bottom.F90:12
mod_atmos_vars::rhot_tp
real(rp), dimension(:,:,:), allocatable, public rhot_tp
Definition: mod_atmos_vars.F90:118
mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_rain
real(rp), dimension(:,:), allocatable, public atmos_phy_cp_sflx_rain
Definition: mod_atmos_phy_cp_vars.F90:70
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_enddef
subroutine, public atmos_phy_ae_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_ae_vars.F90:358
mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow
real(rp), dimension(:,:), allocatable, public atmos_phy_mp_sflx_snow
Definition: mod_atmos_phy_mp_vars.F90:74
mod_atmos_vars::momy_tp
real(rp), dimension(:,:,:), allocatable, public momy_tp
Definition: mod_atmos_vars.F90:124
mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh
real(rp), dimension(:,:), allocatable, public atmos_phy_sf_sflx_sh
Definition: mod_atmos_phy_sf_vars.F90:79
mod_atmos_admin::atmos_sw_phy_bl
logical, public atmos_sw_phy_bl
Definition: mod_atmos_admin.F90:58
mod_atmos_phy_tb_vars::atmos_phy_tb_vars_restart_open
subroutine, public atmos_phy_tb_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_tb_vars.F90:183
scale_atmos_diagnostic::atmos_diagnostic_get_n2
subroutine, public atmos_diagnostic_get_n2(KA, KS, KE, IA, IS, IE, JA, JS, JE, POTT, Rtot, CZ, FZ, F2H, N2)
ATMOS_DIAGNOSTIC_get_n2 N^2.
Definition: scale_atmos_diagnostic.F90:219
scale_atmos_hydrometeor::lhf
real(rp), public lhf
latent heat of fusion for use [J/kg]
Definition: scale_atmos_hydrometeor.F90:128
scale_file::file_get_aggregate
logical function, public file_get_aggregate(fid)
Definition: scale_file.F90:4844
mod_atmos_phy_ae_vars
module ATMOSPHERE / Physics Aerosol Microphysics
Definition: mod_atmos_phy_ae_vars.F90:12
mod_atmos_vars::momx_tp
real(rp), dimension(:,:,:), allocatable, public momx_tp
Definition: mod_atmos_vars.F90:123
scale_atmos_thermodyn
module atmosphere / thermodyn
Definition: scale_atmos_thermodyn.F90:11
scale_file_history::file_history_reg
subroutine, public file_history_reg(name, desc, unit, itemid, standard_name, ndims, dim_type, cell_measures, fill_halo)
Register/Append variable to history file.
Definition: scale_file_history.F90:650
mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_open
subroutine, public atmos_phy_ae_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_ae_vars.F90:243
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_write
subroutine, public atmos_phy_bl_vars_restart_write
Write restart.
Definition: mod_atmos_phy_bl_vars.F90:356
mod_atmos_phy_rd_vars::atmos_phy_rd_vars_restart_enddef
subroutine, public atmos_phy_rd_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_rd_vars.F90:376
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_close
subroutine, public atmos_phy_bl_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_bl_vars.F90:316
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_open
subroutine, public atmos_phy_mp_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_mp_vars.F90:288
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_real_fz
geopotential height [m] (wxy)
Definition: scale_atmos_grid_cartesC_real.F90:42
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_setup
subroutine, public atmos_phy_bl_vars_setup
Setup.
Definition: mod_atmos_phy_bl_vars.F90:97
scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol
real(rp), public atmos_grid_cartesc_real_totvol
total volume (zxy, local) [m3]
Definition: scale_atmos_grid_cartesC_real.F90:87
scale_atmos_grid_cartesc
module atmosphere / grid / cartesC
Definition: scale_atmos_grid_cartesC.F90:12
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_open
subroutine, public atmos_phy_lt_vars_restart_open
Open restart file for read.
Definition: mod_atmos_phy_lt_vars.F90:193
scale_const::const_undef
real(rp), public const_undef
Definition: scale_const.F90:41
scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdy
real(rp), dimension(:), allocatable, public atmos_grid_cartesc_rfdy
reciprocal of face-dy
Definition: scale_atmos_grid_cartesC.F90:68
mod_atmos_phy_sf_vars::atmos_phy_sf_vars_restart_read
subroutine, public atmos_phy_sf_vars_restart_read
Read restart.
Definition: mod_atmos_phy_sf_vars.F90:428
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_close
subroutine, public atmos_phy_mp_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_mp_vars.F90:408
scale_atmos_saturation
module atmosphere / saturation
Definition: scale_atmos_saturation.F90:12
mod_atmos_phy_lt_vars::atmos_phy_lt_vars_restart_close
subroutine, public atmos_phy_lt_vars_restart_close
Close restart file.
Definition: mod_atmos_phy_lt_vars.F90:332
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_history
subroutine, public atmos_phy_mp_vars_history(DENS, TEMP, QTRC)
Definition: mod_atmos_phy_mp_vars.F90:472
mod_atmos_phy_cp_vars
module Atmosphere / Physics Cumulus
Definition: mod_atmos_phy_cp_vars.F90:12
scale_prc_cartesc::prc_twod
logical, public prc_twod
2D experiment
Definition: scale_prc_cartesC.F90:55
scale_file_cartesc
module file / cartesianC
Definition: scale_file_cartesC.F90:11
mod_cpl_admin::cpl_sw
logical, public cpl_sw
Definition: mod_cpl_admin.F90:33
scale_atmos_hydrometeor::n_hyd
integer, parameter, public n_hyd
Definition: scale_atmos_hydrometeor.F90:79
scale_atmos_bottom::atmos_bottom_estimate
subroutine, public atmos_bottom_estimate(KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, PRES, QV, SFC_TEMP, FZ, SFC_DENS, SFC_PRES)
Calc bottom boundary of atmosphere (just above surface)
Definition: scale_atmos_bottom.F90:51
mod_atmos_phy_ch_vars::atmos_phy_ch_vars_restart_read
subroutine, public atmos_phy_ch_vars_restart_read
Read restart.
Definition: mod_atmos_phy_ch_vars.F90:222
mod_atmos_phy_bl_vars::atmos_phy_bl_vars_restart_enddef
subroutine, public atmos_phy_bl_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_bl_vars.F90:302
mod_cpl_admin
module Coupler admin
Definition: mod_cpl_admin.F90:11
mod_atmos_phy_mp_vars::atmos_phy_mp_vars_restart_enddef
subroutine, public atmos_phy_mp_vars_restart_enddef
Exit netCDF define mode.
Definition: mod_atmos_phy_mp_vars.F90:394
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_def_var
subroutine, public atmos_phy_cp_vars_restart_def_var
Write restart.
Definition: mod_atmos_phy_cp_vars.F90:503
mod_atmos_dyn_vars::atmos_dyn_vars_setup
subroutine, public atmos_dyn_vars_setup
Setup.
Definition: mod_atmos_dyn_vars.F90:80
scale_atmos_hydrometeor::i_hg
integer, parameter, public i_hg
graupel
Definition: scale_atmos_hydrometeor.F90:85
scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc
real(rp), dimension(:,:,:), allocatable, public atmos_grid_cartesc_metric_rotc
rotation coefficient
Definition: scale_atmos_grid_cartesC_metric.F90:35
mod_atmos_phy_cp_vars::atmos_phy_cp_vars_restart_read
subroutine, public atmos_phy_cp_vars_restart_read
Read restart.
Definition: mod_atmos_phy_cp_vars.F90:358
mod_atmos_admin::atmos_sw_phy_sf
logical, public atmos_sw_phy_sf
Definition: mod_atmos_admin.F90:56
scale_monitor
module MONITOR
Definition: scale_monitor.F90:12
Generated by   doxygen 1.8.17