****************************** -------------------------------------------- * PALM (git SHA-1): 6b314732 * atmosphere - restart run ****************************** -------------------------------------------- Date: 2023-01-13 Run: ONKELTOMS_06_1d.r006 Time: 19:10:10 Run-No.: 05 Run on host: default Number of PEs: 64 Processor grid (x,y): ( 8, 8) prescribed ------------------------------------------------------------------------------ Numerical Schemes: ----------------- --> Use the 1.5-order turbulence closure (LES mode). --> Use the boussinesq approximation for the model equations. --> Solve perturbation pressure via multigrid method (w-cycle) number of grid levels: 1 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 --> Additional prognostic equation for the specific humidity ------------------------------------------------------------------------------ Run time and time step information: ---------------------------------- Timestep: variable maximum value: 20.000 s CFL-factor: 0.90 Start time: 18000.131 s End time: 43200.000 s Time reached: 21600.010 s CPU-time used: 15825.083 s per timestep: 1.791 s per second of simulated time: 4.396 s Next restart at: 25200.000 s restart every 3600.000 s Computational grid and domain size: ---------------------------------- Grid length: dx = 2.000 m dy = 2.000 m dz(1) = 2.000 m dz(2) = 6.000 m Domain size: x = 1200.000 m y = 1200.000 m z(u) = 349.254 m Vertical stretching starts at height: 99.0 m, Vertical stretching starts at index: 50, Vertical stretching ends at height: 9999999.9 m, Vertical stretching ends at index: 131, Factor used for stretching: 1.010, Number of gridpoints (x,y,z): (0: 599, 0: 599, 0: 131) Subdomain size (x,y,z): ( 75, 75, 132) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 0.00 0.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) = from soil model T. bound.: p(nzt+1) = p(nzt) | 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) q B. bound.: q(0) = from soil model T. bound.: q(nzt+1) = q(nzt) + dq/dz Bottom surface fluxes are used in diffusion terms at k=1 Constant flux layer between bottom surface and first computational u,v-level: z_mo = 1.00 m z0 = 0.1000 m z0h = 0.10000 m kappa = 0.40 zeta value range: -20.00 <= zeta <= 20.00 Predefined surface temperature Predefined surface humidity Lateral boundaries: left/right: cyclic north/south: cyclic Initial profiles: ---------------- Initial wind profiles (u,v) are interpolated from given profiles Characteristic levels of the initial temperature profile: Height: 0.0 100.0 300.0 m Temperature: 285.00 285.00 285.00 K Gradient: ------ 0.00 0.10 K/100m Gridpoint: 0 50 119 Characteristic levels of the initial humidity profile: Height: 0.0 0.0 m Humidity: 0.9E-02 0.9E-02 kg/kg Gradient: -------- -0.1E-02 (kg/kg)/100m Gridpoint: 0 0 List output: ----------- 1D-Profiles: Output every 9999999.90 s Time averaged over 300.00 s Averaging input every 0.00 s Data output: ----------- Time averaged over 300.00 s Averaging input every 0.00 s 1D-Profiles: Output format: parallel netCDF4/HDF5 Profile: wtheta, w"theta", w*theta*, w*2, theta, theta*2, rad_lw_hr, : rad_sw_hr, u, v, w, p, rho, q, rh, Output every 300.00 s Time averaged over 300.00 s Averaging input every 0.00 s 2D-Arrays: Output format: parallel netCDF4/HDF5 XY-cross-section Arrays: ta_2m*, theta_2m*, wspeed_10m*, tsurf*, Output every 300.00 s and at the start Cross sections at k = /0/2/4/6/8/10/12/14/16/18/20/ scalar-coordinates: /0.0/3.0/7.0/11.0/15.0/19.0/23.0/27.0/31.0/35.0/39.0/ m Number of output time levels allowed: 85 3D-Arrays: Output format: parallel netCDF4/HDF5 Arrays: u, v, w, q, theta, rad_lw_hr, rad_lw_in, rad_lw_out, rad_sw_hr, rad_s Output every 300.00 s and at the start Upper output limit at 349.25 m (GP 131) Number of output time levels allowed: 85 Time series: Output format: parallel netCDF4/HDF5 Output every 300.00 s Restart data format(s): Input format: mpi_shared_memory Output format: mpi_shared_memory ------------------------------------------------------------------------------ Physical quantities: ------------------- Geograph. latitude : latitude = 52.4 degr Geograph. longitude : longitude = 13.2 degr Rotation angle : rotation_angle = 0.0 degr Angular velocity : omega = 0.729E-04 rad/s Coriolis parameter : f = 0.000116 1/s f* = 0.000089 1/s Date and time at model start : 2019-06-01 06:00:00 +00 Gravity : g = 9.8 m/s**2 Reference state used in buoyancy terms: initial_profile Cloud physics quantities / methods: ---------------------------------- Humidity is considered, bu no condensation LES / Turbulence quantities: --------------------------- Mixing length is limited close to surfaces Actions during the simulation: ----------------------------- Disturbance impulse (u,v) every : 150.00 s Disturbance amplitude : 0.25 m/s Lower disturbance level : 5.00 m (GP 3) Upper disturbance level : 85.00 m (GP 43) Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2 Random number generator used : random-parallel Land surface model information: ------------------------------ --> Soil bottom is closed (water content is conserved, default) --> Land surface type : read from file --> Soil porosity type : read from file Initial soil temperature and moisture profile: Height: -0.01 -0.02 -0.05 -0.10 -0.20 -0.40 -0.80 -2. m Temperature: 290.15 290.00 290.00 290.00 290.00 290.00 290.00 290. K Moisture: 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0. m**3/m**3 Root fraction: 0.49 0.26 0.25 0.00 0.00 0.00 0.00 0. Grid point: 0 1 2 3 4 5 6 Vegetation canopy (drag) model: ------------------------------ Canopy mode: read_from_file Canopy height: 30.00m ( 15 grid points) Leaf drag coefficient: 0.30 Predefined constant heatflux at the top of the vegetation: 0.00 K m/s Characteristic levels of the leaf area density: Height: 0.0 0.0 m Leaf area density: 0.0000 0.00 m**2/m**3 Gradient: ------ 0.00 m**2/m**4 Gridpoint: 0 0 Radiation model information: ---------------------------- Albedo is set individually for each xy-location, according to given surface type. --> Albedo is fixed during the run Timestep: dt_radiation = 30.00 s ------------------------------------------------------------------------------