source: palm/trunk/TESTS/cases/urban_environment_salsa/MONITORING/urban_environment_salsa_rc.000 @ 4859

Last change on this file since 4859 was 4859, checked in by raasch, 22 months ago

urbane environment test case results updated to match r4858

File size: 18.9 KB
Line 
1
2 ******************************    --------------------------------------------
3 * PALM 6.0  Rev: 4858        *    atmosphere - run without 1D - prerun
4 ******************************    --------------------------------------------
5
6 Date:               2021-01-30    Run:       urban_environment_salsa__gfortran_
7 Time:                 12:00:23    Run-No.:   00
8 Run on host:        gfortran_d
9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
10 ------------------------------------------------------------------------------
11
12 Numerical Schemes:
13 -----------------
14
15 --> Use the 1.5-order 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       zeta value range:     -20.00 <= zeta <= 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.4E-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 close to surfaces
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 Land surface model information:
208 ------------------------------
209
210    --> Soil bottom is closed (water content is conserved, default)
211    --> Land surface type  : read from file
212    --> Soil porosity type : read from file
213
214    Initial soil temperature and moisture profile:
215
216       Height:              -0.01      -0.02      -0.05      -0.10      -0.20      -0.40      -0.80      -2.  m
217       Temperature:        277.15     277.00     277.00     277.00     277.00     277.00     277.00     277.  K
218       Moisture:             0.30       0.30       0.30       0.30       0.30       0.30       0.30       0.  m**3/m**3
219       Root fraction:        0.49       0.26       0.25       0.00       0.00       0.00       0.00       0. 
220       Grid point:              0          1          2          3          4          5          6
221
222
223 Vegetation canopy (drag) model:
224 ------------------------------
225
226 Canopy mode: read_from_file               
227 Canopy height:  10.00m (   5 grid points)
228 Leaf drag coefficient:   0.30
229
230 Predefined constant heatflux at the top of the vegetation:   0.00 K m/s
231
232 Characteristic levels of the leaf area density:
233
234 Height:                 0.0     0.0  m
235 Leaf area density:    0.0000    0.00  m**2/m**3
236 Gradient:            ------    0.00  m**2/m**4
237 Gridpoint:                0       0
238
239
240 Radiation model information:
241 ----------------------------
242
243    --> RRTMG scheme is used
244
245    Albedo is set individually for each xy-location, according to given surface type.
246    Timestep: dt_radiation =  60.00  s
247
248
249 SALSA information:
250 ------------------------------
251
252    Starts at: skip_time_do_salsa =      10.00  s
253
254    Timestep: dt_salsa =   2.00  s
255
256    Array shape (z,y,x,bins):
257       aerosol_number:     60   19   19    8
258
259       aerosol_mass:       60   19   19   32
260       (advect_particle_water = T)
261       salsa_gas:    60   19   19    5
262       (salsa_gases_from_chem = F)
263
264    Aerosol dynamic processes included:
265
266       dry deposition (on vegetation = T and on topography = T)
267
268    Aerosol bin subrange limits (in metres):  3.00E-009 1.00E-008 2.50E-006
269    Number of size bins for each aerosol subrange:   1  7
270    Aerosol bin lower limits (in metres):  3.00E-009 1.00E-008 2.20E-008 4.84E-008 1.07E-007 2.35E-007 5.16E-007 1.14E-006
271
272    Bin geometric mean diameters (in metres):  5.48E-009 1.48E-008 3.26E-008 7.18E-008 1.58E-007 3.48E-007 7.66E-007 1.69E-006
273
274    Number of chemical components used: 3
275       Species:    OC    NO    NH                         
276    Initial relative contribution of each species to particle volume in:
277       a-bins:  0.800 0.100 0.100 0.000 0.000 0.000 0.000
278       b-bins:  0.000 0.000 0.000 0.000 0.000 0.000 0.000
279
280    Number of gaseous tracers used: 5
281    Initial gas concentrations:
282       H2SO4:  5.0000E+012 #/m**3
283       HNO3:   3.0000E+012 #/m**3
284       NH3:    6.0000E+012 #/m**3
285       OCNV:   1.0000E+012 #/m**3
286       OCSV:   1.0000E+012 #/m**3
287
288   Initialising concentrations:
289      Aerosol size distribution: init_aerosol_type = 1
290      Gas concentrations: init_gases_type = 1
291
292      Size distribution read from a file.
293
294   Emissions: salsa_emission_mode = read_from_file     
295
296      (currently all emissions are soluble!)
297
298   Boundary conditions for aerosols:
299      bottom/top:      neumann /  dirichlet
300      north/south:      cyclic /     cyclic
301      left/right:    dirichlet /  dirichlet
302 ------------------------------------------------------------------------------
303
304
305
306Spinup control output:---------------------------------
307
308ITER.   HH:MM:SS    DT   PT(z_MO)---------------------------------
309    1  -00:18:00  120.00   288.14
310    2  -00:16:00  120.00   288.17
311    3  -00:14:00  120.00   288.21
312    4  -00:12:00  120.00   288.24
313    5  -00:10:00  120.00   288.27
314    6  -00:08:00  120.00   288.30
315    7  -00:06:00  120.00   288.33
316    8  -00:04:00  120.00   288.35
317    9  -00:02:00  120.00   288.38
318   10   00:00:00  120.00   288.41
319
320
321
322Run-control output:
323------------------
324
325RUN  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
326---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
327  0      0 00:00:00.00   0.7580A  -2.3731D -1.1942D  0.4177  0.060  0.00  -0.147E+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
328  0      1 00:00:00.75   0.7580A  -2.3732  -1.5382   0.7410  0.084  0.33  -0.778E+00      2.  0.137E+01  0.116E-01  0.268E-02  0.892E-03   60  18  15   16  15   0   20  18   0    0.000    0.000       0
329  0      2 00:00:01.51   0.7580A  -2.3733  -1.6029   0.8136  0.083  0.00  -0.796E+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
330  0      3 00:00:02.27   0.7580A  -2.3733  -1.5856   0.8389  0.082  0.00  -0.803E+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
331  0      4 00:00:03.03   0.7580A  -2.3734  -1.5378   0.8459  0.082  0.12  -0.808E+00     12.  0.137E+01  0.132E-01  0.778E-03  0.213E-03   60  17  13   17  15   0   20  17   0    0.000    0.000       0
332  0      5 00:00:03.79   0.7580A  -2.3735  -1.5013   0.8437  0.081  0.16  -0.811E+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
333  0      6 00:00:04.54   0.7580A  -2.3736  -1.4549   0.8375  0.080  0.17  -0.812E+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
334  0      7 00:00:05.30   0.7580A  -2.3736  -1.4322   0.8301  0.080  0.18  -0.812E+00     12.  0.137E+01  0.140E-01  0.759E-03  0.199E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
335  0      8 00:00:06.06   0.7580A  -2.3737  -1.4026   0.8224  0.079  0.18  -0.813E+00     12.  0.137E+01  0.143E-01  0.758E-03  0.197E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
336  0      9 00:00:06.82   0.7580A  -2.3738  -1.3753   0.8146  0.079  0.19  -0.813E+00     12.  0.137E+01  0.146E-01  0.758E-03  0.196E-03   60  16   9   11  14   0   20  17   0    0.000    0.000       0
337  0     10 00:00:07.58   0.7580A  -2.3739  -1.3531   0.8073  0.078  0.24  -0.813E+00     18.  0.137E+01  0.149E-01  0.759E-03  0.196E-03   60  15   8   11  14   0   20  17   0    0.000    0.000       0
338  0     11 00:00:08.33   0.7580A  -2.3739  -1.3349   0.8010  0.077  0.24  -0.814E+00     18.  0.137E+01  0.152E-01  0.761E-03  0.196E-03   60  15   7   11  14   0   20  17   0    0.000    0.000       0
339  0     12 00:00:09.09   0.7580A  -2.3740  -1.3227   0.7963  0.077  0.24  -0.815E+00     18.  0.136E+01  0.156E-01  0.765E-03  0.195E-03   60  15   6   16  14   0   20  17   0    0.000    0.000       0
340  0     13 00:00:09.85   0.7580A  -2.3741  -1.3176   0.7934  0.076  0.25  -0.816E+00     18.  0.136E+01  0.160E-01  0.769E-03  0.195E-03   60  14   6   16  14   0   20  17   0    0.000    0.000       0
341  0     14 00:00:10.61   0.7580A  -2.3742  -1.3089   0.7920  0.076  0.25  -0.817E+00     18.  0.136E+01  0.163E-01  0.774E-03  0.195E-03   60  14   5   16  14   0   20  17   0    0.000    0.000       0
342  0     15 00:00:11.37   0.7580A  -2.3743  -1.2968   0.7912  0.076  0.25  -0.817E+00     18.  0.136E+01  0.167E-01  0.778E-03  0.194E-03   60  13   4   16  14   0   20  17   0    0.000    0.000       0
343  0     16 00:00:12.12   0.7580A  -2.3743  -1.2813   0.7902  0.075  0.25  -0.816E+00     18.  0.136E+01  0.171E-01  0.781E-03  0.193E-03   60  13   3   16  14   0   20  17   0    0.000    0.000       0
344  0     17 00:00:12.88   0.7580A  -2.3744  -1.2624   0.8043  0.075  0.25  -0.816E+00     18.  0.136E+01  0.175E-01  0.784E-03  0.192E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
345  0     18 00:00:13.64   0.7580A  -2.3745  -1.2557   0.8321  0.074  0.25  -0.815E+00     18.  0.136E+01  0.179E-01  0.787E-03  0.191E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
346  0     19 00:00:14.40   0.7580A  -2.3746  -1.2518   0.8628  0.074  0.26  -0.815E+00     18.  0.136E+01  0.183E-01  0.790E-03  0.189E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
347  0     20 00:00:15.16   0.7580A  -2.3746  -1.2492   0.8922  0.074  0.26  -0.815E+00     18.  0.136E+01  0.187E-01  0.793E-03  0.188E-03   60  12   0   11  14   0   10  14  18    0.000    0.000       0
348  0     21 00:00:15.91   0.7580A  -2.3747  -1.2474   0.9185  0.073  0.26  -0.815E+00     18.  0.136E+01  0.191E-01  0.797E-03  0.186E-03   60  11  19   11  14   0   10  14  18    0.000    0.000       0
349  0     22 00:00:16.67   0.7580A  -2.3748  -1.2464   0.9402  0.073  0.26  -0.815E+00     18.  0.136E+01  0.195E-01  0.800E-03  0.185E-03   60  11  18   11  14   0   10  14  18    0.000    0.000       0
350  0     23 00:00:17.43   0.7580A  -2.3749  -1.2460   0.9566  0.073  0.26  -0.815E+00     18.  0.136E+01  0.199E-01  0.803E-03  0.183E-03   60  10  17   11  14   0   10  14  18    0.000    0.000       0
351  0     24 00:00:18.19   0.7580A  -2.3749  -1.2461   0.9677  0.072  0.26  -0.815E+00     18.  0.136E+01  0.203E-01  0.805E-03  0.182E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
352  0     25 00:00:18.95   0.7580A  -2.3750  -1.2465   0.9738  0.072  0.26  -0.815E+00     18.  0.136E+01  0.207E-01  0.809E-03  0.180E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
353  0     26 00:00:19.70   0.7580A  -2.3751  -1.2469   0.9757  0.072  0.27  -0.814E+00     18.  0.136E+01  0.211E-01  0.812E-03  0.179E-03   60  10  15   11  14   0   10  14  18    0.000    0.000       0
354  0     27 00:00:20.46   0.7580A  -2.3752  -1.2472   0.9748  0.072  0.27  -0.814E+00     18.  0.136E+01  0.215E-01  0.816E-03  0.178E-03   60   9  14   11  14   0   11  14  18    0.000    0.000       0
355  0     28 00:00:21.22   0.7580A  -2.3752  -1.2474   0.9803  0.071  0.38  -0.814E+00     40.  0.136E+01  0.219E-01  0.819E-03  0.177E-03   60   9  13   11  14   0   11  14  18    0.000    0.000       0
356  0     29 00:00:21.98   0.7580A  -2.3753  -1.2474   0.9811  0.071  0.38  -0.814E+00     40.  0.136E+01  0.223E-01  0.822E-03  0.176E-03   60   9  13   11  14   0   11  14  18    0.000    0.000       0
357  0     30 00:00:22.74   0.7580A  -2.3754  -1.2473   0.9775  0.071  0.38  -0.815E+00     40.  0.136E+01  0.226E-01  0.826E-03  0.175E-03   60   8  11   11  14   0   11  14  18    0.000    0.000       0
358  0     31 00:00:23.49   0.7580A  -2.3755  -1.2471   0.9703  0.070  0.39  -0.815E+00     40.  0.136E+01  0.230E-01  0.830E-03  0.174E-03   60   8  11   11  14   0   11  14  18    0.000    0.000       0
359  0     32 00:00:24.25   0.7580A  -2.3755  -1.2467   0.9603  0.070  0.39  -0.815E+00     40.  0.136E+01  0.234E-01  0.834E-03  0.173E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
360  0     33 00:00:25.01   0.7580A  -2.3756  -1.2462   0.9497  0.070  0.39  -0.815E+00     40.  0.136E+01  0.238E-01  0.838E-03  0.172E-03   60   7   9   11  14   0   10  14  18    0.000    0.000       0
361  0     34 00:00:25.77   0.7580A  -2.3757  -1.2454   0.9458  0.070  0.39  -0.815E+00     40.  0.136E+01  0.242E-01  0.842E-03  0.171E-03   60   7   8   11  14   0   10  14  18    0.000    0.000       0
362  0     35 00:00:26.53   0.7580A  -2.3758  -1.2442   0.9418  0.069  0.40  -0.815E+00     40.  0.136E+01  0.246E-01  0.845E-03  0.170E-03   60   7   7   11  14   0   10  14  18    0.000    0.000       0
363  0     36 00:00:27.28   0.7580A  -2.3758  -1.2425   0.9449  0.069  0.40  -0.815E+00     40.  0.136E+01  0.250E-01  0.848E-03  0.169E-03   60   6   6   11  14   0   11  14  17    0.000    0.000       0
364  0     37 00:00:28.04   0.7580A  -2.3759  -1.2401   0.9570  0.069  0.40  -0.816E+00     40.  0.136E+01  0.254E-01  0.851E-03  0.168E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
365  0     38 00:00:28.80   0.7580A  -2.3760  -1.2369   0.9688  0.069  0.40  -0.816E+00     40.  0.136E+01  0.257E-01  0.853E-03  0.167E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
366  0     39 00:00:29.56   0.7580A  -2.3761  -1.2328   0.9791  0.069  0.41  -0.816E+00     40.  0.136E+01  0.261E-01  0.856E-03  0.167E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
367  0     40 00:00:30.32   0.7580A  -2.3761  -1.2305   0.9870  0.069  0.41  -0.816E+00     40.  0.136E+01  0.265E-01  0.859E-03  0.166E-03   60   5   3   11  13   0   11  14  17    0.000    0.000       0
368  0     41 00:00:31.07   0.7580A  -2.3762  -1.2321   0.9920  0.068  0.41  -0.817E+00     40.  0.136E+01  0.269E-01  0.861E-03  0.165E-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.