1D run control output: ------------------------------ ITER. HH:MM:SS DT UMAX VMAX U* ALPHA ENERG. ------------------------------------------------------------- 0 00:00:00 0.01 3.00 0.10 0.100 45.0 4.45 3325 12:00:04 12.30 3.36 0.87 0.127 26.2 4.48 ****************************** -------------------------------------------- * PALM 6.0 Rev: 3527:35 * atmosphere - run with 1D - prerun ****************************** -------------------------------------------- Date: 2018-11-20 Run: rans_tkee__intel_default__4 Time: 16:32:26 Run-No.: 00 Run on host: intel_defa Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated ------------------------------------------------------------------------------ Numerical Schemes: ----------------- --> Use the TKE-e turbulence closure (RANS mode). --> Use the boussinesq approximation for the model equations. --> Solve perturbation pressure via multigrid method (w-cycle) number of grid levels: 4 Gauss-Seidel red/black iterations: 2 fixed number of multigrid cycles: 2 perturbation pressure is calculated at every Runge-Kutta step --> Momentum advection via Wicker-Skamarock-Scheme 5th order --> Scalar advection via Wicker-Skamarock-Scheme 5th order --> Loop optimization method: cache --> Time differencing scheme: runge-kutta-3 ------------------------------------------------------------------------------ Run time and time step information: ---------------------------------- Timestep: variable maximum value: 20.000 s CFL-factor: 0.90 Start time: 0.000 s End time: 50.000 s Computational grid and domain size: ---------------------------------- Grid length: dx = 10.000 m dy = 10.000 m dz(1) = 10.000 m dz(2) = 999.000 m Domain size: x = 320.000 m y = 320.000 m z(u) = 965.000 m Number of gridpoints (x,y,z): (0: 31, 0: 31, 0: 97) Subdomain size (x,y,z): ( 16, 16, 98) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 3.00 3.00 m/s Gradient: ------ 0.00 1/100s Gridpoint: 0 0 Characteristic levels of the geo. wind component vg: Height: 0.0 0.0 m vg: 0.00 0.00 m/s Gradient: ------ 0.00 1/100s Gridpoint: 0 0 Topography information: ---------------------- Topography: read_from_file Topography grid definition convention: cell center (scalar grid points) Boundary conditions: ------------------- p uv pt B. bound.: p(0) = p(1) | uv(0) = -uv(1) | pt(0) = pt(1) T. bound.: p(nzt+1) = 0 | uv(nzt+1) = ug(nzt+1), vg(nzt+1) | pt(nzt+1) = pt(nzt) + dpt/dz_ e B. bound.: e(0) = e(1) T. bound.: e(nzt+1) = e(nzt) = e(nzt-1) Bottom surface fluxes are used in diffusion terms at k=1 Predefined constant heatflux: 0.001000 K m/s Constant flux layer between bottom surface and first computational u,v-level: z_mo = 5.00 m z0 = 0.1000 m z0h = 0.10000 m kappa = 0.40 Rif value range: -20.00 <= rif <= 20.00 Lateral boundaries: left/right: dirichlet/radiation north/south: cyclic use_cmax: T pt damping layer width = 0.00 m, pt damping factor = 0.0000 Initial profiles: ---------------- Characteristic levels of the initial temperature profile: Height: 0.0 0.0 m Temperature: 300.00 300.00 K Gradient: ------ 0.00 K/100m Gridpoint: 0 0 List output: ----------- 1D-Profiles: Output every 9999999.90 s Data output: ----------- 1D-Profiles: Output format: netCDF 64bit offset Profile: u, w"u", w*u*, wu, u*2, v, w"v", w*v*, wv, v*2, w, w*2, : e, e*, km, l, Output every 10.00 s Time averaged over 0.00 s Averaging input every 0.00 s Time series: Output format: netCDF 64bit offset Output every 1.00 s ------------------------------------------------------------------------------ Physical quantities: ------------------- Geograph. latitude : latitude = 52.5 degr Geograph. longitude : longitude = 13.3 degr Angular velocity : omega = 0.729E-04 rad/s Coriolis parameter : f = 0.000116 1/s f* = 0.000089 1/s Day of the year at model start : day_init = 172 UTC time at model start : time_utc_init = 0.0 s Gravity : g = 9.8 m/s**2 Reference state used in buoyancy terms: initial_profile Cloud physics quantities / methods: ---------------------------------- LES / Turbulence quantities: --------------------------- initial TKE is prescribed as 0.10E+00 (m/s)**2 TKE is not allowed to fall below 0.10E-02 (m/s)**2 Mixing length is limited to 1.80 * z RANS-mode constants: c_0 = 0.55000 c_1 = 1.44000 c_2 = 1.92000 c_3 = 1.44000 c_4 = 0.00000 sigma_e = 1.00000 sigma_diss = 1.30000 Actions during the simulation: ----------------------------- 1D-Model parameters: ------------------- Simulation time: 43200.0 s Run-controll output every: 43200.0 s Vertical profile output every: 86400.0 s Mixing length calculation: blackadar Dissipation calculation: prognostic *** no user-defined variables found ------------------------------------------------------------------------------ Run-control output: ------------------ 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 --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 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 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 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 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 0 4 00:00:00.30 0.3200A 3.3845 1.7361 1.3742 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 0 5 00:00:00.62 0.6400A 3.3838 1.7102 1.3740 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 0 6 00:00:01.26 1.2800A 3.3833 1.6639 1.3659 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 0 7 00:00:02.54 2.5600A 3.3827 1.5852 1.3337 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 0 8 00:00:05.09 2.6600A 3.3826 1.4669 1.2462 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 0 9 00:00:07.75 2.6200D 3.3825 1.4246 1.2019 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 0 10 00:00:10.37 2.6600A 3.3822 1.4068 1.2181 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 0 11 00:00:13.03 2.6600A 3.3818 1.3947 1.2421 0.138 0.00 -0.554E-02 0. 0.444E+01 0.372E-02 0.913E-04 0.361E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2 0 12 00:00:15.69 2.6600A 3.3818 1.3891 1.2528 0.137 0.00 -0.573E-02 0. 0.444E+01 0.366E-02 0.869E-04 0.301E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2 0 13 00:00:18.35 2.6600A 3.3820 1.3872 1.2552 0.135 0.00 -0.581E-02 0. 0.444E+01 0.362E-02 0.866E-04 0.283E-04 19 6 23 2 17 0 1 12 2 0.000 0.000 2 0 14 00:00:21.01 2.6600A 3.3822 1.3883 1.2530 0.134 0.00 -0.617E-02 0. 0.444E+01 0.357E-02 0.864E-04 0.278E-04 19 6 23 2 17 0 1 12 2 0.000 0.000 2 0 15 00:00:23.67 2.6600A 3.3824 1.3916 1.2483 0.133 0.00 -0.643E-02 0. 0.444E+01 0.352E-02 0.802E-04 0.251E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2 0 16 00:00:26.33 2.6600A 3.3825 1.3968 1.2433 0.132 0.00 -0.634E-02 0. 0.444E+01 0.348E-02 0.758E-04 0.252E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 17 00:00:29.00 2.6600A 3.3827 1.4033 1.2388 0.131 0.00 -0.649E-02 0. 0.444E+01 0.345E-02 0.737E-04 0.244E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 18 00:00:31.66 2.6600A 3.3829 1.4106 1.2353 0.130 0.00 -0.651E-02 0. 0.444E+01 0.342E-02 0.723E-04 0.234E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 19 00:00:34.32 2.6600A 3.3831 1.4182 1.2329 0.129 0.00 -0.661E-02 0. 0.444E+01 0.339E-02 0.713E-04 0.230E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 20 00:00:36.98 2.6600A 3.3833 1.4254 1.2314 0.128 0.00 -0.648E-02 0. 0.444E+01 0.336E-02 0.697E-04 0.226E-04 19 8 23 2 17 0 1 12 2 0.000 0.000 2 0 21 00:00:39.64 2.6600D 3.3835 1.4320 1.2306 0.128 0.00 -0.649E-02 0. 0.444E+01 0.333E-02 0.680E-04 0.224E-04 19 8 23 2 17 0 1 12 2 0.000 0.000 2 0 22 00:00:42.29 2.5300D 3.3837 1.4377 1.2303 0.127 0.00 -0.659E-02 0. 0.444E+01 0.331E-02 0.667E-04 0.221E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 23 00:00:44.82 2.4200D 3.3839 1.4475 1.2348 0.126 0.00 -0.665E-02 0. 0.444E+01 0.329E-02 0.633E-04 0.210E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 24 00:00:47.25 2.3200D 3.3840 1.4563 1.2402 0.126 0.00 -0.662E-02 0. 0.444E+01 0.327E-02 0.603E-04 0.200E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 25 00:00:49.57 2.2400D 3.3842 1.4642 1.2457 0.126 0.00 -0.657E-02 0. 0.444E+01 0.326E-02 0.576E-04 0.190E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 26 00:00:51.81 2.1600D 3.3844 1.4711 1.2506 0.125 0.00 -0.652E-02 0. 0.444E+01 0.324E-02 0.553E-04 0.183E-04 19 8 25 2 17 0 1 12 2 0.000 0.000 2