source: palm/trunk/TESTS/cases/oceanml/MONITORING/oceanml_rc @ 4405

Last change on this file since 4405 was 4211, checked in by raasch, 5 years ago

ocean mixed layer test case modified; can now be used as a basis to reproduce the Noh et.al. (2004) results

File size: 7.4 KB
Line 
1
2 ******************************    --------------------------------------------
3 * PALM 6.0  Rev: 4071        *    ocean - run without 1D - prerun
4 ******************************    --------------------------------------------
5
6 Date:               2019-09-02    Run:       oceanml                           
7 Time:                 15:36:39    Run-No.:   00
8 Run on host:           default
9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
10 ------------------------------------------------------------------------------
11
12 Numerical Schemes:
13 -----------------
14
15 --> Use the Moeng_Wyngaard turbulence closure (LES mode).
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 ------------------------------------------------------------------------------
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 =       1.250 m    dy =       1.250 m
38                dz(1) =       1.250 m
39
40 Domain size:       x =     80.000 m     y =     80.000 m  z(u) =    -80.000 m
41
42
43
44 Number of gridpoints (x,y,z):  (0:  63, 0:  63, 0:  65)
45 Subdomain size (x,y,z):        (    32,     32,     66)
46
47
48    Characteristic levels of the geo. wind component ug:
49
50       Height:         0.0     0.0  m
51       ug:            0.00    0.00  m/s
52       Gradient:    ------    0.00  1/100s
53       Gridpoint:        0      65
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      65
61
62
63 Topography information:
64 ----------------------
65
66 Topography: flat                                   
67
68
69 Boundary conditions:
70 -------------------
71
72                     p                    uv                                  pt
73
74 B. bound.: p(0)     = p(1)   | uv(0)     = uv(1)                 | pt(0)     = pt(1)           
75 T. bound.: p(nzt+1) = p(nzt) | uv(nzt+1) = uv(nzt)               | pt(nzt+1) = pt(nzt)         
76
77                     e
78
79 B. bound.: e(0)     = e(1)                                                                     
80 T. bound.: e(nzt+1) = e(nzt) = e(nzt-1)                                                         
81
82                     sa
83
84 B. bound.: sa(0)    = sa(1)                                                                     
85 T. bound.: sa(nzt+1) = sa(nzt)                                                                 
86
87 Top surface fluxes are used in diffusion terms at k=nzt
88       Predefined constant momentumflux:  u: -0.000100 m**2/s**2
89                                          v:  0.000000 m**2/s**2
90       Predefined constant heatflux:    0.000102 K m/s
91       Predefined constant salinityflux:    0.000000 psu m/s
92
93
94 Lateral boundaries:
95       left/right:  cyclic             
96       north/south: cyclic             
97
98
99 Initial profiles:
100 ----------------
101
102    Characteristic levels of the initial temperature profile:
103
104       Height:           0.0     0.0  m
105       Temperature:   300.00  300.00  K
106       Gradient:      ------    0.00  K/100m
107       Gridpoint:          0      65
108
109    Characteristic levels of the initial salinity profile:
110
111       Height:        0.0     0.0  m
112       Salinity:    35.00   35.00  psu
113       Gradient:   ------    0.00  psu/100m
114       Gridpoint:       0      65
115
116
117 List output:
118 -----------
119
120    1D-Profiles:
121       Output every             9999999.90 s
122
123
124 Data output:
125 -----------
126
127
128    Time series:
129       Output format: netCDF 64bit offset                     
130
131       Output every             ******** s
132
133 ------------------------------------------------------------------------------
134
135
136 Physical quantities:
137 -------------------
138
139    Geograph. latitude  :   latitude  = 55.0 degr
140    Geograph. longitude :   longitude =  0.0 degr
141    Angular velocity    :   omega  = 0.729E-04 rad/s
142    Coriolis parameter  :   f      =  0.000119 1/s
143                            f*     =  0.000084 1/s
144
145    Day of the year at model start :   day_init      =     172
146    UTC time at model start        :   time_utc_init = 43200.0 s
147
148    Gravity             :   g      =  9.8 m/s**2
149
150    Reference state used in buoyancy terms: initial_profile
151
152
153 Cloud physics quantities / methods:
154 ----------------------------------
155
156
157
158 LES / Turbulence quantities:
159 ---------------------------
160
161    Mixing length is limited to 1.80 * z
162
163
164 Actions during the simulation:
165 -----------------------------
166
167    Disturbance impulse (u,v) every :   ****** s
168    Disturbance amplitude           :     0.00 m/s
169    Lower disturbance level         :   -28.12 m (GP   42)
170    Upper disturbance level         :    -4.38 m (GP   61)
171    Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2
172    Random number generator used    : random-parallel
173
174
175
176 *** dynamic module disabled
177
178
179
180 Ocean settings:
181 ------------------------------------------
182
183    --> Craik-Leibovich vortex force and Stokes drift switched on
184        waveheight:  1.0 m   wavelength:   40.0 m
185    --> surface heat flux is switched off after    900.0 s
186 ------------------------------------------------------------------------------
187
188
189
190Run-control output:
191------------------
192
193RUN  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
194---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
195  0      0 00:00:00.00  20.0000A   0.0002D  0.0002D  0.0001  0.000  0.00   0.000E+00      0.  0.122E-08  0.122E-08  0.270E-04  0.387E-14   53   1   6   58  26  28   47   9   8    0.000    0.000       0
196  0      1 00:00:20.00  20.0000A   0.0016  -0.0002   0.0001  0.000  0.00   0.000E+00      0.  0.401E-07  0.122E-08  0.380E-06  0.423E-16   65  55  48   60  21  57   47   9   8    0.000    0.000       0
197  0      2 00:00:40.00  20.0000A   0.0032  -0.0002   0.0001  0.000  0.00   0.000E+00      0.  0.156E-06  0.122E-08  0.380E-06  0.422E-16   65  56  48   61  22  63   47   9   8    0.000    0.000       0
198  0      3 00:01:00.00  20.0000A   0.0048  -0.0003   0.0001  0.000  0.00   0.000E+00      0.  0.347E-06  0.122E-08  0.380E-06  0.421E-16   65  56  48   64  54  31   47   9   8    0.000    0.000       0
199  0      4 00:01:20.00  20.0000A   0.0063  -0.0004   0.0001  0.000  0.00   0.000E+00      0.  0.611E-06  0.122E-08  0.380E-06  0.419E-16   65  55  49   64  54  31   47   9   8    0.000    0.000       0
200  0      5 00:01:40.00  20.0000A   0.0079  -0.0005   0.0001  0.000  0.00   0.000E+00      0.  0.945E-06  0.122E-08  0.381E-06  0.417E-16   64  55  50   65  54  31   58  16  15    0.000    0.000       0
201  0      6 00:02:00.00  20.0000A   0.0094  -0.0007   0.0001  0.000  0.00   0.000E+00      0.  0.135E-05  0.122E-08  0.381E-06  0.415E-16   65  55  50   64  63  26   58  16  15    0.000    0.000       0
Note: See TracBrowser for help on using the repository browser.