****************************** -------------------------------------------- * PALM 6.0 Rev: 4781 * atmosphere - run without 1D - prerun ****************************** -------------------------------------------- Date: 2020-11-12 Run: urban_environment_salsa Time: 17:46:37 Run-No.: 00 Run on host: default Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated ------------------------------------------------------------------------------ Numerical Schemes: ----------------- --> Use the 1.5-order turbulence closure (LES mode). --> Use the boussinesq approximation for the model equations. --> Solve perturbation pressure via FFT using fftw 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: 1231.000 s Spinup time: 1200.000 s Computational grid and domain size: ---------------------------------- Grid length: dx = 2.000 m dy = 2.000 m dz(1) = 2.000 m Domain size: x = 40.000 m y = 40.000 m z(u) = 121.000 m Number of gridpoints (x,y,z): (0: 19, 0: 19, 0: 61) Subdomain size (x,y,z): ( 10, 10, 62) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 1.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) = 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 = 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: ---------------- Characteristic levels of the initial temperature profile: Height: 0.0 0.0 m Temperature: 277.15 277.12 K Gradient: ------ 0.00 K/100m Gridpoint: 0 0 Characteristic levels of the initial humidity profile: Height: 0.0 0.0 m Humidity: 0.1E-02 0.4E-02 kg/kg Gradient: -------- 0.0E+00 (kg/kg)/100m Gridpoint: 0 0 List output: ----------- 1D-Profiles: Output every 9999999.90 s Time averaged over 30.00 s Averaging input every 0.00 s Data output: ----------- Time averaged over 30.00 s Averaging input every 0.00 s 1D-Profiles: Output format: netCDF 64bit offset Profile: u, u*2, wu, w*u*, w"u", v, v*2, wv, w*v*, w"v", theta, : theta*2, wtheta, w*theta*, w"theta", q, q*2, e, e*, : km, l, t_soil, m_soil, rad_lw_in, rad_lw_out, salsa_Ntot, : salsa_N_UFP, salsa_LDSA, salsa_PM0.1, salsa_PM2.5, salsa_PM10 : Output every 30.00 s Time averaged over 30.00 s Averaging input every 0.00 s Time series: Output format: netCDF 64bit offset Output every 5.00 s Restart data format(s): Output format: fortran_binary ------------------------------------------------------------------------------ Physical quantities: ------------------- Geograph. latitude : latitude = 52.5 degr Geograph. longitude : longitude = 13.3 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-03-06 10: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 : ****** s Disturbance amplitude : 0.25 m/s Lower disturbance level : 5.00 m (GP 3) Upper disturbance level : 39.00 m (GP 20) Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2 Random number generator used : random-parallel *** dynamic module disabled 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: 277.15 277.00 277.00 277.00 277.00 277.00 277.00 277. K Moisture: 0.30 0.30 0.30 0.30 0.30 0.30 0.30 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: 10.00m ( 5 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: ---------------------------- --> RRTMG scheme is used Albedo is set individually for each xy-location, according to given surface type. Timestep: dt_radiation = 60.00 s SALSA information: ------------------------------ Starts at: skip_time_do_salsa = 10.00 s Timestep: dt_salsa = 2.00 s Array shape (z,y,x,bins): aerosol_number: 60 19 19 8 aerosol_mass: 60 19 19 32 (advect_particle_water = T) salsa_gas: 60 19 19 5 (salsa_gases_from_chem = F) Aerosol dynamic processes included: dry deposition (on vegetation = T and on topography = T) Aerosol bin subrange limits (in metres): 3.00E-009 1.00E-008 2.50E-006 Number of size bins for each aerosol subrange: 1 7 Aerosol bin lower limits (in metres): 3.00E-009 1.00E-008 2.20E-008 4.84E-008 1.07E-007 2.35E-007 5.16E-007 1.14E-006 Bin geometric mean diameters (in metres): 5.48E-009 1.48E-008 3.26E-008 7.18E-008 1.58E-007 3.48E-007 7.66E-007 1.69E-006 Number of chemical components used: 3 Species: OC NO NH Initial relative contribution of each species to particle volume in: a-bins: 0.800 0.100 0.100 0.000 0.000 0.000 0.000 b-bins: 0.000 0.000 0.000 0.000 0.000 0.000 0.000 Number of gaseous tracers used: 5 Initial gas concentrations: H2SO4: 5.0000E+012 #/m**3 HNO3: 3.0000E+012 #/m**3 NH3: 6.0000E+012 #/m**3 OCNV: 1.0000E+012 #/m**3 OCSV: 1.0000E+012 #/m**3 Initialising concentrations: Aerosol size distribution: init_aerosol_type = 1 Gas concentrations: init_gases_type = 1 Size distribution read from a file. Emissions: salsa_emission_mode = read_from_file (currently all emissions are soluble!) Boundary conditions for aerosols: bottom/top: neumann / dirichlet north/south: cyclic / cyclic left/right: dirichlet / dirichlet ------------------------------------------------------------------------------ Spinup control output:--------------------------------- ITER. HH:MM:SS DT PT(z_MO)--------------------------------- 1 -00:18:00 120.00 288.14 2 -00:16:00 120.00 288.17 3 -00:14:00 120.00 288.21 4 -00:12:00 120.00 288.24 5 -00:10:00 120.00 288.27 6 -00:08:00 120.00 288.30 7 -00:06:00 120.00 288.33 8 -00:04:00 120.00 288.35 9 -00:02:00 120.00 288.38 10 00:00:00 120.00 288.41 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.7580A -2.3731D -1.1942D 0.4177 0.060 0.00 -0.147E+00 16. 0.137E+01 0.743E-02 0.234E-01 0.446E-02 60 17 18 19 0 14 20 18 0 0.000 0.000 0 0 1 00:00:00.75 0.7580A -2.3732 -1.5382 0.7410 0.084 0.33 -0.778E+00 2. 0.137E+01 0.116E-01 0.268E-02 0.892E-03 60 18 15 16 15 0 20 18 0 0.000 0.000 0 0 2 00:00:01.51 0.7580A -2.3733 -1.6029 0.8136 0.083 0.00 -0.796E+00 2. 0.137E+01 0.125E-01 0.953E-03 0.351E-03 60 17 14 16 15 0 20 17 0 0.000 0.000 0 0 3 00:00:02.27 0.7580A -2.3733 -1.5856 0.8389 0.082 0.00 -0.801E+00 2. 0.137E+01 0.129E-01 0.818E-03 0.240E-03 60 17 13 17 15 0 20 17 0 0.000 0.000 0 0 4 00:00:03.03 0.7580A -2.3734 -1.5378 0.8459 0.082 0.12 -0.806E+00 12. 0.137E+01 0.132E-01 0.778E-03 0.213E-03 60 17 13 17 15 0 20 17 0 0.000 0.000 0 0 5 00:00:03.79 0.7580A -2.3735 -1.5013 0.8437 0.081 0.16 -0.809E+00 12. 0.137E+01 0.134E-01 0.769E-03 0.205E-03 60 17 12 17 14 0 20 17 0 0.000 0.000 0 0 6 00:00:04.54 0.7580A -2.3736 -1.4549 0.8375 0.080 0.17 -0.810E+00 12. 0.137E+01 0.137E-01 0.763E-03 0.201E-03 60 17 11 17 14 0 20 17 0 0.000 0.000 0 0 7 00:00:05.30 0.7580A -2.3736 -1.4322 0.8301 0.080 0.18 -0.811E+00 12. 0.137E+01 0.140E-01 0.759E-03 0.199E-03 60 16 10 11 14 0 20 17 0 0.000 0.000 0 0 8 00:00:06.06 0.7580A -2.3737 -1.4026 0.8224 0.079 0.18 -0.811E+00 12. 0.137E+01 0.143E-01 0.758E-03 0.197E-03 60 16 10 11 14 0 20 17 0 0.000 0.000 0 0 9 00:00:06.82 0.7580A -2.3738 -1.3753 0.8146 0.078 0.19 -0.811E+00 12. 0.137E+01 0.146E-01 0.758E-03 0.196E-03 60 16 9 11 14 0 20 17 0 0.000 0.000 0 0 10 00:00:07.58 0.7580A -2.3739 -1.3531 0.8073 0.078 0.24 -0.812E+00 18. 0.137E+01 0.149E-01 0.759E-03 0.196E-03 60 15 8 11 14 0 20 17 0 0.000 0.000 0 0 11 00:00:08.33 0.7580A -2.3739 -1.3349 0.8010 0.077 0.24 -0.812E+00 18. 0.137E+01 0.152E-01 0.761E-03 0.196E-03 60 15 7 11 14 0 20 17 0 0.000 0.000 0 0 12 00:00:09.09 0.7580A -2.3740 -1.3227 0.7963 0.077 0.24 -0.813E+00 18. 0.136E+01 0.156E-01 0.765E-03 0.195E-03 60 15 6 16 14 0 20 17 0 0.000 0.000 0 0 13 00:00:09.85 0.7580A -2.3741 -1.3176 0.7934 0.076 0.25 -0.814E+00 18. 0.136E+01 0.160E-01 0.769E-03 0.195E-03 60 14 5 16 14 0 20 17 0 0.000 0.000 0 0 14 00:00:10.61 0.7580A -2.3742 -1.3089 0.7919 0.076 0.25 -0.815E+00 18. 0.136E+01 0.163E-01 0.774E-03 0.195E-03 60 14 5 16 14 0 20 17 0 0.000 0.000 0 0 15 00:00:11.37 0.7580A -2.3743 -1.2968 0.7912 0.076 0.25 -0.815E+00 18. 0.136E+01 0.167E-01 0.778E-03 0.194E-03 60 13 4 16 14 0 20 17 0 0.000 0.000 0 0 16 00:00:12.12 0.7580A -2.3743 -1.2813 0.7902 0.075 0.25 -0.814E+00 18. 0.136E+01 0.171E-01 0.781E-03 0.193E-03 60 13 3 16 14 0 20 17 0 0.000 0.000 0 0 17 00:00:12.88 0.7580A -2.3744 -1.2624 0.8043 0.075 0.25 -0.814E+00 18. 0.136E+01 0.175E-01 0.783E-03 0.192E-03 60 12 2 16 14 0 11 14 18 0.000 0.000 0 0 18 00:00:13.64 0.7580A -2.3745 -1.2557 0.8321 0.074 0.25 -0.813E+00 18. 0.136E+01 0.179E-01 0.787E-03 0.191E-03 60 12 1 11 14 0 10 14 18 0.000 0.000 0 0 19 00:00:14.40 0.7580A -2.3746 -1.2519 0.8629 0.074 0.26 -0.813E+00 18. 0.136E+01 0.183E-01 0.790E-03 0.189E-03 60 12 1 11 14 0 10 14 18 0.000 0.000 0 0 20 00:00:15.16 0.7580A -2.3746 -1.2492 0.8923 0.074 0.26 -0.813E+00 18. 0.136E+01 0.187E-01 0.793E-03 0.188E-03 60 12 0 11 14 0 10 14 18 0.000 0.000 0 0 21 00:00:15.91 0.7580A -2.3747 -1.2475 0.9185 0.073 0.26 -0.813E+00 18. 0.136E+01 0.191E-01 0.797E-03 0.186E-03 60 11 19 11 14 0 10 14 18 0.000 0.000 0 0 22 00:00:16.67 0.7580A -2.3748 -1.2465 0.9402 0.073 0.26 -0.813E+00 18. 0.136E+01 0.195E-01 0.800E-03 0.185E-03 60 11 18 11 14 0 10 14 18 0.000 0.000 0 0 23 00:00:17.43 0.7580A -2.3749 -1.2461 0.9566 0.073 0.26 -0.813E+00 18. 0.136E+01 0.199E-01 0.803E-03 0.183E-03 60 10 17 11 14 0 10 14 18 0.000 0.000 0 0 24 00:00:18.19 0.7580A -2.3749 -1.2461 0.9677 0.072 0.26 -0.813E+00 18. 0.136E+01 0.203E-01 0.805E-03 0.182E-03 60 10 16 11 14 0 10 14 18 0.000 0.000 0 0 25 00:00:18.95 0.7580A -2.3750 -1.2465 0.9738 0.072 0.26 -0.813E+00 18. 0.136E+01 0.207E-01 0.809E-03 0.180E-03 60 10 16 11 14 0 10 14 18 0.000 0.000 0 0 26 00:00:19.70 0.7580A -2.3751 -1.2469 0.9757 0.072 0.27 -0.813E+00 18. 0.136E+01 0.211E-01 0.812E-03 0.179E-03 60 10 15 11 14 0 10 14 18 0.000 0.000 0 0 27 00:00:20.46 0.7580A -2.3752 -1.2473 0.9748 0.071 0.27 -0.813E+00 18. 0.136E+01 0.215E-01 0.816E-03 0.178E-03 60 9 14 11 14 0 11 14 18 0.000 0.000 0 0 28 00:00:21.22 0.7580A -2.3752 -1.2474 0.9803 0.071 0.38 -0.813E+00 40. 0.136E+01 0.219E-01 0.819E-03 0.177E-03 60 9 13 11 14 0 11 14 18 0.000 0.000 0 0 29 00:00:21.98 0.7580A -2.3753 -1.2475 0.9810 0.071 0.38 -0.813E+00 40. 0.136E+01 0.223E-01 0.822E-03 0.176E-03 60 9 12 11 14 0 11 14 18 0.000 0.000 0 0 30 00:00:22.74 0.7580A -2.3754 -1.2473 0.9775 0.071 0.38 -0.813E+00 40. 0.136E+01 0.226E-01 0.826E-03 0.175E-03 60 8 11 11 14 0 11 14 18 0.000 0.000 0 0 31 00:00:23.49 0.7580A -2.3755 -1.2471 0.9703 0.070 0.39 -0.813E+00 40. 0.136E+01 0.230E-01 0.830E-03 0.174E-03 60 8 11 11 14 0 11 14 18 0.000 0.000 0 0 32 00:00:24.25 0.7580A -2.3755 -1.2468 0.9603 0.070 0.39 -0.813E+00 40. 0.136E+01 0.234E-01 0.834E-03 0.173E-03 60 8 10 11 14 0 11 14 18 0.000 0.000 0 0 33 00:00:25.01 0.7580A -2.3756 -1.2463 0.9497 0.070 0.39 -0.813E+00 40. 0.136E+01 0.238E-01 0.838E-03 0.172E-03 60 7 9 11 14 0 10 14 18 0.000 0.000 0 0 34 00:00:25.77 0.7580A -2.3757 -1.2455 0.9459 0.070 0.39 -0.814E+00 40. 0.136E+01 0.242E-01 0.841E-03 0.171E-03 60 7 8 11 14 0 10 14 18 0.000 0.000 0 0 35 00:00:26.53 0.7580A -2.3758 -1.2443 0.9418 0.069 0.40 -0.814E+00 40. 0.136E+01 0.246E-01 0.845E-03 0.170E-03 60 7 7 11 14 0 10 14 18 0.000 0.000 0 0 36 00:00:27.28 0.7580A -2.3758 -1.2425 0.9449 0.069 0.40 -0.814E+00 40. 0.136E+01 0.250E-01 0.848E-03 0.169E-03 60 6 6 11 14 0 11 14 17 0.000 0.000 0 0 37 00:00:28.04 0.7580A -2.3759 -1.2401 0.9570 0.069 0.40 -0.814E+00 40. 0.136E+01 0.254E-01 0.850E-03 0.168E-03 60 6 5 11 14 0 11 14 17 0.000 0.000 0 0 38 00:00:28.80 0.7580A -2.3760 -1.2369 0.9688 0.069 0.40 -0.814E+00 40. 0.136E+01 0.257E-01 0.853E-03 0.167E-03 60 6 5 11 14 0 11 14 17 0.000 0.000 0 0 39 00:00:29.56 0.7580A -2.3761 -1.2329 0.9791 0.069 0.41 -0.814E+00 40. 0.136E+01 0.261E-01 0.856E-03 0.167E-03 60 5 3 11 14 0 11 14 17 0.000 0.000 0 0 40 00:00:30.32 0.7580A -2.3761 -1.2305 0.9870 0.068 0.41 -0.815E+00 40. 0.136E+01 0.265E-01 0.858E-03 0.166E-03 60 5 3 11 13 0 11 14 17 0.000 0.000 0 0 41 00:00:31.07 0.7580A -2.3762 -1.2321 0.9920 0.068 0.41 -0.815E+00 40. 0.136E+01 0.269E-01 0.861E-03 0.165E-03 60 5 2 11 13 0 11 14 17 0.000 0.000 0