source: palm/trunk/TESTS/cases/rans_tkee/MONITORING/rans_tkee_rc @ 4101

Last change on this file since 4101 was 4101, checked in by gronemeier, 2 years ago

timestep.f90:

  • consider 2*Km within diffusion criterion as Km is considered twice within the diffusion of e,
  • in RANS mode, instead of considering each wind component individually use the wind speed of 3d wind vector in CFL criterion
  • do not limit the increase of dt based on its previous value in RANS mode

other:

  • remove dt_old
File size: 14.3 KB
Line 
1
2
3
41D run control output:
5------------------------------
6
7ITER.   HH:MM:SS    DT      UMAX   VMAX    U*   ALPHA   ENERG.
8-------------------------------------------------------------
9      0 00:00:00    0.01    3.00   0.10  0.100   45.0     4.45
10   3325 12:00:04   12.30    3.36   0.87  0.127   26.2     4.48
11
12 ******************************    --------------------------------------------
13 * PALM 6.0  Rev: 4100M       *    atmosphere - run with 1D - prerun
14 ******************************    --------------------------------------------
15
16 Date:               2019-07-17    Run:       rans_tkee__gfortran_default__4   
17 Time:                 16:59:43    Run-No.:   00
18 Run on host:        gfortran_d
19 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
20 ------------------------------------------------------------------------------
21
22 Numerical Schemes:
23 -----------------
24
25 --> Use the TKE-e turbulence closure (RANS mode).
26 --> Use the boussinesq approximation for the model equations.
27 --> Solve perturbation pressure via multigrid method (w-cycle)
28     number of grid levels:                    4
29     Gauss-Seidel red/black iterations:        2
30     fixed number of multigrid cycles:           2
31     perturbation pressure is calculated at every Runge-Kutta step
32 --> Momentum advection via Wicker-Skamarock-Scheme 5th order
33 --> Scalar advection via Wicker-Skamarock-Scheme 5th order
34 --> Loop optimization method: cache
35 --> Time differencing scheme: runge-kutta-3       
36 ------------------------------------------------------------------------------
37
38
39 Run time and time step information:
40 ----------------------------------
41
42 Timestep:             variable     maximum value: 20.000 s    CFL-factor: 0.90
43 Start time:              0.000 s
44 End time:               50.000 s
45
46
47 Computational grid and domain size:
48 ----------------------------------
49
50 Grid length:      dx =      10.000 m    dy =      10.000 m
51                dz(1) =      10.000 m
52                dz(2) =     999.000 m
53
54 Domain size:       x =    320.000 m     y =    320.000 m  z(u) =    965.000 m
55
56
57
58 Number of gridpoints (x,y,z):  (0:  31, 0:  31, 0:  97)
59 Subdomain size (x,y,z):        (    16,     16,     98)
60
61
62    Characteristic levels of the geo. wind component ug:
63
64       Height:         0.0     0.0  m
65       ug:            3.00    3.00  m/s
66       Gradient:    ------    0.00  1/100s
67       Gridpoint:        0       0
68
69    Characteristic levels of the geo. wind component vg:
70
71       Height:         0.0     0.0  m
72       vg:            0.00    0.00  m/s
73       Gradient:    ------    0.00  1/100s
74       Gridpoint:        0       0
75
76
77 Topography information:
78 ----------------------
79
80 Topography: read_from_file                         
81 Topography grid definition convention:
82 cell center (scalar grid points)
83
84
85
86 Boundary conditions:
87 -------------------
88
89                     p                    uv                                  pt
90
91 B. bound.: p(0)     = p(1)   | uv(0)     = -uv(1)                | pt(0)     = pt(1)           
92 T. bound.: p(nzt+1) = 0      | uv(nzt+1) = ug(nzt+1), vg(nzt+1)  | pt(nzt+1) = pt(nzt) + dpt/dz_
93
94                     e
95
96 B. bound.: e(0)     = e(1)                                                                     
97 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
98
99 Bottom surface fluxes are used in diffusion terms at k=1
100       Predefined constant heatflux:    0.001000 K m/s
101
102
103    Constant flux layer between bottom surface and first computational u,v-level:
104
105       z_mo =   5.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
106       Rif value range:     -20.00 <= rif <= 20.00
107
108
109 Lateral boundaries:
110       left/right:  dirichlet/radiation
111       north/south: cyclic             
112
113       use_cmax: T
114       pt damping layer width =     0.00 m, pt damping factor = 0.0000
115
116
117 Initial profiles:
118 ----------------
119
120    Characteristic levels of the initial temperature profile:
121
122       Height:           0.0     0.0  m
123       Temperature:   300.00  300.00  K
124       Gradient:      ------    0.00  K/100m
125       Gridpoint:          0       0
126
127
128 List output:
129 -----------
130
131    1D-Profiles:
132       Output every             9999999.90 s
133
134
135 Data output:
136 -----------
137
138
139    1D-Profiles:
140       Output format: netCDF 64bit offset                     
141
142
143       Profile: u, w"u", w*u*, wu, u*2, v, w"v", w*v*, wv, v*2, w, w*2,     
144
145              : e, e*, km, l,                                               
146       Output every                10.00 s
147       Time averaged over           0.00 s
148       Averaging input every        0.00 s
149
150    Time series:
151       Output format: netCDF 64bit offset                     
152
153       Output every                 1.00 s
154
155 ------------------------------------------------------------------------------
156
157
158 Physical quantities:
159 -------------------
160
161    Geograph. latitude  :   latitude  = 52.5 degr
162    Geograph. longitude :   longitude = 13.3 degr
163    Angular velocity    :   omega  = 0.729E-04 rad/s
164    Coriolis parameter  :   f      =  0.000116 1/s
165                            f*     =  0.000089 1/s
166
167    Day of the year at model start :   day_init      =     172
168    UTC time at model start        :   time_utc_init =     0.0 s
169
170    Gravity             :   g      =  9.8 m/s**2
171
172    Reference state used in buoyancy terms: initial_profile
173
174
175 Cloud physics quantities / methods:
176 ----------------------------------
177
178
179
180 LES / Turbulence quantities:
181 ---------------------------
182
183    initial TKE is prescribed as  0.10E+00 (m/s)**2
184    TKE is not allowed to fall below  0.10E-02 (m/s)**2
185    Mixing length is limited to 1.80 * z
186    RANS-mode constants: c_0 =   0.55000
187                         c_1 =   1.44000
188                         c_2 =   1.92000
189                         c_3 =   1.44000
190                         c_4 =   0.00000
191                         sigma_e    =   1.00000
192                         sigma_diss =   1.30000
193
194
195 Actions during the simulation:
196 -----------------------------
197
198
199
200 1D-Model parameters:
201 -------------------
202
203    Simulation time:                    43200.0 s
204    Run-controll output every:          43200.0 s
205    Vertical profile output every:      86400.0 s
206    Mixing length calculation:         blackadar           
207    Dissipation calculation:           prognostic         
208
209
210
211 *** dynamic module disabled
212
213 ------------------------------------------------------------------------------
214
215
216
217Run-control output:
218------------------
219
220RUN  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
221---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
222  0      0 00:00:00.00   2.6700A   3.3592   0.8728   0.0000  0.127  0.00   0.000E+00      0.  0.442E+01  0.442E-02  0.000E+00  0.000E+00   20  29  29    6  29  29    0  29  29    0.000    0.000       0
223  0      1 00:00:02.67   2.6600A   3.3803   1.5756   1.1290  0.151  0.00  -0.406E-02     10.  0.444E+01  0.437E-02  0.202E-03  0.798E-04   19   7  19    2  14  13    2  13   1    0.000    0.000       2
224  0      2 00:00:05.33   0.2370D   3.3829   1.4756   1.2399  0.147  0.00  -0.468E-02     10.  0.444E+01  0.432E-02  0.137E-03  0.457E-04   19   6  21    2  14  13    2  13   1    0.000    0.000       2
225  0      3 00:00:05.56   0.5360D   3.3841   1.5585   1.3626  0.145  0.00  -0.487E-02      0.  0.444E+01  0.422E-02  0.512E-04  0.118E-04   19   6  22    2  17   0    2  13   1    0.000    0.000       2
226  0      4 00:00:06.10   0.7840D   3.3844   1.5478   1.3603  0.145  0.00  -0.497E-02      0.  0.444E+01  0.408E-02  0.412E-04  0.145E-04   19   5  23    2  17   0    2  13   1    0.000    0.000       2
227  0      5 00:00:06.88   1.0200D   3.3843   1.5248   1.3393  0.144  0.00  -0.510E-02      0.  0.444E+01  0.401E-02  0.459E-04  0.160E-04   19   5  25    2  17   0    2  13   1    0.000    0.000       2
228  0      6 00:00:07.90   1.2500D   3.3842   1.5006   1.3115  0.143  0.00  -0.528E-02      0.  0.444E+01  0.397E-02  0.521E-04  0.121E-04   19   3  27    2  17   0    2  13   1    0.000    0.000       2
229  0      7 00:00:09.15   1.2700D   3.3842   1.4774   1.2804  0.141  0.00  -0.554E-02      0.  0.444E+01  0.393E-02  0.643E-04  0.175E-04   19   0   0    2  17   0    2  13   1    0.000    0.000       2
230  0      8 00:00:10.42   1.3200D   3.3840   1.4655   1.2613  0.140  0.00  -0.558E-02      0.  0.444E+01  0.387E-02  0.621E-04  0.164E-04   19   2   0    2  17   0    2  13   1    0.000    0.000       2
231  0      9 00:00:11.74   1.4300D   3.3835   1.4555   1.2678  0.139  0.00  -0.563E-02      0.  0.444E+01  0.380E-02  0.562E-04  0.155E-04   19   1   0    2  17   0    1  12   2    0.000    0.000       2
232  0     10 00:00:13.17   1.5700D   3.3829   1.4445   1.2852  0.138  0.00  -0.566E-02      0.  0.444E+01  0.376E-02  0.568E-04  0.170E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
233  0     11 00:00:14.74   1.7600D   3.3824   1.4336   1.2960  0.137  0.00  -0.573E-02      0.  0.444E+01  0.371E-02  0.587E-04  0.172E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
234  0     12 00:00:16.50   2.0100D   3.3820   1.4223   1.2996  0.136  0.00  -0.590E-02      0.  0.444E+01  0.367E-02  0.625E-04  0.183E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
235  0     13 00:00:18.51   1.9500D   3.3816   1.4105   1.2953  0.135  0.00  -0.602E-02      0.  0.444E+01  0.362E-02  0.674E-04  0.202E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
236  0     14 00:00:20.46   1.8700D   3.3814   1.4124   1.2963  0.134  0.00  -0.644E-02      0.  0.444E+01  0.359E-02  0.635E-04  0.194E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
237  0     15 00:00:22.33   1.8000D   3.3814   1.4174   1.2972  0.133  0.00  -0.626E-02      0.  0.444E+01  0.356E-02  0.594E-04  0.185E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
238  0     16 00:00:24.13   1.7400D   3.3814   1.4235   1.2974  0.133  0.00  -0.622E-02      0.  0.444E+01  0.353E-02  0.562E-04  0.174E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
239  0     17 00:00:25.87   1.6900D   3.3815   1.4301   1.2973  0.132  0.00  -0.630E-02      0.  0.444E+01  0.351E-02  0.536E-04  0.165E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
240  0     18 00:00:27.56   1.6400D   3.3815   1.4371   1.2970  0.131  0.00  -0.639E-02      0.  0.444E+01  0.349E-02  0.513E-04  0.159E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
241  0     19 00:00:29.20   1.6000D   3.3815   1.4448   1.2972  0.131  0.00  -0.659E-02      0.  0.444E+01  0.347E-02  0.491E-04  0.152E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
242  0     20 00:00:30.80   1.5500D   3.3816   1.4526   1.2975  0.130  0.00  -0.661E-02      0.  0.444E+01  0.345E-02  0.471E-04  0.147E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
243  0     21 00:00:32.35   1.4700D   3.3816   1.4611   1.2985  0.129  0.00  -0.654E-02      0.  0.444E+01  0.343E-02  0.451E-04  0.141E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
244  0     22 00:00:33.82   1.4000D   3.3817   1.4711   1.3011  0.129  0.00  -0.649E-02      0.  0.444E+01  0.342E-02  0.426E-04  0.134E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
245  0     23 00:00:35.22   1.3400D   3.3817   1.4808   1.3040  0.129  0.00  -0.646E-02      0.  0.444E+01  0.341E-02  0.403E-04  0.126E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
246  0     24 00:00:36.56   1.2800D   3.3818   1.4899   1.3070  0.128  0.00  -0.645E-02      0.  0.444E+01  0.340E-02  0.383E-04  0.120E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
247  0     25 00:00:37.84   1.2400D   3.3818   1.4988   1.3101  0.128  0.00  -0.647E-02      0.  0.444E+01  0.339E-02  0.363E-04  0.115E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
248  0     26 00:00:39.08   1.1900D   3.3819   1.5066   1.3128  0.127  0.00  -0.650E-02      0.  0.444E+01  0.338E-02  0.349E-04  0.110E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
249  0     27 00:00:40.27   1.1500D   3.3819   1.5146   1.3157  0.127  0.00  -0.654E-02      0.  0.444E+01  0.337E-02  0.333E-04  0.106E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
250  0     28 00:00:41.42   1.1100D   3.3820   1.5219   1.3184  0.127  0.00  -0.659E-02      0.  0.444E+01  0.336E-02  0.319E-04  0.102E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
251  0     29 00:00:42.53   1.0800D   3.3820   1.5289   1.3211  0.127  0.00  -0.663E-02      0.  0.444E+01  0.335E-02  0.306E-04  0.979E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
252  0     30 00:00:43.61   1.0500D   3.3821   1.5352   1.3234  0.126  0.00  -0.667E-02      0.  0.444E+01  0.334E-02  0.296E-04  0.949E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
253  0     31 00:00:44.66   1.0200D   3.3822   1.5412   1.3257  0.126  0.00  -0.667E-02      0.  0.444E+01  0.334E-02  0.286E-04  0.921E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
254  0     32 00:00:45.68   0.9930D   3.3822   1.5470   1.3278  0.126  0.00  -0.664E-02      0.  0.444E+01  0.333E-02  0.276E-04  0.893E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
255  0     33 00:00:46.68   0.9680D   3.3823   1.5525   1.3299  0.126  0.00  -0.661E-02      0.  0.444E+01  0.332E-02  0.267E-04  0.868E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
256  0     34 00:00:47.64   0.9450D   3.3823   1.5576   1.3318  0.126  0.00  -0.659E-02      0.  0.444E+01  0.332E-02  0.259E-04  0.844E-05   19   8  24    2  17   0    1  12   2    0.000    0.000       2
257  0     35 00:00:48.59   0.9240D   3.3824   1.5625   1.3335  0.125  0.00  -0.657E-02      0.  0.444E+01  0.331E-02  0.252E-04  0.823E-05   19   8  24    2  17   0    1  12   2    0.000    0.000       2
258  0     36 00:00:49.51   0.9040D   3.3825   1.5670   1.3351  0.125  0.00  -0.655E-02      0.  0.444E+01  0.331E-02  0.245E-04  0.803E-05   19   9  24    2  17   0    1  12   2    0.000    0.000       2
259  0     37 00:00:50.42   0.8850D   3.3826   1.5714   1.3366  0.125  0.00  -0.654E-02      0.  0.444E+01  0.330E-02  0.239E-04  0.785E-05   19   9  24    2  17   0    1  12   2    0.000    0.000       2
Note: See TracBrowser for help on using the repository browser.