Home

Context Navigation

rans_tkee_rc @ 4102

Last change on this file since 4102 was 4102, checked in by suehring, 5 years ago

Bugfix, set Neumann boundary conditions for the subgrid TKE at vertical walls instead of implicit Dirichlet conditions that always act as a sink term for the subgrid TKE. Therefore, add new data structure for vertical surfaces and revise the setting of the boundary grid point index space. Moreover, slightly revise setting of boundary conditions at upward- and downward facing surfaces. Finally, setting of boundary conditions for subgrid TKE and dissipation (in RANS mode) is now modularized. Update test case results.

File size: 11.9 KB

Line
1
2
3
4	1D run control output:
5	------------------------------
6
7	ITER. 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: 3915M * atmosphere - run with 1D - prerun
14	****************************** --------------------------------------------
15
16	Date: 2019-07-17 Run: rans_tkee__gfortran_default__4
17	Time: 17:23:23 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", wu, wu, u2, v, w"v", wv, wv, v2, 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
217	Run-control output:
218	------------------
219
220	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
221	---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
222	0 0 00:00:00.00 0.0200A 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:00.02 0.0400A 3.3802 1.7430 -1.2281 0.153 0.00 -0.390E-02 10. 0.444E+01 0.460E-02 0.155E-03 0.701E-04 19 7 19 2 14 13 2 26 29 0.000 0.000 2
224	0 2 00:00:00.06 0.0800A 3.3832 1.7620 -1.3655 0.151 0.20 -0.421E-02 10. 0.444E+01 0.493E-02 0.874E-04 0.359E-04 19 6 21 2 14 13 2 26 29 0.000 0.000 2
225	0 3 00:00:00.14 0.1600A 3.3844 1.7515 1.3704 0.150 0.00 -0.428E-02 0. 0.444E+01 0.496E-02 0.977E-04 0.323E-04 19 5 22 2 14 13 2 13 1 0.000 0.000 2
226	0 4 00:00:00.30 0.3200A 3.3845 1.7361 1.3743 0.150 0.00 -0.429E-02 0. 0.444E+01 0.456E-02 0.521E-04 0.687E-05 19 5 24 2 14 13 2 13 1 0.000 0.000 2
227	0 5 00:00:00.62 0.6400A 3.3838 1.7102 1.3741 0.149 0.00 -0.432E-02 0. 0.444E+01 0.437E-02 0.385E-04 0.278E-04 19 4 26 2 14 13 2 13 1 0.000 0.000 2
228	0 6 00:00:01.26 1.2800A 3.3833 1.6639 1.3660 0.148 0.00 -0.437E-02 0. 0.444E+01 0.430E-02 0.704E-04 0.332E-04 19 1 0 2 14 13 2 13 1 0.000 0.000 2
229	0 7 00:00:02.54 2.5600A 3.3827 1.5852 1.3342 0.147 0.39 -0.449E-02 30. 0.444E+01 0.420E-02 0.618E-04 0.165E-04 19 1 0 2 14 13 2 13 1 0.000 0.000 2
230	0 8 00:00:05.09 2.6600A 3.3826 1.4669 1.2489 0.145 0.00 -0.484E-02 20. 0.444E+01 0.421E-02 0.144E-03 0.499E-04 20 3 31 2 14 13 2 13 1 0.000 0.000 2
231	0 9 00:00:07.75 2.6200D 3.3825 1.4248 1.2091 0.142 0.00 -0.540E-02 0. 0.444E+01 0.399E-02 0.137E-03 0.424E-04 20 3 0 2 17 0 2 13 1 0.000 0.000 2
232	0 10 00:00:10.37 2.6600A 3.3822 1.4058 1.2176 0.140 0.00 -0.536E-02 0. 0.444E+01 0.380E-02 0.971E-04 0.324E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2
233	0 11 00:00:13.03 2.6600A 3.3818 1.3920 1.2409 0.138 0.00 -0.554E-02 0. 0.444E+01 0.372E-02 0.912E-04 0.361E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2
234	0 12 00:00:15.69 2.6600A 3.3818 1.3841 1.2506 0.137 0.00 -0.573E-02 0. 0.444E+01 0.366E-02 0.867E-04 0.300E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2
235	0 13 00:00:18.35 2.6600A 3.3820 1.3790 1.2520 0.135 0.00 -0.581E-02 0. 0.444E+01 0.362E-02 0.864E-04 0.282E-04 19 6 23 2 17 0 1 12 2 0.000 0.000 2
236	0 14 00:00:21.01 2.6600A 3.3822 1.3757 1.2484 0.134 0.00 -0.617E-02 0. 0.444E+01 0.358E-02 0.862E-04 0.277E-04 19 6 23 2 17 0 1 12 2 0.000 0.000 2
237	0 15 00:00:23.67 2.6600A 3.3824 1.3733 1.2423 0.133 0.00 -0.644E-02 0. 0.444E+01 0.353E-02 0.800E-04 0.249E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2
238	0 16 00:00:26.33 2.6600A 3.3825 1.3716 1.2356 0.132 0.00 -0.633E-02 0. 0.444E+01 0.348E-02 0.755E-04 0.250E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2
239	0 17 00:00:29.00 2.6600A 3.3827 1.3703 1.2296 0.131 0.00 -0.648E-02 0. 0.444E+01 0.345E-02 0.734E-04 0.242E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2
240	0 18 00:00:31.66 2.6600A 3.3829 1.3692 1.2245 0.130 0.00 -0.650E-02 0. 0.444E+01 0.342E-02 0.719E-04 0.232E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2
241	0 19 00:00:34.32 2.6600A 3.3830 1.3683 1.2206 0.129 0.00 -0.661E-02 0. 0.444E+01 0.339E-02 0.709E-04 0.228E-04 19 7 24 2 17 0 1 12 2 0.000 0.000 2
242	0 20 00:00:36.98 2.6600A 3.3832 1.3675 1.2176 0.128 0.00 -0.648E-02 0. 0.444E+01 0.336E-02 0.692E-04 0.224E-04 19 7 24 2 17 0 1 12 2 0.000 0.000 2
243	0 21 00:00:39.64 2.6600A 3.3834 1.3668 1.2153 0.128 0.00 -0.649E-02 0. 0.444E+01 0.334E-02 0.674E-04 0.221E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2
244	0 22 00:00:42.29 2.6600A 3.3836 1.3662 1.2136 0.127 0.00 -0.658E-02 0. 0.444E+01 0.331E-02 0.661E-04 0.218E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2
245	0 23 00:00:44.95 2.6600A 3.3837 1.3657 1.2122 0.127 0.00 -0.664E-02 0. 0.444E+01 0.329E-02 0.651E-04 0.215E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2
246	0 24 00:00:47.61 2.6600A 3.3839 1.3652 1.2111 0.126 0.00 -0.661E-02 0. 0.444E+01 0.327E-02 0.642E-04 0.212E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2
247	0 25 00:00:50.27 2.6600A 3.3841 1.3647 1.2101 0.126 0.00 -0.654E-02 0. 0.444E+01 0.325E-02 0.634E-04 0.210E-04 19 8 25 2 17 0 1 12 2 0.000 0.000 2

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

Download in other formats:

Original Format