source: palm/trunk/INSTALL/ex_lsm_clearsky_rc @ 2697

Last change on this file since 2697 was 2696, checked in by kanani, 7 years ago

Merge of branch palm4u into trunk

File size: 12.2 KB
RevLine 
[2542]1
2 ******************************    --------------------------------------------
[2696]3 * PALM 4.0  Rev: 2693        *    atmosphere - run without 1D - prerun
[2542]4 ******************************    --------------------------------------------
5
[2696]6 Date:                 13-12-17    Run:       ex_lsm_clearsky     
7 Time:                 18:28:57    Run-No.:   00
8 Run on host:        imukpalm4u
[2556]9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
[2542]10 ------------------------------------------------------------------------------
11
12 Numerical Schemes:
13 -----------------
14
[2696]15 --> Use the Moeng_Wyngaard turbulence closure (LES mode).
[2542]16 --> Use the boussinesq approximation for the model equations.
17 --> Solve perturbation pressure via FFT using temperton-algorithm 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
[2556]32 End time:              120.000 s
[2542]33
34
35 Computational grid and domain size:
36 ----------------------------------
37
38 Grid length:      dx =     50.000 m    dy =     50.000 m    dz =     50.000 m
39 Domain size:       x =   2000.000 m     y =   2000.000 m  z(u) =   2862.511 m
40
41 dz constant up to   1200.000 m (k=  25), then stretched by factor: 1.080
42 maximum dz not to be exceeded is dz_max = ********** m
43
44 Number of gridpoints (x,y,z):  (0:  39, 0:  39, 0:  41)
[2556]45 Subdomain size (x,y,z):        (    20,     20,     42)
[2542]46
47
48    Characteristic levels of the geo. wind component ug:
49
50       Height:         0.0     0.0  m
[2556]51       ug:            5.00    5.00  m/s
[2542]52       Gradient:    ------    0.00  1/100s
53       Gridpoint:        0       0
54
55    Characteristic levels of the geo. wind component vg:
56
57       Height:         0.0     0.0  m
58       vg:            0.00    0.00  m/s
59       Gradient:    ------    0.00  1/100s
60       Gridpoint:        0       0
61
62
63 Topography information:
64 ----------------------
65
66 Topography: flat                                   
67
68
69 Land surface model information:
70 ------------------------------
71
72    --> Soil bottom is closed (water content is conserved, default)
73    --> Land surface type  : crops, mixed farming
74    --> Soil porosity type : medium-fine
75
76    Initial soil temperature and moisture profile:
77
78       Height:              -0.01      -0.02      -0.05      -0.11      -0.22      -0.40      -0.86      -2.  m
79       Temperature:        293.00     293.00     293.00     293.00     293.00     293.00     293.00     293.  K
80       Moisture:             0.35       0.35       0.35       0.35       0.35       0.35       0.35       0.  m**3/m**3
81       Root fraction:        0.03       0.07       0.14       0.14       0.28       0.09       0.22       0. 
82       Grid point:              0          1          2          3          4          5          6
83
84
85 Radiation model information:
86 ----------------------------
87
88    --> Simple radiation scheme for clear sky is used (no clouds, default)
89
90    Albedo is set for land surface type: mixed farming, tall grassland
91    Timestep: dt_radiation =  60.00  s
92
93
94 Boundary conditions:
95 -------------------
96
97                     p                    uv                                  pt
98
99 B. bound.: p(0)     = p(1)   | uv(0)     = -uv(1)                | pt(0)     = from soil model 
100 T. bound.: p(nzt+1) = 0      | uv(nzt+1) = ug(nzt+1), vg(nzt+1)  | pt(nzt+1) = pt(nzt) + dpt/dz_
101
102                     e
103
104 B. bound.: e(0)     = e(1)                                                                     
105 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
106
107                     q
108
109 B. bound.: q(0)     = from soil model                                                           
110 T. bound.: q(nzt+1) = q(nzt) + dq/dz                                                           
111
112 Bottom surface fluxes are used in diffusion terms at k=1
113
114
115    Constant flux layer between bottom surface and first computational u,v-level:
116
117       z_mo =  25.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
118       Rif value range:     -20.00 <= rif <= 20.00
119       Predefined surface temperature
120       Predefined surface humidity
121
122
123 Lateral boundaries:
124       left/right:  cyclic             
125       north/south: cyclic             
126
127
128 Initial profiles:
129 ----------------
130
131    Characteristic levels of the initial temperature profile:
132
133       Height:           0.0     0.0   800.0  m
134       Temperature:   293.00  293.00  293.00  K
135       Gradient:      ------    0.00    1.00  K/100m
136       Gridpoint:          0       0      16
137
138    Characteristic levels of the initial humidity profile:
139
140       Height:           0.0       0.0  m
141       Humidity:     0.1E-01   0.1E-01  kg/kg
142       Gradient:    --------   0.0E+00  (kg/kg)/100m
143       Gridpoint:          0         0
144
145
146 List output:
147 -----------
148
149    1D-Profiles:
150       Output every             ******** s
151       Time averaged over         900.00 s
152       Averaging input every        0.00 s
153
154
155 Data output:
156 -----------
157
158       Time averaged over         900.00 s
159       Averaging input every        0.00 s
160
161    1D-Profiles:
162       Output format: netCDF 64bit offset                     
163
164
165       Profile: u, v, pt, km, kh, l, w, e, e*, p, w"pt", w*pt*, wpt,         
166
167              : w*e*, pt*2, u*2, v*2, w*2, t_soil,                           
168       Output every               900.00 s
169       Time averaged over         900.00 s
170       Averaging input every        0.00 s
171
172    2D-Arrays:
173       Output format: netCDF 64bit offset                     
174
175
176       XY-cross-section  Arrays:  u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l
177       Output every              3600.00 s 
178       Cross sections at k = /0/
179       scalar-coordinates:   /0.0/ m
180
181
182       Number of output time levels allowed: unlimited
183
184
185    2D-Arrays(time-averaged):
186       Output format: netCDF 64bit offset                     
187
188
189       XY-cross-section  Arrays:  u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l
190       Output every              3600.00 s 
191       Time averaged over        1800.00 s
192       Averaging input every        0.00 s
193       Cross sections at k = /0/
194       scalar-coordinates:   /0.0/ m
195
196
197       Number of output time levels allowed: unlimited
198
199
200    3D-Arrays:
201       Output format: netCDF 64bit offset                     
202
203
204       Arrays:  pt, vpt, q, u, v, w, t_soil, m_soil,                                 
205       Output every              3600.00 s 
206       Upper output limit at     2862.51 m  (GP   41)
207
208
209       Number of output time levels allowed: unlimited
210
211
212    3D-Arrays(time-averaged):
213       Output format: netCDF 64bit offset                     
214
215
216       Arrays:  t_soil, m_soil,                                                     
217       Output every              3600.00 s 
218       Time averaged over        1800.00 s
219       Averaging input every        0.00 s
220       Upper output limit at     2862.51 m  (GP   41)
221
222
223       Number of output time levels allowed: unlimited
224
225
226    Time series:
227       Output format: netCDF 64bit offset                     
228
229       Output every                 0.00 s
230
231 ------------------------------------------------------------------------------
232
233
234 Physical quantities:
235 -------------------
236
[2608]237    Geograph. latitude  :   latitude  = 55.0 degr
238    Geograph. longitude :   longitude =  0.0 degr
[2542]239    Angular velocity    :   omega  = 0.729E-04 rad/s
240    Coriolis parameter  :   f      =  0.000119 1/s
241                            f*     =  0.000084 1/s
[2608]242    Day of the year at model start :   day_init = 172
243    UTC time at model start        :   time_utc_init = 43200.0 s
[2542]244
245    Gravity             :   g      =  9.8 m/s**2
246
[2543]247    Reference state used in buoyancy terms: initial_profile
[2542]248
249
250 Cloud physics quantities / methods:
251 ----------------------------------
252
253    Humidity is considered, bu no condensation
254
255
256 LES / Turbulence quantities:
257 ---------------------------
258
259    Mixing length is limited to 1.80 * z
260
261
262 Actions during the simulation:
263 -----------------------------
264
[2556]265    Disturbance impulse (u,v) every :   ****** s
[2542]266    Disturbance amplitude           :     0.25 m/s
267    Lower disturbance level         :   125.00 m (GP    3)
268    Upper disturbance level         :   625.00 m (GP   13)
269    Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2
270    Random number generator used    : random-parallel
271
272
273
274 *** no user-defined variables found
275
276 ------------------------------------------------------------------------------
277
278
279
280Run-control output:
281------------------
282
283RUN  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
284---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
[2696]285  0      0 00:00:00.00   8.6700A   5.1898D -0.2288D  0.1182  0.100  0.00   0.000E+00    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
286  0      1 00:00:08.67   8.6900A   5.1806   0.2159   0.1162  0.004  0.26  -0.882E+01    100.  0.119E+02  0.867E-03  0.953E-05  0.439E-15    9   9  24    4   5  26    5  26   8    0.000    0.000       0
287  0      2 00:00:17.35   8.6900A   5.1762   0.1915   0.1151  0.000  0.04  -0.174E-01    100.  0.119E+02  0.793E-03  0.276E-05  0.403E-15    5  16  32    4   5  27    5  26   9    0.000    0.000       0
288  0      3 00:00:26.04   8.7000A   5.1723   0.1786   0.1113  0.005  0.27  -0.365E+01    100.  0.119E+02  0.750E-03  0.314E-05  0.376E-15    5  16  33    4   5  28    5  26  10    0.000    0.000       0
289  0      4 00:00:34.75   8.7000A   5.1732  -0.1675   0.1061  0.004  0.43  -0.144E+01    250.  0.119E+02  0.721E-03  0.272E-05  0.356E-15    5  16  34    8  29  32    5  26  11    0.000    0.000       0
290  0      5 00:00:43.45   8.7000A   5.1719  -0.1671   0.1028  0.005  0.43  -0.170E+01    250.  0.119E+02  0.698E-03  0.292E-05  0.335E-15    5  16  35    8  29  33    5  26  11    0.000    0.000       0
291  0      6 00:00:52.15   8.7000A   5.1700  -0.1621   0.1039  0.007  0.45  -0.194E+01    250.  0.119E+02  0.680E-03  0.303E-05  0.319E-15    5  16  36    8  29  34    5  26  12    0.000    0.000       0
292  0      7 00:01:00.85   8.7100A   5.1674  -0.1565   0.1044  0.009  0.65  -0.207E+01    550.  0.119E+02  0.664E-03  0.315E-05  0.303E-15    5  16  37    8  29  35    5  26  13    0.000    0.000       0
293  0      8 00:01:09.56   8.7100A   5.1641   0.1515   0.1041  0.011  0.70  -0.181E+01    550.  0.119E+02  0.651E-03  0.325E-05  0.290E-15    5  16  38   11   3  22    5  26  14    0.000    0.000       0
294  0      9 00:01:18.27   8.7200A   5.1634  -0.1506   0.1035  0.013  0.69  -0.138E+01    550.  0.119E+02  0.639E-03  0.332E-05  0.277E-15    5  16  38    8  29  36    5  26  15    0.000    0.000       0
295  0     10 00:01:26.99   8.7200A   5.1633  -0.1516   0.1026  0.016  0.73  -0.141E+01    550.  0.119E+02  0.629E-03  0.361E-05  0.266E-15    5  16  39    8  29  37    5  26  16    0.000    0.000       0
296  0     11 00:01:35.71   8.7200A   5.1628  -0.1516   0.1014  0.018  0.76  -0.136E+01    550.  0.119E+02  0.619E-03  0.379E-05  0.256E-15    5  16   0    8  29  38    5  26  17    0.000    0.000       0
297  0     12 00:01:44.43   8.7200A   5.1621  -0.1506   0.1000  0.020  0.79  -0.132E+01    550.  0.119E+02  0.610E-03  0.403E-05  0.246E-15    5  16   1    8  29  39    5  26  18    0.000    0.000       0
298  0     13 00:01:53.15   8.7200A   5.1610  -0.1488   0.0983  0.023  0.81  -0.129E+01    550.  0.119E+02  0.602E-03  0.430E-05  0.237E-15    5  16   2    8  29   0    5  26  19    0.000    0.000       0
299  0     14 00:02:01.87   8.7200A   5.1597  -0.1462   0.0964  0.025  0.83  -0.126E+01    550.  0.119E+02  0.593E-03  0.459E-05  0.227E-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.