Index

Rundeck: P1SoM40

Created: Wed May 8 03:30:18 EDT 2024

VERTICALRES

File: AtmL40p.F90

Summary: Vertical Resolution file, 40 layers, top at .1 mb but recoded to be able to import notEarth info
Author : Original Development Team
Version:

---

Subroutines:

---

Global Variables:

grav : used from constant

mb2kg : used from constant

planetparams : used from planetparams_mod

lm : integer*4, parameter

= number of dynamical layers
Initial Value = 40
Used by: | :read_aic |

ls1 : integer, parameter

= lowest layer of strtosphere
Initial Value = wt*ls1_lower + (1-wt)

pmtop :

global mean surface, model top pressure (mb)

psf : real*8, parameter

global mean surface, model top pressure (mb)
Initial Value = PlanetParams%psf

ptop :

pressure at interface level sigma/const press coord syst (mb)
Used by: | :read_aic |

pstrat :

pressure due to troposhere,stratosphere

psfmpt :

pressure due to troposhere,stratosphere

plbot :

pressure levels at bottom of layers (mb)

pratio : real*8, parameter, private

Initial Value = PSF/984d0

iii : integer, private

delp : real*8, dimension(lm), parameter, private

nominal pressure thicknesses of layers (mb)
Initial Value = plbot(1:lm)

ptop0 : real*8, parameter, private

Initial Value = PlanetParams%ptop

lev : integer, dimension(lm), parameter, private

Initial Value = (/ (iii,iii=1,lm) /)

pdiff : real*8, parameter, dimension(lm+1), private

Initial Value = 1d6*(plbot-ptop0)

ls1_lower : integer, parameter, private

Initial Value = &

tropomask :

unity from layers 1-ls1-1, zero above

stratmask :

zero from layers 1-ls1-1, unity above

mdrya : real*8, parameter

= dry atmospheric mass (kg/m^2) = 100*PSF/GRAV
Initial Value = psf*mb2kg

mtop : real*8, parameter

= mass above dynamical top (kg/m^2) = 100*PMTOP/GRAV
Initial Value = pmtop*mb2kg

mfixs : real*8, parameter

= summation of MFIX (kg/m^2) = 100*(PTOP-PMTOP)/GRAV
Initial Value = pstrat*mb2kg

mvar :

= spatially and temporally varying column mass (kg/m^2)

msurf :

= MTOP + MFIXs + MVAR (kg/m^2)

am :

(L) = MFIX(L) + MVAR*MFRAC(L) (kg/m^2)

mfix :

(L) = fixed mass in each layer (kg/m^2) = 100*[PLBOT(L)-PLBOT(L+1)]/GRAV

mfrac :

(L) = fraction of variable mass in each layer = DSIG(L)