[2542] | 1 | |
---|
| 2 | ****************************** -------------------------------------------- |
---|
[2608] | 3 | * PALM 4.0 Rev: 2603:26 * atmosphere - run without 1D - prerun |
---|
[2542] | 4 | ****************************** -------------------------------------------- |
---|
| 5 | |
---|
[2608] | 6 | Date: 13-11-17 Run: ex_lsm_clearsky |
---|
| 7 | Time: 14:31:43 Run-No.: 00 |
---|
| 8 | Run on host: lcmuk |
---|
[2556] | 9 | Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated |
---|
[2542] | 10 | ------------------------------------------------------------------------------ |
---|
| 11 | |
---|
| 12 | Numerical Schemes: |
---|
| 13 | ----------------- |
---|
| 14 | |
---|
| 15 | --> Use the boussinesq approximation for the model equations. |
---|
| 16 | --> Solve perturbation pressure via FFT using temperton-algorithm routines |
---|
| 17 | perturbation pressure is calculated at every Runge-Kutta step |
---|
| 18 | --> Momentum advection via Wicker-Skamarock-Scheme 5th order |
---|
| 19 | --> Scalar advection via Wicker-Skamarock-Scheme 5th order |
---|
| 20 | --> Loop optimization method: cache |
---|
| 21 | --> Time differencing scheme: runge-kutta-3 |
---|
| 22 | --> Additional prognostic equation for the specific humidity |
---|
| 23 | ------------------------------------------------------------------------------ |
---|
| 24 | |
---|
| 25 | |
---|
| 26 | Run time and time step information: |
---|
| 27 | ---------------------------------- |
---|
| 28 | |
---|
| 29 | Timestep: variable maximum value: 20.000 s CFL-factor: 0.90 |
---|
| 30 | Start time: 0.000 s |
---|
[2556] | 31 | End time: 120.000 s |
---|
[2542] | 32 | |
---|
| 33 | |
---|
| 34 | Computational grid and domain size: |
---|
| 35 | ---------------------------------- |
---|
| 36 | |
---|
| 37 | Grid length: dx = 50.000 m dy = 50.000 m dz = 50.000 m |
---|
| 38 | Domain size: x = 2000.000 m y = 2000.000 m z(u) = 2862.511 m |
---|
| 39 | |
---|
| 40 | dz constant up to 1200.000 m (k= 25), then stretched by factor: 1.080 |
---|
| 41 | maximum dz not to be exceeded is dz_max = ********** m |
---|
| 42 | |
---|
| 43 | Number of gridpoints (x,y,z): (0: 39, 0: 39, 0: 41) |
---|
[2556] | 44 | Subdomain size (x,y,z): ( 20, 20, 42) |
---|
[2542] | 45 | |
---|
| 46 | |
---|
| 47 | Characteristic levels of the geo. wind component ug: |
---|
| 48 | |
---|
| 49 | Height: 0.0 0.0 m |
---|
[2556] | 50 | ug: 5.00 5.00 m/s |
---|
[2542] | 51 | Gradient: ------ 0.00 1/100s |
---|
| 52 | Gridpoint: 0 0 |
---|
| 53 | |
---|
| 54 | Characteristic levels of the geo. wind component vg: |
---|
| 55 | |
---|
| 56 | Height: 0.0 0.0 m |
---|
| 57 | vg: 0.00 0.00 m/s |
---|
| 58 | Gradient: ------ 0.00 1/100s |
---|
| 59 | Gridpoint: 0 0 |
---|
| 60 | |
---|
| 61 | |
---|
| 62 | Topography information: |
---|
| 63 | ---------------------- |
---|
| 64 | |
---|
| 65 | Topography: flat |
---|
| 66 | |
---|
| 67 | |
---|
| 68 | Land surface model information: |
---|
| 69 | ------------------------------ |
---|
| 70 | |
---|
| 71 | --> Soil bottom is closed (water content is conserved, default) |
---|
| 72 | --> Land surface type : crops, mixed farming |
---|
| 73 | --> Soil porosity type : medium-fine |
---|
| 74 | |
---|
| 75 | Initial soil temperature and moisture profile: |
---|
| 76 | |
---|
| 77 | Height: -0.01 -0.02 -0.05 -0.11 -0.22 -0.40 -0.86 -2. m |
---|
| 78 | Temperature: 293.00 293.00 293.00 293.00 293.00 293.00 293.00 293. K |
---|
| 79 | Moisture: 0.35 0.35 0.35 0.35 0.35 0.35 0.35 0. m**3/m**3 |
---|
| 80 | Root fraction: 0.03 0.07 0.14 0.14 0.28 0.09 0.22 0. |
---|
| 81 | Grid point: 0 1 2 3 4 5 6 |
---|
| 82 | |
---|
| 83 | |
---|
| 84 | Radiation model information: |
---|
| 85 | ---------------------------- |
---|
| 86 | |
---|
| 87 | --> Simple radiation scheme for clear sky is used (no clouds, default) |
---|
| 88 | |
---|
| 89 | Albedo is set for land surface type: mixed farming, tall grassland |
---|
| 90 | Timestep: dt_radiation = 60.00 s |
---|
| 91 | |
---|
| 92 | |
---|
| 93 | Boundary conditions: |
---|
| 94 | ------------------- |
---|
| 95 | |
---|
| 96 | p uv pt |
---|
| 97 | |
---|
| 98 | B. bound.: p(0) = p(1) | uv(0) = -uv(1) | pt(0) = from soil model |
---|
| 99 | T. bound.: p(nzt+1) = 0 | uv(nzt+1) = ug(nzt+1), vg(nzt+1) | pt(nzt+1) = pt(nzt) + dpt/dz_ |
---|
| 100 | |
---|
| 101 | e |
---|
| 102 | |
---|
| 103 | B. bound.: e(0) = e(1) |
---|
| 104 | T. bound.: e(nzt+1) = e(nzt) = e(nzt-1) |
---|
| 105 | |
---|
| 106 | q |
---|
| 107 | |
---|
| 108 | B. bound.: q(0) = from soil model |
---|
| 109 | T. bound.: q(nzt+1) = q(nzt) + dq/dz |
---|
| 110 | |
---|
| 111 | Bottom surface fluxes are used in diffusion terms at k=1 |
---|
| 112 | |
---|
| 113 | |
---|
| 114 | Constant flux layer between bottom surface and first computational u,v-level: |
---|
| 115 | |
---|
| 116 | z_mo = 25.00 m z0 = 0.1000 m z0h = 0.10000 m kappa = 0.40 |
---|
| 117 | Rif value range: -20.00 <= rif <= 20.00 |
---|
| 118 | Predefined surface temperature |
---|
| 119 | Predefined surface humidity |
---|
| 120 | |
---|
| 121 | |
---|
| 122 | Lateral boundaries: |
---|
| 123 | left/right: cyclic |
---|
| 124 | north/south: cyclic |
---|
| 125 | |
---|
| 126 | |
---|
| 127 | Initial profiles: |
---|
| 128 | ---------------- |
---|
| 129 | |
---|
| 130 | Characteristic levels of the initial temperature profile: |
---|
| 131 | |
---|
| 132 | Height: 0.0 0.0 800.0 m |
---|
| 133 | Temperature: 293.00 293.00 293.00 K |
---|
| 134 | Gradient: ------ 0.00 1.00 K/100m |
---|
| 135 | Gridpoint: 0 0 16 |
---|
| 136 | |
---|
| 137 | Characteristic levels of the initial humidity profile: |
---|
| 138 | |
---|
| 139 | Height: 0.0 0.0 m |
---|
| 140 | Humidity: 0.1E-01 0.1E-01 kg/kg |
---|
| 141 | Gradient: -------- 0.0E+00 (kg/kg)/100m |
---|
| 142 | Gridpoint: 0 0 |
---|
| 143 | |
---|
| 144 | |
---|
| 145 | List output: |
---|
| 146 | ----------- |
---|
| 147 | |
---|
| 148 | 1D-Profiles: |
---|
| 149 | Output every ******** s |
---|
| 150 | Time averaged over 900.00 s |
---|
| 151 | Averaging input every 0.00 s |
---|
| 152 | |
---|
| 153 | |
---|
| 154 | Data output: |
---|
| 155 | ----------- |
---|
| 156 | |
---|
| 157 | Time averaged over 900.00 s |
---|
| 158 | Averaging input every 0.00 s |
---|
| 159 | |
---|
| 160 | 1D-Profiles: |
---|
| 161 | Output format: netCDF 64bit offset |
---|
| 162 | |
---|
| 163 | |
---|
| 164 | Profile: u, v, pt, km, kh, l, w, e, e*, p, w"pt", w*pt*, wpt, |
---|
| 165 | |
---|
| 166 | : w*e*, pt*2, u*2, v*2, w*2, t_soil, |
---|
| 167 | Output every 900.00 s |
---|
| 168 | Time averaged over 900.00 s |
---|
| 169 | Averaging input every 0.00 s |
---|
| 170 | |
---|
| 171 | 2D-Arrays: |
---|
| 172 | Output format: netCDF 64bit offset |
---|
| 173 | |
---|
| 174 | |
---|
| 175 | XY-cross-section Arrays: u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l |
---|
| 176 | Output every 3600.00 s |
---|
| 177 | Cross sections at k = /0/ |
---|
| 178 | scalar-coordinates: /0.0/ m |
---|
| 179 | |
---|
| 180 | |
---|
| 181 | Number of output time levels allowed: unlimited |
---|
| 182 | |
---|
| 183 | |
---|
| 184 | 2D-Arrays(time-averaged): |
---|
| 185 | Output format: netCDF 64bit offset |
---|
| 186 | |
---|
| 187 | |
---|
| 188 | XY-cross-section Arrays: u*, t*, r_a*, r_s*, ghf*, shf*, qsws_liq*, qsws_soil*, qsws_veg*, c_l |
---|
| 189 | Output every 3600.00 s |
---|
| 190 | Time averaged over 1800.00 s |
---|
| 191 | Averaging input every 0.00 s |
---|
| 192 | Cross sections at k = /0/ |
---|
| 193 | scalar-coordinates: /0.0/ m |
---|
| 194 | |
---|
| 195 | |
---|
| 196 | Number of output time levels allowed: unlimited |
---|
| 197 | |
---|
| 198 | |
---|
| 199 | 3D-Arrays: |
---|
| 200 | Output format: netCDF 64bit offset |
---|
| 201 | |
---|
| 202 | |
---|
| 203 | Arrays: pt, vpt, q, u, v, w, t_soil, m_soil, |
---|
| 204 | Output every 3600.00 s |
---|
| 205 | Upper output limit at 2862.51 m (GP 41) |
---|
| 206 | |
---|
| 207 | |
---|
| 208 | Number of output time levels allowed: unlimited |
---|
| 209 | |
---|
| 210 | |
---|
| 211 | 3D-Arrays(time-averaged): |
---|
| 212 | Output format: netCDF 64bit offset |
---|
| 213 | |
---|
| 214 | |
---|
| 215 | Arrays: t_soil, m_soil, |
---|
| 216 | Output every 3600.00 s |
---|
| 217 | Time averaged over 1800.00 s |
---|
| 218 | Averaging input every 0.00 s |
---|
| 219 | Upper output limit at 2862.51 m (GP 41) |
---|
| 220 | |
---|
| 221 | |
---|
| 222 | Number of output time levels allowed: unlimited |
---|
| 223 | |
---|
| 224 | |
---|
| 225 | Time series: |
---|
| 226 | Output format: netCDF 64bit offset |
---|
| 227 | |
---|
| 228 | Output every 0.00 s |
---|
| 229 | |
---|
| 230 | ------------------------------------------------------------------------------ |
---|
| 231 | |
---|
| 232 | |
---|
| 233 | Physical quantities: |
---|
| 234 | ------------------- |
---|
| 235 | |
---|
[2608] | 236 | Geograph. latitude : latitude = 55.0 degr |
---|
| 237 | Geograph. longitude : longitude = 0.0 degr |
---|
[2542] | 238 | Angular velocity : omega = 0.729E-04 rad/s |
---|
| 239 | Coriolis parameter : f = 0.000119 1/s |
---|
| 240 | f* = 0.000084 1/s |
---|
[2608] | 241 | Day of the year at model start : day_init = 172 |
---|
| 242 | UTC time at model start : time_utc_init = 43200.0 s |
---|
[2542] | 243 | |
---|
| 244 | Gravity : g = 9.8 m/s**2 |
---|
| 245 | |
---|
[2543] | 246 | Reference state used in buoyancy terms: initial_profile |
---|
[2542] | 247 | |
---|
| 248 | |
---|
| 249 | Cloud physics quantities / methods: |
---|
| 250 | ---------------------------------- |
---|
| 251 | |
---|
| 252 | Humidity is considered, bu no condensation |
---|
| 253 | |
---|
| 254 | |
---|
| 255 | LES / Turbulence quantities: |
---|
| 256 | --------------------------- |
---|
| 257 | |
---|
| 258 | Mixing length is limited to 1.80 * z |
---|
| 259 | |
---|
| 260 | |
---|
| 261 | Actions during the simulation: |
---|
| 262 | ----------------------------- |
---|
| 263 | |
---|
[2556] | 264 | Disturbance impulse (u,v) every : ****** s |
---|
[2542] | 265 | Disturbance amplitude : 0.25 m/s |
---|
| 266 | Lower disturbance level : 125.00 m (GP 3) |
---|
| 267 | Upper disturbance level : 625.00 m (GP 13) |
---|
| 268 | Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2 |
---|
| 269 | Random number generator used : random-parallel |
---|
| 270 | |
---|
| 271 | |
---|
| 272 | |
---|
| 273 | *** no user-defined variables found |
---|
| 274 | |
---|
| 275 | ------------------------------------------------------------------------------ |
---|
| 276 | |
---|
| 277 | |
---|
| 278 | |
---|
| 279 | Run-control output: |
---|
| 280 | ------------------ |
---|
| 281 | |
---|
| 282 | 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 |
---|
| 283 | --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- |
---|
[2556] | 284 | 0 0 00:00:00.00 8.6700A 5.1898D -0.2288D 0.1182 0.100 0.00 0.200E-11 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 |
---|
| 285 | 0 1 00:00:08.67 8.6900A 5.1806 0.2159 0.1162 0.001 0.00 0.000E+00 550. 0.119E+02 0.867E-03 0.287E-05 0.440E-15 9 9 24 4 5 26 5 26 8 0.000 0.000 0 |
---|
[2608] | 286 | 0 2 00:00:17.35 8.6900A 5.1762 0.1915 0.1151 0.021 1.06 -0.495E+02 100. 0.119E+02 0.793E-03 0.744E-05 0.407E-15 5 16 32 4 5 27 5 26 9 0.000 0.000 0 |
---|
| 287 | 0 3 00:00:26.04 8.7000A 5.1723 0.1786 0.1113 0.021 0.70 -0.516E+01 100. 0.119E+02 0.751E-03 0.247E-05 0.385E-15 5 16 33 4 5 28 5 26 10 0.000 0.000 0 |
---|
| 288 | 0 4 00:00:34.75 8.7000A 5.1732 -0.1675 0.1061 0.027 0.69 -0.360E+01 100. 0.119E+02 0.721E-03 0.396E-05 0.358E-15 5 16 34 8 29 32 5 26 11 0.000 0.000 0 |
---|
| 289 | 0 5 00:00:43.45 8.7000A 5.1719 -0.1671 0.1028 0.034 0.93 -0.287E+01 250. 0.119E+02 0.698E-03 0.471E-05 0.337E-15 5 16 35 8 29 33 5 26 11 0.000 0.000 0 |
---|
| 290 | 0 6 00:00:52.15 8.7000A 5.1701 -0.1621 0.1039 0.039 0.93 -0.242E+01 250. 0.119E+02 0.680E-03 0.548E-05 0.321E-15 5 16 36 8 29 34 5 26 12 0.000 0.000 0 |
---|
| 291 | 0 7 00:01:00.85 8.7100A 5.1675 -0.1565 0.1043 0.045 1.21 -0.211E+01 550. 0.119E+02 0.665E-03 0.627E-05 0.306E-15 5 16 37 8 29 35 5 26 13 0.000 0.000 0 |
---|
| 292 | 0 8 00:01:09.56 8.7100A 5.1643 0.1515 0.1040 0.051 1.21 -0.189E+01 550. 0.119E+02 0.652E-03 0.709E-05 0.294E-15 5 16 38 11 3 22 5 26 14 0.000 0.000 0 |
---|
| 293 | 0 9 00:01:18.27 8.7100A 5.1636 -0.1506 0.1034 0.056 1.20 -0.166E+01 550. 0.119E+02 0.640E-03 0.785E-05 0.282E-15 5 16 38 8 29 36 5 26 15 0.000 0.000 0 |
---|
| 294 | 0 10 00:01:26.98 8.7100A 5.1635 -0.1516 0.1024 0.061 1.20 -0.153E+01 550. 0.119E+02 0.630E-03 0.868E-05 0.269E-15 5 16 39 8 29 37 5 26 16 0.000 0.000 0 |
---|
| 295 | 0 11 00:01:35.69 8.7200A 5.1632 -0.1516 0.1013 0.067 1.20 -0.143E+01 550. 0.119E+02 0.620E-03 0.949E-05 0.260E-15 5 16 0 8 29 38 5 26 17 0.000 0.000 0 |
---|
| 296 | 0 12 00:01:44.41 8.7200A 5.1624 -0.1506 0.0999 0.072 1.21 -0.135E+01 550. 0.119E+02 0.612E-03 0.103E-04 0.249E-15 5 16 1 8 29 39 5 26 18 0.000 0.000 0 |
---|
| 297 | 0 13 00:01:53.13 8.7200A 5.1614 -0.1488 0.0983 0.077 1.22 -0.128E+01 550. 0.119E+02 0.604E-03 0.111E-04 0.239E-15 5 16 2 8 29 0 5 26 19 0.000 0.000 0 |
---|
| 298 | 0 14 00:02:01.85 8.7200A 5.1600 -0.1462 0.0966 0.081 1.22 -0.122E+01 550. 0.119E+02 0.596E-03 0.119E-04 0.230E-15 5 16 3 8 29 1 5 26 20 0.000 0.000 0 |
---|