Index | Rundeck: E6TomaF40 | Created: Fri May 17 03:30:39 EDT 2024 |
get_dep_vel | Module: | File: TRDRYDEP.f |
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:
subroutine get_dep_vel(i, j, itype, obk, zhh, ustarr, tempk, dep_vel, stomatal_dep_vel, trnmm )
- 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