source: flexpart.git/src/com_mod.f90 @ 2bec33e

10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bugunivie
Last change on this file since 2bec33e was 2bec33e, checked in by Sabine <sabine.eckhardt@…>, 7 years ago

resolved kernel merge: removed all uskernel

  • Property mode set to 100644
File size: 36.0 KB
Line 
1!*******************************************************************************
2!        Include file for particle diffusion model FLEXPART                    *
3!        This file contains a global common block used by FLEXPART             *
4!                                                                              *
5!        Author: A. Stohl                                                      *
6!                                                                              *
7!        June 1996                                                             *
8!                                                                              *
9!        Last update:15 August 2013 IP                                         *
10!                                                                              *
11!*******************************************************************************
12
13module com_mod
14
15  use par_mod, only: dp, numpath, maxnests, maxageclass, maxspec, ni, &
16       numclass, nymax, nxmax, maxcolumn, maxwf, nzmax, nxmaxn, nymaxn, &
17       maxreceptor, maxpart, maxrand, nwzmax, nuvzmax, numwfmem
18
19  implicit none
20
21  !****************************************************************
22  ! Variables defining where FLEXPART input/output files are stored
23  !****************************************************************
24
25  character :: path(numpath+2*maxnests)*120
26  integer :: length(numpath+2*maxnests)
27  character(len=256) :: pathfile, flexversion, flexversion_major, arg1, arg2
28  character(len=256) :: ohfields_path
29 
30  ! path                    path names needed for trajectory model
31  ! length                  length of path names needed for trajectory model
32  ! pathfile                file where pathnames are stored
33  ! flexversion             version of flexpart (descriptive long string)
34  ! flexversion_major       version of flexpart (major version number)
35  ! arg                     input arguments from launch at command line
36  ! ohfields_path           path to binary files for OH fields
37
38  !********************************************************
39  ! Variables defining the general model run specifications
40  !********************************************************
41
42  integer :: ibdate,ibtime,iedate,ietime
43  real(kind=dp) :: bdate,edate
44
45
46  ! ibdate                  beginning date (YYYYMMDD)
47  ! ibtime                  beginning time (HHMISS)
48  ! iedate                  ending date (YYYYMMDD)
49  ! ietime                  ending time (HHMISS)
50  ! bdate                   beginning date of simulation (julian date)
51  ! edate                   ending date of simulation (julian date)
52
53
54  integer :: ldirect,ideltas
55
56  ! ldirect                 1 for forward, -1 for backward simulation
57  ! ideltas                 length of trajectory loop from beginning to
58  !                    ending date (s)
59
60  integer :: loutstep,loutaver,loutsample,method,lsynctime
61  real :: outstep
62
63  ! loutstep [s]            gridded concentration output every loutstep seconds
64  ! loutaver [s]            concentration output is an average over [s] seconds
65  ! loutsample [s]          sampling interval of gridded concentration output
66  ! lsynctime [s]           synchronisation time of all particles
67  ! method                  indicator which dispersion method is to be used
68  ! outstep = real(abs(loutstep))
69
70  real :: ctl,fine
71  integer :: ifine,iout,ipout,ipin,iflux,mdomainfill
72  integer :: mquasilag,nested_output,ind_source,ind_receptor
73  integer :: ind_rel,ind_samp,ioutputforeachrelease,linit_cond,surf_only
74  logical :: turbswitch
75  integer :: cblflag !added by mc for cbl
76
77  ! ctl      factor, by which time step must be smaller than Lagrangian time scale
78  ! ifine    reduction factor for time step used for vertical wind
79  !     Langevin equation for the vertical wind component
80  ! ioutputforeachrelease Should each release be a seperate output field?
81  ! iflux    flux calculation options: 1 calculation of fluxes, 2 no fluxes
82  ! iout     output options: 1 conc. output (ng/m3), 2 mixing ratio (pptv), 3 both
83  ! ipout    particle dump options: 0 no, 1 every output interval, 2 only at end
84  ! ipin     read in particle positions from dumped file from a previous run
85  ! fine     real(ifine)
86  ! mdomainfill 0: normal run
87  !        1: particles are initialized according to atmospheric mass distribution
88  ! ind_source switches between different units for concentrations at the source
89  !  NOTE that in backward simulations the release of computational particles
90  !  takes place at the "receptor" and the sampling of particles at the "source".
91  !     1= mass units
92  !     2= mass mixing ratio units
93  ! ind_receptor switches between different units for FLEXPART concentration at the receptor
94  !     1= mass units
95  !     2= mass mixing ratio units
96  ! linit_cond  switch on the output of sensitivity to initial conditions for backward runs
97  !     0=no, 1=mass unit, 2=mass mixing ratio unit
98  ! mquasilag 0: normal run
99  !      1: Particle position output is produced in a condensed format and particles are numbered
100  ! surf_only   switch output in grid_time files for surface only or full vertical resolution
101  !      0=no (full vertical resolution), 1=yes (surface only)
102  ! nested_output: 0 no, 1 yes
103  ! turbswitch              determines how the Markov chain is formulated
104
105  ! ind_rel and ind_samp  are used within the code to change between mass and mass-mix (see readcommand.f)
106  ! cblflag !: 1 activate cbl skewed pdf routines with bi-gaussina pdf whan OL<0 added by mc
107
108
109  integer :: mintime,itsplit
110
111  ! mintime                 minimum time step to be used by FLEXPART
112  ! itsplit                 time constant for splitting particles
113
114  integer :: lsubgrid,lconvection,lagespectra
115
116  ! lsubgrid     1 if subgrid topography parameterization switched on, 2 if not
117  ! lconvection  1 if convection parameterization switched on, 0 if not
118  ! lagespectra  1 if age spectra calculation switched on, 2 if not
119
120  integer :: lnetcdfout
121  ! lnetcdfout   1 for netcdf grid output, 0 if not. Set in COMMAND (namelist input)
122
123  integer :: nageclass,lage(maxageclass)
124
125  ! nageclass               number of ageclasses for the age spectra calculation
126  ! lage [s]                ageclasses for the age spectra calculation
127
128
129  logical :: gdomainfill
130
131  ! gdomainfill             .T., if domain-filling is global, .F. if not
132
133!ZHG SEP 2015 wheather or not to read clouds from GRIB
134  logical :: readclouds=.false.
135!ESO DEC 2015 whether or not both clwc and ciwc are present (if so they are summed)
136  logical :: sumclouds=.false.
137
138  logical,dimension(maxnests) :: readclouds_nest, sumclouds_nest
139 
140
141!NIK 16.02.2015
142  integer(selected_int_kind(16)), dimension(maxspec) :: tot_blc_count=0, &
143       &tot_inc_count=0
144
145
146  !*********************************************************************
147  ! Variables defining the release locations, released species and their
148  ! properties, etc.
149  !*********************************************************************
150
151  !change Sabine Eckhardt, only save the first 1000 identifier for releasepoints
152  character :: compoint(1001)*45
153  integer :: numpoint
154  !sec, now dynamically allocated:
155  ! ireleasestart(maxpoint),ireleaseend(maxpoint)
156  !      real xpoint1(maxpoint),ypoint1(maxpoint)
157  !real xpoint2(maxpoint),ypoint2(maxpoint)
158  !real zpoint1(maxpoint),zpoint2(maxpoint)
159  !integer*2 kindz(maxpoint)
160  integer :: specnum(maxspec)
161  !real xmass(maxpoint,maxspec)
162  real :: decay(maxspec)
163  real :: weta_gas(maxspec),wetb_gas(maxspec)
164  real :: crain_aero(maxspec),csnow_aero(maxspec)
165! NIK: 31.01.2013- parameters for in-cloud scavening
166  real :: ccn_aero(maxspec),in_aero(maxspec)
167  real :: reldiff(maxspec),henry(maxspec),f0(maxspec)
168  real :: density(maxspec),dquer(maxspec),dsigma(maxspec)
169  real :: vsetaver(maxspec),cunningham(maxspec),weightmolar(maxspec)
170  real :: vset(maxspec,ni),schmi(maxspec,ni),fract(maxspec,ni)
171  real :: ri(5,numclass),rac(5,numclass),rcl(maxspec,5,numclass)
172  real :: rgs(maxspec,5,numclass),rlu(maxspec,5,numclass)
173  real :: rm(maxspec),dryvel(maxspec),kao(maxspec)
174  real :: ohcconst(maxspec),ohdconst(maxspec),ohnconst(maxspec)
175  ! se  it is possible to associate a species with a second one to make transfer from gas to aerosol
176  integer :: spec_ass(maxspec)
177
178  real :: area_hour(maxspec,24),point_hour(maxspec,24)
179  real :: area_dow(maxspec,7),point_dow(maxspec,7)
180
181  !integer npart(maxpoint)
182  integer :: nspec,maxpointspec_act
183  character(len=10) :: species(maxspec)
184
185
186  ! compoint                comment, also "name" of each starting point
187  ! numpoint                actual number of trajectory starting/ending points
188  ! ireleasestart,ireleaseend [s] starting and ending time of each release
189  ! xmass                   total mass emitted
190  ! xpoint1,ypoint1         lower left coordinates of release area
191  ! xpoint2,ypoint2         upper right coordinates of release area
192  ! zpoint1,zpoint2         min./max. z-coordinates of release points
193  ! kindz                   1: zpoint is in m agl, 2: zpoint is in m asl
194  ! npart                   number of particles per release point
195  ! nspec                   number of different species allowed for one release
196  ! maxpointspec_act        number of releaspoints for which a different output shall be created
197  ! species                 name of species
198  ! decay                   decay constant of radionuclide
199
200  ! WET DEPOSITION
201  ! weta_gas, wetb_gas     parameters for below-cloud wet scavenging coefficients (gasses)
202  ! crain_aero, csnow_aero parameters for below-cloud wet scavenging coefficients (aerosols)
203  ! ccn_aero, cin_aero     parameters for in-cloud wet scavenging coefficients (aerosols)
204
205  ! GAS DEPOSITION
206  ! reldiff                 diffusivitiy of species relative to diff. of H2O
207  ! henry [M/atm]           Henry constant
208  ! f0                      reactivity relative to that of O3
209  ! ri [s/m]                stomatal resistance
210  ! rcl [s/m]               lower canopy resistance
211  ! rgs [s/m]               ground resistance
212  ! rlu [s/m]               leaf cuticular resistance
213  ! rm [s/m]                mesophyll resistance
214  ! dryvel [m/s]            constant dry deposition velocity
215
216  ! PARTICLE DEPOSITION
217  ! density [kg/m3]         density of particles
218  ! dquer [m]               mean diameter of particles
219  ! dsigma                  dsigma=10 or dsigma=0.1 means that 68% of the
220  !                         mass are between 0.1*dquer and 10*dquer
221
222  ! fract                   mass fraction of each diameter interval
223  ! vset [m/s]              gravitational settling velocity in ni intervals
224  ! cunningham              Cunningham slip correction (strictly valid only near surface)
225  ! vsetaver [m/s]          average gravitational settling velocity
226  ! schmi                   Schmidt number**2/3 of each diameter interval
227  ! weightmolar [g/mol]     molecular weight
228
229  ! TIME VARIATION OF EMISSION
230  ! area_hour, point_hour   daily variation of emission strengths for area and point sources
231  ! area_dow, point_dow     day-of-week variation of emission strengths for area and point sources
232
233
234
235  !**********************************************************
236  ! Variables used for domain-filling trajectory calculations
237  !**********************************************************
238
239  integer :: nx_we(2),ny_sn(2)
240  integer :: numcolumn
241  integer :: numcolumn_we(2,0:nymax-1),numcolumn_sn(2,0:nxmax-1)
242  real :: zcolumn_we(2,0:nymax-1,maxcolumn)
243  real :: zcolumn_sn(2,0:nxmax-1,maxcolumn)
244  real :: xmassperparticle
245  real :: acc_mass_we(2,0:nymax-1,maxcolumn)
246  real :: acc_mass_sn(2,0:nxmax-1,maxcolumn)
247
248  ! nx_we(2)                x indices of western and eastern boundary of domain-filling
249  ! ny_sn(2)                y indices of southern and northern boundary of domain-filling
250  ! numcolumn_we            number of particles to be released within one column
251  !                    at the western and eastern boundary surfaces
252  ! numcolumn_sn            same as numcolumn_we, but for southern and northern domain boundary
253  ! numcolumn               maximum number of particles to be released within a single
254  !                    column
255  ! zcolumn_we              altitudes where particles are to be released
256  !                    at the western and eastern boundary surfaces
257  ! zcolumn_sn              same as zcolumn_we, but for southern and northern domain boundary
258  ! xmassperparticle        air mass per particle in the domain-filling traj. option
259  ! acc_mass_we             mass that has accumulated at the western and eastern boundary;
260  !                    if it exceeds xmassperparticle, a particle is released and
261  !                    acc_mass_we is reduced accordingly
262  ! acc_mass_sn             same as acc_mass_we, but for southern and northern domain boundary
263
264
265
266  !******************************************************************************
267  ! Variables associated with the ECMWF meteorological input data ("wind fields")
268  !******************************************************************************
269
270  integer :: numbwf,wftime(maxwf),lwindinterv
271  character(len=255) :: wfname(maxwf),wfspec(maxwf)
272
273  ! lwindinterv [s]         Interval between wind fields currently in memory
274  ! numbwf                  actual number of wind fields
275  ! wftime(maxwf) [s]       times relative to beginning time of wind fields
276  ! wfname(maxwf)           file names of wind fields
277  ! wfspec(maxwf)           specifications of wind field file, e.g. if on hard
278  !                         disc or on tape
279
280  integer :: memtime(numwfmem),memind(3) ! eso: or memind(numwfmem)
281
282  ! memtime [s]             validation times of wind fields in memory
283  ! memind                  pointer to wind field, in order to avoid shuffling
284  !                         of wind fields
285
286
287
288  !****************************************************************************
289  ! Variables defining actual size and geographical location of the wind fields
290  !****************************************************************************
291
292  integer :: nx,ny,nxmin1,nymin1,nxfield,nuvz,nwz,nz,nmixz,nlev_ec
293  real :: dx,dy,xlon0,ylat0,dxconst,dyconst,height(nzmax)
294
295  ! nx,ny,nz                actual dimensions of wind fields in x,y and z
296  !                    direction, respectively
297  ! nxmin1,nymin1           nx-1, ny-1, respectively
298  ! nuvz,nwz                vertical dimension of original ECMWF data
299  ! nxfield                 same as nx for limited area fields,
300  !                    but for global fields nx=nxfield+1
301  ! nmixz                   number of levels up to maximum PBL height (3500 m)
302
303  ! nuvz is used for u,v components
304  ! nwz is used for w components (staggered grid)
305  ! nz is used for the levels in transformed coordinates (terrain-following Cartesian
306  ! coordinates)
307
308  ! nlev_ec  number of levels ECMWF model
309  ! dx                      grid distance in x direction
310  ! dy                      grid distance in y direction
311  ! dxconst,dyconst         auxiliary variables for utransform,vtransform
312  ! height                  heights of all levels
313  ! xlon0                   geographical longitude and
314  ! ylat0                   geographical latitude of lower left grid point
315
316
317
318  !*************************************************
319  ! Variables used for vertical model discretization
320  !*************************************************
321
322  real :: akm(nwzmax),bkm(nwzmax)
323  real :: akz(nuvzmax),bkz(nuvzmax)
324  real :: aknew(nzmax),bknew(nzmax)
325
326  ! akm,bkm: coeffizients which regulate vertical discretization of ecmwf model
327  !     (at the border of model layers)
328  ! akz,bkz: model discretization coeffizients at the centre of the layers
329  ! aknew,bknew model discretization coeffizients at the interpolated levels
330
331
332
333  ! Fixed fields, unchangeable with time
334  !*************************************
335
336  real :: oro(0:nxmax-1,0:nymax-1)
337  real :: excessoro(0:nxmax-1,0:nymax-1)
338  real :: lsm(0:nxmax-1,0:nymax-1)
339  real :: xlanduse(0:nxmax-1,0:nymax-1,numclass)
340
341  ! oro [m]              orography of the ECMWF model
342  ! excessoro            excess orography mother domain
343  ! lsm                  land sea mask of the ECMWF model
344  ! xlanduse [0-1]       area fractions in percent
345
346  ! 3d fields
347  !**********
348
349  real :: uu(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
350  real :: vv(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
351  real :: uupol(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
352  real :: vvpol(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
353  real :: ww(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
354  real :: tt(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
355  real :: qv(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
356!ZHG adding cloud water
357  real :: clwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem)=0.0 !liquid   [kg/kg]
358  real :: ciwc(0:nxmax-1,0:nymax-1,nzmax,numwfmem)=0.0 !ice      [kg/kg]
359  real :: clw(0:nxmax-1,0:nymax-1,nzmax,numwfmem)=0.0  !combined [m3/m3]
360
361  real :: pv(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
362  real :: rho(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
363  real :: drhodz(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
364  real :: tth(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)
365  real :: qvh(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)
366  real :: clwch(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)=0.0
367  real :: ciwch(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)=0.0
368
369  real :: pplev(0:nxmax-1,0:nymax-1,nuvzmax,numwfmem)
370  !scavenging NIK, PS
371  integer(kind=1) :: clouds(0:nxmax-1,0:nymax-1,nzmax,numwfmem)
372  integer :: cloudsh(0:nxmax-1,0:nymax-1,numwfmem)
373
374!ZHG Sep 2015 
375!   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
376   real :: ctwc(0:nxmax-1,0:nymax-1,numwfmem) ! ESO: =icloud_stats(:,:,4,:)
377
378
379  ! uu,vv,ww [m/2]       wind components in x,y and z direction
380  ! uupol,vvpol [m/s]    wind components in polar stereographic projection
381  ! tt [K]               temperature data
382  ! qv                   specific humidity data
383  ! pv (pvu)             potential vorticity
384  ! rho [kg/m3]          air density
385  ! drhodz [kg/m2]       vertical air density gradient
386  ! tth,qvh              tth,qvh on original eta levels
387  ! clouds:   no cloud, no precipitation   0
388  !      cloud, no precipitation      1
389  !      rainout  conv/lsp dominated  2/3
390  !      washout  conv/lsp dominated  4/5
391  ! PS 2013
392  !c icloudbot (m)        cloud bottom height
393  !c icloudthck (m)       cloud thickness     
394
395  ! pplev for the GFS version
396  ! ctwc                 total cloud water content
397
398  ! 2d fields
399  !**********
400
401  real :: ps(0:nxmax-1,0:nymax-1,1,numwfmem)
402  real :: sd(0:nxmax-1,0:nymax-1,1,numwfmem)
403  real :: msl(0:nxmax-1,0:nymax-1,1,numwfmem)
404  real :: tcc(0:nxmax-1,0:nymax-1,1,numwfmem)
405  real :: u10(0:nxmax-1,0:nymax-1,1,numwfmem)
406  real :: v10(0:nxmax-1,0:nymax-1,1,numwfmem)
407  real :: tt2(0:nxmax-1,0:nymax-1,1,numwfmem)
408  real :: td2(0:nxmax-1,0:nymax-1,1,numwfmem)
409  real :: lsprec(0:nxmax-1,0:nymax-1,1,numwfmem)
410  real :: convprec(0:nxmax-1,0:nymax-1,1,numwfmem)
411  real :: sshf(0:nxmax-1,0:nymax-1,1,numwfmem)
412  real :: ssr(0:nxmax-1,0:nymax-1,1,numwfmem)
413  real :: surfstr(0:nxmax-1,0:nymax-1,1,numwfmem)
414  real :: ustar(0:nxmax-1,0:nymax-1,1,numwfmem)
415  real :: wstar(0:nxmax-1,0:nymax-1,1,numwfmem)
416  real :: hmix(0:nxmax-1,0:nymax-1,1,numwfmem)
417  real :: tropopause(0:nxmax-1,0:nymax-1,1,numwfmem)
418  real :: oli(0:nxmax-1,0:nymax-1,1,numwfmem)
419! real :: diffk(0:nxmax-1,0:nymax-1,1,numwfmem) ESO: this is not in use?
420! logical :: beneath_cloud=.true.
421  ! ps                   surface pressure
422  ! sd                   snow depth
423  ! msl                  mean sea level pressure
424  ! tcc                  total cloud cover
425  ! u10                  10 meter u
426  ! v10                  10 meter v
427  ! tt2                  2 meter temperature
428  ! td2                  2 meter dew point
429  ! lsprec [mm/h]        large scale total precipitation
430  ! convprec [mm/h]      convective precipitation
431  ! sshf                 surface sensible heat flux
432  ! ssr                  surface solar radiation
433  ! surfstr              surface stress
434  ! ustar [m/s]          friction velocity
435  ! wstar [m/s]          convective velocity scale
436  ! hmix  [m]            mixing height
437  ! tropopause [m]       altitude of thermal tropopause
438  ! oli [m]              inverse Obukhov length (1/L)
439  ! diffk [m2/s]         diffusion coefficient at reference height
440
441
442  real :: vdep(0:nxmax-1,0:nymax-1,maxspec,numwfmem)
443
444  ! vdep [m/s]           deposition velocities
445
446
447  !********************************************************************
448  ! Variables associated with the ECMWF input data (nested wind fields)
449  !********************************************************************
450
451  ! NOTE: all nested variables have the same name as the variables used
452  ! for the mother domain, except with a 'n' appended at the end
453  !********************************************************************
454
455  integer :: numbnests
456
457  ! numbnests    number of nested grids
458
459  character(len=255) :: wfnamen(maxnests,maxwf)
460  character(len=18) :: wfspecn(maxnests,maxwf)
461
462  ! wfnamen      nested wind field names
463  ! wfspecn      specifications of wind field file, e.g. if on hard
464  !         disc or on tape
465
466
467  !*********************************************************************
468  ! Variables characterizing size and location of the nested wind fields
469  !*********************************************************************
470
471  integer :: nxn(maxnests),nyn(maxnests)
472  real :: dxn(maxnests),dyn(maxnests),xlon0n(maxnests),ylat0n(maxnests)
473
474  ! nxn,nyn      actual dimensions of nested wind fields in x and y direction
475  ! dxn,dyn      grid distances in x,y direction for the nested grids
476  ! xlon0n       geographical longitude of lower left grid point of nested wind fields
477  ! ylat0n       geographical latitude of lower left grid point of nested wind fields
478
479
480  ! Nested fields, unchangeable with time
481  !**************************************
482
483  real :: oron(0:nxmaxn-1,0:nymaxn-1,maxnests)
484  real :: excessoron(0:nxmaxn-1,0:nymaxn-1,maxnests)
485  real :: lsmn(0:nxmaxn-1,0:nymaxn-1,maxnests)
486  real :: xlandusen(0:nxmaxn-1,0:nymaxn-1,numclass,maxnests)
487
488
489  ! 3d nested fields
490  !*****************
491
492  real,allocatable,dimension(:,:,:,:,:) :: uun, vvn, wwn, ttn, qvn, pvn,&
493       & rhon, drhodzn, tthn, qvhn, clwcn, ciwcn, clwn, clwchn, ciwchn
494  real,allocatable,dimension(:,:,:,:) :: ctwcn
495  integer,allocatable,dimension(:,:,:,:) :: cloudshn
496  integer(kind=1),allocatable,dimension(:,:,:,:,:) :: cloudsn
497
498  ! 2d nested fields
499  !*****************
500
501  real :: psn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
502  real :: sdn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
503  real :: msln(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
504  real :: tccn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
505  real :: u10n(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
506  real :: v10n(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
507  real :: tt2n(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
508  real :: td2n(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
509  real :: lsprecn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
510  real :: convprecn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
511  real :: sshfn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
512  real :: ssrn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
513  real :: surfstrn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
514  real :: ustarn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
515  real :: wstarn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
516  real :: hmixn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
517  real :: tropopausen(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
518  real :: olin(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests)
519  ! real :: diffkn(0:nxmaxn-1,0:nymaxn-1,1,numwfmem,maxnests) ! not in use?
520  real :: vdepn(0:nxmaxn-1,0:nymaxn-1,maxspec,numwfmem,maxnests)
521
522
523  !*************************************************
524  ! Certain auxiliary variables needed for the nests
525  !*************************************************
526
527  real :: xresoln(0:maxnests),yresoln(0:maxnests)
528
529  ! xresoln, yresoln   Factors by which the resolutions in the nests
530  !               are enhanced compared to mother grid
531
532  real :: xln(maxnests),yln(maxnests),xrn(maxnests),yrn(maxnests)
533
534  ! xln,yln,xrn,yrn    Corner points of nested grids in grid coordinates
535  !               of mother grid
536
537
538  !******************************************************
539  ! Variables defining the polar stereographic projection
540  !******************************************************
541
542  logical :: xglobal,sglobal,nglobal
543  real :: switchnorthg,switchsouthg
544
545  !xglobal             T for global fields, F for limited area fields
546  !sglobal             T if domain extends towards south pole
547  !nglobal             T if domain extends towards north pole
548  !switchnorthg,switchsouthg   same as parameters switchnorth,
549  !                    switchsouth, but in grid units
550
551  real :: southpolemap(9),northpolemap(9)
552
553  !southpolemap,northpolemap   define stereographic projections
554  !                    at the two poles
555
556
557  !******************
558  ! Landuse inventory
559  ! Sabine Eckhardt Dec 06: change to new landuse inventary - 11 classes, 1200 x 600 global
560  !******************
561
562  integer(kind=1) :: landinvent(1200,600,6)
563  real :: z0(numclass)
564
565  ! landinvent         landuse inventory (numclass=11 classes)
566  ! z0                  roughness length for the landuse classes
567
568
569
570  !**************************************************************************
571  ! Variables characterizing the output grid and containing the model results
572  !**************************************************************************
573
574  integer :: numxgrid,numygrid,numzgrid
575  real :: dxout,dyout,outlon0,outlat0,xoutshift,youtshift
576  integer :: numxgridn,numygridn
577  real :: dxoutn,dyoutn,outlon0n,outlat0n,xoutshiftn,youtshiftn
578  !real outheight(maxzgrid),outheighthalf(maxzgrid)
579
580  logical :: DEP,DRYDEP,DRYDEPSPEC(maxspec),WETDEP,WETDEPSPEC(maxspec),&
581       & OHREA,ASSSPEC
582  logical :: DRYBKDEP,WETBKDEP
583
584  ! numxgrid,numygrid       number of grid points in x,y-direction
585  ! numxgridn,numygridn     number of grid points in x,y-direction for nested output grid
586  ! numzgrid                number of vertical levels of output grid
587  ! dxout,dyout             grid distance of output grid
588  ! dxoutn,dyoutn           grid distance of nested output grid
589  ! outlon0,outlat0         lower left corner of output grid
590  ! outlon0n,outlat0n       lower left corner of nested output grid
591  ! xoutshift,youtshift     xlon0-outlon0, ylat0-outlat0
592  ! xoutshiftn,youtshiftn   xlon0-outlon0n, ylat0-outlat0n
593  ! outheight [m]           upper levels of the output grid
594  ! outheighthalf [m]       half (middle) levels of the output grid cells
595  ! DEP                     .true., if either dry or wet depos. is switched on
596  ! DRYDEP                  .true., if dry deposition is switched on
597  ! DRYDEPSPEC              .true., if dry deposition is switched on for that species
598  ! WETDEP                  .true., if wet deposition is switched on
599  ! WETDEPSPEC              .true., if wet deposition is switched on for that species
600  ! OHREA                   .true., if OH reaction is switched on
601  ! ASSSPEC                 .true., if there are two species asscoiated
602  ! DRYBKDEP,WETBKDEP        .true., for bkwd runs, where mass deposited and source regions is calculated - either for dry or for wet deposition
603  !                    (i.e. transfer of mass between these two occurs
604
605
606
607  !  if output for each releasepoint shall be created maxpointspec=number of releasepoints
608  !  else maxpointspec is 1 -> moved to unc_mod
609  !  the OUTGRID is moved to the module outg_mod
610  !******************************************************************************
611
612  !real gridunc(0:maxxgrid-1,0:maxygrid-1,maxzgrid,maxspec,
613  !    +             maxpointspec_act,nclassunc,maxageclass)
614  !real griduncn(0:maxxgridn-1,0:maxygridn-1,maxzgrid,maxspec,
615  !    +              maxpointspec_act,nclassunc,maxageclass)
616  !real wetgridunc(0:maxxgrid-1,0:maxygrid-1,maxspec,
617  !    +                maxpointspec_act,nclassunc,maxageclass)
618  !real wetgriduncn(0:maxxgridn-1,0:maxygridn-1,maxspec,
619  !    +ct                 maxpointspec,nclassunc,maxageclass)
620  !real drygridunc(0:maxxgrid-1,0:maxygrid-1,maxspec,maxpointspec,
621  !    +                nclassunc,maxageclass)
622  !real drygriduncn(0:maxxgridn-1,0:maxygridn-1,maxspec,
623  !    +                 maxpointspec,nclassunc,maxageclass)
624
625  !real oroout(0:maxxgrid-1,0:maxygrid-1)
626  !real orooutn(0:maxxgridn-1,0:maxygridn-1)
627  !     real area(0:maxxgrid-1,0:maxygrid-1)
628  !real arean(0:maxxgridn-1,0:maxygridn-1)
629  !real volume(0:maxxgrid-1,0:maxygrid-1,maxzgrid)
630  !real volumen(0:maxxgridn-1,0:maxygridn-1,maxzgrid)
631
632  !real areaeast(0:maxxgrid-1,0:maxygrid-1,maxzgrid)
633  !real areanorth(0:maxxgrid-1,0:maxygrid-1,maxzgrid)
634
635
636  ! gridunc,griduncn        uncertainty of outputted concentrations
637  ! wetgridunc,wetgriduncn  uncertainty of accumulated wet deposited mass on output grid
638  ! drygridunc,drygriduncn  uncertainty of accumulated dry deposited mass on output grid
639  ! oroout,orooutn [m]      height of model topography at output grid
640  ! area,arean [m2]         area of each grid cell
641  ! volume,volumen [m3]     volume of each grid cell
642  ! ... field names with n at the end indicate a nested output grid
643
644
645  !***********************************
646  ! Variables defining receptor points
647  !***********************************
648
649  real :: xreceptor(maxreceptor),yreceptor(maxreceptor)
650  real :: receptorarea(maxreceptor)
651  real :: creceptor(maxreceptor,maxspec)
652  character(len=16) :: receptorname(maxreceptor)
653  integer :: numreceptor
654
655  ! xreceptor,yreceptor     receptor position
656  ! creceptor               concentrations at receptor points
657  ! receptorarea            area of 1*1 grid cell at receptor point
658
659
660
661  !***************************************
662  ! Variables characterizing each particle
663  !***************************************
664
665  integer :: numpart=0
666  integer :: numparticlecount
667
668  integer, allocatable, dimension(:) :: itra1, npoint, nclass, idt, itramem, itrasplit
669
670  real(kind=dp), allocatable, dimension(:) :: xtra1, ytra1
671  real, allocatable, dimension(:) :: ztra1
672  real, allocatable, dimension(:,:) :: xmass1
673  real, allocatable, dimension(:,:) :: xscav_frac1
674
675  ! eso: Moved from timemanager
676  real, allocatable, dimension(:) :: uap,ucp,uzp,us,vs,ws
677  integer(kind=2), allocatable, dimension(:) :: cbt
678
679
680  !CGZ-lifetime
681  real, allocatable, dimension(:,:) ::checklifetime, species_lifetime
682  !CGZ-lifetime
683
684  ! numpart                 actual number of particles in memory
685  ! itra1 (maxpart) [s]     temporal positions of the particles
686  ! npoint(maxpart)         indicates the release point of each particle
687  ! nclass (maxpart)        one of nclassunc classes to which the particle is attributed
688  ! itramem (maxpart) [s]   memorized release times of the particles
689  ! itrasplit (maxpart) [s] next time when particle is to be split into two
690  ! idt(maxpart) [s]        time step to be used for next integration
691  ! numparticlecount        counts the total number of particles that have been released
692  ! xtra1,ytra1,ztra1       spatial positions of the particles
693  ! xmass1 [kg]             particle masses
694  ! xscav_frac1             fraction of particle masse which has been scavenged at receptor
695 
696
697
698  !*******************************************************
699  ! Info table on available chemical species/radionuclides
700  !*******************************************************
701
702  !character*10 specname(maxtable)
703  !real decaytime(maxtable),wetscava(maxtable),wetscavb(maxtable)
704  !real drydiff(maxtable),dryhenry(maxtable),dryactiv(maxtable)
705  !real partrho(maxtable),partmean(maxtable),partsig(maxtable)
706  !real dryvelo(maxtable),weightmol(maxtable),ohreact(maxtable)
707
708  ! specname            Name of chemical species/radionuclide
709  ! decaytime           Half time of radionuclides
710  ! wetscava, wetscavb  Parameters for calculating scavenging coefficients
711  ! drydiff             diffusivitiy of species relative to diff. of H2O
712  ! dryhenry [M/atm]    Henry constant
713  ! dryactiv            reactivity relative to that of O3
714  ! partrho [kg/m3]     density of particles
715  ! partmean [m]        mean diameter of particles
716  ! partsig [m]         mean stand. deviation of particle diameter
717  ! dryvelo [cm/s]      constant dry deposition velocity
718  ! weightmol [g/mol]   molecular weight
719  ! ohreact             OH reaction rate
720
721
722  !********************
723  ! Random number field
724  !********************
725
726  real :: rannumb(maxrand)
727
728  ! rannumb                 field of normally distributed random numbers
729 
730  !********************************************************************
731  ! variables to control stability of CBL scheme under variation
732  ! of statistics in time and space
733  !********************************************************************
734  integer :: nan_count,nan_count2,sum_nan_count(3600),maxtl=1200
735  !added by mc , note that for safety sum_nan_count(N) with N>maxtl
736
737  !********************************************************************
738  ! variables to test well-mixed state of CBL scheme not to be included in final release
739  !********************************************************************
740  real :: well_mixed_vector(50),h_well,well_mixed_norm,avg_air_dens(50),avg_ol,avg_wst,avg_h
741  ! modified by mc to test well-mixed for cbl
742
743  !********************
744  ! Verbosity, testing flags, namelist I/O
745  !********************   
746  integer :: verbosity=0
747  integer :: info_flag=0
748  integer :: count_clock, count_clock0,  count_rate, count_max
749  real    :: tins
750  logical, parameter :: nmlout=.true.
751
752  ! These variables are used to avoid having separate versions of
753  ! files in cases where differences with MPI version are minor (eso)
754  !*****************************************************************
755  integer :: mpi_mode=0 ! .gt. 0 if running MPI version
756  logical :: lroot=.true. ! true if serial version, or if MPI .and. root process
757 
758  logical :: interpolhmix=.false. ! true if the hmix shall be interpolated
759  logical :: turboff=.false.       ! true if the turbulence shall be switched off
760 
761 
762contains
763  subroutine com_mod_allocate_part(nmpart)
764  !*******************************************************************************   
765  ! Dynamic allocation of arrays
766  !
767  ! For FLEXPART version 9.2 and earlier these arrays were statically declared
768  ! with size maxpart. This function is introduced so that the MPI version
769  ! can declare these arrays with smaller size ("maxpart per process"), while
770  ! the serial version allocate at run-time with size maxpart
771  !
772  !*******************************************************************************
773    implicit none
774
775    integer, intent(in) :: nmpart ! maximum number of particles (per process)
776   
777! Arrays, previously static of size maxpart
778    allocate(itra1(nmpart),npoint(nmpart),nclass(nmpart),&
779         & idt(nmpart),itramem(nmpart),itrasplit(nmpart),&
780         & xtra1(nmpart),ytra1(nmpart),ztra1(nmpart),&
781         & xmass1(nmpart, maxspec),&
782         & checklifetime(nmpart,maxspec), species_lifetime(maxspec,2))!CGZ-lifetime
783
784
785    allocate(uap(nmpart),ucp(nmpart),uzp(nmpart),us(nmpart),&
786         & vs(nmpart),ws(nmpart),cbt(nmpart))
787   
788  end subroutine com_mod_allocate_part
789
790
791  subroutine com_mod_allocate_nests
792  !*******************************************************************************   
793  ! Dynamic allocation of arrays
794  !
795  ! For nested wind fields.
796  !
797  !*******************************************************************************
798    implicit none
799
800    allocate(uun(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
801    allocate(vvn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
802    allocate(wwn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
803    allocate(ttn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
804    allocate(qvn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
805    allocate(pvn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
806    allocate(clwcn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
807    allocate(ciwcn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
808    allocate(clwn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
809
810    allocate(cloudsn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
811    allocate(cloudshn(0:nxmaxn-1,0:nymaxn-1,numwfmem,numbnests))
812    allocate(rhon(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
813    allocate(drhodzn(0:nxmaxn-1,0:nymaxn-1,nzmax,numwfmem,numbnests))
814    allocate(tthn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
815    allocate(qvhn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
816    allocate(clwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
817    allocate(ciwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
818    allocate(ctwcn(0:nxmax-1,0:nymax-1,numwfmem,numbnests))
819
820    clwcn(:,:,:,:,:)=0.
821    ciwcn(:,:,:,:,:)=0.
822    clwchn(:,:,:,:,:)=0.
823    ciwchn(:,:,:,:,:)=0.
824   
825  end subroutine com_mod_allocate_nests
826   
827
828end module com_mod
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG