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