IndexRundeck: E6TomaF40Created: Fri May 17 03:30:39 EDT 2024

get_dep_vel

Module:
File: TRDRYDEP.f
Summary: get_dep_vel computes the Bulk surface reistance to tracer dry deposition using a resistance-in-series model from a portion of the Harvard CTM dry deposition routine. The deposition velocity is the reciprocal of the Bulk surface reistance...
Author : D.J. Jacob and Y.H. Wang, modularized by G.M. Gardner, adapted for GISS GCM by D. Koch modelEified by G. Faluvegi Harvard version 3.1: 12/17/97)
Version:

Declaration:
Variables:
imaxj : used from geom

tf : used from constant

pi : used from constant

grav : used from constant

gasc : used from constant

cosz1 : used from rad_com

cfrac : used from rad_com

srdn : used from rad_com

tr_wd_type : used from oldtracer_mod

npart : used from oldtracer_mod

trname : used from oldtracer_mod

dodrydep : used from oldtracer_mod

f0_glob : used from oldtracer_mod
Initial Value = >F0
hstar_glob : used from oldtracer_mod
Initial Value = >HSTAR
ntm : used from tracer_com

n_nox : used from tracer_com

pnox : used from trchem_shindell_com

npoly : used from tracers_drydep

ijreg : used from tracers_drydep

ijland : used from tracers_drydep

xylai : used from tracers_drydep

drycoeff : used from tracers_drydep

ijuse : used from tracers_drydep

ntype : used from tracers_drydep

idep : used from tracers_drydep

iri : used from tracers_drydep

irlu : used from tracers_drydep

irac : used from tracers_drydep

irgss : used from tracers_drydep

irgso : used from tracers_drydep

ircls : used from tracers_drydep

irclo : used from tracers_drydep

ivsmax : used from tracers_drydep

xmwh2o : real*8, parameter
molecular weight of water in KG/mole
Initial Value = 18.d-3
n : integer
the current tracer number
ldt : integer
dummy loop variable for land type
iw : integer
dummy loop variables
k : integer
dummy loop variables
ri : real*8, dimension(ntype)
internal resistance (minimum stomatal resistance for
rlu : real*8, dimension(ntype)
cuticular resistance for the bulk canopy
rac : real*8, dimension(ntype)
______ resistance
rgss : real*8, dimension(ntype)
______ resistance
rgso : real*8, dimension(ntype)
______ resistance
rcls : real*8, dimension(ntype)
______ resistance
rclo : real*8, dimension(ntype)
______ resistance
rsurface : real*8, dimension(ntm, ntype)
Bulk surface resistance for species K landtype LDT
ri : real*8, dimension(ntype)
internal resistance (minimum stomatal resistance for
rstomatal : real*8, dimension(ntype)
Stomatal portion of RSURFACE for ozone, landtype LDT
rad0 : real*8
downward solar radiation flux at surface (w/m2)
rt : real*8
correction term for bulk surface resistance for gases?
rix : real*8
______ resistance
czh : real*8
?
rclx : real*8
?
rgsx : real*8
?
gfaci : real*8
?
gfact : real*8
?
rixx : real*8
corrected RIX
rluxx : real*8
corrected RLU
rdc : real*8
aerodynamic resistance to elements in lower part of
dtmp4 : real*8
temp arrays for recipricol resistances
dtmp3 : real*8
temp arrays for recipricol resistances
dtmp2 : real*8
temp arrays for recipricol resistances
dtmp1 : real*8
temp arrays for recipricol resistances
vds : real*8
temp deposition velocity
dummy4 : real*8
dummy temp variables
dummy3 : real*8
dummy temp variables
dummy2 : real*8
dummy temp variables
dummy1 : real*8
dummy temp variables
tempk : real*8, intent(in)
Surface air temperature (Kelvin)
tempc : real*8
Surface air temperature (oC)
bytempc : real*8
1/TEMPC
j : integer, intent(in)
GCM grid box horizontal position
i : integer, intent(in)
GCM grid box horizontal position
itype : integer, intent(in)
GCM surface type 1=ocean; 2=ocean ice; 3=land ice; 4=land
tota : real*8, dimension(ntm)
total landuse area, excluding water and ice dep. types.
problem_point : logical
logical that is true if ITYPE=4, but there are no
obk : real*8, intent(in)
Monin-Obukhov length (m) (now is the PBL lmonin variable)
zhh : real*8, intent(in)
boundary layer depth (m) ("mixing depth") (PBL dbl variable)
ustarr : real*8, intent(in)
Friction velocity (m s-1) (PBL ustar variable)
tr_mm_temp : real*8
temporary variable to hold trnmm(k), etc.
suncos : real*8
Cosine of solar zenith angle
iolson : integer
integer index for olson surface types?
vd : real*8, dimension(ntm)
deposition velocity temp array (s m-1)
svd : real*8
stomatal deposition velocity temp array (s m-1)
hstar : real*8, dimension(ntm)

f0 : real*8, dimension(ntm)

trnmm : real*8, dimension(ntm)

biofit : real*8

diffg : real*8

dep_vel : real*8, intent(out), dimension(ntm)
the deposition velocity = 1/bulk sfc. res. (m/s)
stomatal_dep_vel : real*8, intent(out)

ii : integer

Simplex Website Curator: Igor Aleinov — NASA Official: Gavin A. Schmidt

Contact GISS NASA Privacy PolicyAccessibility