source: palm/trunk/INSTALL/ex_lsm_clearsky_rc @ 2630

Last change on this file since 2630 was 2608, checked in by schwenkel, 7 years ago

Inital revision of diagnostic_quantities_mod allows unified calculation of magnus equation and saturion mixing ratio

File size: 12.1 KB
Line 
1
2 ******************************    --------------------------------------------
3 * PALM 4.0  Rev: 2603:26     *    atmosphere - run without 1D - prerun
4 ******************************    --------------------------------------------
5
6 Date:                 13-11-17    Run:       ex_lsm_clearsky     
7 Time:                 14:31:43    Run-No.:   00
8 Run on host:             lcmuk
9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
10 ------------------------------------------------------------------------------
11
12 Numerical Schemes:
13 -----------------
14
15 --> Use the boussinesq approximation for the model equations.
16 --> Solve perturbation pressure via FFT using temperton-algorithm routines
17     perturbation pressure is calculated at every Runge-Kutta step
18 --> Momentum advection via Wicker-Skamarock-Scheme 5th order
19 --> Scalar advection via Wicker-Skamarock-Scheme 5th order
20 --> Loop optimization method: cache
21 --> Time differencing scheme: runge-kutta-3       
22 --> Additional prognostic equation for the specific humidity
23 ------------------------------------------------------------------------------
24
25
26 Run time and time step information:
27 ----------------------------------
28
29 Timestep:             variable     maximum value: 20.000 s    CFL-factor: 0.90
30 Start time:              0.000 s
31 End time:              120.000 s
32
33
34 Computational grid and domain size:
35 ----------------------------------
36
37 Grid length:      dx =     50.000 m    dy =     50.000 m    dz =     50.000 m
38 Domain size:       x =   2000.000 m     y =   2000.000 m  z(u) =   2862.511 m
39
40 dz constant up to   1200.000 m (k=  25), then stretched by factor: 1.080
41 maximum dz not to be exceeded is dz_max = ********** m
42
43 Number of gridpoints (x,y,z):  (0:  39, 0:  39, 0:  41)
44 Subdomain size (x,y,z):        (    20,     20,     42)
45
46
47    Characteristic levels of the geo. wind component ug:
48
49       Height:         0.0     0.0  m
50       ug:            5.00    5.00  m/s
51       Gradient:    ------    0.00  1/100s
52       Gridpoint:        0       0
53
54    Characteristic levels of the geo. wind component vg:
55
56       Height:         0.0     0.0  m
57       vg:            0.00    0.00  m/s
58       Gradient:    ------    0.00  1/100s
59       Gridpoint:        0       0
60
61
62 Topography information:
63 ----------------------
64
65 Topography: flat                                   
66
67
68 Land surface model information:
69 ------------------------------
70
71    --> Soil bottom is closed (water content is conserved, default)
72    --> Land surface type  : crops, mixed farming
73    --> Soil porosity type : medium-fine
74
75    Initial soil temperature and moisture profile:
76
77       Height:              -0.01      -0.02      -0.05      -0.11      -0.22      -0.40      -0.86      -2.  m
78       Temperature:        293.00     293.00     293.00     293.00     293.00     293.00     293.00     293.  K
79       Moisture:             0.35       0.35       0.35       0.35       0.35       0.35       0.35       0.  m**3/m**3
80       Root fraction:        0.03       0.07       0.14       0.14       0.28       0.09       0.22       0. 
81       Grid point:              0          1          2          3          4          5          6
82
83
84 Radiation model information:
85 ----------------------------
86
87    --> Simple radiation scheme for clear sky is used (no clouds, default)
88
89    Albedo is set for land surface type: mixed farming, tall grassland
90    Timestep: dt_radiation =  60.00  s
91
92
93 Boundary conditions:
94 -------------------
95
96                     p                    uv                                  pt
97
98 B. bound.: p(0)     = p(1)   | uv(0)     = -uv(1)                | pt(0)     = from soil model 
99 T. bound.: p(nzt+1) = 0      | uv(nzt+1) = ug(nzt+1), vg(nzt+1)  | pt(nzt+1) = pt(nzt) + dpt/dz_
100
101                     e
102
103 B. bound.: e(0)     = e(1)                                                                     
104 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
105
106                     q
107
108 B. bound.: q(0)     = from soil model                                                           
109 T. bound.: q(nzt+1) = q(nzt) + dq/dz                                                           
110
111 Bottom surface fluxes are used in diffusion terms at k=1
112
113
114    Constant flux layer between bottom surface and first computational u,v-level:
115
116       z_mo =  25.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
117       Rif value range:     -20.00 <= rif <= 20.00
118       Predefined surface temperature
119       Predefined surface humidity
120
121
122 Lateral boundaries:
123       left/right:  cyclic             
124       north/south: cyclic             
125
126
127 Initial profiles:
128 ----------------
129
130    Characteristic levels of the initial temperature profile:
131
132       Height:           0.0     0.0   800.0  m
133       Temperature:   293.00  293.00  293.00  K
134       Gradient:      ------    0.00    1.00  K/100m
135       Gridpoint:          0       0      16
136
137    Characteristic levels of the initial humidity profile:
138
139       Height:           0.0       0.0  m
140       Humidity:     0.1E-01   0.1E-01  kg/kg
141       Gradient:    --------   0.0E+00  (kg/kg)/100m
142       Gridpoint:          0         0
143
144
145 List output:
146 -----------
147
148    1D-Profiles:
149       Output every             ******** s
150       Time averaged over         900.00 s
151       Averaging input every        0.00 s
152
153
154 Data output:
155 -----------
156
157       Time averaged over         900.00 s
158       Averaging input every        0.00 s
159
160    1D-Profiles:
161       Output format: netCDF 64bit offset                     
162
163
164       Profile: u, v, pt, km, kh, l, w, e, e*, p, w"pt", w*pt*, wpt,         
165
166              : w*e*, pt*2, u*2, v*2, w*2, t_soil,                           
167       Output every               900.00 s
168       Time averaged over         900.00 s
169       Averaging input every        0.00 s
170
171    2D-Arrays:
172       Output format: netCDF 64bit offset                     
173
174
175       XY-cross-section  Arrays:  u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l
176       Output every              3600.00 s 
177       Cross sections at k = /0/
178       scalar-coordinates:   /0.0/ m
179
180
181       Number of output time levels allowed: unlimited
182
183
184    2D-Arrays(time-averaged):
185       Output format: netCDF 64bit offset                     
186
187
188       XY-cross-section  Arrays:  u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l
189       Output every              3600.00 s 
190       Time averaged over        1800.00 s
191       Averaging input every        0.00 s
192       Cross sections at k = /0/
193       scalar-coordinates:   /0.0/ m
194
195
196       Number of output time levels allowed: unlimited
197
198
199    3D-Arrays:
200       Output format: netCDF 64bit offset                     
201
202
203       Arrays:  pt, vpt, q, u, v, w, t_soil, m_soil,                                 
204       Output every              3600.00 s 
205       Upper output limit at     2862.51 m  (GP   41)
206
207
208       Number of output time levels allowed: unlimited
209
210
211    3D-Arrays(time-averaged):
212       Output format: netCDF 64bit offset                     
213
214
215       Arrays:  t_soil, m_soil,                                                     
216       Output every              3600.00 s 
217       Time averaged over        1800.00 s
218       Averaging input every        0.00 s
219       Upper output limit at     2862.51 m  (GP   41)
220
221
222       Number of output time levels allowed: unlimited
223
224
225    Time series:
226       Output format: netCDF 64bit offset                     
227
228       Output every                 0.00 s
229
230 ------------------------------------------------------------------------------
231
232
233 Physical quantities:
234 -------------------
235
236    Geograph. latitude  :   latitude  = 55.0 degr
237    Geograph. longitude :   longitude =  0.0 degr
238    Angular velocity    :   omega  = 0.729E-04 rad/s
239    Coriolis parameter  :   f      =  0.000119 1/s
240                            f*     =  0.000084 1/s
241    Day of the year at model start :   day_init = 172
242    UTC time at model start        :   time_utc_init = 43200.0 s
243
244    Gravity             :   g      =  9.8 m/s**2
245
246    Reference state used in buoyancy terms: initial_profile
247
248
249 Cloud physics quantities / methods:
250 ----------------------------------
251
252    Humidity is considered, bu no condensation
253
254
255 LES / Turbulence quantities:
256 ---------------------------
257
258    Mixing length is limited to 1.80 * z
259
260
261 Actions during the simulation:
262 -----------------------------
263
264    Disturbance impulse (u,v) every :   ****** s
265    Disturbance amplitude           :     0.25 m/s
266    Lower disturbance level         :   125.00 m (GP    3)
267    Upper disturbance level         :   625.00 m (GP   13)
268    Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2
269    Random number generator used    : random-parallel
270
271
272
273 *** no user-defined variables found
274
275 ------------------------------------------------------------------------------
276
277
278
279Run-control output:
280------------------
281
282RUN  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
283---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
284  0      0 00:00:00.00   8.6700A   5.1898D -0.2288D  0.1182  0.100  0.00   0.200E-11    800.  0.119E+02  0.103E-02  0.611E-03  0.888E-13    9   9  23   11   7  19    5  26   7    0.000    0.000       0
285  0      1 00:00:08.67   8.6900A   5.1806   0.2159   0.1162  0.001  0.00   0.000E+00    550.  0.119E+02  0.867E-03  0.287E-05  0.440E-15    9   9  24    4   5  26    5  26   8    0.000    0.000       0
286  0      2 00:00:17.35   8.6900A   5.1762   0.1915   0.1151  0.021  1.06  -0.495E+02    100.  0.119E+02  0.793E-03  0.744E-05  0.407E-15    5  16  32    4   5  27    5  26   9    0.000    0.000       0
287  0      3 00:00:26.04   8.7000A   5.1723   0.1786   0.1113  0.021  0.70  -0.516E+01    100.  0.119E+02  0.751E-03  0.247E-05  0.385E-15    5  16  33    4   5  28    5  26  10    0.000    0.000       0
288  0      4 00:00:34.75   8.7000A   5.1732  -0.1675   0.1061  0.027  0.69  -0.360E+01    100.  0.119E+02  0.721E-03  0.396E-05  0.358E-15    5  16  34    8  29  32    5  26  11    0.000    0.000       0
289  0      5 00:00:43.45   8.7000A   5.1719  -0.1671   0.1028  0.034  0.93  -0.287E+01    250.  0.119E+02  0.698E-03  0.471E-05  0.337E-15    5  16  35    8  29  33    5  26  11    0.000    0.000       0
290  0      6 00:00:52.15   8.7000A   5.1701  -0.1621   0.1039  0.039  0.93  -0.242E+01    250.  0.119E+02  0.680E-03  0.548E-05  0.321E-15    5  16  36    8  29  34    5  26  12    0.000    0.000       0
291  0      7 00:01:00.85   8.7100A   5.1675  -0.1565   0.1043  0.045  1.21  -0.211E+01    550.  0.119E+02  0.665E-03  0.627E-05  0.306E-15    5  16  37    8  29  35    5  26  13    0.000    0.000       0
292  0      8 00:01:09.56   8.7100A   5.1643   0.1515   0.1040  0.051  1.21  -0.189E+01    550.  0.119E+02  0.652E-03  0.709E-05  0.294E-15    5  16  38   11   3  22    5  26  14    0.000    0.000       0
293  0      9 00:01:18.27   8.7100A   5.1636  -0.1506   0.1034  0.056  1.20  -0.166E+01    550.  0.119E+02  0.640E-03  0.785E-05  0.282E-15    5  16  38    8  29  36    5  26  15    0.000    0.000       0
294  0     10 00:01:26.98   8.7100A   5.1635  -0.1516   0.1024  0.061  1.20  -0.153E+01    550.  0.119E+02  0.630E-03  0.868E-05  0.269E-15    5  16  39    8  29  37    5  26  16    0.000    0.000       0
295  0     11 00:01:35.69   8.7200A   5.1632  -0.1516   0.1013  0.067  1.20  -0.143E+01    550.  0.119E+02  0.620E-03  0.949E-05  0.260E-15    5  16   0    8  29  38    5  26  17    0.000    0.000       0
296  0     12 00:01:44.41   8.7200A   5.1624  -0.1506   0.0999  0.072  1.21  -0.135E+01    550.  0.119E+02  0.612E-03  0.103E-04  0.249E-15    5  16   1    8  29  39    5  26  18    0.000    0.000       0
297  0     13 00:01:53.13   8.7200A   5.1614  -0.1488   0.0983  0.077  1.22  -0.128E+01    550.  0.119E+02  0.604E-03  0.111E-04  0.239E-15    5  16   2    8  29   0    5  26  19    0.000    0.000       0
298  0     14 00:02:01.85   8.7200A   5.1600  -0.1462   0.0966  0.081  1.22  -0.122E+01    550.  0.119E+02  0.596E-03  0.119E-04  0.230E-15    5  16   3    8  29   1    5  26  20    0.000    0.000       0
Note: See TracBrowser for help on using the repository browser.