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: 4102M * atmosphere - run with 1D - prerun ****************************** -------------------------------------------- Date: 2019-07-17 Run: rans_tkee__gfortran_default__4 Time: 18:21:01 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 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 0 1 00:00:02.67 2.6600A 3.3803 1.5756 1.1292 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 0 2 00:00:05.33 0.2370D 3.3829 1.4756 1.2412 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 0 3 00:00:05.56 0.5360D 3.3841 1.5597 1.3639 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 0 4 00:00:06.10 0.7840D 3.3844 1.5492 1.3622 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 0 5 00:00:06.88 1.0200D 3.3843 1.5263 1.3422 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 0 6 00:00:07.90 1.2500D 3.3842 1.5019 1.3159 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 0 7 00:00:09.15 1.2700D 3.3842 1.4782 1.2869 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 0 8 00:00:10.42 1.3200D 3.3840 1.4657 1.2702 0.140 0.00 -0.558E-02 0. 0.444E+01 0.387E-02 0.621E-04 0.163E-04 19 2 0 2 17 0 2 13 1 0.000 0.000 2 0 9 00:00:11.74 1.4300D 3.3835 1.4548 1.2674 0.139 0.00 -0.563E-02 0. 0.444E+01 0.380E-02 0.561E-04 0.155E-04 19 1 0 2 17 0 1 12 2 0.000 0.000 2 0 10 00:00:13.17 1.5700D 3.3829 1.4428 1.2843 0.138 0.00 -0.566E-02 0. 0.444E+01 0.376E-02 0.567E-04 0.170E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2 0 11 00:00:14.74 1.7500D 3.3824 1.4305 1.2946 0.137 0.00 -0.573E-02 0. 0.444E+01 0.371E-02 0.586E-04 0.172E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2 0 12 00:00:16.49 2.0100D 3.3820 1.4178 1.2979 0.136 0.00 -0.590E-02 0. 0.444E+01 0.367E-02 0.621E-04 0.182E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2 0 13 00:00:18.50 2.2400D 3.3816 1.4031 1.2924 0.135 0.00 -0.602E-02 0. 0.444E+01 0.362E-02 0.673E-04 0.202E-04 20 2 0 2 17 0 1 12 2 0.000 0.000 2 0 14 00:00:20.74 2.3200D 3.3814 1.3906 1.2814 0.134 0.00 -0.646E-02 0. 0.444E+01 0.358E-02 0.702E-04 0.216E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2 0 15 00:00:23.06 2.3300D 3.3815 1.3849 1.2710 0.133 0.00 -0.625E-02 0. 0.444E+01 0.354E-02 0.697E-04 0.221E-04 19 6 24 2 17 0 1 12 2 0.000 0.000 2 0 16 00:00:25.39 2.3200D 3.3816 1.3828 1.2623 0.132 0.00 -0.627E-02 0. 0.444E+01 0.350E-02 0.681E-04 0.217E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 17 00:00:27.71 2.3000D 3.3817 1.3820 1.2552 0.131 0.00 -0.638E-02 0. 0.444E+01 0.347E-02 0.663E-04 0.211E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 18 00:00:30.01 2.2900D 3.3818 1.3818 1.2496 0.130 0.00 -0.655E-02 0. 0.444E+01 0.344E-02 0.644E-04 0.205E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 19 00:00:32.30 2.2900D 3.3820 1.3815 1.2451 0.130 0.00 -0.652E-02 0. 0.444E+01 0.342E-02 0.627E-04 0.201E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 20 00:00:34.59 2.2900D 3.3821 1.3808 1.2413 0.129 0.00 -0.647E-02 0. 0.444E+01 0.339E-02 0.613E-04 0.198E-04 19 7 23 2 17 0 1 12 2 0.000 0.000 2 0 21 00:00:36.88 2.3000D 3.3823 1.3802 1.2384 0.128 0.00 -0.646E-02 0. 0.444E+01 0.337E-02 0.601E-04 0.195E-04 19 8 23 2 17 0 1 12 2 0.000 0.000 2 0 22 00:00:39.18 2.3100D 3.3825 1.3792 1.2359 0.128 0.00 -0.650E-02 0. 0.444E+01 0.335E-02 0.592E-04 0.193E-04 19 8 23 2 17 0 1 12 2 0.000 0.000 2 0 23 00:00:41.49 2.3100D 3.3826 1.3783 1.2338 0.127 0.00 -0.658E-02 0. 0.444E+01 0.333E-02 0.585E-04 0.191E-04 19 8 23 2 17 0 1 12 2 0.000 0.000 2 0 24 00:00:43.80 2.3200D 3.3828 1.3778 1.2324 0.127 0.00 -0.665E-02 0. 0.444E+01 0.331E-02 0.576E-04 0.189E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 25 00:00:46.12 2.3200D 3.3830 1.3770 1.2310 0.126 0.00 -0.659E-02 0. 0.444E+01 0.329E-02 0.570E-04 0.188E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 26 00:00:48.44 2.3300D 3.3831 1.3766 1.2300 0.126 0.00 -0.653E-02 0. 0.444E+01 0.327E-02 0.563E-04 0.187E-04 19 8 24 2 17 0 1 12 2 0.000 0.000 2 0 27 00:00:50.77 2.3300D 3.3833 1.3759 1.2289 0.125 0.00 -0.649E-02 0. 0.444E+01 0.326E-02 0.558E-04 0.186E-04 19 9 24 2 17 0 1 12 2 0.000 0.000 2