****************************** -------------------------------------------- * PALM 4.0 Rev: 2693 * atmosphere - run without 1D - prerun ****************************** -------------------------------------------- Date: 13-12-17 Run: ex_lsm_clearsky Time: 18:28:57 Run-No.: 00 Run on host: imukpalm4u Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated ------------------------------------------------------------------------------ Numerical Schemes: ----------------- --> Use the Moeng_Wyngaard turbulence closure (LES mode). --> Use the boussinesq approximation for the model equations. --> Solve perturbation pressure via FFT using temperton-algorithm routines 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: 0.000 s End time: 120.000 s Computational grid and domain size: ---------------------------------- Grid length: dx = 50.000 m dy = 50.000 m dz = 50.000 m Domain size: x = 2000.000 m y = 2000.000 m z(u) = 2862.511 m dz constant up to 1200.000 m (k= 25), then stretched by factor: 1.080 maximum dz not to be exceeded is dz_max = ********** m Number of gridpoints (x,y,z): (0: 39, 0: 39, 0: 41) Subdomain size (x,y,z): ( 20, 20, 42) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 5.00 5.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: flat Land surface model information: ------------------------------ --> Soil bottom is closed (water content is conserved, default) --> Land surface type : crops, mixed farming --> Soil porosity type : medium-fine Initial soil temperature and moisture profile: Height: -0.01 -0.02 -0.05 -0.11 -0.22 -0.40 -0.86 -2. m Temperature: 293.00 293.00 293.00 293.00 293.00 293.00 293.00 293. K Moisture: 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0. m**3/m**3 Root fraction: 0.03 0.07 0.14 0.14 0.28 0.09 0.22 0. Grid point: 0 1 2 3 4 5 6 Radiation model information: ---------------------------- --> Simple radiation scheme for clear sky is used (no clouds, default) Albedo is set for land surface type: mixed farming, tall grassland Timestep: dt_radiation = 60.00 s 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) = 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) 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 = 25.00 m z0 = 0.1000 m z0h = 0.10000 m kappa = 0.40 Rif value range: -20.00 <= rif <= 20.00 Predefined surface temperature Predefined surface humidity Lateral boundaries: left/right: cyclic north/south: cyclic Initial profiles: ---------------- Characteristic levels of the initial temperature profile: Height: 0.0 0.0 800.0 m Temperature: 293.00 293.00 293.00 K Gradient: ------ 0.00 1.00 K/100m Gridpoint: 0 0 16 Characteristic levels of the initial humidity profile: Height: 0.0 0.0 m Humidity: 0.1E-01 0.1E-01 kg/kg Gradient: -------- 0.0E+00 (kg/kg)/100m Gridpoint: 0 0 List output: ----------- 1D-Profiles: Output every ******** s Time averaged over 900.00 s Averaging input every 0.00 s Data output: ----------- Time averaged over 900.00 s Averaging input every 0.00 s 1D-Profiles: Output format: netCDF 64bit offset Profile: u, v, pt, km, kh, l, w, e, e*, p, w"pt", w*pt*, wpt, : w*e*, pt*2, u*2, v*2, w*2, t_soil, Output every 900.00 s Time averaged over 900.00 s Averaging input every 0.00 s 2D-Arrays: Output format: netCDF 64bit offset XY-cross-section Arrays: u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l Output every 3600.00 s Cross sections at k = /0/ scalar-coordinates: /0.0/ m Number of output time levels allowed: unlimited 2D-Arrays(time-averaged): Output format: netCDF 64bit offset XY-cross-section Arrays: u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l Output every 3600.00 s Time averaged over 1800.00 s Averaging input every 0.00 s Cross sections at k = /0/ scalar-coordinates: /0.0/ m Number of output time levels allowed: unlimited 3D-Arrays: Output format: netCDF 64bit offset Arrays: pt, vpt, q, u, v, w, t_soil, m_soil, Output every 3600.00 s Upper output limit at 2862.51 m (GP 41) Number of output time levels allowed: unlimited 3D-Arrays(time-averaged): Output format: netCDF 64bit offset Arrays: t_soil, m_soil, Output every 3600.00 s Time averaged over 1800.00 s Averaging input every 0.00 s Upper output limit at 2862.51 m (GP 41) Number of output time levels allowed: unlimited Time series: Output format: netCDF 64bit offset Output every 0.00 s ------------------------------------------------------------------------------ Physical quantities: ------------------- Geograph. latitude : latitude = 55.0 degr Geograph. longitude : longitude = 0.0 degr Angular velocity : omega = 0.729E-04 rad/s Coriolis parameter : f = 0.000119 1/s f* = 0.000084 1/s Day of the year at model start : day_init = 172 UTC time at model start : time_utc_init = 43200.0 s 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 to 1.80 * z Actions during the simulation: ----------------------------- Disturbance impulse (u,v) every : ****** s Disturbance amplitude : 0.25 m/s Lower disturbance level : 125.00 m (GP 3) Upper disturbance level : 625.00 m (GP 13) Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2 Random number generator used : random-parallel *** 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 8.6700A 5.1898D -0.2288D 0.1182 0.100 0.00 0.000E+00 800. 0.119E+02 0.103E-02 0.611E-03 0.888E-13 9 9 23 11 7 19 5 26 7 0.000 0.000 0 0 1 00:00:08.67 8.6900A 5.1806 0.2159 0.1162 0.004 0.26 -0.882E+01 100. 0.119E+02 0.867E-03 0.953E-05 0.439E-15 9 9 24 4 5 26 5 26 8 0.000 0.000 0 0 2 00:00:17.35 8.6900A 5.1762 0.1915 0.1151 0.000 0.04 -0.174E-01 100. 0.119E+02 0.793E-03 0.276E-05 0.403E-15 5 16 32 4 5 27 5 26 9 0.000 0.000 0 0 3 00:00:26.04 8.7000A 5.1723 0.1786 0.1113 0.005 0.27 -0.365E+01 100. 0.119E+02 0.750E-03 0.314E-05 0.376E-15 5 16 33 4 5 28 5 26 10 0.000 0.000 0 0 4 00:00:34.75 8.7000A 5.1732 -0.1675 0.1061 0.004 0.43 -0.144E+01 250. 0.119E+02 0.721E-03 0.272E-05 0.356E-15 5 16 34 8 29 32 5 26 11 0.000 0.000 0 0 5 00:00:43.45 8.7000A 5.1719 -0.1671 0.1028 0.005 0.43 -0.170E+01 250. 0.119E+02 0.698E-03 0.292E-05 0.335E-15 5 16 35 8 29 33 5 26 11 0.000 0.000 0 0 6 00:00:52.15 8.7000A 5.1700 -0.1621 0.1039 0.007 0.45 -0.194E+01 250. 0.119E+02 0.680E-03 0.303E-05 0.319E-15 5 16 36 8 29 34 5 26 12 0.000 0.000 0 0 7 00:01:00.85 8.7100A 5.1674 -0.1565 0.1044 0.009 0.65 -0.207E+01 550. 0.119E+02 0.664E-03 0.315E-05 0.303E-15 5 16 37 8 29 35 5 26 13 0.000 0.000 0 0 8 00:01:09.56 8.7100A 5.1641 0.1515 0.1041 0.011 0.70 -0.181E+01 550. 0.119E+02 0.651E-03 0.325E-05 0.290E-15 5 16 38 11 3 22 5 26 14 0.000 0.000 0 0 9 00:01:18.27 8.7200A 5.1634 -0.1506 0.1035 0.013 0.69 -0.138E+01 550. 0.119E+02 0.639E-03 0.332E-05 0.277E-15 5 16 38 8 29 36 5 26 15 0.000 0.000 0 0 10 00:01:26.99 8.7200A 5.1633 -0.1516 0.1026 0.016 0.73 -0.141E+01 550. 0.119E+02 0.629E-03 0.361E-05 0.266E-15 5 16 39 8 29 37 5 26 16 0.000 0.000 0 0 11 00:01:35.71 8.7200A 5.1628 -0.1516 0.1014 0.018 0.76 -0.136E+01 550. 0.119E+02 0.619E-03 0.379E-05 0.256E-15 5 16 0 8 29 38 5 26 17 0.000 0.000 0 0 12 00:01:44.43 8.7200A 5.1621 -0.1506 0.1000 0.020 0.79 -0.132E+01 550. 0.119E+02 0.610E-03 0.403E-05 0.246E-15 5 16 1 8 29 39 5 26 18 0.000 0.000 0 0 13 00:01:53.15 8.7200A 5.1610 -0.1488 0.0983 0.023 0.81 -0.129E+01 550. 0.119E+02 0.602E-03 0.430E-05 0.237E-15 5 16 2 8 29 0 5 26 19 0.000 0.000 0 0 14 00:02:01.87 8.7200A 5.1597 -0.1462 0.0964 0.025 0.83 -0.126E+01 550. 0.119E+02 0.593E-03 0.459E-05 0.227E-15 5 16 3 8 29 1 5 26 20 0.000 0.000 0