|
integer, parameter | nvegtype = 25 |
|
real(kind_phys), dimension(nvegtype), parameter | kpart = (/500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500., 500. /) |
|
real(kind_phys), parameter | max_dep_vel = 0.005 |
|
real(kind_phys), parameter | dep_ref_hgt = 2.0 |
|
real(kind_phys), parameter | pi = 3.1415926536 |
|
real(kind_phys), parameter | threepi =3.0*pi |
|
real(kind_phys), parameter | gravity = 9.81 |
|
real(kind_phys), parameter | grav =9.80622 |
|
real(kind_phys), parameter | boltzmann = 1.3807e-16 |
|
real(kind_phys), parameter | rgasuniv =8.314510 |
|
real, parameter | avo =6.0221367E23 |
|
real(kind_phys), parameter | boltz =rgasuniv/avo |
|
real(kind_phys), parameter | cb = 2. |
|
real(kind_phys), parameter | cim = 0.4 |
|
real(kind_phys), parameter | alpha = 0.8 |
|
real(kind_phys), parameter | cin = 2.5 |
|
real(kind_phys), parameter | vv = 0.8 |
|
real(kind_phys), parameter | a_for = 0.1 |
|
real(kind_phys), parameter | a_grs = 0.2 |
|
real(kind_phys), parameter | a_wat = 100. |
|
real(kind_phys), parameter | eps0_for = 0.8*0.01 |
|
real(kind_phys), parameter | eps0_grs = 0.4*0.01 |
|
real(kind_phys), parameter | eps0_wat = 0.6*0.01 |
|
real(kind_phys), parameter | one3 =1.0/3.0 |
|
real(kind_phys), parameter | two3 =2.0/3.0 |
|
real(kind_phys), parameter | sqrt2 =1.4142135623731 |
|
real(kind_phys), parameter | sqrtpi =1.7724539 |
|
real(kind_phys) | karman = 0.4 |
|
real(kind_phys), parameter | conmin = 1.E-16 |
|
real(kind_phys), parameter | pirs =3.14159265358979324 |
|
real(kind_phys), parameter | f6dpi =6.0/pirs |
|
real(kind_phys), parameter | f6dpim9 =1.0E-9*f6dpi |
|
real(kind_phys), parameter | rhosmoke = 1.4E3 |
|
real(kind_phys), parameter | rhodust = 2.6E3 |
|
real(kind_phys), parameter | smokefac =f6dpim9/rhosmoke |
|
real(kind_phys), parameter | dustfac =f6dpim9/rhodust |
|
real(kind_phys), parameter | tss0 =288.15 |
|
real(kind_phys), parameter | sigma1 = 1.8 |
|
real(kind_phys), parameter | mean_diameter1 = 4.e-8 |
|
real(kind_phys), parameter | fact_wfa = 1.e-9*6.0/pirs*exp(4.5*log(sigma1)**2)/mean_diameter1**3 |
|
real(kind_phys), parameter | sginia =2.00 |
|
real(kind_phys), parameter | sginin =1.70 |
|
real(kind_phys), parameter | sginic =2.5 |
|
real(kind_phys), parameter | pss0 =101325.0 |
|
real(kind_phys), parameter | dgmin =1.0E-09 |
|
real(kind_phys), parameter | densmin =1.0E03 |
|
integer, parameter | vdnnuc =1 |
|
integer, parameter | vdnacc =2 |
|
integer, parameter | vdncor =3 |
|
integer, parameter | vdmnuc =4 |
|
integer, parameter | vdmacc =5 |
|
integer, parameter | vdmfine =6 |
|
integer, parameter | vdmcor =7 |
|
integer, parameter | vsnnuc =1 |
|
integer, parameter | vsnacc =2 |
|
integer, parameter | vsncor =3 |
|
integer, parameter | vsmnuc =4 |
|
integer, parameter | vsmacc =5 |
|
integer, parameter | vsmcor =6 |
|
real(kind_phys) | esn04 |
|
real(kind_phys) | esa04 |
|
real(kind_phys) | esc04 |
|
real(kind_phys) | esn05 |
|
real(kind_phys) | esa05 |
|
real(kind_phys) | esn08 |
|
real(kind_phys) | esa08 |
|
real(kind_phys) | esc08 |
|
real(kind_phys) | esn09 |
|
real(kind_phys) | esa09 |
|
real(kind_phys) | esn12 |
|
real(kind_phys) | esa12 |
|
real(kind_phys) | esc12 |
|
real(kind_phys) | esn16 |
|
real(kind_phys) | esa16 |
|
real(kind_phys) | esc16 |
|
real(kind_phys) | esn20 |
|
real(kind_phys) | esa20 |
|
real(kind_phys) | esc20 |
|
real(kind_phys) | esn25 |
|
real(kind_phys) | esa25 |
|
real(kind_phys) | esn24 |
|
real(kind_phys) | esa24 |
|
real(kind_phys) | esc24 |
|
real(kind_phys) | esn28 |
|
real(kind_phys) | esa28 |
|
real(kind_phys) | esc28 |
|
real(kind_phys) | esn32 |
|
real(kind_phys) | esa32 |
|
real(kind_phys) | esc32 |
|
real(kind_phys) | esn36 |
|
real(kind_phys) | esa36 |
|
real(kind_phys) | esc36 |
|
real(kind_phys) | esn49 |
|
real(kind_phys) | esa49 |
|
real(kind_phys) | esn52 |
|
real(kind_phys) | esa52 |
|
real(kind_phys) | esn64 |
|
real(kind_phys) | esa64 |
|
real(kind_phys) | esc64 |
|
real(kind_phys) | esn100 |
|
real(kind_phys) | esnm20 |
|
real(kind_phys) | esam20 |
|
real(kind_phys) | escm20 |
|
real(kind_phys) | esnm32 |
|
real(kind_phys) | esam32 |
|
real(kind_phys) | escm32 |
|
integer, parameter | ngausdv = 7 |
|
real(kind_phys) | xxlsgn |
|
real(kind_phys) | xxlsga |
|
real(kind_phys) | xxlsgc |
|
real(kind_phys), dimension(ngausdv) | y_gq |
|
real(kind_phys), dimension(ngausdv) | wgaus |
|