source: flexpart.git/src/par_mod.f90 @ 79e0349

devrelease-10
Last change on this file since 79e0349 was 79e0349, checked in by Sabine <sabine.eckhardt@…>, 9 months ago

Doing the fit for scavenging when using metfields with no clouds

  • Property mode set to 100644
File size: 13.2 KB
Line 
1!**********************************************************************
2! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
3! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
4! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann   *
5!                                                                     *
6! This file is part of FLEXPART.                                      *
7!                                                                     *
8! FLEXPART is free software: you can redistribute it and/or modify    *
9! it under the terms of the GNU General Public License as published by*
10! the Free Software Foundation, either version 3 of the License, or   *
11! (at your option) any later version.                                 *
12!                                                                     *
13! FLEXPART is distributed in the hope that it will be useful,         *
14! but WITHOUT ANY WARRANTY; without even the implied warranty of      *
15! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       *
16! GNU General Public License for more details.                        *
17!                                                                     *
18! You should have received a copy of the GNU General Public License   *
19! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
20!**********************************************************************
21
22!*******************************************************************************
23!   Include file for calculation of particle trajectories (Program FLEXPART)   *
24!        This file contains the parameter statements used in FLEXPART          *
25!                                                                              *
26!        Author: A. Stohl                                                      *
27!                                                                              *
28!        1997                                                                  *
29!                                                                              *
30!        Update 15 August 2013 IP                                              *
31!                                                                              *
32!        ESO 2016:                                                             *
33!          GFS specific parameters moved to gfs_mod.f90                        *
34!          ECMWF specific parameters moved to ecmwf_mod.f90                    *
35!                                                                              *
36!*******************************************************************************
37
38module par_mod
39
40  implicit none
41
42  !****************************************************************
43  ! Parameters defining KIND parameter for double/single precision
44  !****************************************************************
45
46  integer,parameter :: dp=selected_real_kind(P=15)
47  integer,parameter :: sp=selected_real_kind(6)
48
49  !****************************************************************
50  ! dep_prec sets the precision for deposition calculations (sp or
51  ! dp). sp is default, dp can be used for increased precision.
52  !****************************************************************
53
54  integer,parameter :: dep_prec=sp
55
56  !****************************************************************
57  ! Set to F to disable use of kernel for concentrations/deposition
58  !****************************************************************
59
60  logical, parameter :: lusekerneloutput=.true.
61
62  !*********************************************************************
63  ! Set to T to change output units to number of particles per grid cell
64  !*********************************************************************
65  logical, parameter :: lparticlecountoutput=.false.
66
67  !***********************************************************
68  ! number of directories/files used for FLEXPART input/output
69  !***********************************************************
70
71  integer,parameter :: numpath=4
72
73  ! numpath                 Number of different pathnames for input/output files
74
75
76  !*****************************
77  ! Physical and other constants
78  !*****************************
79
80  real,parameter :: pi=3.14159265, r_earth=6.371e6, r_air=287.05, ga=9.81
81  real,parameter :: cpa=1004.6, kappa=0.286, pi180=pi/180., vonkarman=0.4
82
83  ! pi                      number "pi"
84  ! pi180                   pi/180.
85  ! r_earth                 radius of earth [m]
86  ! r_air                   individual gas constant for dry air [J/kg/K]
87  ! ga                      gravity acceleration of earth [m/s**2]
88  ! cpa                     specific heat for dry air
89  ! kappa                   exponent of formula for potential temperature
90  ! vonkarman               von Karman constant
91
92  real,parameter :: karman=0.40, href=15., convke=2.0
93  real,parameter :: hmixmin=100., hmixmax=4500., turbmesoscale=0.16
94  real,parameter :: d_trop=50., d_strat=0.1
95  real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3]
96  !ZHG MAR2016
97  real,parameter :: incloud_ratio=6.2
98
99  ! karman                  Karman's constant
100  ! href [m]                Reference height for dry deposition
101  ! konvke                  Relative share of kinetic energy used for parcel lifting
102  ! hmixmin,hmixmax         Minimum and maximum allowed PBL height
103  ! turbmesoscale           the factor by which standard deviations of winds at grid
104  !                    points surrounding the particle positions are scaled to
105  !                    yield the scales for the mesoscale wind velocity fluctuations
106  ! d_trop [m2/s]           Turbulent diffusivity for horizontal components in the troposphere
107  ! d_strat [m2/s]          Turbulent diffusivity for vertical component in the stratosphere
108
109  real,parameter :: xmwml=18.016/28.960
110
111  ! xmwml   ratio of molar weights of water vapor and dry air
112  !****************************************************
113  ! Constants related to the stratospheric ozone tracer
114  !****************************************************
115
116  real,parameter :: ozonescale=60., pvcrit=2.0
117
118  ! ozonescale              ppbv O3 per PV unit
119  ! pvcrit                  PV level of the tropopause
120
121
122
123  !********************
124  ! Some time constants
125  !********************
126
127  integer,parameter :: idiffnorm=10800, idiffmax=2*idiffnorm, minstep=1
128
129  ! idiffnorm [s]           normal time interval between two wind fields
130  ! idiffmax [s]            maximum time interval between two wind fields
131  ! minstep [s]             minimum time step to be used within FLEXPART
132
133
134  !*****************************************************************
135  ! Parameters for polar stereographic projection close to the poles
136  !*****************************************************************
137
138  real,parameter :: switchnorth=75., switchsouth=-75.
139
140  ! switchnorth    use polar stereographic grid north of switchnorth
141  ! switchsouth    use polar stereographic grid south of switchsouth
142
143
144  !*********************************************
145  ! Maximum dimensions of the input mother grids
146  !*********************************************
147 
148  ! ECMWF
149! integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=359 ! 1.0 degree 92 level
150!   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 1.0 degree 138 level
151! integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359  ! 0.5 degree 138 level
152!  integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0  ! CERA 2.0 degree 92 level
153
154! GFS
155   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
156
157
158  !*********************************************
159  ! Maximum dimensions of the nested input grids
160  !*********************************************
161
162  integer,parameter :: maxnests=0,nxmaxn=0,nymaxn=0
163
164  ! nxmax,nymax        maximum dimension of wind fields in x and y
165  !                    direction, respectively
166  ! nuvzmax,nwzmax     maximum dimension of (u,v) and (w) wind fields in z
167  !                    direction (for fields on eta levels)
168  ! nzmax              maximum dimension of wind fields in z direction
169  !                    for the transformed Cartesian coordinates
170  ! nxshift            for global grids (in x), the grid can be shifted by
171  !                    nxshift grid points, in order to accomodate nested
172  !                    grids, and output grids overlapping the domain "boundary"
173  !                    nxshift must not be negative; "normal" setting would be 0
174
175 
176  integer,parameter :: nconvlevmax = nuvzmax-1
177  integer,parameter :: na = nconvlevmax+1
178
179  ! ntracermax         maximum number of tracer species in convection
180  ! nconvlevmax        maximum number of levels for convection
181  ! na                 parameter used in Emanuel's convect subroutine
182
183
184  !*********************************
185  ! Parmaters for GRIB file decoding
186  !*********************************
187
188  integer,parameter :: jpack=4*nxmax*nymax, jpunp=4*jpack
189
190  ! jpack,jpunp             maximum dimensions needed for GRIB file decoding
191
192
193  !**************************************
194  ! Maximum dimensions of the output grid
195  !**************************************
196
197  !integer,parameter :: maxageclass=1,maxzgrid=10,nclassunc=1
198  integer,parameter :: maxageclass=1,nclassunc=1
199
200  ! nclassunc               number of classes used to calculate the uncertainty
201  !                         of the output
202  ! maxageclass             maximum number of age classes used for output
203
204  ! Sabine Eckhardt, June, 2008
205  ! the dimensions of the OUTGRID are now set dynamically during runtime
206  ! maxxgrid,maxygrid,maxzgrid    maximum dimensions in x,y,z direction
207  ! maxxgridn,maxygridn           maximum dimension of the nested grid
208  !integer maxxgrid,maxygrid,maxzgrid,maxxgridn,maxygridn
209  !integer,parameter :: maxxgrid=361,maxygrid=181,maxxgridn=0,maxygridn=0)
210
211  integer,parameter :: maxreceptor=20
212
213  ! maxreceptor             maximum number of receptor points
214
215
216  !**************************************************
217  ! Maximum number of particles, species, and similar
218  !**************************************************
219
220  integer,parameter :: maxpart=100000
221  integer,parameter :: maxspec=1
222
223  real,parameter :: minmass=0.0001
224
225  ! maxpart                 Maximum number of particles
226  ! maxspec                 Maximum number of chemical species per release
227  ! minmass                 Terminate particles carrying less mass
228
229  ! maxpoint is also set dynamically during runtime
230  ! maxpoint                Maximum number of release locations
231
232  ! ---------
233  ! Sabine Eckhardt: change of landuse inventary numclass=13
234  ! ---------
235  integer,parameter :: maxwf=50000, maxtable=1000, numclass=13, ni=11
236  integer,parameter :: numwfmem=2 ! Serial version/MPI with 2 fields
237  !integer,parameter :: numwfmem=3 ! MPI with 3 fields
238
239  ! maxwf                   maximum number of wind fields to be used for simulation
240  ! maxtable                Maximum number of chemical species that can be
241  !                         tabulated for FLEXPART
242  ! numclass                Number of landuse classes available to FLEXPART
243  ! ni                      Number of diameter classes of particles
244  ! numwfmem                Number of windfields kept in memory. 2 for serial
245  !                         version, 2 or 3 for MPI version
246
247  !**************************************************************************
248  ! dimension of the OH field
249  !**************************************************************************
250  integer,parameter :: maxxOH=72, maxyOH=46, maxzOH=7
251
252  !**************************************************************************
253  ! Maximum number of particles to be released in a single atmospheric column
254  ! for the domain-filling trajectories option
255  !**************************************************************************
256
257  integer,parameter :: maxcolumn=3000
258
259
260  !*********************************
261  ! Dimension of random number field
262  !*********************************
263
264  integer,parameter :: maxrand=1000000
265
266  ! maxrand                 number of random numbers used
267 
268
269  !*****************************************************
270  ! Number of clusters to be used for plume trajectories
271  !*****************************************************
272
273  integer,parameter :: ncluster=5
274
275  !************************************
276  ! Unit numbers for input/output files
277  !************************************
278
279  integer,parameter :: unitpath=1, unitcommand=1, unitageclasses=1, unitgrid=1
280  integer,parameter :: unitavailab=1, unitreleases=88, unitpartout=93
281  integer,parameter :: unitpartin=93, unitflux=98, unitouttraj=96
282  integer,parameter :: unitvert=1, unitoro=1, unitpoin=1, unitreceptor=1
283  integer,parameter :: unitoutgrid=97, unitoutgridppt=99, unitoutinfo=1
284  integer,parameter :: unitspecies=1, unitoutrecept=91, unitoutreceptppt=92
285  integer,parameter :: unitlsm=1, unitsurfdata=1, unitland=1, unitwesely=1
286  integer,parameter :: unitOH=1
287  integer,parameter :: unitdates=94, unitheader=90,unitheader_txt=100, unitshortpart=95, unitprecip=101
288  integer,parameter :: unitboundcond=89
289  integer,parameter :: unittmp=101
290
291!******************************************************
292! integer code for missing values, used in wet scavenging (PS, 2012)
293!******************************************************
294
295  integer,parameter ::  icmv=-9999
296
297
298end module par_mod
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG