source: palm/trunk/TESTS/cases/urban_environment_salsa/MONITORING/urban_environment_salsa_rc @ 4256

Last change on this file since 4256 was 4256, checked in by monakurppa, 2 years ago

update test case urban_environment_salsa

File size: 19.1 KB
Line 
1
2 ******************************    --------------------------------------------
3 * PALM 6.0  Rev: 4250M       *    atmosphere - run without 1D - prerun
4 ******************************    --------------------------------------------
5
6 Date:               2019-10-07    Run:       urban_environment_salsa           
7 Time:                 12:41:32    Run-No.:   00
8 Run on host:             salsa
9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
10 ------------------------------------------------------------------------------
11
12 Numerical Schemes:
13 -----------------
14
15 --> Use the Moeng_Wyngaard turbulence closure (LES mode).
16 --> Use the boussinesq approximation for the model equations.
17 --> Solve perturbation pressure via FFT using fftw routines
18     perturbation pressure is calculated at every Runge-Kutta step
19 --> Momentum advection via Wicker-Skamarock-Scheme 5th order
20 --> Scalar advection via Wicker-Skamarock-Scheme 5th order
21 --> Loop optimization method: cache
22 --> Time differencing scheme: runge-kutta-3       
23 --> Additional prognostic equation for the specific humidity
24 ------------------------------------------------------------------------------
25
26
27 Run time and time step information:
28 ----------------------------------
29
30 Timestep:             variable     maximum value: 20.000 s    CFL-factor: 0.90
31 Start time:              0.000 s
32 End time:             1231.000 s
33 Spinup time:          1200.000 s
34
35
36 Computational grid and domain size:
37 ----------------------------------
38
39 Grid length:      dx =       2.000 m    dy =       2.000 m
40                dz(1) =       2.000 m
41
42 Domain size:       x =     40.000 m     y =     40.000 m  z(u) =    121.000 m
43
44
45
46 Number of gridpoints (x,y,z):  (0:  19, 0:  19, 0:  61)
47 Subdomain size (x,y,z):        (    10,     10,     62)
48
49
50    Characteristic levels of the geo. wind component ug:
51
52       Height:         0.0     0.0  m
53       ug:            1.00    0.00  m/s
54       Gradient:    ------    0.00  1/100s
55       Gridpoint:        0       0
56
57    Characteristic levels of the geo. wind component vg:
58
59       Height:         0.0     0.0  m
60       vg:            0.00    0.00  m/s
61       Gradient:    ------    0.00  1/100s
62       Gridpoint:        0       0
63
64
65 Topography information:
66 ----------------------
67
68 Topography: read_from_file                         
69 Topography grid definition convention:
70 cell center (scalar grid points)
71
72
73
74 Boundary conditions:
75 -------------------
76
77                     p                    uv                                  pt
78
79 B. bound.: p(0)     = p(1)   | uv(0)     = -uv(1)                | pt(0)     = from soil model 
80 T. bound.: p(nzt+1) = 0      | uv(nzt+1) = ug(nzt+1), vg(nzt+1)  | pt(nzt+1) = pt(nzt) + dpt/dz_
81
82                     e
83
84 B. bound.: e(0)     = e(1)                                                                     
85 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
86
87                     q
88
89 B. bound.: q(0)     = from soil model                                                           
90 T. bound.: q(nzt+1) = q(nzt) + dq/dz                                                           
91
92 Bottom surface fluxes are used in diffusion terms at k=1
93
94
95    Constant flux layer between bottom surface and first computational u,v-level:
96
97       z_mo =   1.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
98       Rif value range:     -20.00 <= rif <= 20.00
99       Predefined surface temperature
100       Predefined surface humidity
101
102
103 Lateral boundaries:
104       left/right:  cyclic             
105       north/south: cyclic             
106
107
108 Initial profiles:
109 ----------------
110
111    Characteristic levels of the initial temperature profile:
112
113       Height:           0.0     0.0  m
114       Temperature:   277.15  277.12  K
115       Gradient:      ------    0.00  K/100m
116       Gridpoint:          0       0
117
118    Characteristic levels of the initial humidity profile:
119
120       Height:           0.0       0.0  m
121       Humidity:     0.1E-02   0.7E-02  kg/kg
122       Gradient:    --------   0.0E+00  (kg/kg)/100m
123       Gridpoint:          0         0
124
125
126 List output:
127 -----------
128
129    1D-Profiles:
130       Output every             9999999.90 s
131       Time averaged over          30.00 s
132       Averaging input every        0.00 s
133
134
135 Data output:
136 -----------
137
138       Time averaged over          30.00 s
139       Averaging input every        0.00 s
140
141    1D-Profiles:
142       Output format: netCDF 64bit offset                     
143
144
145       Profile: u, u*2, wu, w*u*, w"u", v, v*2, wv, w*v*, w"v", theta,       
146
147              : theta*2, wtheta, w*theta*, w"theta", q, q*2, e, e*,         
148
149              : km, l, t_soil, m_soil, rad_lw_in, rad_lw_out, salsa_Ntot,   
150
151              : salsa_N_UFP, salsa_LDSA, salsa_PM0.1, salsa_PM2.5, salsa_PM10
152
153              :                                                             
154       Output every                30.00 s
155       Time averaged over          30.00 s
156       Averaging input every        0.00 s
157
158    Time series:
159       Output format: netCDF 64bit offset                     
160
161       Output every                 5.00 s
162
163 ------------------------------------------------------------------------------
164
165
166 Physical quantities:
167 -------------------
168
169    Geograph. latitude  :   latitude  =  52.5 degr
170    Geograph. longitude :   longitude =  13.3 degr
171    Rotation angle      :   rotation_angle =   0.0 degr
172    Angular velocity    :   omega  = 0.729E-04 rad/s
173    Coriolis parameter  :   f      =  0.000116 1/s
174                            f*     =  0.000089 1/s
175
176    Date and time at model start : 2019-03-06 10:00:00 +00
177
178    Gravity             :   g      =  9.8 m/s**2
179
180    Reference state used in buoyancy terms: initial_profile
181
182
183 Cloud physics quantities / methods:
184 ----------------------------------
185
186    Humidity is considered, bu no condensation
187
188
189 LES / Turbulence quantities:
190 ---------------------------
191
192    Mixing length is limited to 1.80 * z
193
194
195 Actions during the simulation:
196 -----------------------------
197
198    Disturbance impulse (u,v) every :   ****** s
199    Disturbance amplitude           :     0.25 m/s
200    Lower disturbance level         :     5.00 m (GP    3)
201    Upper disturbance level         :    39.00 m (GP   20)
202    Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2
203    Random number generator used    : random-parallel
204
205
206
207 *** dynamic module disabled
208
209
210
211 Land surface model information:
212 ------------------------------
213
214    --> Soil bottom is closed (water content is conserved, default)
215    --> Land surface type  : read from file
216    --> Soil porosity type : read from file
217
218    Initial soil temperature and moisture profile:
219
220       Height:              -0.01      -0.02      -0.05      -0.10      -0.20      -0.40      -0.80      -2.  m
221       Temperature:        277.15     277.00     277.00     277.00     277.00     277.00     277.00     277.  K
222       Moisture:             0.30       0.30       0.30       0.30       0.30       0.30       0.30       0.  m**3/m**3
223       Root fraction:        0.49       0.26       0.25       0.00       0.00       0.00       0.00       0. 
224       Grid point:              0          1          2          3          4          5          6
225
226
227 Vegetation canopy (drag) model:
228 ------------------------------
229
230 Canopy mode: read_from_file_3d             
231 Canopy height:  10.00m (   5 grid points)
232 Leaf drag coefficient:   0.30
233
234 Predefined constant heatflux at the top of the vegetation:   0.00 K m/s
235
236 Characteristic levels of the leaf area density:
237
238 Height:                 0.0     0.0  m
239 Leaf area density:    0.0000    0.00  m**2/m**3
240 Gradient:            ------    0.00  m**2/m**4
241 Gridpoint:                0       0
242
243
244 Radiation model information:
245 ----------------------------
246
247    --> RRTMG scheme is used
248
249    Albedo is set individually for each xy-location, according to given surface type.
250    Timestep: dt_radiation =  60.00  s
251
252
253 SALSA information:
254 ------------------------------
255
256    Starts at: skip_time_do_salsa =      10.00  s
257
258    Timestep: dt_salsa =   2.00  s
259
260    Array shape (z,y,x,bins):
261       aerosol_number:   60 19 19  8
262
263       aerosol_mass:     60 19 19 32
264       (advect_particle_water = T)
265       salsa_gas:  60 19 19  5
266       (salsa_gases_from_chem = F)
267
268    Aerosol dynamic processes included:
269
270       coagulation
271
272       condensation (of precursor gases = T and water vapour = F)
273
274       dry deposition (on vegetation = T and on topography = T)
275
276    Aerosol bin subrange limits (in metres):  3.00E-009 1.00E-008 2.50E-006
277    Number of size bins for each aerosol subrange:   1  7
278    Aerosol bin limits (in metres):  0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000
279    Initial number concentration in bins at the lowest level (#/m**3): 1.19E+008 9.63E+008 5.30E+008 1.96E+008 6.13E+007 6.01E+006 1.95E+005 1.25E+004
280
281    Number of chemical components used: 3
282       Species:    OC    NO    NH                         
283    Initial relative contribution of each species to particle volume in:
284       a-bins:  0.800 0.100 0.100 0.000 0.000 0.000 0.000
285       b-bins:  0.000 0.000 0.000 0.000 0.000 0.000 0.000
286
287    Number of gaseous tracers used: 5
288    Initial gas concentrations:
289       H2SO4:  5.0000E+012 #/m**3
290       HNO3:   3.0000E+012 #/m**3
291       NH3:    6.0000E+012 #/m**3
292       OCNV:   1.0000E+012 #/m**3
293       OCSV:   1.0000E+012 #/m**3
294
295   Initialising concentrations:
296      Aerosol size distribution: init_aerosol_type = 0
297      Gas concentrations: init_gases_type = 0
298      Mode diametres: dpg(nmod) =   0.000  0.000  0.000  0.000  0.000  0.000  0.000 (m)
299      Standard deviation: sigmag(nmod) =    1.80   2.16   2.21   2.00   2.00   2.00   2.00
300      Number concentration: n_lognorm(nmod) =  1.4300E+009 4.4500E+008 7.0000E+004 0.0000E+000 0.0000E+000 0.0000E+000 0.0000E+000 (#/m3)
301
302   Emissions: salsa_emission_mode = read_from_file     
303
304      (currently all emissions are soluble!)
305 ------------------------------------------------------------------------------
306
307
308
309Spinup control output:
310---------------------------------
311
312ITER.   HH:MM:SS    DT   PT(z_MO)
313---------------------------------
314    1  -00:18:00  120.00   288.14
315    2  -00:16:00  120.00   288.17
316    3  -00:14:00  120.00   288.21
317    4  -00:12:00  120.00   288.24
318    5  -00:10:00  120.00   288.27
319    6  -00:08:00  120.00   288.30
320    7  -00:06:00  120.00   288.33
321    8  -00:04:00  120.00   288.35
322    9  -00:02:00  120.00   288.38
323   10   00:00:00  120.00   288.41
324
325
326
327Run-control output:
328------------------
329
330RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
331---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
332  0      0 00:00:00.00   0.7580A  -2.3731D -1.1942D  0.4177  0.064  0.00  -0.164E+00     16.  0.137E+01  0.743E-02  0.234E-01  0.446E-02   60  17  18   19   0  14   20  18   0    0.000    0.000       0
333  0      1 00:00:00.75   0.7580A  -2.3732  -1.5382   0.7410  0.085  0.33  -0.816E+00      2.  0.137E+01  0.116E-01  0.269E-02  0.892E-03   60  18  15   16  15   0   20  18   0    0.000    0.000       0
334  0      2 00:00:01.51   0.7580A  -2.3733  -1.6029   0.8137  0.083  0.00  -0.827E+00      2.  0.137E+01  0.125E-01  0.953E-03  0.351E-03   60  17  14   16  15   0   20  17   0    0.000    0.000       0
335  0      3 00:00:02.27   0.7580A  -2.3733  -1.5856   0.8391  0.083  0.00  -0.842E+00      2.  0.137E+01  0.129E-01  0.818E-03  0.240E-03   60  17  13   17  15   0   20  17   0    0.000    0.000       0
336  0      4 00:00:03.03   0.7580A  -2.3734  -1.5378   0.8462  0.082  0.13  -0.845E+00     12.  0.137E+01  0.132E-01  0.779E-03  0.213E-03   60  17  13   17  15   0   20  17   0    0.000    0.000       0
337  0      5 00:00:03.79   0.7580A  -2.3735  -1.5013   0.8441  0.081  0.17  -0.848E+00     12.  0.137E+01  0.134E-01  0.769E-03  0.205E-03   60  17  12   17  14   0   20  17   0    0.000    0.000       0
338  0      6 00:00:04.54   0.7580A  -2.3736  -1.4550   0.8382  0.081  0.18  -0.849E+00     12.  0.137E+01  0.137E-01  0.763E-03  0.201E-03   60  17  11   17  14   0   20  17   0    0.000    0.000       0
339  0      7 00:00:05.30   0.7580A  -2.3736  -1.4339   0.8311  0.080  0.24  -0.849E+00     18.  0.137E+01  0.140E-01  0.760E-03  0.199E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
340  0      8 00:00:06.06   0.7580A  -2.3737  -1.4045   0.8236  0.080  0.24  -0.849E+00     18.  0.137E+01  0.143E-01  0.759E-03  0.197E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
341  0      9 00:00:06.82   0.7580A  -2.3738  -1.3772   0.8161  0.079  0.25  -0.849E+00     18.  0.137E+01  0.146E-01  0.759E-03  0.196E-03   60  16   9   11  14   0   20  17   0    0.000    0.000       0
342  0     10 00:00:07.58   0.7580A  -2.3739  -1.3551   0.8090  0.078  0.25  -0.849E+00     18.  0.137E+01  0.149E-01  0.760E-03  0.196E-03   60  15   8   11  14   0   20  17   0    0.000    0.000       0
343  0     11 00:00:08.33   0.7580A  -2.3739  -1.3371   0.8029  0.078  0.25  -0.849E+00     18.  0.137E+01  0.153E-01  0.763E-03  0.196E-03   60  15   7   11  14   0   20  17   0    0.000    0.000       0
344  0     12 00:00:09.09   0.7580A  -2.3740  -1.3233   0.7984  0.077  0.26  -0.849E+00     18.  0.136E+01  0.156E-01  0.767E-03  0.195E-03   60  15   6   16  14   0   20  17   0    0.000    0.000       0
345  0     13 00:00:09.85   0.7580A  -2.3741  -1.3182   0.7956  0.077  0.26  -0.849E+00     18.  0.136E+01  0.160E-01  0.772E-03  0.195E-03   60  14   6   16  14   0   20  17   0    0.000    0.000       0
346  0     14 00:00:10.61   0.7580A  -2.3742  -1.3097   0.7943  0.076  0.26  -0.850E+00     18.  0.136E+01  0.164E-01  0.777E-03  0.195E-03   60  14   5   16  14   0   20  17   0    0.000    0.000       0
347  0     15 00:00:11.37   0.7580A  -2.3743  -1.2976   0.7936  0.076  0.26  -0.850E+00     18.  0.136E+01  0.168E-01  0.781E-03  0.194E-03   60  13   4   16  14   0   20  17   0    0.000    0.000       0
348  0     16 00:00:12.12   0.7580A  -2.3743  -1.2822   0.8042  0.076  0.27  -0.849E+00     18.  0.136E+01  0.172E-01  0.784E-03  0.193E-03   60  13   3   16  14   0   11  14  18    0.000    0.000       0
349  0     17 00:00:12.88   0.7580A  -2.3744  -1.2633   0.8294  0.075  0.27  -0.848E+00     18.  0.136E+01  0.176E-01  0.786E-03  0.192E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
350  0     18 00:00:13.64   0.7580A  -2.3745  -1.2575   0.8603  0.075  0.27  -0.847E+00     18.  0.136E+01  0.180E-01  0.790E-03  0.191E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
351  0     19 00:00:14.40   0.7580A  -2.3746  -1.2532   0.8935  0.075  0.27  -0.846E+00     18.  0.136E+01  0.184E-01  0.793E-03  0.190E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
352  0     20 00:00:15.16   0.7580A  -2.3746  -1.2499   0.9249  0.074  0.28  -0.846E+00     18.  0.136E+01  0.188E-01  0.797E-03  0.188E-03   60  12   0   11  14   0   10  14  18    0.000    0.000       0
353  0     21 00:00:15.91   0.7580A  -2.3747  -1.2473   0.9529  0.074  0.28  -0.845E+00     18.  0.136E+01  0.192E-01  0.800E-03  0.187E-03   60  11  19   11  14   0   10  14  18    0.000    0.000       0
354  0     22 00:00:16.67   0.7580A  -2.3748  -1.2454   0.9759  0.074  0.40  -0.845E+00     40.  0.136E+01  0.196E-01  0.804E-03  0.185E-03   60  11  18   11  14   0   10  14  18    0.000    0.000       0
355  0     23 00:00:17.43   0.7580A  -2.3749  -1.2440   0.9934  0.073  0.40  -0.845E+00     40.  0.136E+01  0.200E-01  0.806E-03  0.184E-03   60  10  17   11  14   0   10  14  18    0.000    0.000       0
356  0     24 00:00:18.19   0.7580A  -2.3749  -1.2431   1.0052  0.073  0.40  -0.844E+00     40.  0.136E+01  0.205E-01  0.809E-03  0.182E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
357  0     25 00:00:18.95   0.7580A  -2.3750  -1.2426   1.0120  0.073  0.41  -0.844E+00     40.  0.136E+01  0.209E-01  0.813E-03  0.181E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
358  0     26 00:00:19.70   0.7580A  -2.3751  -1.2423   1.0159  0.072  0.41  -0.843E+00     40.  0.136E+01  0.213E-01  0.816E-03  0.180E-03   60  10  15   11  14   0   11  14  18    0.000    0.000       0
359  0     27 00:00:20.46   0.7580A  -2.3752  -1.2421   1.0277  0.072  0.41  -0.843E+00     40.  0.136E+01  0.217E-01  0.820E-03  0.179E-03   60   9  14   11  14   0   11  14  18    0.000    0.000       0
360  0     28 00:00:21.22   0.7580A  -2.3752  -1.2419   1.0340  0.072  0.41  -0.843E+00     40.  0.136E+01  0.222E-01  0.824E-03  0.177E-03   60   9  13   11  14   0   11  14  18    0.000    0.000       0
361  0     29 00:00:21.98   0.7580A  -2.3753  -1.2417   1.0351  0.072  0.42  -0.843E+00     40.  0.136E+01  0.226E-01  0.828E-03  0.177E-03   60   9  12   11  14   0   11  14  18    0.000    0.000       0
362  0     30 00:00:22.74   0.7580A  -2.3754  -1.2416   1.0316  0.071  0.42  -0.843E+00     40.  0.136E+01  0.230E-01  0.832E-03  0.176E-03   60   8  11   11  14   0   11  14  18    0.000    0.000       0
363  0     31 00:00:23.49   0.7580A  -2.3755  -1.2416   1.0243  0.071  0.42  -0.843E+00     40.  0.136E+01  0.235E-01  0.836E-03  0.175E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
364  0     32 00:00:24.25   0.7580A  -2.3755  -1.2416   1.0146  0.071  0.43  -0.843E+00     40.  0.136E+01  0.239E-01  0.840E-03  0.174E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
365  0     33 00:00:25.01   0.7580A  -2.3756  -1.2415   1.0040  0.071  0.43  -0.843E+00     40.  0.136E+01  0.243E-01  0.845E-03  0.173E-03   60   7   9   11  14   0   11  14  18    0.000    0.000       0
366  0     34 00:00:25.77   0.7580A  -2.3757  -1.2412   0.9939  0.070  0.43  -0.843E+00     40.  0.136E+01  0.247E-01  0.849E-03  0.172E-03   60   7   8   11  14   0   10  14  18    0.000    0.000       0
367  0     35 00:00:26.53   0.7580A  -2.3758  -1.2405   1.0099  0.070  0.44  -0.843E+00     40.  0.136E+01  0.252E-01  0.853E-03  0.172E-03   60   7   7   11  14   0   11  14  17    0.000    0.000       0
368  0     36 00:00:27.28   0.7580A  -2.3759  -1.2391   1.0275  0.070  0.44  -0.843E+00     40.  0.136E+01  0.256E-01  0.857E-03  0.171E-03   60   6   6   11  14   0   11  14  17    0.000    0.000       0
369  0     37 00:00:28.04   0.7580A  -2.3759  -1.2367   1.0435  0.070  0.44  -0.843E+00     40.  0.136E+01  0.260E-01  0.860E-03  0.170E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
370  0     38 00:00:28.80   0.7580A  -2.3760  -1.2333   1.0561  0.070  0.44  -0.843E+00     40.  0.136E+01  0.265E-01  0.863E-03  0.170E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
371  0     39 00:00:29.56   0.7580A  -2.3761  -1.2285   1.0642  0.070  0.45  -0.843E+00     40.  0.136E+01  0.269E-01  0.867E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
372  0     40 00:00:30.32   0.7580A  -2.3762  -1.2224   1.0670  0.069  0.45  -0.844E+00     40.  0.136E+01  0.274E-01  0.871E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
373  0     41 00:00:31.07   0.7580A  -2.3762  -1.2209   1.0647  0.069  0.45  -0.844E+00     40.  0.136E+01  0.278E-01  0.874E-03  0.168E-03   60   5   2   11  13   0   11  14  17    0.000    0.000       0
Note: See TracBrowser for help on using the repository browser.