| Index | Rundeck: P1SoM40 | Created: Thu May 9 03:30:17 EDT 2024 |
SUBDAILY | File: DIAG.f |
Author : Gavin Schmidt
Version:
Subroutines:
- accsubdd
- accSubdd accumulates variables for subdaily diagnostics
- close_subdd
- close_subdd closes sub daily diag files
- get_subdd
- get_SUBDD saves variables at sub-daily frequency
every ABS(NSUBDD)
Note that TMIN,TMAX,{ ,c,t,ct}AOD, are only output once/day.
If there is a choice between outputting pressure levels or
model levels, use lower case for the model levels:
Current options: SLP, PS, SAT, PREC, QS, LCLD, MCLD, HCLD, PTRO
QLAT, QSEN, SWD, SWU, LWD, LWU, LWT, SWT, STX, STY,
ICEF, SNOWD, TCLD, SST, SIT, US, VS, TMIN, TMAX
MCP, SNOWC, RS, GT1, GTD, GW0, GWD, GI0, GID,
GTALL, GWALL, GIALL (on soil levels)
{L,M,H}CLDI,CTPI,TAUI (ISCCP quantities)
LGTN, c2gLGTN (lightning flashs/cloud-to-ground)
TRP*, TRE* (water tracers only)
Z*, R*, T*, Q* (on any fixed pressure level)
z*, r*, t*, q* (on any model level)
U*, V*, W*, C* (on any model level only)
O*, X*, M*, N* (Ox,NOx,CO,NO2 on fixed pres lvl)
o*, x*, m*, n* (Ox,NOx,CO,NO2 on any model lvl)
oAVG (SFC Ox time-average ppbv)
nxAVG (SFC NOx time-average ppbv)
cAVG (SFC CO time-average ppbv)
oAVG1,nAVG1 (L=1 Ox and NO2 time-average ppbv)
PM2p5, PM10 (SFC time-average PM2.5 and PM10 ppmm)
PM2p51,PM101(L=1 time-average PM2.5 and PM10 ppmm)
cPM2p5,cPM10 (SFC time-average PM2.5, PM10 kg/m3)
NO2col NO2 column amount, instant., (kg/m2)
D* (HDO on any model level)
B* (BE7 on any model level)
SO4, RAPR
7BEW, 7BED, BE7ATM
LWC,IWC,TLH,SLH,DLH,LLH,OLR,PW,EVP,PBH ! mjo_subdd
SWH,LWH,TDRY,DDRY,SDRY,LDRY ! mjo_subdd
LWP,IWP,TTRO,PBLH ! additional etc_subdd
CTEM,CD3D,CI3D,CL3D,CDN3D,CRE3D,CLWP ! aerosol
TAUSS,TAUMC,CLDSS,CLDMC,MCCTP,CLDTOT
SO4_d1,SO4_d2,SO4_d3, ! het. chem
Clay, Silt1, Silt2, Silt3 ! dust
TrSMIXR surface mixing ratio for all tracers [kg/kg]
TrSCONC surface concentration for all tracers [kg/m^3]
DuEMIS soil dust aerosol emission flux [kg/m^2/s]
DuEMIS2 soil dust aerosol emission flux [kg/m^2/s]
from cubed wind speed (only diagnostic)
DuDEPTURB turbulent depo of soil dust aer [kg/m^2/s]
DuDEPGRAV grav settling of soil dust aerosols [kg/m^2/s]
DuDEPWET wet deposition of soil dust aerosols [kg/m^2/s]
DuLOAD soil dust aer load of atmospheric column [kg/m^2]
DuCONC three-dimensional soil dust concentrations [kg/m^3]
DuSMIXR surface mix ratio of soil dust aerosols [kg/kg]
DuSCONC surface conc of soil dust aerosols [kg/m^3]
DuAOD dust aer opt depth daily avg [1]
DuCSAOD clear sky dust aer opt depth daily avg [1]
AOD aer opt dep (1,nTracerRadiaActive in rad code) daily avg
tAOD aer opt dep (sum 1,nTracerRadiaActive) daily avg
ctAOD and cAOD are clr-sky versions of tAOD/AOD
ictAOD clr-sky sum AOD, 'instantaneous', 3D
itAAOD all-sky sum AOD ext-scat band6, instant., 3D
FRAC land fractions over 6 types
More options can be added as extra cases in this routine
- init_subdd
- init_subdd initialise sub daily diags and position files
- open_subdd
- open_subdd opens sub daily diag files
- reset_subdd
- reset_subdd resets sub daily diag files
- write_data
- write out subdd data array with optional pole fix
Global Variables:
- getdomainbounds : used from domain_decomp_atm
- grid : used from domain_decomp_atm
- am_i_root : used from domain_decomp_atm
- im : used from resolution
- jm : used from resolution
- lm : used from resolution
- modeleclock : used from model_com
- calendar : used from model_com
- itime : used from model_com
- itime0 : used from model_com
- nday : used from model_com
- iyear1 : used from model_com
- dtsrc : used from model_com
- xlabel : used from model_com
- lrunid : used from model_com
- int_days_per_year : used from timeconstants_mod
- openunit : used from filemanager
- closeunit : used from filemanager
- nameunit : used from filemanager
- gdeep : used from ghy_com
- gsavel : used from ghy_com
- ngm : used from ghy_com
- kgz_max : used from diag_com
- pmname : used from diag_com
- p_acc : used from diag_com
- pm_acc : used from diag_com
- kddmax :
integer, parameter
- maximum number of sub-daily diags output files
Initial Value= 55 - kdd :
integer
- total number of sub-daily diags
- not_open :
integer, parameter
- Initial Value
= -999 - kddunit :
integer
- total number of sub-daily files
- namedd :
character*10, dimension(kddmax)
- array of names of sub-daily diags
- iu_subdd :
integer, dimension(kddmax)
- array of unit numbers for sub-daily diags output
Initial Value= not_open - subddt :
character*320
- = subdd + subdd1,2,3 = all variables for sub-daily diags
Initial Value= " " - subdd :
character*64
- string contains variables to save for sub-daily diags
Initial Value= "SLP" - subdd1 :
character*64
- additional string of variables for sub-daily diags
Initial Value= " " - subdd2 :
character*64
- additional string of variables for sub-daily diags
Initial Value= " " - subdd3 :
character*64
- additional string of variables for sub-daily diags
Initial Value= " " - subdd4 :
character*64
- additional string of variables for sub-daily diags
Initial Value= " " - nsubdd :
integer
- : DT_save_SUBDD = Nsubdd*DTsrc sub-daily diag freq.
Initial Value= 0
Used by: | :atm_phase2 |
- lmaxsubdd :
integer
- : the max L when writing "ALL" levels
Initial Value= LM - lst :
character*2, dimension(lm)
- level strings
- lmaxsubdd_array :
real(kind=8), allocatable, dimension(:, :, :)
- array for three-dimensional fields for subdd diagnostics
- ngm_array :
real(kind=8), allocatable, dimension(:, :, :)
- three-dimensional array for subdd ground diagnostics
- adate_sv :
character*14, private
- kgz_max_suffixes :
character(len=8), allocatable, dimension(:)
- array of names for subdd-diagnostic on pressure levels
- kgz_max_array :
real(kind=8), allocatable, dimension(:, :, :)
- three-dimensional array for diagnostics on pressure levels