Home

Context Navigation

urban_environment_salsa_rc.000 @ 4494

Last change on this file since 4494 was 4494, checked in by raasch, 4 years ago
some test case results updated to match last commit
File size: 18.8 KB

Line
1
2	****************************** --------------------------------------------
3	* PALM 6.0 Rev: 4493 * atmosphere - run without 1D - prerun
4	****************************** --------------------------------------------
5
6	Date: 2020-04-11 Run: urban_environment_salsa__gfortran_
7	Time: 12:57:13 Run-No.: 00
8	Run on host: gfortran_d
9	Number of PEs: 2 Processor grid (x,y): ( 1, 2) calculated
10	A 1d-decomposition along y is used
11	------------------------------------------------------------------------------
12
13	Numerical Schemes:
14	-----------------
15
16	--> Use the 1.5-order turbulence closure (LES mode).
17	--> Use the boussinesq approximation for the model equations.
18	--> Solve perturbation pressure via FFT using fftw routines
19	perturbation pressure is calculated at every Runge-Kutta step
20	--> Momentum advection via Wicker-Skamarock-Scheme 5th order
21	--> Scalar advection via Wicker-Skamarock-Scheme 5th order
22	--> Loop optimization method: cache
23	--> Time differencing scheme: runge-kutta-3
24	--> Additional prognostic equation for the specific humidity
25	------------------------------------------------------------------------------
26
27
28	Run time and time step information:
29	----------------------------------
30
31	Timestep: variable maximum value: 20.000 s CFL-factor: 0.90
32	Start time: 0.000 s
33	End time: 1231.000 s
34	Spinup time: 1200.000 s
35
36
37	Computational grid and domain size:
38	----------------------------------
39
40	Grid length: dx = 2.000 m dy = 2.000 m
41	dz(1) = 2.000 m
42
43	Domain size: x = 40.000 m y = 40.000 m z(u) = 121.000 m
44
45
46
47	Number of gridpoints (x,y,z): (0: 19, 0: 19, 0: 61)
48	Subdomain size (x,y,z): ( 20, 10, 62)
49
50
51	Characteristic levels of the geo. wind component ug:
52
53	Height: 0.0 0.0 m
54	ug: 1.00 0.00 m/s
55	Gradient: ------ 0.00 1/100s
56	Gridpoint: 0 0
57
58	Characteristic levels of the geo. wind component vg:
59
60	Height: 0.0 0.0 m
61	vg: 0.00 0.00 m/s
62	Gradient: ------ 0.00 1/100s
63	Gridpoint: 0 0
64
65
66	Topography information:
67	----------------------
68
69	Topography: read_from_file
70	Topography grid definition convention:
71	cell center (scalar grid points)
72
73
74
75	Boundary conditions:
76	-------------------
77
78	p uv pt
79
80	B. bound.: p(0) = p(1) \| uv(0) = -uv(1) \| pt(0) = from soil model
81	T. bound.: p(nzt+1) = 0 \| uv(nzt+1) = ug(nzt+1), vg(nzt+1) \| pt(nzt+1) = pt(nzt) + dpt/dz_
82
83	e
84
85	B. bound.: e(0) = e(1)
86	T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)
87
88	q
89
90	B. bound.: q(0) = from soil model
91	T. bound.: q(nzt+1) = q(nzt) + dq/dz
92
93	Bottom surface fluxes are used in diffusion terms at k=1
94
95
96	Constant flux layer between bottom surface and first computational u,v-level:
97
98	z_mo = 1.00 m z0 = 0.1000 m z0h = 0.10000 m kappa = 0.40
99	Rif value range: -20.00 <= rif <= 20.00
100	Predefined surface temperature
101	Predefined surface humidity
102
103
104	Lateral boundaries:
105	left/right: cyclic
106	north/south: cyclic
107
108
109	Initial profiles:
110	----------------
111
112	Characteristic levels of the initial temperature profile:
113
114	Height: 0.0 0.0 m
115	Temperature: 277.15 277.12 K
116	Gradient: ------ 0.00 K/100m
117	Gridpoint: 0 0
118
119	Characteristic levels of the initial humidity profile:
120
121	Height: 0.0 0.0 m
122	Humidity: 0.1E-02 0.4E-02 kg/kg
123	Gradient: -------- 0.0E+00 (kg/kg)/100m
124	Gridpoint: 0 0
125
126
127	List output:
128	-----------
129
130	1D-Profiles:
131	Output every 9999999.90 s
132	Time averaged over 30.00 s
133	Averaging input every 0.00 s
134
135
136	Data output:
137	-----------
138
139	Time averaged over 30.00 s
140	Averaging input every 0.00 s
141
142	1D-Profiles:
143	Output format: netCDF 64bit offset
144
145
146	Profile: u, u2, wu, wu, w"u", v, v2, wv, wv, w"v", theta,
147
148	: theta2, wtheta, wtheta, w"theta", q, q2, e, e*,
149
150	: km, l, t_soil, m_soil, rad_lw_in, rad_lw_out, salsa_Ntot,
151
152	: salsa_N_UFP, salsa_LDSA, salsa_PM0.1, salsa_PM2.5, salsa_PM10
153
154	:
155	Output every 30.00 s
156	Time averaged over 30.00 s
157	Averaging input every 0.00 s
158
159	Time series:
160	Output format: netCDF 64bit offset
161
162	Output every 5.00 s
163
164	------------------------------------------------------------------------------
165
166
167	Physical quantities:
168	-------------------
169
170	Geograph. latitude : latitude = 52.5 degr
171	Geograph. longitude : longitude = 13.3 degr
172	Rotation angle : rotation_angle = 0.0 degr
173	Angular velocity : omega = 0.729E-04 rad/s
174	Coriolis parameter : f = 0.000116 1/s
175	f* = 0.000089 1/s
176
177	Date and time at model start : 2019-03-06 10:00:00 +00
178
179	Gravity : g = 9.8 m/s**2
180
181	Reference state used in buoyancy terms: initial_profile
182
183
184	Cloud physics quantities / methods:
185	----------------------------------
186
187	Humidity is considered, bu no condensation
188
189
190	LES / Turbulence quantities:
191	---------------------------
192
193	Mixing length is limited close to surfaces
194
195
196	Actions during the simulation:
197	-----------------------------
198
199	Disturbance impulse (u,v) every : ****** s
200	Disturbance amplitude : 0.25 m/s
201	Lower disturbance level : 5.00 m (GP 3)
202	Upper disturbance level : 39.00 m (GP 20)
203	Disturbances cease as soon as the disturbance energy exceeds 0.010 m2/s2
204	Random number generator used : random-parallel
205
206
207
208	*** dynamic module disabled
209
210
211
212	Land surface model information:
213	------------------------------
214
215	--> Soil bottom is closed (water content is conserved, default)
216	--> Land surface type : read from file
217	--> Soil porosity type : read from file
218
219	Initial soil temperature and moisture profile:
220
221	Height: -0.01 -0.02 -0.05 -0.10 -0.20 -0.40 -0.80 -2. m
222	Temperature: 277.15 277.00 277.00 277.00 277.00 277.00 277.00 277. K
223	Moisture: 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0. m3/m3
224	Root fraction: 0.49 0.26 0.25 0.00 0.00 0.00 0.00 0.
225	Grid point: 0 1 2 3 4 5 6
226
227
228	Vegetation canopy (drag) model:
229	------------------------------
230
231	Canopy mode: read_from_file
232	Canopy height: 10.00m ( 5 grid points)
233	Leaf drag coefficient: 0.30
234
235	Predefined constant heatflux at the top of the vegetation: 0.00 K m/s
236
237	Characteristic levels of the leaf area density:
238
239	Height: 0.0 0.0 m
240	Leaf area density: 0.0000 0.00 m2/m3
241	Gradient: ------ 0.00 m2/m4
242	Gridpoint: 0 0
243
244
245	Radiation model information:
246	----------------------------
247
248	--> RRTMG scheme is used
249
250	Albedo is set individually for each xy-location, according to given surface type.
251	Timestep: dt_radiation = 60.00 s
252
253
254	SALSA information:
255	------------------------------
256
257	Starts at: skip_time_do_salsa = 10.00 s
258
259	Timestep: dt_salsa = 2.00 s
260
261	Array shape (z,y,x,bins):
262	aerosol_number: 60 19 19 8
263
264	aerosol_mass: 60 19 19 32
265	(advect_particle_water = T)
266	salsa_gas: 60 19 19 5
267	(salsa_gases_from_chem = F)
268
269	Aerosol dynamic processes included:
270
271	dry deposition (on vegetation = T and on topography = T)
272
273	Aerosol bin subrange limits (in metres): 3.00E-009 1.00E-008 2.50E-006
274	Number of size bins for each aerosol subrange: 1 7
275	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
276
277	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
278
279	Number of chemical components used: 3
280	Species: OC NO NH
281	Initial relative contribution of each species to particle volume in:
282	a-bins: 0.800 0.100 0.100 0.000 0.000 0.000 0.000
283	b-bins: 0.000 0.000 0.000 0.000 0.000 0.000 0.000
284
285	Number of gaseous tracers used: 5
286	Initial gas concentrations:
287	H2SO4: 5.0000E+012 #/m**3
288	HNO3: 3.0000E+012 #/m**3
289	NH3: 6.0000E+012 #/m**3
290	OCNV: 1.0000E+012 #/m**3
291	OCSV: 1.0000E+012 #/m**3
292
293	Initialising concentrations:
294	Aerosol size distribution: init_aerosol_type = 1
295	Gas concentrations: init_gases_type = 1
296
297	Size distribution read from a file.
298
299	Emissions: salsa_emission_mode = read_from_file
300
301	(currently all emissions are soluble!)
302	------------------------------------------------------------------------------
303
304
305
306	Spinup control output:
307	---------------------------------
308
309	ITER. HH:MM:SS DT PT(z_MO)
310	---------------------------------
311	1 -00:18:00 120.00 288.14
312	2 -00:16:00 120.00 288.17
313	3 -00:14:00 120.00 288.21
314	4 -00:12:00 120.00 288.24
315	5 -00:10:00 120.00 288.27
316	6 -00:08:00 120.00 288.30
317	7 -00:06:00 120.00 288.33
318	8 -00:04:00 120.00 288.35
319	9 -00:02:00 120.00 288.38
320	10 00:00:00 120.00 288.41
321
322
323
324	Run-control output:
325	------------------
326
327	RUN 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
328	---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
329	0 0 00:00:00.00 0.7580A -2.3731D -1.1942D 0.4177 0.064 0.00 -0.207E+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
330	0 1 00:00:00.75 0.7580A -2.3732 -1.5382 0.7410 0.085 0.33 -0.907E+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
331	0 2 00:00:01.51 0.7580A -2.3733 -1.6029 0.8137 0.084 0.00 -0.928E+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
332	0 3 00:00:02.27 0.7580A -2.3733 -1.5857 0.8392 0.083 0.00 -0.950E+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
333	0 4 00:00:03.03 0.7580A -2.3734 -1.5379 0.8466 0.083 0.13 -0.960E+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
334	0 5 00:00:03.79 0.7580A -2.3735 -1.5015 0.8447 0.082 0.17 -0.967E+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
335	0 6 00:00:04.54 0.7580A -2.3736 -1.4553 0.8392 0.081 0.23 -0.971E+00 18. 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
336	0 7 00:00:05.30 0.7580A -2.3736 -1.4292 0.8324 0.081 0.24 -0.974E+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
337	0 8 00:00:06.06 0.7580A -2.3737 -1.3992 0.8253 0.080 0.25 -0.975E+00 18. 0.137E+01 0.143E-01 0.759E-03 0.198E-03 60 16 10 11 14 0 20 17 0 0.000 0.000 0
338	0 9 00:00:06.82 0.7580A -2.3738 -1.3712 0.8181 0.080 0.26 -0.976E+00 18. 0.137E+01 0.146E-01 0.760E-03 0.197E-03 60 16 9 11 14 0 20 17 0 0.000 0.000 0
339	0 10 00:00:07.58 0.7580A -2.3739 -1.3483 0.8112 0.079 0.26 -0.977E+00 18. 0.137E+01 0.150E-01 0.761E-03 0.196E-03 60 15 8 11 14 0 20 17 0 0.000 0.000 0
340	0 11 00:00:08.33 0.7580A -2.3739 -1.3295 0.8054 0.079 0.26 -0.978E+00 18. 0.137E+01 0.153E-01 0.764E-03 0.196E-03 60 15 7 11 14 0 20 17 0 0.000 0.000 0
341	0 12 00:00:09.09 0.7580A -2.3740 -1.3242 0.8010 0.078 0.27 -0.979E+00 18. 0.136E+01 0.157E-01 0.768E-03 0.196E-03 60 14 6 16 14 0 20 17 0 0.000 0.000 0
342	0 13 00:00:09.85 0.7580A -2.3741 -1.3192 0.7984 0.078 0.27 -0.980E+00 18. 0.136E+01 0.161E-01 0.773E-03 0.196E-03 60 14 6 16 14 0 20 17 0 0.000 0.000 0
343	0 14 00:00:10.61 0.7580A -2.3742 -1.3107 0.7972 0.077 0.28 -0.980E+00 18. 0.136E+01 0.165E-01 0.778E-03 0.195E-03 60 14 5 16 14 0 20 17 0 0.000 0.000 0
344	0 15 00:00:11.37 0.7580A -2.3743 -1.2987 0.8183 0.077 0.28 -0.981E+00 18. 0.136E+01 0.169E-01 0.783E-03 0.195E-03 60 13 4 16 14 0 11 14 18 0.000 0.000 0
345	0 16 00:00:12.12 0.7580A -2.3743 -1.2833 0.8548 0.076 0.35 -0.980E+00 28. 0.136E+01 0.173E-01 0.785E-03 0.194E-03 60 13 3 16 14 0 11 14 18 0.000 0.000 0
346	0 17 00:00:12.88 0.7580A -2.3744 -1.2644 0.8840 0.076 0.35 -0.979E+00 28. 0.136E+01 0.177E-01 0.788E-03 0.193E-03 60 12 2 16 14 0 11 14 18 0.000 0.000 0
347	0 18 00:00:13.64 0.7580A -2.3745 -1.2476 0.9195 0.076 0.36 -0.978E+00 28. 0.136E+01 0.181E-01 0.792E-03 0.192E-03 60 12 1 11 14 0 10 14 18 0.000 0.000 0
348	0 19 00:00:14.40 0.7580A -2.3746 -1.2430 0.9579 0.075 0.42 -0.978E+00 40. 0.136E+01 0.186E-01 0.796E-03 0.190E-03 60 12 1 11 14 0 10 14 18 0.000 0.000 0
349	0 20 00:00:15.16 0.7580A -2.3746 -1.2392 0.9938 0.075 0.42 -0.977E+00 40. 0.136E+01 0.190E-01 0.799E-03 0.189E-03 60 12 0 11 14 0 10 14 18 0.000 0.000 0
350	0 21 00:00:15.91 0.7580A -2.3747 -1.2361 1.0254 0.075 0.43 -0.976E+00 40. 0.136E+01 0.195E-01 0.802E-03 0.187E-03 60 11 19 11 14 0 10 14 18 0.000 0.000 0
351	0 22 00:00:16.67 0.7580A -2.3748 -1.2334 1.0514 0.074 0.43 -0.975E+00 40. 0.136E+01 0.200E-01 0.806E-03 0.186E-03 60 11 18 11 14 0 10 14 18 0.000 0.000 0
352	0 23 00:00:17.43 0.7580A -2.3749 -1.2312 1.0713 0.074 0.44 -0.975E+00 40. 0.136E+01 0.204E-01 0.809E-03 0.185E-03 60 10 17 11 14 0 10 14 18 0.000 0.000 0
353	0 24 00:00:18.19 0.7580A -2.3749 -1.2294 1.0868 0.074 0.44 -0.974E+00 40. 0.136E+01 0.209E-01 0.812E-03 0.183E-03 60 10 16 11 14 0 11 14 18 0.000 0.000 0
354	0 25 00:00:18.95 0.7580A -2.3750 -1.2278 1.1142 0.074 0.45 -0.973E+00 40. 0.136E+01 0.214E-01 0.817E-03 0.182E-03 60 10 16 11 14 0 11 14 18 0.000 0.000 0
355	0 26 00:00:19.70 0.7580A -2.3751 -1.2266 1.1359 0.073 0.45 -0.973E+00 40. 0.136E+01 0.219E-01 0.821E-03 0.181E-03 60 10 15 11 14 0 11 14 18 0.000 0.000 0
356	0 27 00:00:20.46 0.7580A -2.3752 -1.2256 1.1510 0.073 0.46 -0.972E+00 40. 0.136E+01 0.224E-01 0.825E-03 0.180E-03 60 9 14 11 14 0 11 14 18 0.000 0.000 0
357	0 28 00:00:21.22 0.7580A -2.3752 -1.2247 1.1596 0.073 0.46 -0.972E+00 40. 0.136E+01 0.229E-01 0.829E-03 0.179E-03 60 9 13 11 14 0 11 14 18 0.000 0.000 0
358	0 29 00:00:21.98 0.7580A -2.3753 -1.2237 1.1624 0.073 0.46 -0.972E+00 40. 0.136E+01 0.234E-01 0.834E-03 0.179E-03 60 9 12 11 14 0 11 14 18 0.000 0.000 0
359	0 30 00:00:22.74 0.7580A -2.3754 -1.2224 1.1607 0.072 0.47 -0.971E+00 40. 0.136E+01 0.239E-01 0.840E-03 0.178E-03 60 8 11 11 14 0 11 14 18 0.000 0.000 0
360	0 31 00:00:23.49 0.7580A -2.3755 -1.2207 1.1562 0.072 0.47 -0.971E+00 40. 0.136E+01 0.244E-01 0.846E-03 0.178E-03 60 8 10 11 14 0 11 14 18 0.000 0.000 0
361	0 32 00:00:24.25 0.7580A -2.3756 -1.2184 1.1525 0.072 0.47 -0.971E+00 40. 0.136E+01 0.250E-01 0.853E-03 0.177E-03 60 8 10 11 14 0 13 14 18 0.000 0.000 0
362	0 33 00:00:25.01 0.7580A -2.3756 -1.2151 1.1588 0.072 0.48 -0.971E+00 40. 0.136E+01 0.255E-01 0.860E-03 0.177E-03 60 7 9 11 14 0 13 14 18 0.000 0.000 0
363	0 34 00:00:25.77 0.7580A -2.3757 -1.2105 1.1879 0.072 0.48 -0.970E+00 40. 0.136E+01 0.261E-01 0.866E-03 0.177E-03 60 7 8 11 14 0 11 14 17 0.000 0.000 0
364	0 35 00:00:26.53 0.7580A -2.3758 -1.2043 1.2127 0.071 0.48 -0.970E+00 40. 0.136E+01 0.266E-01 0.874E-03 0.177E-03 60 7 7 11 14 0 11 14 17 0.000 0.000 0
365	0 36 00:00:27.28 0.7580A -2.3759 -1.1962 1.2272 0.071 0.49 -0.970E+00 40. 0.136E+01 0.272E-01 0.881E-03 0.177E-03 60 6 6 11 14 0 11 14 17 0.000 0.000 0
366	0 37 00:00:28.04 0.7580A -2.3759 -1.1860 1.2410 0.071 0.49 -0.970E+00 40. 0.136E+01 0.278E-01 0.889E-03 0.177E-03 60 6 5 11 14 0 12 14 17 0.000 0.000 0
367	0 38 00:00:28.80 0.7580A -2.3760 -1.1733 1.2634 0.071 0.49 -0.970E+00 40. 0.136E+01 0.284E-01 0.896E-03 0.178E-03 60 6 4 11 14 0 14 14 16 0.000 0.000 0
368	0 39 00:00:29.56 0.7580A -2.3761 -1.1618 1.2906 0.071 0.50 -0.969E+00 40. 0.136E+01 0.290E-01 0.905E-03 0.178E-03 60 5 3 11 13 0 14 14 16 0.000 0.000 0
369	0 40 00:00:30.32 0.7580A -2.3762 -1.1510 1.3236 0.071 0.50 -0.969E+00 40. 0.136E+01 0.296E-01 0.914E-03 0.179E-03 60 5 3 11 13 0 12 14 16 0.000 0.000 0
370	0 41 00:00:31.07 0.7580A -2.3762 -1.1369 1.3552 0.070 0.50 -0.969E+00 40. 0.136E+01 0.302E-01 0.922E-03 0.179E-03 60 5 2 11 13 0 13 14 16 0.000 0.000 0

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format