Index | Rundeck: E6F40 | Created: Fri May 17 03:30:27 EDT 2024 |
ATM_DRV.f
Summary: Author :
Version:
Modules:
Global Subroutines:
- alloc_drv_atm
- atm_exports_phasesrf
- atm_phase1
- atm_phase1_exports
- atm_phase2
- checkt
- CHECKT Checks arrays for NaN/INF and reasonablness
- daily_atm
- def_rsf_atmvars
- def_rsf_atmvars defines atm prognostic array structure in rsf
- finalize_atm
- get_atm_layer1
- input_atm
- new_io_atmvars
- read_aic
- read_AIC for a cold start, read the atmospheric IC file.
Three input options are currently recognized
(1) The input file has already been remapped to the model layering
and contains the variables traditionally expected by the model
(winds, temperature, specific humidity, surface pressure).
(2) Winds, temperature, geopotential height, and RH are available
on constant-pressure levels. Surface pressure is obtained via
the hydrostatic assumption. Remapping from constant-pressure
levels to the model layering is peformed, and RH is converted
to specific humidity. The initial coding for this option was
imported from init_cond/AIC.D771201.f and generalized considerably.
If tropopause pressure and temperature are available, they
are used in the vertical remapping. The coding for the
assumed vertical structure of RH above the upper troposphere
will be made more configurable as AIC files with better-quality
RH data (than the historical AIC) become available.
(3) Surface pressure, winds, temperature, humidity on "standard-hybrid"
layers, where the pressure of the kth input layer at horizontal
location i,j is equal to hyam(k) + hybm(k)*srfp(i,j) and srfp is
the surface pressure. Remapping to model layering is performed.
Condensate species are also (optionally) read on other branches,
but not here (yet), for reasons including the inability of liquid
and ice to coexist under the cloud microphysics on this branch.
Logic will be added to handle other possible input layerings and
combinations of available fields.
- read_pzrh_input
- relayer_hybridlayer_input
- relayer_hybridlayer_input read and relayer atm state on standard hybrid layers
to whatever the model layering is. Relayering currently performed by vint_logp
for historical continuity with hindcasting runs that performed the relayering
with external codes; this means that the relayering is not conservative.
Other continuity choices:
- reading u,v at B-grid locations (and aux B-grid surface pressure for relayering)
- vertically interpolate log(q) rather than q.
- vint_logp
- vint_logp linear interpolation in log(p)-space
- vntrp1
Depends on the following files:
- ATMDYN.f
- ATMDYN_COM.F90
- ATM_COM.f
- AtmLayering.F90
- AtmRes.F90
- DIAG.f
- DIAG_COM.f
- DIAG_ZONAL.f
- FLUXES.f
- GEOM_B.f
- GHY_DRV.f
- LAKES_COM.f
- MODEL_COM.f
- MOMEN2ND.f
- RAD_COM.f
- SEAICE.f
- ATMDYN.f
Used by the following files: