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: 3915M * atmosphere - run with 1D - prerun ****************************** -------------------------------------------- Date: 2019-07-17 Run: rans_tkee__gfortran_default__4 Time: 17:23:23 Run-No.: 00 Run on host: gfortran_d 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 *** dynamic module disabled ------------------------------------------------------------------------------ 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.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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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