par_mod Module


Parameters defining KIND parameter for double/single precision



Variables

Type Visibility Attributes Name Initial
integer, public, parameter :: dp = selected_real_kind(P=15)
integer, public, parameter :: sp = selected_real_kind(6)

dep_prec sets the precision for deposition calculations (sp or dp). sp is default, dp can be used for increased precision.


integer, public, parameter :: dep_prec = dp

Set to F to disable use of kernel for concentrations/deposition


logical, public, parameter :: lusekerneloutput = .true.

Set to T to change output units to number of particles per grid cell


logical, public, parameter :: lparticlecountoutput = .false.

number of directories/files used for FLEXPART input/output


integer, public, parameter :: numpath = 4

Physical and other constants


real, public, parameter :: pi = 3.14159265
real, public, parameter :: r_earth = 6.371e6
real, public, parameter :: r_air = 287.05
real, public, parameter :: ga = 9.81
real, public, parameter :: cpa = 1004.6
real, public, parameter :: kappa = 0.286
real, public, parameter :: pi180 = pi/180.
real, public, parameter :: vonkarman = 0.4
real, public, parameter :: rgas = 8.31447
real, public, parameter :: r_water = 461.495
real, public, parameter :: karman = 0.40
real, public, parameter :: href = 15.
real, public, parameter :: convke = 2.0
real, public, parameter :: hmixmin = 100.
real, public, parameter :: hmixmax = 4500.
real, public :: d_trop = 50.
real, public :: d_strat = 0.1
real, public :: fturbmeso = 0.16
real, public, parameter :: rho_water = 1000.
real, public, parameter :: ratio_incloud = 0.005
real, public, parameter :: wet_a = 1.e-5
real, public, parameter :: wet_b = 0.8
real, public, parameter :: xmwml = 18.016/28.960

Constants related to the stratospheric ozone tracer


real, public, parameter :: ozonescale = 60.

Shape parameters from Bagheri & Bonadonna, 2016


real, public, parameter :: pvcrit = 2.

Shape parameters from Bagheri & Bonadonna, 2016


real, public, parameter :: alpha2 = 0.77

Some time constants


real, public, parameter :: beta2 = 0.63

Some time constants


integer, public, parameter :: idiffnorm = 10800

Parameters for polar stereographic projection close to the poles


integer, public, parameter :: idiffmax = 2*idiffnorm

Parameters for polar stereographic projection close to the poles


integer, public, parameter :: minstep = 1

Parameters for polar stereographic projection close to the poles


real, public, parameter :: switchnorth = 75.

Parmaters for GRIB file decoding



Maximum dimensions of the nested input grids


real, public, parameter :: switchsouth = -75.

Parmaters for GRIB file decoding



Maximum dimensions of the nested input grids


integer, public, parameter :: maxnests = 5

Maximum dimensions of the output grid


integer, public, parameter :: nclassunc = 1

Maximum number of particles, species, and similar


real, public, parameter :: minmassfrac = 0.0
integer, public, parameter :: maxtable = 1000
integer, public, parameter :: numclass = 13
integer, public, parameter :: numpf = 1
integer, public, parameter :: numwfmem = 2
integer, public, parameter :: maxndia = 1

Variables for LCM


integer, public, parameter :: maxreagent = 5
integer, public, parameter :: maxrecsample = 2000

dimension of the OH field


integer, public, parameter :: maxxOH = 72

aerosol below-cloud scavenging removal polynomial constants for rain & snow


integer, public, parameter :: maxyOH = 46

aerosol below-cloud scavenging removal polynomial constants for rain & snow


integer, public, parameter :: maxzOH = 7

aerosol below-cloud scavenging removal polynomial constants for rain & snow


real, public, parameter :: bclr_a(4) = (/-6.2609, 0.682, 0.8676, 0.1282/)
real, public, parameter :: bclr_b(7) = (/-14.707, 51.043, -97.306, 97.946, -53.923, 15.311, -1.751/)
real, public, parameter :: bclr_c(2) = (/0.723, 0.0303/)
real, public, parameter :: bclr_e(7) = (/-0.6492, 9.3483, -21.929, 25.317, -15.395, 4.7242, -0.5766/)
real, public, parameter :: bcls_a(7) = (/-4.426, 1.394, -1.202, -3.2942, -1.9521, -0.4904, -0.0457/)
real, public, parameter :: bcls_b(7) = (/-4.3521, -0.7828, 12.768, -19.864, 13.618, -4.4350, 0.5551/)
real, public, parameter :: bcls_c(7) = (/0.5664, 0.0085, -0.1948, -0.6532, -0.5462, -0.1778, -0.0201/)
real, public, parameter :: bcls_e(7) = (/0.5689, -0.0923, 0.0402, 1.4523, -2.078, 1.05, -0.1821/)
integer, public, parameter :: max_cloudthck = 19000
integer, public, parameter :: min_cloudthck = 50
integer, public, parameter :: conv_clrange(2) = (/3000, 6000/)
integer, public, parameter :: highconvp_clrange(2) = (/0, 10000/)
integer, public, parameter :: lowconvp_clrange(2) = (/500, 8000/)
real, public, parameter :: rhmin = 0.90

Maximum number of particles to be released in a single atmospheric column for the domain-filling trajectories option


integer, public, parameter :: maxcolumn = 3000

Dimension of random number field


integer, public, parameter :: maxrand = 6000000

Number of clusters to be used for plume trajectories


integer, public, parameter :: ncluster = 5

Unit numbers for input/output files


integer, public, parameter :: unitpath = 1
integer, public, parameter :: unitcommand = 1
integer, public, parameter :: unitageclasses = 1
integer, public, parameter :: unitgrid = 1
integer, public, parameter :: unitavailab = 1
integer, public, parameter :: unitreleases = 88
integer, public, parameter :: unitpartout = 93
integer, public, parameter :: unitpartout_average = 105
integer, public, parameter :: unitpartoptions = 106
integer, public, parameter :: unitrestart = 106
integer, public, parameter :: unitheightlevels = 107
integer, public, parameter :: unitpartin = 93
integer, public, parameter :: unitflux = 98
integer, public, parameter :: unitouttraj = 96
integer, public, parameter :: unitvert = 1
integer, public, parameter :: unitoro = 1
integer, public, parameter :: unitpoin = 1
integer, public, parameter :: unitreceptor = 1
integer, public, parameter :: unitreceptorout = 2
integer, public, parameter :: unitoutgrid = 97
integer, public, parameter :: unitoutgridppt = 99
integer, public, parameter :: unitoutinfo = 1
integer, public, parameter :: unitspecies = 1
integer, public, parameter :: unitoutrecept = 91
integer, public, parameter :: unitoutreceptppt = 92
integer, public, parameter :: unitoutsatellite = 108
integer, public, parameter :: unitlsm = 1
integer, public, parameter :: unitsfcdata = 1
integer, public, parameter :: unitland = 1
integer, public, parameter :: unitwesely = 1
integer, public, parameter :: unitOH = 1
integer, public, parameter :: unitdates = 94
integer, public, parameter :: unitheader = 90
integer, public, parameter :: unitheader_txt = 100
integer, public, parameter :: unitshortpart = 95
integer, public, parameter :: unitprecip = 101
integer, public, parameter :: unitboundcond = 89
integer, public, parameter :: unittmp = 101
integer, public, parameter :: unitoutfactor = 102

integer code for missing values, used in wet scavenging (PS, 2012)


integer, public, parameter :: icmv = -9999.
logical, public, parameter :: lgridfraction = .false.
logical, public, parameter :: lpartoutputperfield = .false.
integer, public, parameter :: ispeed = 1