****************************** -------------------------------------------- * PALM 6.0 Rev: 3915M * atmosphere - run without 1D - prerun ****************************** -------------------------------------------- Date: 2019-09-26 Run: urban_environment_salsa__gfortran_ Time: 10:42:08 Run-No.: 00 Run on host: gfortran_d 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 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: 1270.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) = 45.000 m Number of gridpoints (x,y,z): (0: 19, 0: 19, 0: 23) Subdomain size (x,y,z): ( 10, 10, 24) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 1.00 1.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 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 m Temperature: 277.15 277.15 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.1E-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 60.00 s Averaging input every 0.00 s Data output: ----------- Time averaged over 60.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, salsa_Ntot, : salsa_N_UFP, salsa_LDSA, salsa_PM0.1, salsa_PM2.5, salsa_PM10 : Output every 60.00 s Time averaged over 60.00 s Averaging input every 0.00 s Time series: Output format: netCDF 64bit offset Output every 10.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 to 1.80 * z 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 : 13.00 m (GP 7) 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_3d 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: ---------------------------- --> Simple radiation scheme for clear sky is used (no clouds, default) 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: 22 19 19 8 aerosol_mass: 24 16 16 32 (advect_particle_water = T) salsa_gas: 24 16 16 5 (salsa_gases_from_chem = F) Aerosol dynamic processes included: coagulation condensation (of precursor gases = T and water vapour = F) 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 limits (in metres): 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 0.00E+000 Initial number concentration in bins at the lowest level (#/m**3): 1.19E+008 9.63E+008 5.30E+008 1.96E+008 6.13E+007 6.01E+006 1.95E+005 1.25E+004 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 = 0 Gas concentrations: init_gases_type = 0 Mode diametres: dpg(nmod) = 0.000 0.000 0.000 0.000 0.000 0.000 0.000 (m) Standard deviation: sigmag(nmod) = 1.80 2.16 2.21 2.00 2.00 2.00 2.00 Number concentration: n_lognorm(nmod) = 1.4300E+009 4.4500E+008 7.0000E+004 0.0000E+000 0.0000E+000 0.0000E+000 0.0000E+000 (#/m3) Emissions: salsa_emission_mode = read_from_file (currently all emissions are soluble!) ------------------------------------------------------------------------------ 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 1.3300A 1.3569D 0.6524D 0.4803 0.070 0.00 0.861E-01 0. 0.366E+00 0.309E-01 0.378E-01 0.160E-01 16 0 18 6 15 5 20 17 14 0.000 0.000 0 0 1 00:00:01.32 1.0500A 1.7173 -0.9917 0.9157 0.099 0.00 -0.442E+00 2. 0.349E+00 0.513E-01 0.144E-01 0.399E-02 16 0 0 6 15 14 15 1 14 0.000 0.000 0 0 2 00:00:02.37 1.0700A 1.6851 -1.0410 1.0211 0.097 0.00 -0.452E+00 2. 0.345E+00 0.558E-01 0.602E-02 0.189E-02 16 0 0 6 15 14 15 1 14 0.000 0.000 0 0 3 00:00:03.44 1.0900A 1.6547 1.0544 1.0394 0.095 0.28 -0.458E+00 30. 0.341E+00 0.574E-01 0.533E-02 0.156E-02 16 0 1 6 5 14 15 1 14 0.000 0.000 0 0 4 00:00:04.53 1.1200A 1.6140 1.0493 1.0340 0.093 0.00 -0.472E+00 8. 0.338E+00 0.585E-01 0.515E-02 0.147E-02 16 0 1 6 5 14 15 1 14 0.000 0.000 0 0 5 00:00:05.65 1.1400A 1.5788 1.0336 1.0182 0.091 0.00 -0.495E+00 4. 0.335E+00 0.595E-01 0.514E-02 0.145E-02 16 0 2 6 5 14 15 1 14 0.000 0.000 0 0 6 00:00:06.79 1.1900A 1.5169 1.0134 0.9996 0.089 0.00 -0.525E+00 4. 0.332E+00 0.606E-01 0.512E-02 0.144E-02 16 0 2 6 5 14 15 1 14 0.000 0.000 0 0 7 00:00:07.98 1.2200A 1.4729 0.9879 0.9791 0.088 0.00 -0.546E+00 4. 0.329E+00 0.616E-01 0.522E-02 0.145E-02 16 1 15 6 5 14 15 1 14 0.000 0.000 0 0 8 00:00:09.20 1.2400A 1.4512 0.9609 0.9619 0.087 0.00 -0.537E+00 4. 0.326E+00 0.628E-01 0.527E-02 0.146E-02 16 1 15 6 5 14 15 2 14 0.000 0.000 0 0 9 00:00:10.44 1.2600A 1.4301 0.9376 0.9485 0.085 0.00 -0.507E+00 4. 0.323E+00 0.640E-01 0.529E-02 0.146E-02 16 1 15 10 5 14 15 2 14 0.000 0.000 0 0 10 00:00:11.70 1.2800A 1.4106 0.9225 0.9363 0.084 0.00 -0.494E+00 4. 0.321E+00 0.652E-01 0.534E-02 0.146E-02 16 1 15 10 5 14 15 2 14 0.000 0.000 0 0 11 00:00:12.98 1.2900A 1.3930 0.9100 0.9248 0.083 0.00 -0.491E+00 4. 0.318E+00 0.664E-01 0.542E-02 0.146E-02 16 1 15 3 5 14 15 2 14 0.000 0.000 0 0 12 00:00:14.27 1.3100A 1.3776 0.8963 0.9141 0.083 0.00 -0.484E+00 4. 0.316E+00 0.677E-01 0.545E-02 0.145E-02 16 1 15 3 5 14 15 2 14 0.000 0.000 0 0 13 00:00:15.58 1.3200A 1.3623 0.8811 0.9030 0.082 0.00 -0.484E+00 4. 0.313E+00 0.690E-01 0.550E-02 0.146E-02 16 1 15 8 5 14 15 2 14 0.000 0.000 0 0 14 00:00:16.90 1.3400A 1.3479 0.8688 0.8921 0.081 0.00 -0.480E+00 4. 0.311E+00 0.704E-01 0.552E-02 0.145E-02 16 1 15 8 5 14 15 2 14 0.000 0.000 0 0 15 00:00:18.24 1.3500A 1.3331 0.8550 0.8804 0.080 0.00 -0.476E+00 4. 0.309E+00 0.717E-01 0.556E-02 0.145E-02 16 1 15 8 5 14 15 2 14 0.000 0.000 0 0 16 00:00:19.59 1.3600A 1.3189 -0.8433 0.8687 0.079 0.00 -0.475E+00 4. 0.306E+00 0.731E-01 0.556E-02 0.143E-02 16 1 15 7 3 5 15 2 14 0.000 0.000 0 0 17 00:00:20.95 1.3800A 1.3050 -0.8519 0.8567 0.079 0.00 -0.474E+00 4. 0.304E+00 0.744E-01 0.554E-02 0.141E-02 16 1 15 7 3 5 15 2 14 0.000 0.000 0 0 18 00:00:22.33 1.3900A 1.2904 -0.8451 0.8441 0.078 0.00 -0.477E+00 4. 0.302E+00 0.758E-01 0.555E-02 0.140E-02 16 1 15 7 3 5 15 2 14 0.000 0.000 0 0 19 00:00:23.72 1.4100A 1.2763 -0.8250 0.8317 0.077 0.00 -0.477E+00 4. 0.300E+00 0.772E-01 0.553E-02 0.137E-02 16 1 15 7 3 6 15 2 14 0.000 0.000 0 0 20 00:00:25.13 1.4300A 1.2618 -0.8139 0.8189 0.077 0.00 -0.481E+00 4. 0.298E+00 0.785E-01 0.552E-02 0.135E-02 16 1 15 7 3 6 15 2 14 0.000 0.000 0 0 21 00:00:26.56 1.4400A 1.2475 0.7696 0.8062 0.076 0.00 -0.486E+00 4. 0.296E+00 0.799E-01 0.551E-02 0.133E-02 16 1 15 11 5 14 15 2 14 0.000 0.000 0 0 22 00:00:28.00 1.4500A 1.2416 0.7569 0.7939 0.076 0.00 -0.479E+00 4. 0.294E+00 0.813E-01 0.546E-02 0.129E-02 11 7 2 12 5 14 15 2 14 0.000 0.000 0 0 23 00:00:29.45 1.4500A 1.2384 0.7470 0.7814 0.076 0.00 -0.471E+00 4. 0.293E+00 0.827E-01 0.540E-02 0.126E-02 10 7 3 12 5 14 15 2 14 0.000 0.000 0 0 24 00:00:30.90 1.4600A 1.2345 0.7384 0.7685 0.075 0.00 -0.470E+00 6. 0.291E+00 0.841E-01 0.530E-02 0.122E-02 11 7 3 12 5 14 15 2 14 0.000 0.000 0 0 25 00:00:32.36 1.4600A 1.2301 0.7299 0.7760 0.075 0.00 -0.470E+00 6. 0.290E+00 0.855E-01 0.523E-02 0.118E-02 11 7 3 12 5 14 10 14 4 0.000 0.000 0 0 26 00:00:33.82 1.4600A 1.2308 0.7216 0.8160 0.075 0.00 -0.472E+00 6. 0.288E+00 0.868E-01 0.513E-02 0.113E-02 18 2 5 12 5 14 6 2 5 0.000 0.000 0 0 27 00:00:35.28 1.4500A 1.2382 0.7121 0.8576 0.074 0.00 -0.470E+00 6. 0.287E+00 0.882E-01 0.503E-02 0.109E-02 18 2 5 12 5 14 10 14 4 0.000 0.000 0 0 28 00:00:36.74 1.4500A 1.2442 0.7006 0.9018 0.074 0.00 -0.471E+00 6. 0.285E+00 0.896E-01 0.492E-02 0.104E-02 18 2 5 12 5 14 10 14 4 0.000 0.000 0 0 29 00:00:38.19 1.4400A 1.2487 0.6854 0.9549 0.074 0.00 -0.470E+00 6. 0.284E+00 0.909E-01 0.484E-02 0.100E-02 18 2 5 12 5 14 11 14 4 0.000 0.000 0 0 30 00:00:39.63 1.4400A 1.2512 0.6695 1.0069 0.074 0.00 -0.470E+00 6. 0.283E+00 0.923E-01 0.473E-02 0.955E-03 18 2 5 11 5 14 11 14 4 0.000 0.000 0 0 31 00:00:41.07 1.4400A 1.2532 0.6498 1.0521 0.074 0.00 -0.471E+00 6. 0.282E+00 0.937E-01 0.466E-02 0.916E-03 18 2 6 11 5 14 11 14 4 0.000 0.000 0 0 32 00:00:42.51 1.4300A 1.2566 0.6250 1.0928 0.073 0.00 -0.471E+00 6. 0.281E+00 0.950E-01 0.458E-02 0.878E-03 18 2 6 11 5 14 12 14 4 0.000 0.000 0 0 33 00:00:43.94 1.4400A 1.2540 -0.6223 1.1443 0.073 0.00 -0.472E+00 6. 0.281E+00 0.964E-01 0.448E-02 0.836E-03 18 2 6 6 4 9 12 14 4 0.000 0.000 0 0 34 00:00:45.38 1.4500A 1.2448 -0.6272 1.1867 0.073 0.00 -0.474E+00 6. 0.280E+00 0.977E-01 0.443E-02 0.800E-03 18 2 6 6 3 10 11 14 3 0.000 0.000 0 0 35 00:00:46.83 1.4600A 1.2295 -0.6604 1.2293 0.073 0.00 -0.478E+00 6. 0.279E+00 0.991E-01 0.440E-02 0.762E-03 18 2 6 6 3 10 13 14 4 0.000 0.000 0 0 36 00:00:48.29 1.4200A 1.2284 -0.6898 1.2691 0.073 0.00 -0.491E+00 6. 0.279E+00 0.100E+00 0.437E-02 0.722E-03 18 3 6 6 3 10 13 14 4 0.000 0.000 0 0 37 00:00:49.71 1.3900A 1.2317 -0.7102 1.2976 0.073 0.00 -0.483E+00 6. 0.278E+00 0.102E+00 0.421E-02 0.664E-03 18 3 6 6 3 10 14 14 4 0.000 0.000 0 0 38 00:00:51.10 1.3400A 1.2352 -0.7193 1.3406 0.072 0.00 -0.484E+00 6. 0.278E+00 0.103E+00 0.409E-02 0.610E-03 18 3 7 5 4 10 14 14 4 0.000 0.000 0 0 39 00:00:52.44 1.3300A 1.2408 -0.7301 1.3578 0.072 0.00 -0.486E+00 10. 0.278E+00 0.104E+00 0.394E-02 0.555E-03 18 3 7 5 4 10 14 14 4 0.000 0.000 0 0 40 00:00:53.77 1.3100A 1.2409 -0.7383 1.3752 0.072 0.00 -0.489E+00 10. 0.277E+00 0.106E+00 0.390E-02 0.518E-03 18 3 7 5 4 10 15 14 4 0.000 0.000 0 0 41 00:00:55.08 1.2800A 1.2351 -0.7513 1.4071 0.072 0.00 -0.491E+00 10. 0.277E+00 0.107E+00 0.384E-02 0.482E-03 18 3 7 6 3 11 15 14 4 0.000 0.000 0 0 42 00:00:56.36 1.2400A 1.2236 -0.7700 1.4495 0.072 0.00 -0.495E+00 10. 0.277E+00 0.108E+00 0.376E-02 0.447E-03 18 3 7 6 3 11 13 14 4 0.000 0.000 0 0 43 00:00:57.60 1.2000A 1.2192 -0.8028 1.5040 0.072 0.00 -0.499E+00 10. 0.277E+00 0.109E+00 0.367E-02 0.412E-03 18 3 8 5 4 11 13 14 4 0.000 0.000 0 0 44 00:00:58.80 1.1700A 1.2093 -0.8350 1.5334 0.072 0.00 -0.504E+00 10. 0.277E+00 0.110E+00 0.357E-02 0.381E-03 18 3 8 5 4 11 13 14 4 0.000 0.000 0 0 45 00:00:59.97 1.1600A 1.1987 -0.8605 1.5491 0.071 0.00 -0.506E+00 10. 0.277E+00 0.111E+00 0.350E-02 0.357E-03 19 2 8 5 4 11 14 14 4 0.000 0.000 0 0 46 00:01:01.13 1.1700A 1.1976 -0.8763 1.5334 0.071 0.00 -0.506E+00 10. 0.276E+00 0.111E+00 0.349E-02 0.341E-03 6 6 10 5 4 11 14 14 4 0.000 0.000 0 0 47 00:01:02.30 1.2200A 1.2002 -0.8796 1.4801 0.071 0.00 -0.505E+00 10. 0.276E+00 0.112E+00 0.355E-02 0.333E-03 6 6 10 5 4 11 14 14 4 0.000 0.000 0 0 48 00:01:03.52 1.2800A 1.1985 -0.8848 1.4105 0.071 0.00 -0.506E+00 10. 0.276E+00 0.113E+00 0.374E-02 0.340E-03 6 6 10 4 4 11 11 14 5 0.000 0.000 0 0 49 00:01:04.80 1.2100A 1.1920 -0.8819 1.4918 0.071 0.00 -0.507E+00 10. 0.276E+00 0.114E+00 0.400E-02 0.356E-03 6 6 10 4 4 11 12 14 5 0.000 0.000 0 0 50 00:01:06.01 1.1600A 1.1878 -0.8670 1.5556 0.071 0.00 -0.506E+00 10. 0.276E+00 0.114E+00 0.389E-02 0.348E-03 13 11 9 4 4 11 13 14 5 0.000 0.000 0 0 51 00:01:07.17 1.1400A 1.1898 -0.8896 1.5802 0.071 0.07 -0.505E+00 10. 0.276E+00 0.115E+00 0.382E-02 0.342E-03 13 11 10 5 3 12 14 14 5 0.000 0.000 0 0 52 00:01:08.31 1.1400A 1.1965 -0.9231 1.5795 0.071 0.09 -0.505E+00 10. 0.275E+00 0.115E+00 0.383E-02 0.342E-03 13 11 10 5 3 12 14 14 5 0.000 0.000 0 0 53 00:01:09.45 1.1800A 1.2005 -0.9448 1.5249 0.071 0.10 -0.504E+00 10. 0.275E+00 0.116E+00 0.391E-02 0.350E-03 13 11 10 4 4 12 15 14 5 0.000 0.000 0 0 54 00:01:10.63 1.2500A 1.2111 -0.9540 1.4388 0.071 0.10 -0.504E+00 10. 0.275E+00 0.116E+00 0.413E-02 0.370E-03 13 12 11 4 4 12 15 14 5 0.000 0.000 0