****************************** -------------------------------------------- * PALM 6.0 Rev: 4746M * atmosphere - run without 1D - prerun ****************************** -------------------------------------------- Date: 2020-10-16 Run: urban_environment_salsa__gfortran_ Time: 09:29:27 Run-No.: 00 Run on host: gfortran_d 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 ------------------------------------------------------------------------------ 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.065 0.00 -0.212E+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.085 0.33 -0.917E+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.8137 0.084 0.00 -0.940E+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.5857 0.8393 0.083 0.00 -0.956E+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.5379 0.8466 0.083 0.13 -0.968E+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.5016 0.8447 0.082 0.17 -0.975E+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.4553 0.8392 0.081 0.23 -0.979E+00 18. 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.4291 0.8325 0.081 0.24 -0.982E+00 18. 0.137E+01 0.140E-01 0.760E-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.3990 0.8253 0.080 0.25 -0.983E+00 18. 0.137E+01 0.143E-01 0.759E-03 0.198E-03 60 16 9 11 14 0 20 17 0 0.000 0.000 0 0 9 00:00:06.82 0.7580A -2.3738 -1.3711 0.8181 0.080 0.26 -0.985E+00 18. 0.137E+01 0.146E-01 0.760E-03 0.197E-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.3481 0.8113 0.079 0.26 -0.986E+00 18. 0.137E+01 0.150E-01 0.761E-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.3293 0.8055 0.079 0.27 -0.987E+00 18. 0.137E+01 0.153E-01 0.764E-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.3242 0.8011 0.078 0.27 -0.988E+00 18. 0.136E+01 0.157E-01 0.768E-03 0.196E-03 60 14 6 16 14 0 20 17 0 0.000 0.000 0 0 13 00:00:09.85 0.7580A -2.3741 -1.3192 0.7985 0.078 0.27 -0.988E+00 18. 0.136E+01 0.161E-01 0.773E-03 0.196E-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.3107 0.7973 0.077 0.28 -0.989E+00 18. 0.136E+01 0.165E-01 0.778E-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.2988 0.8191 0.077 0.28 -0.989E+00 18. 0.136E+01 0.169E-01 0.783E-03 0.195E-03 60 13 4 16 14 0 11 14 18 0.000 0.000 0 0 16 00:00:12.12 0.7580A -2.3743 -1.2833 0.8557 0.076 0.35 -0.989E+00 28. 0.136E+01 0.173E-01 0.786E-03 0.194E-03 60 13 3 16 14 0 11 14 18 0.000 0.000 0 0 17 00:00:12.88 0.7580A -2.3744 -1.2644 0.8849 0.076 0.35 -0.988E+00 28. 0.136E+01 0.177E-01 0.788E-03 0.193E-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.2474 0.9204 0.076 0.42 -0.987E+00 40. 0.136E+01 0.181E-01 0.792E-03 0.192E-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.2427 0.9588 0.075 0.42 -0.986E+00 40. 0.136E+01 0.186E-01 0.796E-03 0.190E-03 60 12 0 11 14 0 10 14 18 0.000 0.000 0 0 20 00:00:15.16 0.7580A -2.3746 -1.2390 0.9948 0.075 0.43 -0.985E+00 40. 0.136E+01 0.190E-01 0.799E-03 0.189E-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.2358 1.0265 0.075 0.43 -0.985E+00 40. 0.136E+01 0.195E-01 0.803E-03 0.187E-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.2331 1.0525 0.075 0.43 -0.984E+00 40. 0.136E+01 0.200E-01 0.806E-03 0.186E-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.2309 1.0724 0.074 0.44 -0.983E+00 40. 0.136E+01 0.204E-01 0.809E-03 0.185E-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.2290 1.0883 0.074 0.44 -0.983E+00 40. 0.136E+01 0.209E-01 0.812E-03 0.183E-03 60 10 16 11 14 0 11 14 18 0.000 0.000 0 0 25 00:00:18.95 0.7580A -2.3750 -1.2275 1.1158 0.074 0.45 -0.982E+00 40. 0.136E+01 0.214E-01 0.817E-03 0.182E-03 60 10 15 11 14 0 11 14 18 0.000 0.000 0 0 26 00:00:19.70 0.7580A -2.3751 -1.2262 1.1375 0.073 0.45 -0.981E+00 40. 0.136E+01 0.219E-01 0.821E-03 0.181E-03 60 10 15 11 14 0 11 14 18 0.000 0.000 0 0 27 00:00:20.46 0.7580A -2.3752 -1.2252 1.1526 0.073 0.46 -0.981E+00 40. 0.136E+01 0.224E-01 0.825E-03 0.180E-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.2243 1.1612 0.073 0.46 -0.980E+00 40. 0.136E+01 0.229E-01 0.829E-03 0.179E-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.2232 1.1641 0.073 0.46 -0.980E+00 40. 0.136E+01 0.234E-01 0.835E-03 0.179E-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.2219 1.1623 0.072 0.47 -0.980E+00 40. 0.136E+01 0.239E-01 0.840E-03 0.178E-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.2202 1.1578 0.072 0.47 -0.980E+00 40. 0.136E+01 0.245E-01 0.846E-03 0.178E-03 60 8 10 11 14 0 11 14 18 0.000 0.000 0 0 32 00:00:24.25 0.7580A -2.3756 -1.2178 1.1545 0.072 0.47 -0.979E+00 40. 0.136E+01 0.250E-01 0.853E-03 0.177E-03 60 8 10 11 14 0 13 14 18 0.000 0.000 0 0 33 00:00:25.01 0.7580A -2.3756 -1.2145 1.1607 0.072 0.48 -0.979E+00 40. 0.136E+01 0.255E-01 0.860E-03 0.177E-03 60 7 8 11 14 0 13 14 18 0.000 0.000 0 0 34 00:00:25.77 0.7580A -2.3757 -1.2098 1.1901 0.072 0.48 -0.979E+00 40. 0.136E+01 0.261E-01 0.867E-03 0.177E-03 60 7 8 11 14 0 11 14 17 0.000 0.000 0 0 35 00:00:26.53 0.7580A -2.3758 -1.2035 1.2147 0.071 0.48 -0.979E+00 40. 0.136E+01 0.267E-01 0.874E-03 0.177E-03 60 7 7 11 14 0 11 14 17 0.000 0.000 0 0 36 00:00:27.28 0.7580A -2.3759 -1.1953 1.2288 0.071 0.49 -0.979E+00 40. 0.136E+01 0.272E-01 0.882E-03 0.177E-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.1850 1.2434 0.071 0.49 -0.978E+00 40. 0.136E+01 0.278E-01 0.889E-03 0.177E-03 60 6 5 11 14 0 12 14 17 0.000 0.000 0 0 38 00:00:28.80 0.7580A -2.3760 -1.1722 1.2661 0.071 0.49 -0.978E+00 40. 0.136E+01 0.284E-01 0.897E-03 0.178E-03 60 6 4 11 14 0 14 14 16 0.000 0.000 0 0 39 00:00:29.56 0.7580A -2.3761 -1.1605 1.2932 0.071 0.50 -0.978E+00 40. 0.136E+01 0.290E-01 0.906E-03 0.178E-03 60 5 3 11 13 0 14 14 16 0.000 0.000 0 0 40 00:00:30.32 0.7580A -2.3762 -1.1495 1.3257 0.071 0.50 -0.978E+00 40. 0.136E+01 0.296E-01 0.915E-03 0.179E-03 60 5 2 11 13 0 12 14 16 0.000 0.000 0 0 41 00:00:31.07 0.7580A -2.3762 -1.1352 1.3579 0.070 0.50 -0.978E+00 40. 0.136E+01 0.302E-01 0.923E-03 0.179E-03 60 5 2 11 13 0 13 14 16 0.000 0.000 0