- Timestamp:
- Jun 12, 2018 7:03:02 AM (6 years ago)
- Location:
- palm/trunk
- Files:
-
- 22 edited
-
INSTALL/example_cbl_rc (modified) (9 diffs)
-
SOURCE/check_parameters.f90 (modified) (3 diffs)
-
SOURCE/header.f90 (modified) (3 diffs)
-
SOURCE/init_grid.f90 (modified) (9 diffs)
-
SOURCE/init_masks.f90 (modified) (5 diffs)
-
SOURCE/lpm_advec.f90 (modified) (11 diffs)
-
SOURCE/lpm_exchange_horiz.f90 (modified) (2 diffs)
-
SOURCE/lpm_init.f90 (modified) (4 diffs)
-
SOURCE/lpm_set_attributes.f90 (modified) (5 diffs)
-
SOURCE/modules.f90 (modified) (7 diffs)
-
SOURCE/parin.f90 (modified) (3 diffs)
-
SOURCE/plant_canopy_model_mod.f90 (modified) (5 diffs)
-
SOURCE/pmc_interface_mod.f90 (modified) (3 diffs)
-
SOURCE/radiation_model_mod.f90 (modified) (14 diffs)
-
SOURCE/read_restart_data_mod.f90 (modified) (3 diffs)
-
SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (4 diffs)
-
SOURCE/urban_surface_mod.f90 (modified) (6 diffs)
-
SOURCE/user_init_grid.f90 (modified) (2 diffs)
-
SOURCE/vertical_nesting_mod.f90 (modified) (11 diffs)
-
SOURCE/virtual_flight_mod.f90 (modified) (3 diffs)
-
SOURCE/wind_turbine_model_mod.f90 (modified) (22 diffs)
-
SOURCE/write_restart_data_mod.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/INSTALL/example_cbl_rc
r2696 r3065 1 1 2 2 ****************************** -------------------------------------------- 3 * PALM 4.0 Rev: 2683* atmosphere - run without 1D - prerun3 * PALM 5.0 Rev: 3058M * atmosphere - run without 1D - prerun 4 4 ****************************** -------------------------------------------- 5 5 6 Date: 12-12-17 Run: example_cbl 7 Time: 10:31:35 Run-No.: 00 8 Run on host: imuk 9 Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated 6 Date: 11-06-18 Run: example_cbl 7 Time: 15:43:48 Run-No.: 00 8 Run on host: imuk_defau 9 Number of PEs: 1 Processor grid (x,y): ( 1, 1) calculated 10 A 1d-decomposition along x is used 10 11 ------------------------------------------------------------------------------ 11 12 … … 35 36 ---------------------------------- 36 37 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 38 Grid length: dx = 50.000 m dy = 50.000 m 39 dz(1) = 50.000 m 40 dz(2) = 1000.000 m 41 42 Domain size: x = 2000.000 m y = 2000.000 m z(u) = 2997.512 m 43 44 Vertical stretching starts at height: 1175.0 m, 45 Vertical stretching starts at index: 24, 46 Vertical stretching ends at height: 9999999.9 m, 47 Vertical stretching ends at index: 41, 48 Factor used for stretching: 1.080, 42 49 43 50 Number of gridpoints (x,y,z): (0: 39, 0: 39, 0: 41) 44 Subdomain size (x,y,z): ( 20, 20, 42)51 Subdomain size (x,y,z): ( 40, 40, 42) 45 52 46 53 … … 109 116 110 117 1D-Profiles: 111 Output every ********s118 Output every 9999999.90 s 112 119 Time averaged over 600.00 s 113 120 Averaging input every 10.00 s … … 227 234 0 1 00:00:20.00 20.0000D -0.2125 -0.2273 0.1172 0.009 0.69 -0.132E+02 100. 0.103E-02 0.102E-02 0.570E-05 0.411E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 228 235 0 2 00:00:40.00 20.0000D -0.2117 -0.2252 0.1162 0.009 0.69 -0.133E+02 100. 0.102E-02 0.102E-02 0.935E-05 0.416E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 229 0 3 00:01:00.00 20.0000D -0.2106 -0.2227 0.1152 0.009 0.69 -0.134E+02 100. 0.101E-02 0.101E-02 0.130E-04 0.42 1E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0236 0 3 00:01:00.00 20.0000D -0.2106 -0.2227 0.1152 0.009 0.69 -0.134E+02 100. 0.101E-02 0.101E-02 0.130E-04 0.420E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 230 237 0 4 00:01:20.00 20.0000D -0.2093 -0.2197 0.1141 0.009 0.69 -0.134E+02 100. 0.100E-02 0.997E-03 0.166E-04 0.424E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 231 238 0 5 00:01:40.00 20.0000D -0.2077 -0.2164 0.1129 0.009 0.93 -0.134E+02 250. 0.992E-03 0.988E-03 0.203E-04 0.427E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 232 239 0 6 00:02:00.00 20.0000D -0.2058 -0.2126 0.1117 0.009 0.93 -0.133E+02 250. 0.983E-03 0.979E-03 0.239E-04 0.429E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 233 0 7 00:02:20.00 20.0000D -0.2037 -0.2086 0.1103 0.009 1.22 -0.131E+02 550. 0.974E-03 0.970E-03 0.275E-04 0.43 1E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0240 0 7 00:02:20.00 20.0000D -0.2037 -0.2086 0.1103 0.009 1.22 -0.131E+02 550. 0.974E-03 0.970E-03 0.275E-04 0.430E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 234 241 0 8 00:02:40.00 20.0000D 0.3348D 0.3126D 0.1636 0.011 1.22 -0.105E+02 550. 0.202E-02 0.201E-02 0.601E-03 0.856E-13 10 3 39 10 8 27 5 36 28 0.000 0.000 0 235 0 9 00:03:00.00 20.0000D 0.3330 0.3039 0.1619 0.011 1.22 -0.105E+02 550. 0.199E-02 0.198E-02 0.374E-04 0.8 30E-15 10 3 39 10 8 27 5 36 28 0.000 0.000 0242 0 9 00:03:00.00 20.0000D 0.3330 0.3039 0.1619 0.011 1.22 -0.105E+02 550. 0.199E-02 0.198E-02 0.374E-04 0.829E-15 10 3 39 10 8 27 5 36 28 0.000 0.000 0 236 243 0 10 00:03:20.00 20.0000D 0.3296 -0.2989 0.1601 0.011 1.22 -0.103E+02 550. 0.196E-02 0.195E-02 0.409E-04 0.827E-15 10 3 39 8 4 6 5 36 28 0.000 0.000 0 237 244 0 11 00:03:40.00 20.0000D 0.3249 -0.2986 0.1580 0.011 1.22 -0.101E+02 550. 0.194E-02 0.193E-02 0.444E-04 0.828E-15 10 3 39 8 4 6 5 36 28 0.000 0.000 0 238 245 0 12 00:04:00.00 20.0000D 0.3189 -0.2989 0.1558 0.012 1.22 -0.975E+01 550. 0.191E-02 0.190E-02 0.478E-04 0.831E-15 10 3 39 8 0 7 5 36 28 0.000 0.000 0 239 246 0 13 00:04:20.00 20.0000D 0.3119 -0.3005 0.1535 0.012 1.22 -0.940E+01 550. 0.189E-02 0.188E-02 0.513E-04 0.835E-15 10 3 39 8 0 7 5 36 28 0.000 0.000 0 240 0 14 00:04:40.00 20.0000D 0.3056 -0.3012 0.1512 0.013 1.22 -0.905E+01 550. 0.187E-02 0.185E-02 0.546E-04 0.84 2E-15 10 3 0 8 0 7 5 36 28 0.000 0.000 0247 0 14 00:04:40.00 20.0000D 0.3056 -0.3012 0.1512 0.013 1.22 -0.905E+01 550. 0.187E-02 0.185E-02 0.546E-04 0.841E-15 10 3 0 8 0 7 5 36 28 0.000 0.000 0 241 248 0 15 00:05:00.00 20.0000D 0.4151D -0.3964D -0.1812 0.013 1.22 -0.851E+01 550. 0.292E-02 0.289E-02 0.625E-03 0.860E-13 8 17 5 8 3 6 8 19 3 0.000 0.000 0 242 249 0 16 00:05:20.00 20.0000D 0.4183 -0.4076 -0.1810 0.014 1.22 -0.811E+01 550. 0.287E-02 0.285E-02 0.640E-04 0.129E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 … … 262 269 0 36 00:12:00.00 20.0000A -0.5522 -0.6460 1.4350 0.052 1.18 -0.218E+01 500. 0.920E-02 0.917E-02 0.112E-03 0.522E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 263 270 0 37 00:12:20.00 20.0000A 0.6173 -0.7139 1.6131 0.055 1.18 -0.205E+01 500. 0.109E-01 0.108E-01 0.112E-03 0.591E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 264 0 38 00:12:40.00 20.0000A 0.6866 -0.7795 1.7862 0.059 1.18 -0.193E+01 500. 0.130E-01 0.130E-01 0.112E-03 0.67 3E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0271 0 38 00:12:40.00 20.0000A 0.6866 -0.7795 1.7862 0.059 1.18 -0.193E+01 500. 0.130E-01 0.130E-01 0.112E-03 0.672E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 265 272 0 39 00:13:00.00 20.0000A 0.7574 -0.8398 1.9133 0.062 1.18 -0.183E+01 500. 0.158E-01 0.158E-01 0.112E-03 0.767E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 266 273 0 40 00:13:20.00 20.0000A 0.8276 -0.8933 2.1010 0.066 1.18 -0.173E+01 500. 0.193E-01 0.193E-01 0.114E-03 0.877E-14 1 35 11 1 29 17 5 6 12 0.000 0.000 0 … … 284 291 0 58 00:17:41.39 10.4000A -1.4206 -1.3889 4.3431 0.094 1.38 -0.122E+01 800. 0.162E+00 0.162E+00 0.144E-03 0.204E-13 5 29 17 1 28 18 11 34 13 0.000 0.000 0 285 292 0 59 00:17:51.79 10.2000A -1.4342 1.3908 4.3989 0.095 1.38 -0.121E+01 800. 0.170E+00 0.170E+00 0.146E-03 0.208E-13 5 29 17 7 37 19 11 34 13 0.000 0.000 0 286 0 60 00:18:02.00 10.1000A -1.4413 1.490 24.4530 0.095 1.38 -0.120E+01 800. 0.177E+00 0.177E+00 0.148E-03 0.212E-13 4 32 19 7 37 19 12 34 13 0.000 0.000 0293 0 60 00:18:02.00 10.1000A -1.4413 1.4901 4.4530 0.095 1.38 -0.120E+01 800. 0.177E+00 0.177E+00 0.148E-03 0.212E-13 4 32 19 7 37 19 12 34 13 0.000 0.000 0 287 294 0 61 00:18:12.09 10.0000A -1.4694 -1.5081 4.4813 0.096 1.38 -0.120E+01 800. 0.185E+00 0.185E+00 0.151E-03 0.217E-13 4 32 19 8 7 16 12 34 13 0.000 0.000 0 288 295 0 62 00:18:22.09 10.0000A -1.4854 -1.6266 4.4842 0.096 1.38 -0.119E+01 800. 0.193E+00 0.192E+00 0.154E-03 0.222E-13 4 32 19 8 7 16 13 34 13 0.000 0.000 0 289 296 0 63 00:18:32.09 10.0000A -1.4915 -1.7186 4.4785 0.097 1.38 -0.118E+01 800. 0.200E+00 0.200E+00 0.158E-03 0.229E-13 4 32 19 8 7 16 12 34 13 0.000 0.000 0 290 297 0 64 00:18:42.09 10.1000A -1.5053 -1.7631 4.4346 0.097 1.38 -0.116E+01 800. 0.207E+00 0.207E+00 0.163E-03 0.235E-13 4 32 18 8 7 16 14 34 13 0.000 0.000 0 291 0 65 00:18:52.19 10.2000A -1.5591 -1.7454 4.404 60.098 1.41 -0.116E+01 850. 0.214E+00 0.214E+00 0.168E-03 0.244E-13 4 32 18 8 7 16 13 34 13 0.000 0.000 0298 0 65 00:18:52.19 10.2000A -1.5591 -1.7454 4.4047 0.098 1.41 -0.116E+01 850. 0.214E+00 0.214E+00 0.168E-03 0.244E-13 4 32 18 8 7 16 13 34 13 0.000 0.000 0 292 299 0 66 00:19:02.39 10.4000A -1.5902 -1.6634 4.3471 0.099 1.41 -0.115E+01 850. 0.221E+00 0.221E+00 0.175E-03 0.253E-13 4 32 18 8 7 16 13 34 13 0.000 0.000 0 293 300 0 67 00:19:12.79 10.5000A -1.5977 -1.6328 4.2947 0.099 1.41 -0.115E+01 850. 0.228E+00 0.228E+00 0.182E-03 0.264E-13 4 32 18 9 7 16 15 35 12 0.000 0.000 0 … … 297 304 0 71 00:19:55.39 11.0000A 1.6043 -1.8306 4.1022 0.102 1.41 -0.112E+01 850. 0.253E+00 0.252E+00 0.204E-03 0.295E-13 7 1 37 10 7 17 12 23 18 0.000 0.000 0 298 305 0 72 00:20:06.39 10.9000A 1.7066 -1.8092 4.1388 0.103 1.41 -0.111E+01 850. 0.258E+00 0.258E+00 0.214E-03 0.308E-13 7 24 30 10 7 17 12 23 18 0.000 0.000 0 299 0 73 00:20:17.29 10.9000A 1.839 61.7475 4.1429 0.104 1.41 -0.109E+01 850. 0.262E+00 0.262E+00 0.215E-03 0.309E-13 7 24 30 14 18 28 13 23 18 0.000 0.000 0306 0 73 00:20:17.29 10.9000A 1.8395 1.7475 4.1429 0.104 1.41 -0.109E+01 850. 0.262E+00 0.262E+00 0.215E-03 0.309E-13 7 24 30 14 18 28 13 23 18 0.000 0.000 0 300 307 0 74 00:20:28.20 10.7000A 1.9139 1.8281 4.1963 0.105 1.41 -0.107E+01 850. 0.267E+00 0.266E+00 0.219E-03 0.313E-13 7 24 30 14 18 28 13 23 18 0.000 0.000 0 301 308 0 75 00:20:38.90 10.7000A 1.9324 -1.8944 4.1945 0.106 1.41 -0.106E+01 850. 0.270E+00 0.270E+00 0.218E-03 0.311E-13 7 24 30 9 4 0 13 23 17 0.000 0.000 0 … … 305 312 0 79 00:21:21.30 10.4000A 2.0495 -2.1041 4.3252 0.110 1.41 -0.102E+01 850. 0.279E+00 0.279E+00 0.223E-03 0.316E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 306 313 0 80 00:21:31.70 10.4000A 2.0215 -2.0580 4.3215 0.111 1.41 -0.101E+01 850. 0.281E+00 0.280E+00 0.222E-03 0.315E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 307 0 81 00:21:42.10 10.5000A 1.9801 -1.9716 4.299 70.112 1.41 -0.100E+01 850. 0.282E+00 0.281E+00 0.224E-03 0.316E-13 7 24 31 9 4 0 15 23 17 0.000 0.000 0314 0 81 00:21:42.10 10.5000A 1.9801 -1.9716 4.2996 0.112 1.41 -0.100E+01 850. 0.282E+00 0.281E+00 0.224E-03 0.316E-13 7 24 31 9 4 0 15 23 17 0.000 0.000 0 308 315 0 82 00:21:52.60 10.4000A 1.9327 -1.9463 4.3286 0.113 1.41 -0.982E+00 850. 0.283E+00 0.282E+00 0.227E-03 0.321E-13 7 24 31 10 4 0 15 23 17 0.000 0.000 0 309 316 0 83 00:22:03.00 10.5000A 1.8937 -2.0053 4.3023 0.114 1.41 -0.972E+00 850. 0.283E+00 0.283E+00 0.226E-03 0.320E-13 14 4 0 10 4 0 15 23 17 0.000 0.000 0 … … 312 319 0 86 00:22:35.00 11.0000A 1.8894 -1.9584 4.0813 0.117 1.38 -0.945E+00 800. 0.283E+00 0.283E+00 0.237E-03 0.336E-13 14 4 0 10 4 0 16 23 17 0.000 0.000 0 313 320 0 87 00:22:46.00 11.4000A 1.8295 -1.8632 3.9326 0.118 1.38 -0.941E+00 800. 0.283E+00 0.283E+00 0.242E-03 0.344E-13 14 4 0 10 3 0 16 23 17 0.000 0.000 0 314 0 88 00:22:57.40 12.0000A 1.899 2-1.8015 3.7404 0.119 1.38 -0.929E+00 800. 0.283E+00 0.283E+00 0.252E-03 0.357E-13 15 4 1 7 24 4 16 23 17 0.000 0.000 0321 0 88 00:22:57.40 12.0000A 1.8991 -1.8015 3.7404 0.119 1.38 -0.929E+00 800. 0.283E+00 0.283E+00 0.252E-03 0.357E-13 15 4 1 7 24 4 16 23 17 0.000 0.000 0 315 322 0 89 00:23:09.40 12.6000A 1.9542 -1.8235 3.5820 0.119 1.38 -0.923E+00 800. 0.283E+00 0.282E+00 0.265E-03 0.376E-13 15 4 1 8 23 4 17 23 17 0.000 0.000 0 316 323 0 90 00:23:22.00 12.9000A 1.9250 -1.9310 3.4822 0.120 1.38 -0.920E+00 800. 0.282E+00 0.282E+00 0.278E-03 0.395E-13 15 4 1 8 23 4 16 12 29 0.000 0.000 0 317 0 91 00:23:34.90 13.2000A 1.849 8-2.0264 3.4174 0.121 1.38 -0.928E+00 800. 0.282E+00 0.281E+00 0.284E-03 0.404E-13 15 4 1 8 23 4 11 22 39 0.000 0.000 0324 0 91 00:23:34.90 13.2000A 1.8497 -2.0264 3.4174 0.121 1.38 -0.928E+00 800. 0.282E+00 0.281E+00 0.284E-03 0.404E-13 15 4 1 8 23 4 11 22 39 0.000 0.000 0 318 325 0 92 00:23:48.10 13.4000A 1.7617 -2.1071 3.3530 0.121 1.38 -0.909E+00 800. 0.281E+00 0.281E+00 0.290E-03 0.412E-13 15 4 1 8 23 4 12 22 39 0.000 0.000 0 319 0 93 00:24:01.50 13.7000A 1.7141 -2.1684 3.296 50.122 1.35 -0.903E+00 750. 0.280E+00 0.280E+00 0.293E-03 0.418E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0320 0 94 00:24:15.20 13.9000A 1.7172 -2.2027 3.2386 0.122 1.35 -0.901E+00 750. 0.280E+00 0.279E+00 0.29 8E-03 0.426E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0321 0 95 00:24:29.10 14.1000A 1.7197 -2.1962 3.187 4 0.123 1.35 -0.896E+00 750. 0.279E+00 0.279E+00 0.300E-03 0.431E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0326 0 93 00:24:01.50 13.7000A 1.7141 -2.1684 3.2966 0.122 1.35 -0.903E+00 750. 0.280E+00 0.280E+00 0.293E-03 0.418E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0 327 0 94 00:24:15.20 13.9000A 1.7172 -2.2027 3.2386 0.122 1.35 -0.901E+00 750. 0.280E+00 0.279E+00 0.297E-03 0.426E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0 328 0 95 00:24:29.10 14.1000A 1.7197 -2.1962 3.1875 0.123 1.35 -0.896E+00 750. 0.279E+00 0.279E+00 0.300E-03 0.430E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0 322 329 0 96 00:24:43.20 14.1000A 1.7148 -2.1716 3.1978 0.123 1.35 -0.891E+00 750. 0.279E+00 0.278E+00 0.303E-03 0.435E-13 3 22 38 7 23 5 10 23 1 0.000 0.000 0 323 0 97 00:24:57.30 13.8000A 1.6954 -2.1302 3.253 20.123 1.35 -0.890E+00 750. 0.279E+00 0.278E+00 0.302E-03 0.432E-13 3 22 38 7 23 5 10 23 1 0.000 0.000 0324 0 98 00:25:11.10 13.8000A -1.706 0 -2.1049 3.26910.123 1.35 -0.891E+00 750. 0.278E+00 0.278E+00 0.295E-03 0.420E-13 8 25 1 8 23 5 11 23 1 0.000 0.000 0325 0 99 00:25:24.90 13.8000A 1.7235 -2.1802 3.2581 0.123 1.35 -0.892E+00 750. 0.278E+00 0.278E+00 0.29 5E-03 0.416E-13 10 26 31 8 23 5 11 23 1 0.000 0.000 0326 0 100 00:25:38.70 13.8000A 1.7646 -2.1973 3.255 1 0.124 1.41 -0.893E+00 850. 0.278E+00 0.277E+00 0.294E-03 0.413E-13 10 26 31 8 23 5 12 23 1 0.000 0.000 0327 0 101 00:25:52.50 14.1000A 1.78 50 -2.1490 3.2018 0.124 1.38 -0.899E+00 800. 0.278E+00 0.277E+00 0.293E-03 0.409E-13 10 26 31 8 23 5 12 23 1 0.000 0.000 0328 0 102 00:26:06.60 14.1000A 1.7815 -2.092 23.1823 0.124 1.38 -0.897E+00 800. 0.278E+00 0.277E+00 0.299E-03 0.413E-13 10 26 31 8 22 5 13 23 1 0.000 0.000 0329 0 103 00:26:20.70 14.4000A -1.80 30-2.1453 3.1344 0.124 1.38 -0.895E+00 800. 0.277E+00 0.277E+00 0.298E-03 0.409E-13 16 22 38 8 22 5 12 23 19 0.000 0.000 0330 0 104 00:26:35.10 14.4000A 1.7946 -2.1085 3.1234 0.124 1.38 -0.894E+00 800. 0.277E+00 0.276E+00 0.30 3E-03 0.413E-13 1 36 26 8 22 5 12 23 19 0.000 0.000 0331 0 105 00:26:49.50 14.4000A 1.8514 -1.9889 3.1209 0.124 1.38 -0.894E+00 800. 0.276E+00 0.27 6E+00 0.301E-03 0.410E-13 1 37 26 8 22 5 9 38 26 0.000 0.000 0332 0 106 00:27:03.90 14.1000A 1.8658 1.864 63.1965 0.125 1.38 -0.898E+00 800. 0.276E+00 0.275E+00 0.299E-03 0.406E-13 1 37 26 8 35 24 10 38 26 0.000 0.000 0333 0 107 00:27:18.00 13.8000A 1.8273 1.8719 3.254 4 0.125 1.38 -0.900E+00 800. 0.275E+00 0.274E+00 0.292E-03 0.396E-13 1 37 26 8 35 24 10 38 26 0.000 0.000 0334 0 108 00:27:31.80 13.8000A 1.8485 1.899 8 3.2650 0.125 1.38 -0.899E+00 800. 0.274E+00 0.274E+00 0.285E-03 0.385E-13 1 16 38 8 35 23 10 38 26 0.000 0.000 0335 0 109 00:27:45.60 13.9000A 1.9161 1.9 299 3.2455 0.125 1.38 -0.900E+00 800. 0.274E+00 0.273E+00 0.284E-03 0.384E-13 1 16 39 8 35 23 10 38 26 0.000 0.000 0336 0 110 00:27:59.50 13.9000A 1.9930 1.948 8 3.2423 0.126 1.38 -0.902E+00 800. 0.273E+00 0.272E+00 0.285E-03 0.387E-13 1 16 39 8 36 23 11 38 26 0.000 0.000 0337 0 111 00:28:13.40 13.9000A 2.0631 2.0386 3.231 0 0.126 1.38 -0.897E+00 800. 0.272E+00 0.271E+00 0.285E-03 0.386E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0338 0 112 00:28:27.30 14.0000A 2.1131 2.0896 3.2197 0.126 1.38 -0.896E+00 800. 0.271E+00 0.270E+00 0.28 5E-03 0.386E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0339 0 113 00:28:41.30 14.0000A 2.1239 2.107 7 3.2062 0.126 1.41 -0.895E+00 850. 0.270E+00 0.269E+00 0.287E-03 0.388E-13 1 16 39 9 36 23 13 37 37 0.000 0.000 0340 0 114 00:28:55.30 14.0000A 2.1610 2.1093 3.2184 0.127 1.41 -0.893E+00 850. 0.269E+00 0.268E+00 0.286E-03 0.38 8E-13 1 16 0 9 36 23 13 37 37 0.000 0.000 0341 0 115 00:29:09.30 13.9000A 2.1914 2.143 0 3.2457 0.127 1.41 -0.889E+00 850. 0.268E+00 0.267E+00 0.286E-03 0.387E-13 1 16 0 9 37 23 14 37 37 0.000 0.000 0342 0 116 00:29:23.20 13.9000A 2.162 9 2.1318 3.2336 0.127 1.41 -0.888E+00 850. 0.266E+00 0.266E+00 0.284E-03 0.384E-13 1 16 0 9 37 23 14 37 37 0.000 0.000 0343 0 117 00:29:37.10 14.2000A 2.108 8 2.0451 3.1754 0.127 1.38 -0.891E+00 800. 0.265E+00 0.264E+00 0.285E-03 0.382E-13 1 16 1 9 37 23 14 37 37 0.000 0.000 0344 0 118 00:29:51.30 14.3000A 2.1133 1.903 0 3.1387 0.128 1.38 -0.904E+00 800. 0.263E+00 0.262E+00 0.291E-03 0.389E-13 1 16 1 9 37 23 14 9 10 0.000 0.000 0345 0 119 00:30:05.60 14.4000A 2.0602 -1.907 1 3.1173 0.128 1.38 -0.888E+00 800. 0.261E+00 0.260E+00 0.293E-03 0.390E-13 1 16 1 10 6 21 13 6 16 0.000 0.000 0346 0 120 00:30:20.00 14.1000A 2.0168 -2.016 8 3.1893 0.128 1.38 -0.885E+00 800. 0.259E+00 0.258E+00 0.294E-03 0.391E-13 1 17 1 10 6 21 7 25 19 0.000 0.000 0347 0 121 00:30:34.10 13.8000A -2.0339 -2.0424 3.255 4 0.128 1.38 -0.882E+00 800. 0.256E+00 0.256E+00 0.288E-03 0.381E-13 9 10 18 10 6 21 8 25 19 0.000 0.000 0348 0 122 00:30:47.90 13.4000A 1.9253 -2.10 19 3.3500 0.129 1.41 -0.884E+00 850. 0.254E+00 0.253E+00 0.281E-03 0.371E-13 1 17 1 10 5 21 8 25 19 0.000 0.000 0349 0 123 00:31:01.30 13.2000A -1.915 1 -2.1226 3.3987 0.129 1.41 -0.879E+00 850. 0.252E+00 0.251E+00 0.272E-03 0.358E-13 10 8 19 10 5 21 9 25 19 0.000 0.000 0350 0 124 00:31:14.50 13.1000A -1.960 6 -2.0641 3.4475 0.129 1.41 -0.879E+00 850. 0.249E+00 0.249E+00 0.268E-03 0.350E-13 10 9 18 10 5 21 9 25 19 0.000 0.000 0351 0 125 00:31:27.60 13.1000A -1.8503 -1.9343 3.438 6 0.128 1.41 -0.885E+00 850. 0.247E+00 0.246E+00 0.265E-03 0.345E-13 10 9 18 10 5 21 10 25 19 0.000 0.000 0352 0 126 00:31:40.70 13.1000A -1.915 3 -1.8985 3.4354 0.128 1.41 -0.884E+00 850. 0.244E+00 0.244E+00 0.264E-03 0.343E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0353 0 127 00:31:53.80 13.4000A -1.93 19 -1.8848 3.3614 0.128 1.41 -0.890E+00 850. 0.242E+00 0.241E+00 0.263E-03 0.340E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0354 0 128 00:32:07.20 13.7000A -1.961 0 -1.8736 3.2883 0.128 1.41 -0.894E+00 850. 0.239E+00 0.239E+00 0.268E-03 0.346E-13 13 26 19 1 13 16 9 9 19 0.000 0.000 0355 0 129 00:32:20.90 13.9000A -2.053 1 -1.8801 3.2398 0.127 1.41 -0.899E+00 850. 0.237E+00 0.236E+00 0.272E-03 0.351E-13 13 26 19 1 13 16 9 9 19 0.000 0.000 0356 0 130 00:32:34.80 14.1000A -2.092 5 -1.8885 3.1975 0.127 1.41 -0.901E+00 850. 0.234E+00 0.233E+00 0.275E-03 0.353E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0357 0 131 00:32:48.90 14.2000A -2.118 1 -1.8974 3.1694 0.127 1.41 -0.903E+00 850. 0.232E+00 0.231E+00 0.277E-03 0.355E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0358 0 132 00:33:03.10 14.6000A -2.1477 -1.90 40 3.0733 0.126 1.41 -0.903E+00 850. 0.229E+00 0.228E+00 0.278E-03 0.354E-13 14 26 19 1 13 16 10 8 19 0.000 0.000 0359 0 133 00:33:17.70 14.7000A -2.235 3 -1.9052 3.0618 0.126 1.41 -0.908E+00 850. 0.227E+00 0.226E+00 0.285E-03 0.360E-13 14 26 18 1 13 16 10 8 19 0.000 0.000 0360 0 134 00:33:32.40 14.7000A -2.323 1 -1.8976 3.0628 0.125 1.41 -0.911E+00 850. 0.225E+00 0.224E+00 0.286E-03 0.359E-13 14 26 18 1 13 16 11 8 19 0.000 0.000 0361 0 135 00:33:47.10 14.6000A -2.4176 -1.878 4 3.0724 0.125 1.41 -0.913E+00 850. 0.223E+00 0.222E+00 0.284E-03 0.355E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0362 0 136 00:34:01.70 14.9000A -2.5034 -1.877 7 3.0298 0.124 1.38 -0.912E+00 800. 0.221E+00 0.220E+00 0.281E-03 0.349E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0363 0 137 00:34:16.60 15.2000A -2.542 8 -1.8641 2.9686 0.124 1.38 -0.911E+00 800. 0.219E+00 0.218E+00 0.286E-03 0.353E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0364 0 138 00:34:31.80 15.3000A -2.493 3 -1.8350 2.9454 0.123 1.38 -0.912E+00 800. 0.218E+00 0.217E+00 0.290E-03 0.357E-13 14 26 18 1 12 16 12 8 19 0.000 0.000 0365 0 139 00:34:47.10 15.4000A -2.359 6 -1.7912 2.9241 0.122 1.38 -0.917E+00 800. 0.216E+00 0.215E+00 0.292E-03 0.356E-13 14 27 18 1 12 16 12 8 19 0.000 0.000 0366 0 140 00:35:02.50 15.5000A -2.391 9 -1.7717 2.8943 0.122 1.38 -0.917E+00 800. 0.215E+00 0.214E+00 0.293E-03 0.355E-13 15 26 18 1 11 15 12 8 19 0.000 0.000 0367 0 141 00:35:18.00 15.8000A -2.4 704 -1.7378 2.8456 0.121 1.38 -0.922E+00 800. 0.214E+00 0.213E+00 0.294E-03 0.354E-13 15 26 17 1 11 15 12 8 19 0.000 0.000 0368 0 142 00:35:33.80 15.8000A -2.46 81 1.7976 2.8489 0.121 1.38 -0.929E+00 800. 0.213E+00 0.212E+00 0.299E-03 0.357E-13 15 26 17 8 34 23 14 26 24 0.000 0.000 0369 0 143 00:35:49.60 15.7000A -2.36 24 -1.7673 2.8671 0.120 1.38 -0.927E+00 800. 0.213E+00 0.212E+00 0.299E-03 0.354E-13 15 26 17 1 8 35 6 38 32 0.000 0.000 0370 0 144 00:36:05.30 15.2000A -2.24 71 -1.7490 2.9525 0.119 1.38 -0.930E+00 800. 0.212E+00 0.211E+00 0.297E-03 0.349E-13 15 26 16 1 7 34 6 38 32 0.000 0.000 0371 0 145 00:36:20.50 14.7000A -2.09 89 -1.8137 3.0686 0.119 1.38 -0.939E+00 800. 0.212E+00 0.211E+00 0.288E-03 0.337E-13 15 26 16 1 6 34 7 38 32 0.000 0.000 0372 0 146 00:36:35.20 14.2000A -2.033 7 -1.8432 3.1685 0.118 1.38 -0.941E+00 800. 0.212E+00 0.211E+00 0.278E-03 0.325E-13 11 11 13 1 6 34 7 38 32 0.000 0.000 0373 0 147 00:36:49.40 13.9000A -2.049 5 -1.8095 3.2329 0.118 1.38 -0.945E+00 800. 0.211E+00 0.211E+00 0.269E-03 0.313E-13 10 12 14 1 6 34 7 38 32 0.000 0.000 0374 0 148 00:37:03.30 13.6000A -2.1261 -1.887 4 3.2988 0.117 1.38 -0.949E+00 800. 0.211E+00 0.210E+00 0.264E-03 0.306E-13 10 12 14 1 5 34 8 38 32 0.000 0.000 0375 0 149 00:37:16.90 13.4000A -2.1436 -1.8857 3.3518 0.117 1.38 -0.952E+00 800. 0.211E+00 0.210E+00 0.25 8E-03 0.299E-13 11 12 13 1 5 34 8 38 32 0.000 0.000 0376 0 150 00:37:30.30 13.4000A -2.1473 -1.8654 3.3558 0.116 1.38 -0.957E+00 800. 0.211E+00 0.210E+00 0.25 5E-03 0.294E-13 11 12 13 1 6 35 8 38 32 0.000 0.000 0377 0 151 00:37:43.70 13.5000A -2.115 0 -1.8966 3.3332 0.116 1.38 -0.961E+00 800. 0.211E+00 0.210E+00 0.255E-03 0.293E-13 12 12 12 1 5 35 9 38 32 0.000 0.000 0378 0 152 00:37:57.20 13.6000A -2.05 18 -1.9679 3.3052 0.116 1.38 -0.964E+00 800. 0.211E+00 0.210E+00 0.257E-03 0.293E-13 12 12 12 1 5 35 9 38 32 0.000 0.000 0379 0 153 00:38:10.80 13.9000A -1.97 07 -1.9302 3.2273 0.115 1.41 -0.968E+00 850. 0.211E+00 0.210E+00 0.259E-03 0.294E-13 13 12 12 1 5 35 9 38 32 0.000 0.000 0380 0 154 00:38:24.70 14.2000A 1.9577 -2.034 7 3.1608 0.115 1.41 -0.973E+00 850. 0.211E+00 0.210E+00 0.265E-03 0.298E-13 7 28 21 1 4 35 8 10 7 0.000 0.000 0381 0 155 00:38:38.90 13.4000A 1.988 5 -2.0495 3.3484 0.115 1.41 -0.974E+00 850. 0.211E+00 0.210E+00 0.271E-03 0.303E-13 7 28 22 1 4 35 8 9 7 0.000 0.000 0382 0 156 00:38:52.30 12.8000A 2.00 49 -1.9578 3.5085 0.114 1.41 -0.976E+00 850. 0.211E+00 0.210E+00 0.256E-03 0.286E-13 7 28 22 1 4 35 9 9 7 0.000 0.000 0383 0 157 00:39:05.10 12.6000A 1.949 5 -2.0192 3.5708 0.114 1.41 -0.980E+00 850. 0.211E+00 0.210E+00 0.246E-03 0.273E-13 7 28 22 1 3 35 9 9 7 0.000 0.000 0384 0 158 00:39:17.70 12.5000A -1.899 4 -1.9976 3.5936 0.114 1.41 -0.982E+00 850. 0.211E+00 0.210E+00 0.243E-03 0.268E-13 14 13 11 1 3 35 10 9 7 0.000 0.000 0385 0 159 00:39:30.20 12.8000A -1.834 3 -1.8951 3.5263 0.114 1.41 -0.982E+00 850. 0.211E+00 0.210E+00 0.241E-03 0.266E-13 14 13 11 1 3 35 10 8 6 0.000 0.000 0386 0 160 00:39:43.00 13.0000A -1.842 3 1.8539 3.4660 0.114 1.41 -0.988E+00 850. 0.212E+00 0.211E+00 0.248E-03 0.272E-13 8 13 14 17 15 11 10 8 6 0.000 0.000 0387 0 161 00:39:56.00 13.1000A -1.873 4 1.9441 3.4252 0.113 1.41 -0.987E+00 850. 0.212E+00 0.211E+00 0.252E-03 0.276E-13 8 13 14 17 15 11 11 8 6 0.000 0.000 0388 0 162 00:40:09.10 13.1000A -1.885 7 1.9860 3.4418 0.113 1.41 -0.989E+00 850. 0.212E+00 0.211E+00 0.255E-03 0.278E-13 8 13 14 17 15 11 10 9 8 0.000 0.000 0389 0 163 00:40:22.20 12.9000A -1.894 6 1.9867 3.4953 0.113 1.41 -0.992E+00 850. 0.212E+00 0.211E+00 0.256E-03 0.279E-13 9 13 14 17 16 11 10 9 8 0.000 0.000 0390 0 164 00:40:35.10 13.0000A -1.8988 2.10 83 3.4484 0.113 1.41 -0.994E+00 850. 0.213E+00 0.211E+00 0.253E-03 0.275E-13 9 13 14 17 16 11 11 9 8 0.000 0.000 0391 0 165 00:40:48.10 13.3000A -1.9469 2.1 590 3.3863 0.113 1.41 -0.994E+00 850. 0.213E+00 0.212E+00 0.256E-03 0.278E-13 18 13 9 17 16 11 11 8 7 0.000 0.000 0392 0 166 00:41:01.40 13.5000A -1.976 8 2.1184 3.3307 0.113 1.41 -0.996E+00 850. 0.213E+00 0.212E+00 0.263E-03 0.285E-13 18 13 9 17 16 11 12 8 7 0.000 0.000 0393 0 167 00:41:14.90 13.7000A -1.983 7 2.1241 3.2823 0.113 1.41 -0.100E+01 850. 0.213E+00 0.212E+00 0.268E-03 0.290E-13 18 14 9 17 17 11 11 8 8 0.000 0.000 0394 0 168 00:41:28.60 13.6000A -1.991 3 2.1833 3.3064 0.113 1.41 -0.100E+01 850. 0.214E+00 0.212E+00 0.274E-03 0.295E-13 18 14 9 17 17 11 12 8 8 0.000 0.000 0395 0 169 00:41:42.20 13.4000A -1.98 13 2.1737 3.3472 0.113 1.41 -0.101E+01 850. 0.214E+00 0.213E+00 0.273E-03 0.294E-13 18 14 9 18 18 9 12 8 8 0.000 0.000 0396 0 170 00:41:55.60 13.5000A -1.994 4 2.2409 3.3230 0.113 1.38 -0.101E+01 800. 0.214E+00 0.213E+00 0.270E-03 0.291E-13 18 14 8 17 18 10 12 8 8 0.000 0.000 0397 0 171 00:42:09.10 13.5000A -2.00 90 2.2335 3.3380 0.113 1.38 -0.101E+01 800. 0.215E+00 0.213E+00 0.272E-03 0.293E-13 18 14 8 17 18 10 13 8 8 0.000 0.000 0398 0 172 00:42:22.60 13.6000A -1.99 63 2.1836 3.3089 0.113 1.38 -0.101E+01 800. 0.215E+00 0.214E+00 0.273E-03 0.294E-13 18 14 7 18 19 8 13 8 8 0.000 0.000 0399 0 173 00:42:36.20 13.9000A -1.98 47 2.2215 3.2261 0.113 1.38 -0.101E+01 800. 0.216E+00 0.214E+00 0.276E-03 0.297E-13 18 14 7 18 19 8 13 8 8 0.000 0.000 0400 0 174 00:42:50.10 14.1000A -2.00 20 2.1920 3.1998 0.113 1.38 -0.102E+01 800. 0.217E+00 0.215E+00 0.282E-03 0.304E-13 18 15 7 17 19 10 14 8 8 0.000 0.000 0401 0 175 00:43:04.20 14.3000A -1.96 99 2.2250 3.1385 0.113 1.38 -0.101E+01 800. 0.218E+00 0.216E+00 0.286E-03 0.308E-13 18 15 7 17 20 9 14 8 8 0.000 0.000 0402 0 176 00:43:18.50 14.8000A -1.921 5 2.2241 3.0482 0.114 1.38 -0.100E+01 800. 0.219E+00 0.217E+00 0.290E-03 0.313E-13 17 14 4 17 20 9 14 8 8 0.000 0.000 0403 0 177 00:43:33.30 15.2000A -1.95 80 2.1654 2.9561 0.114 1.38 -0.100E+01 800. 0.220E+00 0.218E+00 0.300E-03 0.324E-13 17 14 4 17 20 9 14 8 8 0.000 0.000 0404 0 178 00:43:48.50 15.5000A -1.90 76 2.1079 2.9114 0.114 1.38 -0.994E+00 800. 0.221E+00 0.219E+00 0.308E-03 0.332E-13 17 14 4 18 18 12 10 9 18 0.000 0.000 0405 0 179 00:44:04.00 15.0000A -1.87 69 2.1355 2.9998 0.114 1.38 -0.100E+01 800. 0.222E+00 0.221E+00 0.314E-03 0.339E-13 17 15 4 18 19 11 10 9 18 0.000 0.000 0406 0 180 00:44:19.00 14.7000A -1.814 0 2.1157 3.0702 0.115 1.38 -0.992E+00 800. 0.224E+00 0.222E+00 0.305E-03 0.328E-13 16 15 3 18 17 13 11 9 18 0.000 0.000 0407 0 181 00:44:33.70 14.6000A -1.767 8 2.1454 3.0915 0.115 1.43 -0.991E+00 900. 0.225E+00 0.223E+00 0.300E-03 0.323E-13 16 15 3 18 17 13 13 11 16 0.000 0.000 0408 0 182 00:44:48.30 14.2000A 1.756 0 2.1502 3.1645 0.115 1.43 -0.995E+00 900. 0.226E+00 0.225E+00 0.299E-03 0.321E-13 9 27 23 18 17 13 13 11 16 0.000 0.000 0409 0 183 00:45:02.50 13.9000A 1.785 4 2.1719 3.2326 0.116 1.43 -0.990E+00 900. 0.228E+00 0.226E+00 0.292E-03 0.314E-13 9 27 23 18 18 13 14 11 16 0.000 0.000 0410 0 184 00:45:16.40 13.6000A 1.797 7 2.1976 3.3144 0.116 1.43 -0.992E+00 900. 0.230E+00 0.228E+00 0.287E-03 0.308E-13 9 27 23 18 18 13 14 11 16 0.000 0.000 0411 0 185 00:45:30.00 13.3000A 1.81 09 2.2067 3.3797 0.116 1.43 -0.100E+01 900. 0.231E+00 0.230E+00 0.281E-03 0.302E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0412 0 186 00:45:43.30 13.2000A 1.858 1 2.2022 3.4134 0.117 1.43 -0.991E+00 900. 0.233E+00 0.231E+00 0.276E-03 0.297E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0413 0 187 00:45:56.50 13.1000A 1.892 0 2.1844 3.4439 0.117 1.43 -0.987E+00 900. 0.234E+00 0.233E+00 0.275E-03 0.295E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0414 0 188 00:46:09.60 1 2.9000A 1.9095 2.1525 3.4749 0.117 1.43 -0.985E+00 900. 0.236E+00 0.235E+00 0.274E-03 0.294E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0415 0 189 00:46:22. 50 13.0000A 1.9088 2.1037 3.4713 0.118 1.43 -0.984E+00 900. 0.238E+00 0.237E+00 0.271E-03 0.291E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0416 0 190 00:46:35. 50 13.0000A 1.8892 2.0727 3.4501 0.118 1.43 -0.984E+00 900. 0.240E+00 0.238E+00 0.274E-03 0.294E-13 9 27 24 18 19 12 16 11 16 0.000 0.000 0417 0 191 00:46:48. 50 13.0000A 1.8879 2.0544 3.4586 0.118 1.43 -0.986E+00 900. 0.242E+00 0.240E+00 0.275E-03 0.295E-13 9 26 25 18 19 12 16 11 16 0.000 0.000 0418 0 192 00:47:01. 50 13.1000A 1.8897 2.0539 3.4431 0.119 1.41 -0.981E+00 850. 0.243E+00 0.242E+00 0.276E-03 0.296E-13 9 27 25 8 4 0 16 11 16 0.000 0.000 0419 0 193 00:47:14. 60 13.2000A 1.9011 2.0297 3.4101 0.119 1.41 -0.978E+00 850. 0.245E+00 0.244E+00 0.279E-03 0.300E-13 9 27 25 8 4 0 13 10 17 0.000 0.000 0420 0 194 00:47:27. 80 13.3000A 1.8908 2.0100 3.3868 0.119 1.41 -0.974E+00 850. 0.247E+00 0.246E+00 0.283E-03 0.305E-13 9 27 25 19 18 13 14 10 17 0.000 0.000 0421 0 195 00:47:41. 10 13.2000A 1.8585 2.0275 3.4008 0.119 1.41 -0.970E+00 850. 0.249E+00 0.248E+00 0.286E-03 0.309E-13 9 27 25 19 18 1414 10 17 0.000 0.000 0422 0 196 00:47:54. 30 13.2000A -1.8883 -1.9919 3.4148 0.119 1.41 -0.968E+00 850. 0.251E+00 0.250E+00 0.286E-03 0.310E-13 6 3 113 38 15 11 9 18 0.000 0.000 0423 0 197 00:48:07.5 0 12.9000A -2.0109 -1.9287 3.4810 0.119 1.41 -0.973E+00 850. 0.253E+00 0.252E+00 0.288E-03 0.313E-13 6 3 1 13 38 15 12 9 18 0.000 0.000 0424 0 198 00:48:20. 40 12.8000A -2.0676 1.9262 3.5058 0.119 1.41 -0.971E+00 850. 0.255E+00 0.253E+00 0.283E-03 0.309E-13 6 3 1 19 19 13 12 9 18 0.000 0.000 0425 0 199 00:48:33. 20 12.8000A -2.0940 1.8888 3.5269 0.119 1.41 -0.972E+00 850. 0.257E+00 0.255E+00 0.283E-03 0.309E-13 6 4 0 19 19 13 13 9 18 0.000 0.000 0426 0 200 00:48:46. 00 12.7000A 2.1500 1.8285 3.5411 0.119 1.41 -0.975E+00 850. 0.259E+00 0.257E+00 0.285E-03 0.312E-13 7 2 10 12 4 39 13 9 18 0.000 0.000 0427 0 201 00:48:58. 70 12.7000A 2.1955 1.8433 3.5368 0.119 1.41 -0.974E+00 850. 0.261E+00 0.259E+00 0.284E-03 0.313E-13 7 2 10 12 4 39 14 9 18 0.000 0.000 0428 0 202 00:49:11. 40 12.5000A -2.2761 -1.8582 3.5983 0.119 1.41 -0.976E+00 850. 0.262E+00 0.261E+00 0.286E-03 0.316E-13 7 3 0 3 8 2 9 39 1 0.000 0.000 0429 0 203 00:49:2 3.90 12.9000A -2.3204 -1.9224 3.4935 0.119 1.41 -0.973E+00 850. 0.264E+00 0.262E+00 0.284E-03 0.315E-13 7 3 0 3 8 2 15 9 18 0.000 0.000 0430 0 204 00:49:36. 80 12.8000A -2.2110 -1.9685 3.5230 0.119 1.41 -0.974E+00 850. 0.266E+00 0.264E+00 0.295E-03 0.327E-13 8 2 0 3 8 2 9 0 0 0.000 0.000 0431 0 205 00:49:49. 60 12.6000A -2.3818 -2.0373 3.5648 0.119 1.41 -0.975E+00 850. 0.268E+00 0.266E+00 0.294E-03 0.328E-13 8 1 0 3 7 2 9 7 17 0.000 0.000 0432 0 206 00:50:02. 20 12.4000A -2.4161 -2.1012 3.6437 0.118 1.41 -0.973E+00 850. 0.269E+00 0.267E+00 0.291E-03 0.325E-13 8 1 0 3 7 2 9 7 17 0.000 0.000 0433 0 207 00:50:14. 60 12.2000A -2.3757 -2.1220 3.6851 0.118 1.41 -0.972E+00 850. 0.271E+00 0.269E+00 0.287E-03 0.322E-13 8 1 0 3 7 2 10 7 17 0.000 0.000 0434 0 208 00:50:26. 80 12.0000A -2.4802 -2.0799 3.7439 0.118 1.43 -0.972E+00 900. 0.272E+00 0.270E+00 0.284E-03 0.319E-13 9 1 39 3 7 2 10 7 17 0.000 0.000 0435 0 209 00:50:38. 80 12.0000A -2.5304 -2.0332 3.7519 0.118 1.43 -0.970E+00 900. 0.273E+00 0.271E+00 0.280E-03 0.315E-13 9 1 39 4 7 3 10 7 17 0.000 0.000 0436 0 210 00:50:50. 80 12.1000A -2.4939 -2.0329 3.7144 0.118 1.43 -0.968E+00 900. 0.275E+00 0.272E+00 0.281E-03 0.317E-13 9 1 39 3 6 2 10 7 17 0.000 0.000 0437 0 211 00:51:0 2.90 12.0000A -2.4209 -1.9740 3.7581 0.118 1.43 -0.968E+00 900. 0.276E+00 0.274E+00 0.284E-03 0.321E-13 9 1 39 3 6 2 14 6 18 0.000 0.000 0438 0 212 00:51:1 4.90 11.7000A -2.3846 -2.0341 3.8348 0.118 1.43 -0.970E+00 900. 0.277E+00 0.275E+00 0.283E-03 0.320E-13 9 1 38 15 33 38 14 6 18 0.000 0.000 0439 0 213 00:51:26. 59 11.7000A -2.3506 -2.1307 3.8576 0.117 1.43 -0.975E+00 900. 0.278E+00 0.276E+00 0.277E-03 0.314E-13 10 1 38 15 33 38 14 6 18 0.000 0.000 0440 0 214 00:51:38. 29 11.8000A -2.3205 -2.1303 3.8242 0.117 1.43 -0.977E+00 900. 0.279E+00 0.277E+00 0.277E-03 0.315E-13 10 1 38 15 33 38 14 6 18 0.000 0.000 0441 0 215 00:51:50. 09 12.0000A -2.2789 -2.0978 3.7524 0.117 1.43 -0.981E+00 900. 0.280E+00 0.277E+00 0.281E-03 0.319E-13 11 0 38 15 32 38 15 5 18 0.000 0.000 0442 0 216 00:52:02. 09 12.0000A -2.2487 -2.0444 3.7616 0.117 1.43 -0.984E+00 900. 0.281E+00 0.278E+00 0.287E-03 0.326E-13 11 0 38 15 32 38 16 5 18 0.000 0.000 0443 0 217 00:52:14. 09 12.0000A -2.2217 -1.9955 3.7539 0.117 1.43 -0.980E+00 900. 0.281E+00 0.279E+00 0.288E-03 0.327E-13 11 0 37 15 32 37 16 5 18 0.000 0.000 0444 0 218 00:52:26. 09 12.2000A -2.1807 -2.0609 3.6997 0.117 1.43 -0.978E+00 900. 0.282E+00 0.279E+00 0.289E-03 0.329E-13 11 0 37 15 31 37 16 5 18 0.000 0.000 0445 0 219 00:52:38. 29 12.6000A -2.0993 -2.1447 3.5845 0.117 1.43 -0.978E+00 900. 0.282E+00 0.280E+00 0.294E-03 0.335E-13 11 0 37 15 31 37 16 5 18 0.000 0.000 0446 0 220 00:52:50. 89 13.2000A -1.9941 -2.1586 3.4126 0.116 1.43 -0.979E+00 900. 0.283E+00 0.280E+00 0.304E-03 0.346E-13 13 23 4 15 31 37 16 5 18 0.000 0.000 0447 0 221 00:53:04. 09 13.6000A -1.9882 -2.0853 3.3005 0.116 1.43 -0.983E+00 900. 0.283E+00 0.280E+00 0.318E-03 0.362E-13 13 23 4 15 31 37 17 4 18 0.000 0.000 0448 0 222 00:53:17. 69 13.6000A -1.9653 -2.0588 3.2981 0.116 1.43 -0.993E+00 900. 0.283E+00 0.280E+00 0.328E-03 0.372E-13 13 23 4 15 30 37 8 4 4 0.000 0.000 0449 0 223 00:53:31. 29 13.6000A -1.9332 2.0510 3.3020 0.116 1.43 -0.994E+00 900. 0.283E+00 0.280E+00 0.328E-03 0.370E-13 14 243 17 17 10 8 35 5 0.000 0.000 0450 0 224 00:53:44. 89 13.2000A -1.9953 2.0070 3.4031 0.116 1.43 -0.987E+00 900. 0.283E+00 0.280E+00 0.327E-03 0.369E-13 9 6 4 17 17 10 9 35 5 0.000 0.000 0451 0 225 00:53:58. 09 13.1000A -2.0501 1.9133 3.4417 0.116 1.43 -0.989E+00 900. 0.282E+00 0.280E+00 0.318E-03 0.357E-13 9 6 4 17 17 10 9 35 5 0.000 0.000 0452 0 226 00:54:11. 19 13.1000A -2.0938 1.8867 3.4369 0.115 1.43 -0.993E+00 900. 0.282E+00 0.279E+00 0.315E-03 0.352E-13 10 6 4 13 8 4 10 35 5 0.000 0.000 0453 0 227 00:54:24. 29 13.0000A -2.1923 1.9127 3.4680 0.115 1.43 -0.100E+01 900. 0.282E+00 0.279E+00 0.315E-03 0.351E-13 10 6 4 13 8 4 10 35 5 0.000 0.000 0454 0 228 00:54:37. 29 13.1000A -2.2403 1.9712 3.4314 0.115 1.40 -0.995E+00 850. 0.282E+00 0.279E+00 0.312E-03 0.347E-13 10 6 4 14 9 5 10 35 5 0.000 0.000 0455 0 229 00:54:50. 39 12.8000A -2.2360 2.0200 3.5122 0.115 1.40 -0.994E+00 850. 0.282E+00 0.279E+00 0.315E-03 0.349E-13 10 6 4 14 9 5 10 36 6 0.000 0.000 0456 0 230 00:55:03. 19 12.5000A -2.2380 2.0059 3.5910 0.115 1.40 -0.996E+00 850. 0.282E+00 0.279E+00 0.308E-03 0.340E-13 11 6 4 14 9 5 10 36 6 0.000 0.000 0457 0 231 00:55:15. 69 12.4000A -2.2664 1.9510 3.6299 0.115 1.40 -0.990E+00 850. 0.281E+00 0.279E+00 0.301E-03 0.332E-13 11 6 4 14 9 5 10 36 6 0.000 0.000 0458 0 232 00:55:28. 09 12.3000A -2.2644 1.9653 3.6637 0.115 1.40 -0.991E+00 850. 0.281E+00 0.279E+00 0.300E-03 0.329E-13 11 6 4 17 10 6 11 36 6 0.000 0.000 0459 0 233 00:55:40. 39 12.1000A -2.2635 1.9614 3.7037 0.115 1.40 -0.992E+00 850. 0.281E+00 0.279E+00 0.298E-03 0.326E-13 11 6 3 17 10 6 11 36 6 0.000 0.000 0460 0 234 00:55:52. 49 12.1000A -2.2125 1.9119 3.7274 0.115 1.40 -0.989E+00 850. 0.281E+00 0.278E+00 0.294E-03 0.321E-13 11 6 3 12 3 37 11 36 6 0.000 0.000 0461 0 235 00:56:04. 59 12.1000A -2.1711 1.8977 3.7116 0.115 1.40 -0.987E+00 850. 0.281E+00 0.278E+00 0.295E-03 0.321E-13 12 6 3 13 2 38 11 36 6 0.000 0.000 0462 0 236 00:56:16. 69 12.3000A -2.1905 1.9584 3.6493 0.115 1.40 -0.985E+00 850. 0.281E+00 0.278E+00 0.295E-03 0.321E-13 12 6 3 13 3 36 11 36 6 0.000 0.000 0463 0 237 00:56:2 8.99 12.3000A -2.1734 1.9387 3.6577 0.115 1.40 -0.985E+00 850. 0.280E+00 0.278E+00 0.300E-03 0.327E-13 12 6 3 13 3 3612 36 7 0.000 0.000 0464 0 238 00:56:41. 29 12.4000A -2.1181 1.9538 3.6195 0.115 1.40 -0.987E+00 850. 0.280E+00 0.278E+00 0.301E-03 0.327E-13 12 6 3 13 4 36 12 36 7 0.000 0.000 0465 0 239 00:56:53. 69 12.5000A -2.0250 1.9152 3.5962 0.115 1.38 -0.988E+00 800. 0.280E+00 0.277E+00 0.304E-03 0.330E-13 12 6 3 13 4 37 13 36 7 0.000 0.000 0466 0 240 00:57:06. 19 12.7000A -1.9906 1.9357 3.5442 0.115 1.38 -0.988E+00 800. 0.280E+00 0.277E+00 0.307E-03 0.334E-13 13 6 3 14 4 37 13 36 7 0.000 0.000 0467 0 241 00:57:18. 89 12.8000A 2.0433 1.9496 3.5178 0.115 1.38 -0.990E+00 800. 0.279E+00 0.277E+00 0.313E-03 0.340E-13 12 37 9 14 5 37 14 35 6 0.000 0.000 0468 0 242 00:57:31. 69 12.6000A 2.1147 1.9760 3.5706 0.115 1.40 -0.993E+00 850. 0.279E+00 0.277E+00 0.316E-03 0.343E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0469 0 243 00:57:44. 29 12.6000A 2.1795 1.9496 3.5723 0.115 1.40 -0.998E+00 850. 0.279E+00 0.276E+00 0.312E-03 0.339E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0470 0 244 00:57:56. 89 12.7000A 2.1966 1.9318 3.5502 0.114 1.40 -0.996E+00 850. 0.279E+00 0.276E+00 0.312E-03 0.340E-13 13 37 9 15 5 37 15 35 6 0.000 0.000 0471 0 245 00:58:09. 59 12.7000A 2.1739 1.9579 3.5513 0.114 1.40 -0.996E+00 850. 0.279E+00 0.276E+00 0.315E-03 0.343E-13 13 37 10 15 5 37 15 35 6 0.000 0.000 0472 0 246 00:58:22. 29 12.9000A 2.2250 1.9814 3.4946 0.114 1.38 -0.998E+00 800. 0.279E+00 0.276E+00 0.314E-03 0.344E-13 14 37 10 15 6 37 15 35 6 0.000 0.000 0473 0 247 00:58:35. 19 13.3000A 2.2578 1.9773 3.3900 0.114 1.38 -0.997E+00 800. 0.279E+00 0.276E+00 0.319E-03 0.350E-13 14 37 10 15 6 37 15 35 6 0.000 0.000 0474 0 248 00:58:48. 49 13.6000A 2.2207 1.9381 3.3184 0.114 1.38 -0.998E+00 800. 0.279E+00 0.276E+00 0.329E-03 0.361E-13 14 37 10 15 6 37 15 34 3 0.000 0.000 0475 0 249 00:59:02. 09 13.6000A 2.2069 1.9631 3.3056 0.114 1.38 -0.100E+01 800. 0.278E+00 0.276E+00 0.336E-03 0.369E-13 14 37 11 16 6 37 12 37 60.000 0.000 0476 0 250 00:59:15. 69 13.5000A 2.1185 1.9916 3.3423 0.114 1.38 -0.100E+01 800. 0.278E+00 0.276E+00 0.336E-03 0.369E-13 14 37 11 16 7 37 9 37 5 0.000 0.000 0477 0 251 00:59:29. 19 13.4000A 2.0242 2.0192 3.3527 0.113 1.43 -0.101E+01 900. 0.278E+00 0.276E+00 0.333E-03 0.365E-13 16 36 11 16 7 3710 37 5 0.000 0.000 0478 0 252 00:59:42. 59 13.2000A 1.9364 2.0137 3.4014 0.113 1.43 -0.101E+01 900. 0.278E+00 0.276E+00 0.331E-03 0.362E-13 15 37 12 16 7 37 10 37 5 0.000 0.000 0479 0 253 00:59:55. 79 13.3000A 1.9199 2.0227 3.3908 0.113 1.43 -0.102E+01 900. 0.278E+00 0.276E+00 0.325E-03 0.355E-13 16 36 12 16 7 36 10 37 5 0.000 0.000 0480 0 254 01:00:09. 09 13.2000A 1.9010 2.0017 3.4011 0.113 1.43 -0.102E+01 900. 0.278E+00 0.276E+00 0.328E-03 0.356E-13 16 36 12 16 7 36 11 37 5 0.000 0.000 0330 0 97 00:24:57.30 13.8000A 1.6954 -2.1302 3.2533 0.123 1.35 -0.890E+00 750. 0.279E+00 0.278E+00 0.302E-03 0.432E-13 3 22 38 7 23 5 10 23 1 0.000 0.000 0 331 0 98 00:25:11.10 13.8000A -1.7061 -2.1049 3.2692 0.123 1.35 -0.891E+00 750. 0.278E+00 0.278E+00 0.295E-03 0.420E-13 8 25 1 8 23 5 11 23 1 0.000 0.000 0 332 0 99 00:25:24.90 13.8000A 1.7235 -2.1802 3.2581 0.123 1.35 -0.892E+00 750. 0.278E+00 0.278E+00 0.294E-03 0.416E-13 10 26 31 8 23 5 11 23 1 0.000 0.000 0 333 0 100 00:25:38.70 13.8000A 1.7646 -2.1973 3.2552 0.124 1.41 -0.893E+00 850. 0.278E+00 0.277E+00 0.294E-03 0.412E-13 10 26 31 8 23 5 12 23 1 0.000 0.000 0 334 0 101 00:25:52.50 14.1000A 1.7849 -2.1490 3.2018 0.124 1.38 -0.899E+00 800. 0.278E+00 0.277E+00 0.293E-03 0.408E-13 10 26 31 8 23 5 12 23 1 0.000 0.000 0 335 0 102 00:26:06.60 14.1000A 1.7815 -2.0923 3.1823 0.124 1.38 -0.897E+00 800. 0.278E+00 0.277E+00 0.299E-03 0.413E-13 10 26 31 8 22 5 13 23 1 0.000 0.000 0 336 0 103 00:26:20.70 14.4000A -1.8029 -2.1453 3.1344 0.124 1.38 -0.895E+00 800. 0.277E+00 0.277E+00 0.298E-03 0.409E-13 16 22 38 8 22 5 12 23 19 0.000 0.000 0 337 0 104 00:26:35.10 14.4000A 1.7946 -2.1085 3.1234 0.124 1.38 -0.894E+00 800. 0.277E+00 0.276E+00 0.302E-03 0.413E-13 1 36 26 8 22 5 12 23 19 0.000 0.000 0 338 0 105 00:26:49.50 14.4000A 1.8514 -1.9889 3.1209 0.124 1.38 -0.894E+00 800. 0.276E+00 0.275E+00 0.300E-03 0.409E-13 1 37 26 8 22 5 9 38 26 0.000 0.000 0 339 0 106 00:27:03.90 14.1000A 1.8658 1.8647 3.1965 0.125 1.38 -0.898E+00 800. 0.276E+00 0.275E+00 0.299E-03 0.406E-13 1 37 26 8 35 24 10 38 26 0.000 0.000 0 340 0 107 00:27:18.00 13.8000A 1.8273 1.8719 3.2545 0.125 1.38 -0.900E+00 800. 0.275E+00 0.274E+00 0.291E-03 0.395E-13 1 37 26 8 35 24 10 38 26 0.000 0.000 0 341 0 108 00:27:31.80 13.8000A 1.8485 1.8999 3.2651 0.125 1.38 -0.899E+00 800. 0.274E+00 0.274E+00 0.284E-03 0.385E-13 1 16 38 8 35 23 10 38 26 0.000 0.000 0 342 0 109 00:27:45.60 13.9000A 1.9161 1.9300 3.2455 0.125 1.38 -0.900E+00 800. 0.274E+00 0.273E+00 0.283E-03 0.384E-13 1 16 39 8 35 23 10 38 26 0.000 0.000 0 343 0 110 00:27:59.50 13.9000A 1.9930 1.9489 3.2424 0.126 1.38 -0.902E+00 800. 0.273E+00 0.272E+00 0.285E-03 0.386E-13 1 16 39 8 36 23 11 38 26 0.000 0.000 0 344 0 111 00:28:13.40 13.9000A 2.0631 2.0386 3.2311 0.126 1.38 -0.897E+00 800. 0.272E+00 0.271E+00 0.284E-03 0.386E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0 345 0 112 00:28:27.30 14.0000A 2.1131 2.0896 3.2197 0.126 1.38 -0.896E+00 800. 0.271E+00 0.270E+00 0.284E-03 0.385E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0 346 0 113 00:28:41.30 14.0000A 2.1239 2.1078 3.2063 0.126 1.41 -0.895E+00 850. 0.270E+00 0.269E+00 0.286E-03 0.388E-13 1 16 39 9 36 23 13 37 37 0.000 0.000 0 347 0 114 00:28:55.30 14.0000A 2.1610 2.1093 3.2184 0.127 1.41 -0.893E+00 850. 0.269E+00 0.268E+00 0.286E-03 0.387E-13 1 16 0 9 36 23 13 37 37 0.000 0.000 0 348 0 115 00:29:09.30 13.9000A 2.1914 2.1431 3.2457 0.127 1.41 -0.889E+00 850. 0.268E+00 0.267E+00 0.286E-03 0.386E-13 1 16 0 9 37 23 14 37 37 0.000 0.000 0 349 0 116 00:29:23.20 13.9000A 2.1628 2.1319 3.2337 0.127 1.41 -0.888E+00 850. 0.266E+00 0.265E+00 0.284E-03 0.383E-13 1 16 0 9 37 23 14 37 37 0.000 0.000 0 350 0 117 00:29:37.10 14.2000A 2.1087 2.0452 3.1755 0.127 1.38 -0.891E+00 800. 0.265E+00 0.264E+00 0.284E-03 0.382E-13 1 16 1 9 37 23 14 37 37 0.000 0.000 0 351 0 118 00:29:51.30 14.3000A 2.1133 1.9031 3.1386 0.128 1.38 -0.904E+00 800. 0.263E+00 0.262E+00 0.290E-03 0.388E-13 1 16 1 9 37 23 14 9 10 0.000 0.000 0 352 0 119 00:30:05.60 14.4000A 2.0602 -1.9072 3.1172 0.128 1.38 -0.888E+00 800. 0.261E+00 0.260E+00 0.292E-03 0.390E-13 1 16 1 10 6 21 13 6 16 0.000 0.000 0 353 0 120 00:30:20.00 14.1000A 2.0168 -2.0169 3.1893 0.128 1.38 -0.885E+00 800. 0.259E+00 0.258E+00 0.293E-03 0.390E-13 1 17 1 10 6 21 7 25 19 0.000 0.000 0 354 0 121 00:30:34.10 13.8000A -2.0339 -2.0424 3.2555 0.128 1.38 -0.882E+00 800. 0.256E+00 0.256E+00 0.287E-03 0.380E-13 9 10 18 10 6 21 8 25 19 0.000 0.000 0 355 0 122 00:30:47.90 13.4000A 1.9253 -2.1020 3.3501 0.129 1.41 -0.884E+00 850. 0.254E+00 0.253E+00 0.280E-03 0.370E-13 1 17 1 10 5 21 8 25 19 0.000 0.000 0 356 0 123 00:31:01.30 13.2000A -1.9150 -2.1226 3.3988 0.129 1.41 -0.879E+00 850. 0.252E+00 0.251E+00 0.272E-03 0.357E-13 10 8 19 10 5 21 9 25 19 0.000 0.000 0 357 0 124 00:31:14.50 13.1000A -1.9605 -2.0641 3.4475 0.129 1.41 -0.879E+00 850. 0.249E+00 0.249E+00 0.267E-03 0.349E-13 10 9 18 10 5 21 9 25 19 0.000 0.000 0 358 0 125 00:31:27.60 13.1000A -1.8503 -1.9343 3.4387 0.128 1.41 -0.885E+00 850. 0.247E+00 0.246E+00 0.264E-03 0.344E-13 10 9 18 10 5 21 10 25 19 0.000 0.000 0 359 0 126 00:31:40.70 13.1000A -1.9154 -1.8984 3.4355 0.128 1.41 -0.884E+00 850. 0.244E+00 0.244E+00 0.263E-03 0.342E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0 360 0 127 00:31:53.80 13.4000A -1.9320 -1.8848 3.3615 0.128 1.41 -0.890E+00 850. 0.242E+00 0.241E+00 0.262E-03 0.339E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0 361 0 128 00:32:07.20 13.7000A -1.9612 -1.8735 3.2882 0.128 1.41 -0.894E+00 850. 0.239E+00 0.238E+00 0.267E-03 0.344E-13 13 26 19 1 13 16 9 9 19 0.000 0.000 0 362 0 129 00:32:20.90 13.9000A -2.0533 -1.8800 3.2398 0.127 1.41 -0.899E+00 850. 0.237E+00 0.236E+00 0.271E-03 0.349E-13 13 26 19 1 13 16 9 9 19 0.000 0.000 0 363 0 130 00:32:34.80 14.1000A -2.0927 -1.8885 3.1974 0.127 1.41 -0.901E+00 850. 0.234E+00 0.233E+00 0.274E-03 0.352E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0 364 0 131 00:32:48.90 14.2000A -2.1182 -1.8973 3.1694 0.127 1.41 -0.903E+00 850. 0.231E+00 0.231E+00 0.276E-03 0.354E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0 365 0 132 00:33:03.10 14.6000A -2.1477 -1.9039 3.0732 0.126 1.41 -0.903E+00 850. 0.229E+00 0.228E+00 0.277E-03 0.353E-13 14 26 19 1 13 16 10 8 19 0.000 0.000 0 366 0 133 00:33:17.70 14.7000A -2.2354 -1.9052 3.0618 0.126 1.41 -0.908E+00 850. 0.227E+00 0.226E+00 0.283E-03 0.359E-13 14 26 18 1 13 16 10 8 19 0.000 0.000 0 367 0 134 00:33:32.40 14.7000A -2.3232 -1.8975 3.0628 0.125 1.41 -0.911E+00 850. 0.224E+00 0.224E+00 0.284E-03 0.358E-13 14 26 18 1 13 16 11 8 19 0.000 0.000 0 368 0 135 00:33:47.10 14.6000A -2.4176 -1.8783 3.0724 0.125 1.41 -0.913E+00 850. 0.222E+00 0.222E+00 0.283E-03 0.354E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0 369 0 136 00:34:01.70 14.9000A -2.5034 -1.8776 3.0299 0.124 1.38 -0.912E+00 800. 0.221E+00 0.220E+00 0.280E-03 0.348E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0 370 0 137 00:34:16.60 15.2000A -2.5429 -1.8640 2.9687 0.124 1.38 -0.911E+00 800. 0.219E+00 0.218E+00 0.284E-03 0.351E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0 371 0 138 00:34:31.80 15.3000A -2.4934 -1.8349 2.9455 0.123 1.38 -0.912E+00 800. 0.217E+00 0.216E+00 0.289E-03 0.355E-13 14 26 18 1 12 16 12 8 19 0.000 0.000 0 372 0 139 00:34:47.10 15.4000A -2.3593 -1.7912 2.9244 0.122 1.38 -0.917E+00 800. 0.216E+00 0.215E+00 0.290E-03 0.354E-13 14 27 18 1 12 16 12 8 19 0.000 0.000 0 373 0 140 00:35:02.50 15.5000A -2.3911 -1.7717 2.8945 0.122 1.38 -0.917E+00 800. 0.215E+00 0.214E+00 0.291E-03 0.353E-13 15 26 18 1 11 15 12 8 19 0.000 0.000 0 374 0 141 00:35:18.00 15.8000A -2.4694 -1.7378 2.8459 0.121 1.38 -0.922E+00 800. 0.214E+00 0.213E+00 0.292E-03 0.352E-13 15 26 17 1 11 15 12 8 19 0.000 0.000 0 375 0 142 00:35:33.80 15.8000A -2.4673 1.7977 2.8489 0.121 1.38 -0.929E+00 800. 0.213E+00 0.212E+00 0.297E-03 0.355E-13 15 26 17 8 34 23 14 26 24 0.000 0.000 0 376 0 143 00:35:49.60 15.7000A -2.3617 -1.7673 2.8672 0.120 1.38 -0.927E+00 800. 0.212E+00 0.212E+00 0.296E-03 0.352E-13 15 26 17 1 8 35 6 38 32 0.000 0.000 0 377 0 144 00:36:05.30 15.2000A -2.2459 -1.7490 2.9526 0.119 1.38 -0.930E+00 800. 0.212E+00 0.211E+00 0.295E-03 0.347E-13 15 26 16 1 7 34 6 38 32 0.000 0.000 0 378 0 145 00:36:20.50 14.7000A -2.0976 -1.8137 3.0687 0.119 1.38 -0.939E+00 800. 0.212E+00 0.211E+00 0.286E-03 0.335E-13 15 26 16 1 6 34 7 38 32 0.000 0.000 0 379 0 146 00:36:35.20 14.2000A -2.0336 -1.8432 3.1686 0.118 1.38 -0.941E+00 800. 0.211E+00 0.210E+00 0.276E-03 0.323E-13 11 11 13 1 6 34 7 38 32 0.000 0.000 0 380 0 147 00:36:49.40 13.9000A -2.0496 -1.8095 3.2330 0.118 1.38 -0.945E+00 800. 0.211E+00 0.210E+00 0.267E-03 0.311E-13 10 12 14 1 6 34 7 38 32 0.000 0.000 0 381 0 148 00:37:03.30 13.6000A -2.1261 -1.8873 3.2989 0.117 1.38 -0.949E+00 800. 0.211E+00 0.210E+00 0.262E-03 0.304E-13 10 12 14 1 5 34 8 38 32 0.000 0.000 0 382 0 149 00:37:16.90 13.4000A -2.1436 -1.8857 3.3518 0.117 1.38 -0.952E+00 800. 0.211E+00 0.210E+00 0.256E-03 0.297E-13 11 12 13 1 5 34 8 38 32 0.000 0.000 0 383 0 150 00:37:30.30 13.4000A -2.1473 -1.8654 3.3558 0.116 1.38 -0.957E+00 800. 0.211E+00 0.210E+00 0.253E-03 0.292E-13 11 12 13 1 6 35 8 38 32 0.000 0.000 0 384 0 151 00:37:43.70 13.5000A -2.1151 -1.8964 3.3332 0.116 1.38 -0.961E+00 800. 0.211E+00 0.210E+00 0.253E-03 0.291E-13 12 12 12 1 5 35 9 38 32 0.000 0.000 0 385 0 152 00:37:57.20 13.6000A -2.0520 -1.9678 3.3051 0.116 1.38 -0.964E+00 800. 0.211E+00 0.210E+00 0.255E-03 0.292E-13 12 12 12 1 5 35 9 38 32 0.000 0.000 0 386 0 153 00:38:10.80 13.9000A -1.9710 -1.9302 3.2271 0.115 1.41 -0.968E+00 850. 0.211E+00 0.210E+00 0.257E-03 0.292E-13 13 12 12 1 5 35 9 38 32 0.000 0.000 0 387 0 154 00:38:24.70 14.2000A 1.9577 -2.0346 3.1608 0.115 1.41 -0.973E+00 850. 0.211E+00 0.210E+00 0.262E-03 0.297E-13 7 28 21 1 4 35 8 10 7 0.000 0.000 0 388 0 155 00:38:38.90 13.4000A 1.9887 -2.0495 3.3484 0.115 1.41 -0.974E+00 850. 0.211E+00 0.210E+00 0.268E-03 0.301E-13 7 28 22 1 4 35 8 9 7 0.000 0.000 0 389 0 156 00:38:52.30 12.8000A 2.0050 -1.9579 3.5087 0.114 1.41 -0.976E+00 850. 0.211E+00 0.210E+00 0.254E-03 0.284E-13 7 28 22 1 4 35 9 9 7 0.000 0.000 0 390 0 157 00:39:05.10 12.6000A 1.9496 -2.0192 3.5707 0.114 1.41 -0.980E+00 850. 0.211E+00 0.210E+00 0.244E-03 0.271E-13 7 28 22 1 3 35 9 9 7 0.000 0.000 0 391 0 158 00:39:17.70 12.5000A -1.8998 -1.9977 3.5937 0.114 1.41 -0.982E+00 850. 0.211E+00 0.210E+00 0.240E-03 0.266E-13 14 13 11 1 3 35 10 9 7 0.000 0.000 0 392 0 159 00:39:30.20 12.8000A -1.8349 -1.8953 3.5266 0.114 1.41 -0.982E+00 850. 0.211E+00 0.210E+00 0.239E-03 0.264E-13 14 13 11 1 3 35 10 8 6 0.000 0.000 0 393 0 160 00:39:43.00 13.0000A -1.8421 1.8536 3.4659 0.114 1.41 -0.988E+00 850. 0.211E+00 0.210E+00 0.245E-03 0.270E-13 8 13 14 17 15 11 10 8 6 0.000 0.000 0 394 0 161 00:39:56.00 13.1000A -1.8732 1.9441 3.4251 0.113 1.41 -0.987E+00 850. 0.212E+00 0.210E+00 0.250E-03 0.274E-13 8 13 14 17 15 11 11 8 6 0.000 0.000 0 395 0 162 00:40:09.10 13.1000A -1.8856 1.9865 3.4420 0.113 1.41 -0.989E+00 850. 0.212E+00 0.211E+00 0.253E-03 0.276E-13 8 13 14 17 15 11 10 9 8 0.000 0.000 0 396 0 163 00:40:22.20 12.9000A -1.8945 1.9871 3.4953 0.113 1.41 -0.992E+00 850. 0.212E+00 0.211E+00 0.254E-03 0.277E-13 9 13 14 17 16 11 10 9 8 0.000 0.000 0 397 0 164 00:40:35.10 13.0000A -1.8988 2.1091 3.4483 0.113 1.41 -0.994E+00 850. 0.212E+00 0.211E+00 0.251E-03 0.273E-13 9 13 14 17 16 11 11 9 8 0.000 0.000 0 398 0 165 00:40:48.10 13.3000A -1.9469 2.1604 3.3864 0.113 1.41 -0.994E+00 850. 0.212E+00 0.211E+00 0.254E-03 0.276E-13 18 13 9 17 16 11 11 8 7 0.000 0.000 0 399 0 166 00:41:01.40 13.5000A -1.9766 2.1203 3.3305 0.113 1.41 -0.996E+00 850. 0.213E+00 0.211E+00 0.261E-03 0.283E-13 18 13 9 17 16 11 12 8 7 0.000 0.000 0 400 0 167 00:41:14.90 13.7000A -1.9835 2.1254 3.2826 0.113 1.41 -0.100E+01 850. 0.213E+00 0.212E+00 0.266E-03 0.288E-13 18 14 9 17 17 11 11 8 8 0.000 0.000 0 401 0 168 00:41:28.60 13.6000A -1.9910 2.1849 3.3065 0.113 1.41 -0.100E+01 850. 0.213E+00 0.212E+00 0.271E-03 0.293E-13 18 14 9 17 17 11 12 8 8 0.000 0.000 0 402 0 169 00:41:42.20 13.4000A -1.9808 2.1753 3.3472 0.113 1.41 -0.101E+01 850. 0.214E+00 0.212E+00 0.270E-03 0.292E-13 18 14 9 18 18 9 12 8 8 0.000 0.000 0 403 0 170 00:41:55.60 13.5000A -1.9940 2.2417 3.3229 0.113 1.38 -0.101E+01 800. 0.214E+00 0.213E+00 0.267E-03 0.288E-13 18 14 8 17 18 10 12 8 8 0.000 0.000 0 404 0 171 00:42:09.10 13.5000A -2.0081 2.2342 3.3377 0.113 1.38 -0.101E+01 800. 0.214E+00 0.213E+00 0.270E-03 0.291E-13 18 14 8 17 18 10 13 8 8 0.000 0.000 0 405 0 172 00:42:22.60 13.6000A -1.9959 2.1839 3.3086 0.113 1.38 -0.101E+01 800. 0.215E+00 0.213E+00 0.270E-03 0.292E-13 18 14 7 18 19 8 13 8 8 0.000 0.000 0 406 0 173 00:42:36.20 13.9000A -1.9833 2.2213 3.2258 0.113 1.38 -0.101E+01 800. 0.216E+00 0.214E+00 0.273E-03 0.295E-13 18 14 7 18 19 8 13 8 8 0.000 0.000 0 407 0 174 00:42:50.10 14.1000A -2.0015 2.1933 3.1995 0.113 1.38 -0.102E+01 800. 0.216E+00 0.215E+00 0.279E-03 0.301E-13 18 15 7 17 19 10 14 8 8 0.000 0.000 0 408 0 175 00:43:04.20 14.3000A -1.9689 2.2253 3.1383 0.113 1.38 -0.101E+01 800. 0.217E+00 0.216E+00 0.283E-03 0.306E-13 18 15 7 17 20 9 14 8 8 0.000 0.000 0 409 0 176 00:43:18.50 14.8000A -1.9216 2.2244 3.0482 0.114 1.38 -0.100E+01 800. 0.218E+00 0.217E+00 0.287E-03 0.310E-13 17 14 4 17 20 9 14 8 8 0.000 0.000 0 410 0 177 00:43:33.30 15.2000A -1.9577 2.1655 2.9562 0.114 1.38 -0.100E+01 800. 0.219E+00 0.218E+00 0.297E-03 0.321E-13 17 14 4 17 20 9 14 8 8 0.000 0.000 0 411 0 178 00:43:48.50 15.5000A -1.9068 2.1078 2.9113 0.114 1.38 -0.994E+00 800. 0.221E+00 0.219E+00 0.305E-03 0.330E-13 17 14 4 18 18 12 10 9 18 0.000 0.000 0 412 0 179 00:44:04.00 15.0000A -1.8774 2.1347 2.9999 0.114 1.38 -0.100E+01 800. 0.222E+00 0.220E+00 0.311E-03 0.336E-13 17 15 4 18 19 11 10 9 18 0.000 0.000 0 413 0 180 00:44:19.00 14.7000A -1.8143 2.1133 3.0704 0.115 1.38 -0.992E+00 800. 0.223E+00 0.222E+00 0.302E-03 0.326E-13 16 15 3 18 17 13 11 9 18 0.000 0.000 0 414 0 181 00:44:33.70 14.6000A -1.7675 2.1428 3.0917 0.115 1.43 -0.991E+00 900. 0.225E+00 0.223E+00 0.297E-03 0.320E-13 16 15 3 18 17 13 13 11 16 0.000 0.000 0 415 0 182 00:44:48.30 14.2000A 1.7562 2.1477 3.1647 0.115 1.43 -0.994E+00 900. 0.226E+00 0.225E+00 0.296E-03 0.319E-13 9 27 23 18 17 13 13 11 16 0.000 0.000 0 416 0 183 00:45:02.50 13.9000A 1.7855 2.1698 3.2331 0.116 1.43 -0.990E+00 900. 0.228E+00 0.226E+00 0.288E-03 0.311E-13 9 27 23 18 18 13 14 11 16 0.000 0.000 0 417 0 184 00:45:16.40 13.6000A 1.7978 2.1955 3.3149 0.116 1.43 -0.992E+00 900. 0.229E+00 0.228E+00 0.283E-03 0.306E-13 9 27 23 18 18 13 14 11 16 0.000 0.000 0 418 0 185 00:45:30.00 13.3000A 1.8111 2.2046 3.3800 0.116 1.43 -0.100E+01 900. 0.231E+00 0.229E+00 0.278E-03 0.300E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0 419 0 186 00:45:43.30 13.2000A 1.8583 2.1999 3.4136 0.117 1.43 -0.991E+00 900. 0.232E+00 0.231E+00 0.273E-03 0.294E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0 420 0 187 00:45:56.50 13.1000A 1.8922 2.1820 3.4441 0.117 1.43 -0.987E+00 900. 0.234E+00 0.233E+00 0.272E-03 0.293E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 421 0 188 00:46:09.60 13.0000A 1.9096 2.1503 3.4749 0.117 1.43 -0.985E+00 900. 0.236E+00 0.235E+00 0.271E-03 0.292E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 422 0 189 00:46:22.60 13.0000A 1.9088 2.1015 3.4708 0.118 1.43 -0.984E+00 900. 0.238E+00 0.236E+00 0.270E-03 0.291E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 423 0 190 00:46:35.60 13.0000A 1.8889 2.0704 3.4500 0.118 1.43 -0.984E+00 900. 0.239E+00 0.238E+00 0.271E-03 0.292E-13 9 27 24 18 19 12 16 11 16 0.000 0.000 0 424 0 191 00:46:48.60 13.0000A 1.8877 2.0526 3.4578 0.118 1.43 -0.985E+00 900. 0.241E+00 0.240E+00 0.272E-03 0.293E-13 9 26 25 18 19 12 16 11 16 0.000 0.000 0 425 0 192 00:47:01.60 13.1000A 1.8905 2.0533 3.4417 0.119 1.41 -0.981E+00 850. 0.243E+00 0.242E+00 0.273E-03 0.294E-13 9 27 25 8 4 0 16 11 16 0.000 0.000 0 426 0 193 00:47:14.70 13.2000A 1.9015 2.0285 3.4092 0.119 1.41 -0.979E+00 850. 0.245E+00 0.244E+00 0.276E-03 0.298E-13 9 27 25 8 4 0 13 10 17 0.000 0.000 0 427 0 194 00:47:27.90 13.3000A 1.8909 2.0073 3.3854 0.119 1.41 -0.974E+00 850. 0.247E+00 0.246E+00 0.280E-03 0.303E-13 9 27 25 19 18 13 14 10 17 0.000 0.000 0 428 0 195 00:47:41.20 13.2000A 1.8588 -2.0241 3.3991 0.119 1.41 -0.970E+00 850. 0.249E+00 0.247E+00 0.284E-03 0.307E-13 9 27 26 13 38 15 14 10 17 0.000 0.000 0 429 0 196 00:47:54.40 13.2000A 1.8886 -1.9930 3.4144 0.119 1.41 -0.969E+00 850. 0.251E+00 0.249E+00 0.283E-03 0.308E-13 9 27 26 13 38 15 11 9 18 0.000 0.000 0 430 0 197 00:48:07.59 12.9000A -2.0108 -1.9302 3.4806 0.119 1.41 -0.973E+00 850. 0.253E+00 0.251E+00 0.285E-03 0.311E-13 6 3 1 13 38 15 12 9 18 0.000 0.000 0 431 0 198 00:48:20.50 12.8000A -2.0666 1.9227 3.5051 0.119 1.41 -0.971E+00 850. 0.254E+00 0.253E+00 0.280E-03 0.307E-13 6 3 1 19 19 13 12 9 18 0.000 0.000 0 432 0 199 00:48:33.30 12.8000A -2.0933 1.8878 3.5261 0.119 1.41 -0.972E+00 850. 0.256E+00 0.255E+00 0.280E-03 0.307E-13 6 4 0 19 19 13 13 9 18 0.000 0.000 0 433 0 200 00:48:46.10 12.7000A 2.1503 1.8281 3.5401 0.119 1.41 -0.975E+00 850. 0.258E+00 0.257E+00 0.282E-03 0.310E-13 7 2 10 12 4 39 13 9 18 0.000 0.000 0 434 0 201 00:48:58.80 12.7000A 2.1949 1.8433 3.5356 0.119 1.41 -0.974E+00 850. 0.260E+00 0.258E+00 0.281E-03 0.311E-13 7 2 10 12 4 39 14 9 18 0.000 0.000 0 435 0 202 00:49:11.50 12.5000A -2.2768 -1.8583 3.5983 0.119 1.41 -0.976E+00 850. 0.262E+00 0.260E+00 0.283E-03 0.314E-13 7 3 0 3 8 2 9 39 1 0.000 0.000 0 436 0 203 00:49:24.00 12.9000A -2.3194 -1.9224 3.4924 0.119 1.41 -0.973E+00 850. 0.264E+00 0.262E+00 0.281E-03 0.312E-13 7 3 0 3 8 2 15 9 18 0.000 0.000 0 437 0 204 00:49:36.90 12.8000A -2.2131 -1.9682 3.5215 0.119 1.41 -0.975E+00 850. 0.266E+00 0.264E+00 0.291E-03 0.325E-13 8 2 0 3 8 2 9 0 0 0.000 0.000 0 438 0 205 00:49:49.70 12.6000A -2.3826 -2.0375 3.5652 0.119 1.41 -0.975E+00 850. 0.267E+00 0.265E+00 0.291E-03 0.325E-13 8 1 0 3 7 2 9 7 17 0.000 0.000 0 439 0 206 00:50:02.30 12.4000A -2.4158 -2.1010 3.6437 0.118 1.41 -0.973E+00 850. 0.269E+00 0.267E+00 0.287E-03 0.322E-13 8 1 0 3 7 2 9 7 17 0.000 0.000 0 440 0 207 00:50:14.70 12.2000A -2.3751 -2.1213 3.6855 0.118 1.41 -0.972E+00 850. 0.270E+00 0.268E+00 0.284E-03 0.319E-13 8 1 0 3 7 2 10 7 17 0.000 0.000 0 441 0 208 00:50:26.90 12.0000A -2.4818 -2.0787 3.7439 0.118 1.43 -0.972E+00 900. 0.272E+00 0.270E+00 0.281E-03 0.316E-13 9 1 39 3 7 2 10 7 17 0.000 0.000 0 442 0 209 00:50:38.90 12.0000A -2.5308 -2.0324 3.7515 0.118 1.43 -0.970E+00 900. 0.273E+00 0.271E+00 0.277E-03 0.313E-13 9 1 39 4 7 3 10 7 17 0.000 0.000 0 443 0 210 00:50:50.90 12.1000A -2.4937 -2.0321 3.7136 0.118 1.43 -0.968E+00 900. 0.274E+00 0.272E+00 0.278E-03 0.315E-13 9 1 39 3 6 2 10 7 17 0.000 0.000 0 444 0 211 00:51:03.00 12.0000A -2.4204 -1.9727 3.7591 0.118 1.43 -0.968E+00 900. 0.275E+00 0.273E+00 0.281E-03 0.319E-13 9 1 39 3 6 2 14 6 18 0.000 0.000 0 445 0 212 00:51:15.00 11.7000A -2.3845 -2.0374 3.8360 0.118 1.43 -0.970E+00 900. 0.277E+00 0.274E+00 0.280E-03 0.318E-13 9 1 38 15 33 38 14 6 18 0.000 0.000 0 446 0 213 00:51:26.69 11.7000A -2.3517 -2.1313 3.8587 0.117 1.43 -0.975E+00 900. 0.278E+00 0.275E+00 0.274E-03 0.312E-13 10 1 38 15 33 38 14 6 18 0.000 0.000 0 447 0 214 00:51:38.39 11.8000A -2.3208 -2.1277 3.8253 0.117 1.43 -0.977E+00 900. 0.279E+00 0.276E+00 0.275E-03 0.313E-13 10 1 38 15 33 38 14 6 18 0.000 0.000 0 448 0 215 00:51:50.19 12.0000A -2.2790 -2.0966 3.7533 0.117 1.43 -0.981E+00 900. 0.279E+00 0.277E+00 0.278E-03 0.317E-13 11 0 38 15 32 38 15 5 18 0.000 0.000 0 449 0 216 00:52:02.19 12.0000A -2.2486 -2.0403 3.7622 0.117 1.43 -0.984E+00 900. 0.280E+00 0.278E+00 0.284E-03 0.324E-13 11 0 38 15 32 38 16 5 18 0.000 0.000 0 450 0 217 00:52:14.19 12.0000A -2.2214 -1.9990 3.7541 0.117 1.43 -0.980E+00 900. 0.281E+00 0.278E+00 0.285E-03 0.326E-13 11 0 37 15 32 37 16 5 18 0.000 0.000 0 451 0 218 00:52:26.19 12.2000A -2.1797 -2.0654 3.6993 0.117 1.43 -0.978E+00 900. 0.281E+00 0.279E+00 0.286E-03 0.327E-13 11 0 37 15 31 37 16 5 18 0.000 0.000 0 452 0 219 00:52:38.39 12.6000A -2.0979 -2.1467 3.5835 0.117 1.43 -0.978E+00 900. 0.282E+00 0.279E+00 0.291E-03 0.333E-13 11 0 37 15 31 37 16 5 18 0.000 0.000 0 453 0 220 00:52:50.99 13.2000A -1.9948 -2.1573 3.4108 0.116 1.43 -0.979E+00 900. 0.282E+00 0.280E+00 0.301E-03 0.344E-13 13 23 4 15 31 37 16 5 18 0.000 0.000 0 454 0 221 00:53:04.19 13.6000A -1.9884 -2.0801 3.2998 0.116 1.43 -0.983E+00 900. 0.282E+00 0.280E+00 0.315E-03 0.360E-13 13 23 4 15 31 37 17 4 18 0.000 0.000 0 455 0 222 00:53:17.79 13.6000A -1.9649 -2.0529 3.2994 0.116 1.43 -0.992E+00 900. 0.282E+00 0.280E+00 0.324E-03 0.369E-13 13 23 4 15 30 37 8 4 4 0.000 0.000 0 456 0 223 00:53:31.39 13.6000A -1.9354 2.0467 3.3034 0.116 1.43 -0.993E+00 900. 0.282E+00 0.280E+00 0.324E-03 0.368E-13 13 23 3 17 17 10 8 35 5 0.000 0.000 0 457 0 224 00:53:44.99 13.2000A -1.9955 2.0016 3.4048 0.116 1.43 -0.987E+00 900. 0.282E+00 0.279E+00 0.323E-03 0.366E-13 9 6 4 17 17 10 9 35 5 0.000 0.000 0 458 0 225 00:53:58.19 13.1000A -2.0497 1.9075 3.4432 0.116 1.43 -0.988E+00 900. 0.282E+00 0.279E+00 0.314E-03 0.354E-13 9 6 4 17 17 10 9 35 5 0.000 0.000 0 459 0 226 00:54:11.29 13.1000A -2.0943 1.8859 3.4378 0.115 1.43 -0.993E+00 900. 0.282E+00 0.279E+00 0.311E-03 0.350E-13 10 6 4 13 8 4 10 35 5 0.000 0.000 0 460 0 227 00:54:24.39 13.0000A -2.1924 1.9114 3.4684 0.115 1.43 -0.101E+01 900. 0.282E+00 0.279E+00 0.311E-03 0.348E-13 10 6 4 13 8 4 10 35 5 0.000 0.000 0 461 0 228 00:54:37.39 13.1000A -2.2396 1.9699 3.4311 0.115 1.43 -0.995E+00 900. 0.281E+00 0.279E+00 0.309E-03 0.344E-13 10 6 4 14 9 5 10 35 5 0.000 0.000 0 462 0 229 00:54:50.49 12.8000A -2.2345 2.0185 3.5139 0.115 1.40 -0.994E+00 850. 0.281E+00 0.279E+00 0.311E-03 0.346E-13 10 6 4 14 9 5 10 36 6 0.000 0.000 0 463 0 230 00:55:03.29 12.5000A -2.2386 2.0042 3.5923 0.115 1.40 -0.996E+00 850. 0.281E+00 0.278E+00 0.305E-03 0.338E-13 11 6 4 14 9 5 10 36 6 0.000 0.000 0 464 0 231 00:55:15.79 12.4000A -2.2667 1.9491 3.6309 0.115 1.40 -0.990E+00 850. 0.281E+00 0.278E+00 0.298E-03 0.330E-13 11 6 4 14 9 5 10 36 6 0.000 0.000 0 465 0 232 00:55:28.19 12.3000A -2.2639 1.9620 3.6661 0.115 1.40 -0.990E+00 850. 0.281E+00 0.278E+00 0.297E-03 0.327E-13 11 6 4 17 10 6 11 36 6 0.000 0.000 0 466 0 233 00:55:40.49 12.1000A -2.2609 1.9583 3.7059 0.115 1.40 -0.992E+00 850. 0.281E+00 0.278E+00 0.295E-03 0.324E-13 11 6 3 17 10 6 11 36 6 0.000 0.000 0 467 0 234 00:55:52.59 12.1000A -2.2079 1.9095 3.7294 0.115 1.40 -0.989E+00 850. 0.280E+00 0.278E+00 0.291E-03 0.319E-13 11 6 3 12 3 37 11 36 6 0.000 0.000 0 468 0 235 00:56:04.69 12.1000A -2.1707 1.8977 3.7131 0.115 1.40 -0.986E+00 850. 0.280E+00 0.278E+00 0.292E-03 0.319E-13 12 6 3 13 2 38 11 36 6 0.000 0.000 0 469 0 236 00:56:16.79 12.3000A -2.1897 1.9603 3.6503 0.115 1.40 -0.985E+00 850. 0.280E+00 0.277E+00 0.292E-03 0.320E-13 12 6 3 13 3 36 11 36 6 0.000 0.000 0 470 0 237 00:56:29.09 12.3000A -2.1717 1.9372 3.6538 0.115 1.40 -0.985E+00 850. 0.280E+00 0.277E+00 0.297E-03 0.325E-13 12 6 3 13 3 37 12 36 7 0.000 0.000 0 471 0 238 00:56:41.39 12.4000A -2.1149 1.9536 3.6151 0.115 1.40 -0.987E+00 850. 0.280E+00 0.277E+00 0.298E-03 0.325E-13 12 6 3 13 4 36 12 36 7 0.000 0.000 0 472 0 239 00:56:53.79 12.5000A -2.0201 1.9144 3.5932 0.115 1.40 -0.988E+00 850. 0.279E+00 0.277E+00 0.301E-03 0.328E-13 12 6 3 13 4 37 13 36 7 0.000 0.000 0 473 0 240 00:57:06.29 12.7000A -1.9931 1.9369 3.5414 0.115 1.38 -0.988E+00 800. 0.279E+00 0.277E+00 0.304E-03 0.332E-13 13 6 3 14 4 37 13 36 7 0.000 0.000 0 474 0 241 00:57:18.99 12.8000A 2.0384 1.9510 3.5139 0.115 1.38 -0.991E+00 800. 0.279E+00 0.276E+00 0.310E-03 0.338E-13 12 37 9 14 5 37 14 35 6 0.000 0.000 0 475 0 242 00:57:31.79 12.6000A 2.1168 1.9781 3.5652 0.115 1.40 -0.993E+00 850. 0.279E+00 0.276E+00 0.313E-03 0.341E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0 476 0 243 00:57:44.39 12.6000A 2.1807 1.9523 3.5655 0.114 1.40 -0.998E+00 850. 0.279E+00 0.276E+00 0.308E-03 0.337E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0 477 0 244 00:57:56.99 12.7000A 2.1971 1.9347 3.5448 0.114 1.40 -0.997E+00 850. 0.279E+00 0.276E+00 0.308E-03 0.337E-13 13 37 9 15 5 37 15 35 6 0.000 0.000 0 478 0 245 00:58:09.69 12.7000A 2.1710 1.9616 3.5440 0.114 1.40 -0.997E+00 850. 0.279E+00 0.276E+00 0.311E-03 0.341E-13 13 37 10 15 5 37 15 35 6 0.000 0.000 0 479 0 246 00:58:22.39 12.9000A 2.2295 1.9868 3.4858 0.114 1.38 -0.998E+00 800. 0.278E+00 0.276E+00 0.311E-03 0.341E-13 14 37 10 15 6 37 15 35 6 0.000 0.000 0 480 0 247 00:58:35.29 13.3000A 2.2611 1.9834 3.3800 0.114 1.38 -0.997E+00 800. 0.278E+00 0.276E+00 0.315E-03 0.347E-13 14 37 10 15 6 37 15 35 6 0.000 0.000 0 481 0 248 00:58:48.59 13.6000A 2.2233 1.9450 3.3102 0.114 1.38 -0.999E+00 800. 0.278E+00 0.276E+00 0.325E-03 0.358E-13 14 37 10 15 6 37 15 34 3 0.000 0.000 0 482 0 249 00:59:02.19 13.6000A 2.2127 1.9673 3.3017 0.114 1.38 -0.100E+01 800. 0.278E+00 0.276E+00 0.333E-03 0.366E-13 14 37 11 16 6 37 9 37 5 0.000 0.000 0 483 0 250 00:59:15.79 13.5000A 2.1257 1.9989 3.3397 0.113 1.38 -0.100E+01 800. 0.278E+00 0.275E+00 0.333E-03 0.366E-13 14 37 11 16 7 37 9 37 5 0.000 0.000 0 484 0 251 00:59:29.29 13.4000A 2.0291 2.0298 3.3500 0.113 1.38 -0.101E+01 800. 0.278E+00 0.275E+00 0.330E-03 0.363E-13 16 36 11 16 6 36 10 37 5 0.000 0.000 0 485 0 252 00:59:42.69 13.2000A 1.9435 2.0214 3.3980 0.113 1.43 -0.101E+01 900. 0.278E+00 0.275E+00 0.328E-03 0.360E-13 15 37 12 16 7 37 10 37 5 0.000 0.000 0 486 0 253 00:59:55.89 13.3000A 1.9283 2.0368 3.3876 0.113 1.43 -0.102E+01 900. 0.278E+00 0.275E+00 0.322E-03 0.353E-13 16 36 12 16 7 36 10 37 5 0.000 0.000 0 487 0 254 01:00:09.19 13.2000A 1.9061 2.0144 3.3966 0.113 1.43 -0.102E+01 900. 0.277E+00 0.275E+00 0.325E-03 0.355E-13 16 36 12 16 7 36 11 37 5 0.000 0.000 0 -
palm/trunk/SOURCE/check_parameters.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dz(1), error message revised 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! add variable description 28 31 ! … … 4178 4181 IF ( ( constant_flux_layer .OR. & 4179 4182 INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) & 4180 .AND. roughness_length >= 0.5 * dz ) THEN4183 .AND. roughness_length >= 0.5 * dz(1) ) THEN 4181 4184 message_string = 'roughness_length must be smaller than dz/2' 4182 4185 CALL message( 'check_parameters', 'PA0424', 1, 2, 0, 6, 0 ) … … 4198 4201 !-- Check if vertical grid stretching is switched off in case of complex 4199 4202 !-- terrain simulations 4200 IF ( complex_terrain .AND. dz_stretch_level < 100000.0_wp ) THEN 4203 IF ( complex_terrain .AND. & 4204 ANY( dz_stretch_level_start /= -9999999.9_wp ) ) THEN 4201 4205 message_string = 'Vertical grid stretching is not allowed for ' // & 4202 4206 'complex_terrain = .T.' -
palm/trunk/SOURCE/header.f90
r3045 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Header output concerning stretching revised 28 ! 29 ! 3045 2018-05-28 07:55:41Z Giersch 27 30 ! Error messages revised 28 31 ! … … 800 803 !-- Computational grid 801 804 IF ( .NOT. ocean ) THEN 802 WRITE ( io, 250 ) dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(nzt+1) 803 IF ( dz_stretch_level_index < nzt+1 ) THEN 804 WRITE ( io, 252 ) dz_stretch_level, dz_stretch_level_index, & 805 dz_stretch_factor, dz_max 806 ENDIF 805 WRITE ( io, 250 ) dx, dy 806 807 DO i = 1, number_stretch_level_start+1 808 WRITE ( io, 253 ) i, dz(i) 809 ENDDO 810 811 WRITE( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(nzt+1) 812 813 IF ( ANY( dz_stretch_level_start_index < nzt+1 ) ) THEN 814 WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:' 815 DO i = 1, number_stretch_level_start 816 WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_start(i), ' m,' 817 ENDDO 818 WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: ' 819 DO i = 1, number_stretch_level_start 820 WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_start_index(i), ',' 821 ENDDO 822 WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height: ' 823 DO i = 1, number_stretch_level_start 824 WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_end(i), ' m,' 825 ENDDO 826 WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index: ' 827 DO i = 1, number_stretch_level_start 828 WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_end_index(i), ',' 829 ENDDO 830 WRITE( io, '(/,A)', advance='no') ' Factor used for stretching: ' 831 DO i = 1, number_stretch_level_start 832 WRITE ( io, '(F12.3,A1)', advance='no' ) dz_stretch_factor_array(i), ',' 833 ENDDO 834 ENDIF 835 807 836 ELSE 808 WRITE ( io, 250 ) dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(0) 809 IF ( dz_stretch_level_index > 0 ) THEN 810 WRITE ( io, 252 ) dz_stretch_level, dz_stretch_level_index, & 811 dz_stretch_factor, dz_max 812 ENDIF 813 ENDIF 814 WRITE ( io, 254 ) nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), & 837 WRITE ( io, 250 ) dx, dy 838 DO i = 1, number_stretch_level_start+1 839 WRITE ( io, 253 ) i, dz(i) 840 ENDDO 841 842 WRITE ( io, 251 ) (nx+1)*dx, (ny+1)*dy, zu(0) 843 844 IF ( ANY( dz_stretch_level_start_index > 0 ) ) THEN 845 WRITE( io, '(A)', advance='no') ' Vertical stretching starts at height:' 846 DO i = 1, number_stretch_level_start 847 WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_start(i), ' m,' 848 ENDDO 849 WRITE( io, '(/,A)', advance='no') ' Vertical stretching starts at index: ' 850 DO i = 1, number_stretch_level_start 851 WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_start_index(i), ',' 852 ENDDO 853 WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at height: ' 854 DO i = 1, number_stretch_level_start 855 WRITE ( io, '(F10.1,A3)', advance='no' ) dz_stretch_level_end(i), ' m,' 856 ENDDO 857 WRITE( io, '(/,A)', advance='no') ' Vertical stretching ends at index: ' 858 DO i = 1, number_stretch_level_start 859 WRITE ( io, '(I12,A1)', advance='no' ) dz_stretch_level_end_index(i), ',' 860 ENDDO 861 WRITE( io, '(/,A)', advance='no') ' Factor used for stretching: ' 862 DO i = 1, number_stretch_level_start 863 WRITE ( io, '(F12.3,A1)', advance='no' ) dz_stretch_factor_array(i), ',' 864 ENDDO 865 ENDIF 866 ENDIF 867 WRITE ( io, 254 ) nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), & 815 868 MIN( nnz+2, nzt+2 ) 816 869 IF ( sloping_surface ) WRITE ( io, 260 ) alpha_surface … … 2083 2136 250 FORMAT (//' Computational grid and domain size:'/ & 2084 2137 ' ----------------------------------'// & 2085 ' Grid length: dx = ',F7.3,' m dy = ',F7.3, & 2086 ' m dz = ',F7.3,' m'/ & 2087 ' Domain size: x = ',F10.3,' m y = ',F10.3, & 2138 ' Grid length: dx = ',F8.3,' m dy = ',F8.3, ' m') 2139 251 FORMAT ( /' Domain size: x = ',F10.3,' m y = ',F10.3, & 2088 2140 ' m z(u) = ',F10.3,' m'/) 2089 252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', & 2090 ' factor:',F6.3/ & 2091 ' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/) 2092 254 FORMAT (' Number of gridpoints (x,y,z): (0:',I4,', 0:',I4,', 0:',I4,')'/ & 2141 253 FORMAT ( ' dz(',I1,') = ', F8.3, ' m') 2142 254 FORMAT (//' Number of gridpoints (x,y,z): (0:',I4,', 0:',I4,', 0:',I4,')'/ & 2093 2143 ' Subdomain size (x,y,z): ( ',I4,', ',I4,', ',I4,')'/) 2094 2144 260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,& -
palm/trunk/SOURCE/init_grid.f90
r3051 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! New vertical stretching mechanism introduced 28 ! 29 ! 3051 2018-05-30 17:43:55Z suehring 27 30 ! Minor bugfix concerning mapping 3D buildings on top of terrain 28 31 ! … … 325 328 canyon_height, canyon_wall_left, canyon_wall_south, & 326 329 canyon_width_x, canyon_width_y, constant_flux_layer, & 327 dp_level_ind_b, dz, dz_max, dz_stretch_factor, & 328 dz_stretch_level, dz_stretch_level_index, grid_level, & 330 dp_level_ind_b, dz, dz_max, dz_stretch_factor, & 331 dz_stretch_factor_array, dz_stretch_level, dz_stretch_level_end,& 332 dz_stretch_level_end_index, dz_stretch_level_start_index, & 333 dz_stretch_level_start, grid_level, & 329 334 force_bound_l, force_bound_r, force_bound_n, force_bound_s, & 330 335 ibc_uv_b, inflow_l, inflow_n, inflow_r, inflow_s, & 331 336 masking_method, maximum_grid_level, message_string, & 332 momentum_advec, nest_domain, nest_bound_l, nest_bound_n, & 333 nest_bound_r, nest_bound_s, ocean, outflow_l, outflow_n, & 334 outflow_r, outflow_s, psolver, scalar_advec, topography, & 335 topography_grid_convention, tunnel_height, tunnel_length, & 336 tunnel_width_x, tunnel_width_y, tunnel_wall_depth, & 337 use_surface_fluxes, use_top_fluxes, wall_adjustment_factor 337 momentum_advec, nest_domain, nest_bound_l, & 338 nest_bound_n, nest_bound_r, nest_bound_s, & 339 number_stretch_level_end, number_stretch_level_start, ocean, & 340 outflow_l, outflow_n, outflow_r, outflow_s, psolver, & 341 scalar_advec, topography, topography_grid_convention, & 342 tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y, & 343 tunnel_wall_depth, use_surface_fluxes, use_top_fluxes, & 344 wall_adjustment_factor 338 345 339 346 USE grid_variables, & … … 361 368 IMPLICIT NONE 362 369 363 INTEGER(iwp) :: i !< index variable along x 364 INTEGER(iwp) :: j !< index variable along y 365 INTEGER(iwp) :: k !< index variable along z 366 INTEGER(iwp) :: k_top !< topography top index on local PE 367 INTEGER(iwp) :: l !< loop variable 368 INTEGER(iwp) :: nzb_local_max !< vertical grid index of maximum topography height 369 INTEGER(iwp) :: nzb_local_min !< vertical grid index of minimum topography height 370 INTEGER(iwp) :: i !< index variable along x 371 INTEGER(iwp) :: j !< index variable along y 372 INTEGER(iwp) :: k !< index variable along z 373 INTEGER(iwp) :: k_top !< topography top index on local PE 374 INTEGER(iwp) :: n !< loop variable for stretching 375 INTEGER(iwp) :: number_dz !< number of user-specified dz values 376 INTEGER(iwp) :: nzb_local_max !< vertical grid index of maximum topography height 377 INTEGER(iwp) :: nzb_local_min !< vertical grid index of minimum topography height 370 378 371 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_local 372 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_tmp 379 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_local !< index for topography top at cell-center 380 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: nzb_tmp !< dummy to calculate topography indices on u- and v-grid 373 381 374 382 INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE :: topo !< input array for 3D topography and dummy array for setting "outer"-flags 375 383 384 REAL(wp) :: dz_level_end !< distance between calculated height level for u/v-grid and user-specified end level for stretching 376 385 REAL(wp) :: dz_stretched !< stretched vertical grid spacing 386 387 REAL(wp), DIMENSION(:), ALLOCATABLE :: min_dz_stretch_level_end !< Array that contains all minimum heights where the stretching can end 377 388 378 389 … … 391 402 ! 392 403 !-- Compute height of u-levels from constant grid length and dz stretch factors 393 IF ( dz == -1.0_wp ) THEN404 IF ( dz(1) == -1.0_wp ) THEN 394 405 message_string = 'missing dz' 395 406 CALL message( 'init_grid', 'PA0200', 1, 2, 0, 6, 0 ) 396 ELSEIF ( dz <= 0.0_wp ) THEN397 WRITE( message_string, * ) 'dz=',dz ,' <= 0.0'407 ELSEIF ( dz(1) <= 0.0_wp ) THEN 408 WRITE( message_string, * ) 'dz=',dz(1),' <= 0.0' 398 409 CALL message( 'init_grid', 'PA0201', 1, 2, 0, 6, 0 ) 399 410 ENDIF 400 411 401 412 ! 413 !-- Initialize dz_stretch_level_start with the value of dz_stretch_level 414 !-- if it was set by the user 415 IF ( dz_stretch_level /= -9999999.9_wp ) THEN 416 dz_stretch_level_start(1) = dz_stretch_level 417 ENDIF 418 419 ! 420 !-- Determine number of dz values and stretching levels specified by the 421 !-- user to allow right controlling of the stretching mechanism and to 422 !-- perform error checks 423 number_dz = COUNT( dz /= -1.0_wp ) 424 number_stretch_level_start = COUNT( dz_stretch_level_start /= & 425 -9999999.9_wp ) 426 number_stretch_level_end = COUNT( dz_stretch_level_end /= & 427 9999999.9_wp ) 428 429 ! 430 !-- The number of specified end levels +1 has to be the same than the number 431 !-- of specified dz values 432 IF ( number_dz /= number_stretch_level_end + 1 ) THEN 433 WRITE( message_string, * ) 'The number of values for dz = ', & 434 number_dz, 'has to be the same than& ', & 435 'the number of values for ', & 436 'dz_stretch_level_end + 1 = ', & 437 number_stretch_level_end+1 438 CALL message( 'init_grid', 'PA0156', 1, 2, 0, 6, 0 ) 439 ENDIF 440 441 ! 442 !-- The number of specified start levels has to be the same or one less than 443 !-- the number of specified dz values 444 IF ( number_dz /= number_stretch_level_start + 1 .AND. & 445 number_dz /= number_stretch_level_start ) THEN 446 WRITE( message_string, * ) 'The number of values for dz = ', & 447 number_dz, 'has to be the same or one ', & 448 'more than& the number of values for ', & 449 'dz_stretch_level_start = ', & 450 number_stretch_level_start 451 CALL message( 'init_grid', 'PA0211', 1, 2, 0, 6, 0 ) 452 ENDIF 453 454 !-- The number of specified start levels has to be the same or one more than 455 !-- the number of specified end levels 456 IF ( number_stretch_level_start /= number_stretch_level_end + 1 .AND. & 457 number_stretch_level_start /= number_stretch_level_end ) THEN 458 WRITE( message_string, * ) 'The number of values for ', & 459 'dz_stretch_level_start = ', & 460 dz_stretch_level_start, 'has to be the ',& 461 'same or one more than& the number of ', & 462 'values for dz_stretch_level_end = ', & 463 number_stretch_level_end 464 CALL message( 'init_grid', 'PA0216', 1, 2, 0, 6, 0 ) 465 ENDIF 466 467 ! 468 !-- Initialize dz for the free atmosphere with the value of dz_max 469 IF ( dz(number_stretch_level_start+1) == -1.0_wp .AND. & 470 number_stretch_level_start /= 0 ) THEN 471 dz(number_stretch_level_start+1) = dz_max 472 ENDIF 473 474 ! 475 !-- Initialize the stretching factor if (infinitely) stretching in the free 476 !-- atmosphere is desired (dz_stretch_level_end was not specified for the 477 !-- free atmosphere) 478 IF ( number_stretch_level_start == number_stretch_level_end + 1 ) THEN 479 dz_stretch_factor_array(number_stretch_level_start) = & 480 dz_stretch_factor 481 ENDIF 482 483 ! 484 !-- Allocation of arrays for stretching 485 ALLOCATE( min_dz_stretch_level_end(number_stretch_level_start) ) 486 487 ! 402 488 !-- Define the vertical grid levels 403 489 IF ( .NOT. ocean ) THEN 490 491 ! 492 !-- The stretching region has to be large enough to allow for a smooth 493 !-- transition between two different grid spacings 494 DO n = 1, number_stretch_level_start 495 min_dz_stretch_level_end(n) = dz_stretch_level_start(n) + & 496 4 * MAX( dz(n),dz(n+1) ) 497 ENDDO 498 499 IF ( ANY( min_dz_stretch_level_end > dz_stretch_level_end ) ) THEN 500 message_string= 'Eeach dz_stretch_level_end has to be larger ' // & 501 'than its corresponding value for &' // & 502 'dz_stretch_level_start + 4*MAX(dz(n),dz(n+1)) '//& 503 'to allow for smooth grid stretching' 504 CALL message( 'init_grid', 'PA0224', 1, 2, 0, 6, 0 ) 505 ENDIF 506 507 ! 508 !-- Stretching must not be applied within the prandtl_layer 509 !-- (first two grid points). For the default case dz_stretch_level_start 510 !-- is negative. Therefore the absolut value is checked here. 511 IF ( ANY( ABS( dz_stretch_level_start ) < dz(1) * 1.5_wp ) ) THEN 512 WRITE( message_string, * ) 'Eeach dz_stretch_level_start has to be ',& 513 'larger than ', dz(1) * 1.5 514 CALL message( 'init_grid', 'PA0226', 1, 2, 0, 6, 0 ) 515 ENDIF 516 517 ! 518 !-- The stretching has to start and end on a grid level. Therefore 519 !-- user-specified values have to ''interpolate'' to the next lowest level 520 IF ( number_stretch_level_start /= 0 ) THEN 521 dz_stretch_level_start(1) = INT( (dz_stretch_level_start(1) - & 522 dz(1)/2.0) / dz(1) ) & 523 * dz(1) + dz(1)/2.0 524 ENDIF 525 526 IF ( number_stretch_level_start > 1 ) THEN 527 DO n = 2, number_stretch_level_start 528 dz_stretch_level_start(n) = INT( dz_stretch_level_start(n) / & 529 dz(n) ) * dz(n) 530 ENDDO 531 ENDIF 532 533 IF ( number_stretch_level_end /= 0 ) THEN 534 DO n = 1, number_stretch_level_end 535 dz_stretch_level_end(n) = INT( dz_stretch_level_end(n) / & 536 dz(n+1) ) * dz(n+1) 537 ENDDO 538 ENDIF 539 540 ! 541 !-- Determine stretching factor if necessary 542 IF ( number_stretch_level_end >= 1 ) THEN 543 CALL calculate_stretching_factor( number_stretch_level_end ) 544 ENDIF 545 404 546 ! 405 547 !-- Grid for atmosphere with surface at z=0 (k=0, w-grid). 548 !-- First compute the u- and v-levels. In case of dirichlet bc for u and v 549 !-- the first u/v- and w-level (k=0) are defined at same height (z=0). 406 550 !-- The second u-level (k=1) corresponds to the top of the 407 551 !-- Prandtl-layer. 408 409 552 IF ( ibc_uv_b == 0 .OR. ibc_uv_b == 2 ) THEN 410 553 zu(0) = 0.0_wp 411 ! zu(0) = - dz * 0.5_wp412 554 ELSE 413 zu(0) = - dz * 0.5_wp 414 ENDIF 415 zu(1) = dz * 0.5_wp 416 417 dz_stretch_level_index = nzt+1 418 dz_stretched = dz 555 zu(0) = - dz(1) * 0.5_wp 556 ENDIF 557 558 zu(1) = dz(1) * 0.5_wp 559 560 ! 561 !-- Determine u and v height levels considering the possibility of grid 562 !-- stretching in several heights. 563 n = 1 564 dz_stretch_level_start_index = nzt+1 565 dz_stretch_level_end_index = nzt+1 566 dz_stretched = dz(1) 567 568 !-- The default value of dz_stretch_level_start is negative, thus the first 569 !-- condition is always true. Hence, the second condition is necessary. 419 570 DO k = 2, nzt+1 420 IF ( dz_stretch_level <= zu(k-1) .AND. dz_stretched < dz_max ) THEN 421 dz_stretched = dz_stretched * dz_stretch_factor 422 dz_stretched = MIN( dz_stretched, dz_max ) 423 IF ( dz_stretch_level_index == nzt+1 ) dz_stretch_level_index = k-1 424 ENDIF 571 IF ( dz_stretch_level_start(n) <= zu(k-1) .AND. & 572 dz_stretch_level_start(n) /= -9999999.9_wp ) THEN 573 dz_stretched = dz_stretched * dz_stretch_factor_array(n) 574 575 IF ( dz(n) > dz(n+1) ) THEN 576 dz_stretched = MAX( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (higher) dz 577 ELSE 578 dz_stretched = MIN( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (lower) dz 579 ENDIF 580 581 IF ( dz_stretch_level_start_index(n) == nzt+1 ) & 582 dz_stretch_level_start_index(n) = k-1 583 584 ENDIF 585 425 586 zu(k) = zu(k-1) + dz_stretched 587 588 ! 589 !-- Make sure that the stretching ends exactly at dz_stretch_level_end 590 dz_level_end = ABS( zu(k) - dz_stretch_level_end(n) ) 591 592 IF ( dz_level_end < dz(n+1)/3.0 ) THEN 593 zu(k) = dz_stretch_level_end(n) 594 dz_stretched = dz(n+1) 595 dz_stretch_level_end_index(n) = k 596 n = n + 1 597 ENDIF 426 598 ENDDO 427 599 … … 438 610 439 611 ELSE 612 613 ! 614 !-- The stretching region has to be large enough to allow for a smooth 615 !-- transition between two different grid spacings 616 DO n = 1, number_stretch_level_start 617 min_dz_stretch_level_end(n) = dz_stretch_level_start(n) - & 618 4 * MAX( dz(n),dz(n+1) ) 619 ENDDO 620 621 IF ( ANY( min_dz_stretch_level_end < dz_stretch_level_end ) ) THEN 622 message_string= 'Eeach dz_stretch_level_end has to be less ' // & 623 'than its corresponding value for &' // & 624 'dz_stretch_level_start - 4*MAX(dz(n),dz(n+1)) '//& 625 'to allow for smooth grid stretching' 626 CALL message( 'init_grid', 'PA0224', 1, 2, 0, 6, 0 ) 627 ENDIF 628 629 ! 630 !-- Stretching must not be applied within the prandtl_layer 631 !-- (first two grid points). For the default case dz_stretch_level_start 632 !-- is negative. Therefore the absolut value is checked here. 633 IF ( ANY( dz_stretch_level_start > dz(1) * 1.5_wp ) ) THEN 634 WRITE( message_string, * ) 'Eeach dz_stretch_level_start has to be ',& 635 'less than ', dz(1) * 1.5 636 CALL message( 'init_grid', 'PA0226', 1, 2, 0, 6, 0 ) 637 ENDIF 638 639 ! 640 !-- The stretching has to start and end on a grid level. Therefore 641 !-- user-specified values have to ''interpolate'' to the next highest level 642 IF ( number_stretch_level_start /= 0 ) THEN 643 dz_stretch_level_start(1) = INT( (dz_stretch_level_start(1) + & 644 dz(1)/2.0) / dz(1) ) & 645 * dz(1) - dz(1)/2.0 646 ENDIF 647 648 IF ( number_stretch_level_start > 1 ) THEN 649 DO n = 2, number_stretch_level_start 650 dz_stretch_level_start(n) = INT( dz_stretch_level_start(n) / & 651 dz(n) ) * dz(n) 652 ENDDO 653 ENDIF 654 655 IF ( number_stretch_level_end /= 0 ) THEN 656 DO n = 1, number_stretch_level_end 657 dz_stretch_level_end(n) = INT( dz_stretch_level_end(n) / & 658 dz(n+1) ) * dz(n+1) 659 ENDDO 660 ENDIF 661 662 ! 663 !-- Determine stretching factor if necessary 664 IF ( number_stretch_level_end >= 1 ) THEN 665 CALL calculate_stretching_factor( number_stretch_level_end ) 666 ENDIF 667 440 668 ! 441 669 !-- Grid for ocean with free water surface is at k=nzt (w-grid). … … 444 672 !-- w-level are defined at same height, but staggered from the second level. 445 673 !-- The second u-level (k=1) corresponds to the top of the Prandtl-layer. 446 zu(nzt+1) = dz * 0.5_wp 447 zu(nzt) = - dz * 0.5_wp 448 449 dz_stretch_level_index = 0 450 dz_stretched = dz 674 !-- z values are negative starting from z=0 (surface) 675 zu(nzt+1) = dz(1) * 0.5_wp 676 zu(nzt) = - dz(1) * 0.5_wp 677 678 ! 679 !-- Determine u and v height levels considering the possibility of grid 680 !-- stretching in several heights. 681 n = 1 682 dz_stretch_level_start_index = 0 683 dz_stretch_level_end_index = 0 684 dz_stretched = dz(1) 685 451 686 DO k = nzt-1, 0, -1 452 ! 453 !-- The default value of dz_stretch_level is positive, thus the first 454 !-- condition is always true. Hence, the second condition is necessary. 455 IF ( dz_stretch_level >= zu(k+1) .AND. dz_stretch_level <= 0.0 & 456 .AND. dz_stretched < dz_max ) THEN 457 dz_stretched = dz_stretched * dz_stretch_factor 458 dz_stretched = MIN( dz_stretched, dz_max ) 459 IF ( dz_stretch_level_index == 0 ) dz_stretch_level_index = k+1 460 ENDIF 687 688 IF ( dz_stretch_level_start(n) >= zu(k+1) ) THEN 689 dz_stretched = dz_stretched * dz_stretch_factor_array(n) 690 691 IF ( dz(n) > dz(n+1) ) THEN 692 dz_stretched = MAX( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (higher) dz 693 ELSE 694 dz_stretched = MIN( dz_stretched, dz(n+1) ) !Restrict dz_stretched to the user-specified (lower) dz 695 ENDIF 696 697 IF ( dz_stretch_level_start_index(n) == 0 ) & 698 dz_stretch_level_start_index(n) = k+1 699 700 ENDIF 701 461 702 zu(k) = zu(k+1) - dz_stretched 703 704 ! 705 !-- Make sure that the stretching ends exactly at dz_stretch_level_end 706 dz_level_end = ABS( zu(k) - dz_stretch_level_end(n) ) 707 708 IF ( dz_level_end < dz(n+1)/3.0 ) THEN 709 zu(k) = dz_stretch_level_end(n) 710 dz_stretched = dz(n+1) 711 dz_stretch_level_end_index(n) = k 712 n = n + 1 713 ENDIF 462 714 ENDDO 463 715 464 716 ! 465 717 !-- Compute the w-levels. They are always staggered half-way between the … … 468 720 !-- same height. The top w-level (nzt+1) is not used but set for 469 721 !-- consistency, since w and all scalar variables are defined up tp nzt+1. 470 zw(nzt+1) = dz 722 zw(nzt+1) = dz(1) 471 723 zw(nzt) = 0.0_wp 472 724 DO k = 0, nzt … … 855 1107 END SUBROUTINE init_grid 856 1108 1109 1110 ! Description: 1111 ! -----------------------------------------------------------------------------! 1112 !> Calculation of the stretching factor through an iterative method. Ideas were 1113 !> taken from the paper "Regional stretched grid generation and its application 1114 !> to the NCAR RegCM (1999)". Normally, no analytic solution exists because the 1115 !> system of equations has two variables (r,l) but four requirements 1116 !> (l=integer, r=[0,88;1,2], Eq(6), Eq(5) starting from index j=1) which 1117 !> results into an overdetermined system. 1118 !------------------------------------------------------------------------------! 1119 SUBROUTINE calculate_stretching_factor( number_end ) 1120 1121 USE control_parameters, & 1122 ONLY: dz, dz_stretch_factor, dz_stretch_factor_array, & 1123 dz_stretch_level_end, dz_stretch_level_start, message_string 1124 1125 USE kinds 1126 1127 IMPLICIT NONE 1128 1129 INTEGER(iwp) :: iterations !< number of iterations until stretch_factor_lower/upper_limit is reached 1130 INTEGER(iwp) :: l_rounded !< after l_rounded grid levels dz(n) is strechted to dz(n+1) with stretch_factor_2 1131 INTEGER(iwp) :: n !< loop variable for stretching 1132 1133 INTEGER(iwp), INTENT(IN) :: number_end !< number of user-specified end levels for stretching 1134 1135 REAL(wp) :: delta_l !< absolute difference between l and l_rounded 1136 REAL(wp) :: delta_stretch_factor !< absolute difference between stretch_factor_1 and stretch_factor_2 1137 REAL(wp) :: delta_total_new !< sum of delta_l and delta_stretch_factor for the next iteration (should be as small as possible) 1138 REAL(wp) :: delta_total_old !< sum of delta_l and delta_stretch_factor for the last iteration 1139 REAL(wp) :: distance !< distance between dz_stretch_level_start and dz_stretch_level_end (stretching region) 1140 REAL(wp) :: l !< value that fulfil Eq. (5) in the paper mentioned above together with stretch_factor_1 exactly 1141 REAL(wp) :: numerator !< numerator of the quotient 1142 REAL(wp) :: stretch_factor_1 !< stretching factor that fulfil Eq. (5) togehter with l exactly 1143 REAL(wp) :: stretch_factor_2 !< stretching factor that fulfil Eq. (6) togehter with l_rounded exactly 1144 1145 REAL(wp), PARAMETER :: stretch_factor_interval = 1.0E-06 !< interval for sampling possible stretching factors 1146 REAL(wp), PARAMETER :: stretch_factor_lower_limit = 0.88 !< lowest possible stretching factor 1147 REAL(wp), PARAMETER :: stretch_factor_upper_limit = 1.12 !< highest possible stretching factor 1148 1149 1150 l = 0 1151 DO n = 1, number_end 1152 1153 iterations = 1 1154 stretch_factor_1 = 1.0 1155 stretch_factor_2 = 1.0 1156 delta_total_old = 1.0 1157 1158 IF ( dz(n) > dz(n+1) ) THEN 1159 DO WHILE ( stretch_factor_1 >= stretch_factor_lower_limit ) 1160 1161 stretch_factor_1 = 1.0 - iterations * stretch_factor_interval 1162 distance = ABS( dz_stretch_level_end(n) - & 1163 dz_stretch_level_start(n) ) 1164 numerator = distance*stretch_factor_1/dz(n) + & 1165 stretch_factor_1 - distance/dz(n) 1166 1167 IF ( numerator > 0.0 ) THEN 1168 l = LOG( numerator ) / LOG( stretch_factor_1 ) - 1.0 1169 l_rounded = NINT( l ) 1170 delta_l = ABS( l_rounded - l ) / l 1171 ENDIF 1172 1173 stretch_factor_2 = EXP( LOG( dz(n+1)/dz(n) ) / (l_rounded) ) 1174 1175 delta_stretch_factor = ABS( stretch_factor_1 - & 1176 stretch_factor_2 ) / & 1177 stretch_factor_2 1178 1179 delta_total_new = delta_l + delta_stretch_factor 1180 1181 ! 1182 !-- stretch_factor_1 is taken to guarantee that the stretching 1183 !-- procedure ends as close as possible to dz_stretch_level_end. 1184 !-- stretch_factor_2 would guarantee that the stretched dz(n) is 1185 !-- equal to dz(n+1) after l_rounded grid levels. 1186 IF (delta_total_new < delta_total_old) THEN 1187 dz_stretch_factor_array(n) = stretch_factor_1 1188 delta_total_old = delta_total_new 1189 ENDIF 1190 1191 iterations = iterations + 1 1192 1193 ENDDO 1194 1195 ELSEIF ( dz(n) < dz(n+1) ) THEN 1196 DO WHILE ( stretch_factor_1 <= stretch_factor_upper_limit ) 1197 1198 stretch_factor_1 = 1.0 + iterations * stretch_factor_interval 1199 distance = ABS( dz_stretch_level_end(n) - & 1200 dz_stretch_level_start(n) ) 1201 numerator = distance*stretch_factor_1/dz(n) + & 1202 stretch_factor_1 - distance/dz(n) 1203 1204 l = LOG( numerator ) / LOG( stretch_factor_1 ) - 1.0 1205 l_rounded = NINT( l ) 1206 delta_l = ABS( l_rounded - l ) / l 1207 1208 stretch_factor_2 = EXP( LOG( dz(n+1)/dz(n) ) / (l_rounded) ) 1209 1210 delta_stretch_factor = ABS( stretch_factor_1 - & 1211 stretch_factor_2 ) / & 1212 stretch_factor_2 1213 1214 delta_total_new = delta_l + delta_stretch_factor 1215 1216 ! 1217 !-- stretch_factor_1 is taken to guarantee that the stretching 1218 !-- procedure ends as close as possible to dz_stretch_level_end. 1219 !-- stretch_factor_2 would guarantee that the stretched dz(n) is 1220 !-- equal to dz(n+1) after l_rounded grid levels. 1221 IF (delta_total_new < delta_total_old) THEN 1222 dz_stretch_factor_array(n) = stretch_factor_1 1223 delta_total_old = delta_total_new 1224 ENDIF 1225 1226 iterations = iterations + 1 1227 ENDDO 1228 1229 ELSE 1230 message_string= 'Two adjacent values of dz must be different' 1231 CALL message( 'init_grid', 'PA0228', 1, 2, 0, 6, 0 ) 1232 1233 ENDIF 1234 ENDDO 1235 1236 END SUBROUTINE calculate_stretching_factor 1237 1238 857 1239 ! Description: 858 1240 ! -----------------------------------------------------------------------------! … … 1703 2085 !-- Tunnel-wall depth 1704 2086 IF ( tunnel_wall_depth == 9999999.9_wp ) THEN 1705 td = MAX ( dx, dy, dz )2087 td = MAX ( dx, dy, dz(1) ) 1706 2088 ELSE 1707 2089 td = tunnel_wall_depth -
palm/trunk/SOURCE/init_masks.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz_stretch_level was replaced by dz_stretch_level_start 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 133 136 ONLY: constant_diffusion, cloud_droplets, cloud_physics, & 134 137 data_output_masks, data_output_masks_user, & 135 doav, doav_n, domask, domask_no, dz, dz_stretch_level , humidity,&136 mask, masks, mask_scale, mask_i,&138 doav, doav_n, domask, domask_no, dz, dz_stretch_level_start, & 139 humidity, mask, masks, mask_scale, mask_i, & 137 140 mask_i_global, mask_j, mask_j_global, mask_k, mask_k_global, & 138 141 mask_loop, mask_size, mask_size_l, mask_start_l, mask_x, & … … 141 144 microphysics_morrison, microphysics_seifert, passive_scalar, & 142 145 ocean, varnamelength 146 143 147 144 148 USE grid_variables, & … … 490 494 CALL set_mask_locations( 1, dx, 'dx', nx, 'nx', nxl, nxr ) 491 495 CALL set_mask_locations( 2, dy, 'dy', ny, 'ny', nys, nyn ) 492 CALL set_mask_locations( 3, dz , 'dz', nz, 'nz', nzb, nzt )496 CALL set_mask_locations( 3, dz(1), 'dz', nz, 'nz', nzb, nzt ) 493 497 ! 494 498 !-- Set global masks along all three dimensions (required by … … 727 731 !-- The following line assumes a constant vertical grid spacing within 728 732 !-- the vertical mask range; it fails for vertical grid stretching. 729 !-- Maybe revise later. Issue warning but continue execution. 733 !-- Maybe revise later. Issue warning but continue execution. ABS(...) 734 !-- within the IF statement is necessary because the default value of 735 !-- dz_stretch_level_start is -9999999.9_wp. 730 736 loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) * ddxyz ) 731 737 732 IF ( mask_loop(mid,dim,2) * mask_scale(dim) > dz_stretch_level )&733 THEN738 IF ( mask_loop(mid,dim,2) * mask_scale(dim) > & 739 ABS( dz_stretch_level_start(1) ) ) THEN 734 740 WRITE ( message_string, '(A,I3,A,I1,A,F9.3,A,F8.2,3A)' ) & 735 741 'mask_loop(',mid,',',dim,',2)=', mask_loop(mid,dim,2), & 736 ' exceeds dz_stretch_level=',dz_stretch_level ,&742 ' exceeds dz_stretch_level=',dz_stretch_level_start(1), & 737 743 '.&Vertical mask locations will not ', & 738 744 'match the desired heights within the stretching ', & -
palm/trunk/SOURCE/lpm_advec.f90
r2969 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz values were replaced by dzw or dz(1) to allow for right vertical stretching 28 ! 29 ! 2969 2018-04-13 11:55:09Z thiele 27 30 ! Bugfix in Interpolation indices. 28 31 ! … … 158 161 159 162 USE arrays_3d, & 160 ONLY: de_dx, de_dy, de_dz, diss, e, km, u, v, w, zu, zw163 ONLY: de_dx, de_dy, de_dz, diss, dzw, e, km, u, v, w, zu, zw 161 164 162 165 USE cpulog … … 409 412 + ( gg-cc ) * u(k+1,j+1,i) + ( gg-dd ) * & 410 413 u(k+1,j+1,i+1) ) / ( 3.0_wp * gg ) - u_gtrans 411 u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *&414 u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 412 415 ( u_int_u - u_int_l ) 413 416 ENDIF … … 503 506 + ( gg-cc ) * v(k+1,j+1,i) + ( gg-dd ) * v(k+1,j+1,i+1) & 504 507 ) / ( 3.0_wp * gg ) - v_gtrans 505 v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz *&508 v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 506 509 ( v_int_u - v_int_l ) 507 510 ENDIF … … 540 543 ( gg-dd ) * w(k+1,j+1,i+1) & 541 544 ) / ( 3.0_wp * gg ) 542 w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dz *&545 w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dzw(k) * & 543 546 ( w_int_u - w_int_l ) 544 547 ENDIF … … 611 614 ( gg - dd ) * e(k+1,j+1,i+1) & 612 615 ) / ( 3.0_wp * gg ) 613 e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dz *&616 e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 614 617 ( e_int_u - e_int_l ) 615 618 ENDIF … … 637 640 ( gg - dd ) * de_dx(k+1,j+1,i+1) & 638 641 ) / ( 3.0_wp * gg ) 639 de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dz *&642 de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 640 643 ( de_dx_int_u - de_dx_int_l ) 641 644 ENDIF … … 655 658 ( gg - dd ) * de_dy(k+1,j+1,i+1) & 656 659 ) / ( 3.0_wp * gg ) 657 de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dz * &660 de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 658 661 ( de_dy_int_u - de_dy_int_l ) 659 662 ENDIF … … 661 664 ! 662 665 !-- Interpolate the TKE gradient along z 663 IF ( zv(n) < 0.5_wp * dz ) THEN666 IF ( zv(n) < 0.5_wp * dz(1) ) THEN 664 667 de_dz_int(n) = 0.0_wp 665 668 ELSE … … 678 681 ( gg - dd ) * de_dz(k+1,j+1,i+1) & 679 682 ) / ( 3.0_wp * gg ) 680 de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dz * &683 de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 681 684 ( de_dz_int_u - de_dz_int_l ) 682 685 ENDIF … … 699 702 ( gg - dd ) * diss(k+1,j+1,i+1) & 700 703 ) / ( 3.0_wp * gg ) 701 diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dz *&704 diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 702 705 ( diss_int_u - diss_int_l ) 703 706 ENDIF -
palm/trunk/SOURCE/lpm_exchange_horiz.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dz(1) to allow for right vertical stretching 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 910 913 ip = particle_array(n)%x * ddx 911 914 jp = particle_array(n)%y * ddy 912 kp = particle_array(n)%z / dz + 1 + offset_ocean_nzt915 kp = particle_array(n)%z / dz(1) + 1 + offset_ocean_nzt 913 916 ! 914 917 !-- In case of grid stretching a particle might be above or below the -
palm/trunk/SOURCE/lpm_init.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dzw or dz(1) to allow for right vertical stretching 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 794 797 ip = tmp_particle%x * ddx 795 798 jp = tmp_particle%y * ddy 796 kp = tmp_particle%z / dz + 1 + offset_ocean_nzt799 kp = tmp_particle%z / dz(1) + 1 + offset_ocean_nzt 797 800 DO WHILE( zw(kp) < tmp_particle%z ) 798 801 kp = kp + 1 … … 969 972 pdz(particles(n)%group) 970 973 particles(n)%z = particles(n)%z + & 971 MERGE( rand_contr, SIGN( dz , rand_contr ),&972 ABS( rand_contr ) < dz 974 MERGE( rand_contr, SIGN( dzw(kp), rand_contr ), & 975 ABS( rand_contr ) < dzw(kp) & 973 976 ) 974 977 ENDIF … … 987 990 i = particles(n)%x * ddx 988 991 j = particles(n)%y * ddy 989 k = particles(n)%z / dz + 1 + offset_ocean_nzt992 k = particles(n)%z / dz(1) + 1 + offset_ocean_nzt 990 993 DO WHILE( zw(k) < particles(n)%z ) 991 994 k = k + 1 -
palm/trunk/SOURCE/lpm_set_attributes.f90
r2718 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dzw to allow for right vertical stretching 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 80 83 81 84 USE arrays_3d, & 82 ONLY: pt, u, v, w, zu, zw85 ONLY: dzw, pt, u, v, w, zu, zw 83 86 84 87 USE control_parameters, & 85 ONLY: u_gtrans, v_gtrans , dz88 ONLY: u_gtrans, v_gtrans 86 89 87 90 USE cpulog, & … … 207 210 ( gg - dd ) * u(k+1,j+1,i+1) & 208 211 ) / ( 3.0_wp * gg ) - u_gtrans 209 u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *&212 u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 210 213 ( u_int_u - u_int_l ) 211 214 ENDIF … … 240 243 ( gg - dd ) * v(k+1,j+1,i+1) & 241 244 ) / ( 3.0_wp * gg ) - v_gtrans 242 v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz *&245 v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 243 246 ( v_int_u - v_int_l ) 244 247 ENDIF … … 344 347 ) / ( 3.0_wp * gg ) - sums(k,4) 345 348 346 pt_int = pt_int_l + ( zv(n) - zu(k) ) / dz *&349 pt_int = pt_int_l + ( zv(n) - zu(k) ) / dzw(k) * & 347 350 ( pt_int_u - pt_int_l ) 348 351 -
palm/trunk/SOURCE/modules.f90
r3045 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Variables concerning stretching introduced or revised 28 ! 29 ! 3045 2018-05-28 07:55:41Z Giersch 27 30 ! z_max_do2d removed 28 31 ! … … 1148 1151 INTEGER(iwp) :: dp_level_ind_b = 0 !< lowest grid index for external pressure gradient forcing 1149 1152 INTEGER(iwp) :: dvrp_filecount = 0 !< parameter for dvr visualization software 1150 INTEGER(iwp) :: dz_stretch_level_index !< vertical grid level index above which the vertical grid spacing is stretched1151 1153 INTEGER(iwp) :: ensemble_member_nr = 0 !< namelist parameter 1152 1154 INTEGER(iwp) :: gamma_mg !< switch for steering the multigrid cycle: 1: v-cycle, 2: w-cycle … … 1190 1192 INTEGER(iwp) :: num_var_fl !< number of sampling/output variables in virtual flight measurements 1191 1193 INTEGER(iwp) :: num_var_fl_user=0 !< number of user-defined sampling/output variables in virtual flight measurements 1194 INTEGER(iwp) :: number_stretch_level_start !< number of user-specified start levels for stretching 1195 INTEGER(iwp) :: number_stretch_level_end !< number of user-specified end levels for stretching 1192 1196 INTEGER(iwp) :: nz_do3d = -9999 !< namelist parameter 1193 1197 INTEGER(iwp) :: prt_time_count = 0 !< number of output intervals for particle data output … … 1199 1203 INTEGER(iwp) :: timestep_count = 0 !< number of timesteps carried out since the beginning of the initial run 1200 1204 INTEGER(iwp) :: y_shift = 0 !< namelist parameter 1205 1201 1206 INTEGER(iwp) :: dist_nxl(0:1) !< left boundary of disturbance region 1202 1207 INTEGER(iwp) :: dist_nxr(0:1) !< right boundary of disturbance region … … 1211 1216 INTEGER(iwp) :: domask_no(max_masks,0:1) = 0 !< number of masked output quantities 1212 1217 INTEGER(iwp) :: domask_time_count(max_masks,0:1) !< number of output intervals for masked data 1218 INTEGER(iwp) :: dz_stretch_level_end_index(9) !< vertical grid level index until which the vertical grid spacing is stretched 1219 INTEGER(iwp) :: dz_stretch_level_start_index(9) !< vertical grid level index above which the vertical grid spacing is stretched 1213 1220 INTEGER(iwp) :: mask_size(max_masks,3) = -1 !< size of mask array per mask and dimension (for netcdf output) 1214 1221 INTEGER(iwp) :: mask_size_l(max_masks,3) = -1 !< subdomain size of mask array per mask and dimension (for netcdf output) … … 1416 1423 REAL(wp) :: dt_spinup = 60.0_wp !< namelist parameter 1417 1424 REAL(wp) :: dt_3d = 1.0_wp !< time step 1418 REAL(wp) :: dz = -1.0_wp !< namelist parameter 1419 REAL(wp) :: dz_max = 9999999.9_wp !< namelist parameter 1425 REAL(wp) :: dz_max = 1000.0_wp !< namelist parameter 1420 1426 REAL(wp) :: dz_stretch_factor = 1.08_wp !< namelist parameter 1421 REAL(wp) :: dz_stretch_level = 100000.0_wp!< namelist parameter1427 REAL(wp) :: dz_stretch_level = -9999999.9_wp !< namelist parameter 1422 1428 REAL(wp) :: e_init = 0.0_wp !< namelist parameter 1423 1429 REAL(wp) :: e_min = 0.0_wp !< namelist parameter … … 1529 1535 REAL(wp) :: dpdxy(1:2) = 0.0_wp !< namelist parameter 1530 1536 REAL(wp) :: dt_domask(max_masks) = 9999999.9_wp !< namelist parameter 1537 REAL(wp) :: dz(10) = -1.0_wp !< namelist parameter 1538 REAL(wp) :: dz_stretch_level_start(9) = -9999999.9_wp !< namelist parameter 1539 REAL(wp) :: dz_stretch_level_end(9) = 9999999.9_wp !< namelist parameter 1540 REAL(wp) :: dz_stretch_factor_array(9) = 1.08_wp !< namelist parameter 1531 1541 REAL(wp) :: mask_scale(3) !< collective array for mask_scale_x/y/z 1532 1542 REAL(wp) :: pt_vertical_gradient(10) = 0.0_wp !< namelist parameter -
palm/trunk/SOURCE/parin.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! New initialization parameters added 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 509 512 dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius, & 510 513 dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max, & 511 dz_stretch_factor, dz_stretch_level, end_time_1d,&512 ensemble_member_nr, e_init, e_min, fft_method,&513 flux_input_mode, flux_output_mode, forcing,&514 dz_stretch_factor, dz_stretch_level, dz_stretch_level_start, & 515 dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init, & 516 e_min, fft_method, flux_input_mode, flux_output_mode, forcing, & 514 517 galilei_transformation, humidity, & 515 518 inflow_damping_height, inflow_damping_width, & … … 579 582 dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius, & 580 583 dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max, & 581 dz_stretch_factor, dz_stretch_level, end_time_1d,&582 ensemble_member_nr, e_init, e_min, fft_method,&583 flux_input_mode, flux_output_mode, forcing,&584 dz_stretch_factor, dz_stretch_level, dz_stretch_level_start, & 585 dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init, & 586 e_min, fft_method, flux_input_mode, flux_output_mode, forcing, & 584 587 galilei_transformation, humidity, & 585 588 inflow_damping_height, inflow_damping_width, & -
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by the help of zw to allow for vertical stretching 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 550 553 551 554 USE control_parameters, & 552 ONLY: dz,passive_scalar555 ONLY: passive_scalar 553 556 554 557 … … 568 571 REAL(wp) :: canopy_height !< canopy height (in m) 569 572 570 canopy_height = pch_index * dz573 canopy_height = zw(pch_index) 571 574 572 575 WRITE ( io, 1 ) canopy_mode, canopy_height, pch_index, & … … 668 671 669 672 USE control_parameters, & 670 ONLY: dz, humidity, io_blocks, io_group, message_string, ocean,&673 ONLY: humidity, io_blocks, io_group, message_string, ocean, & 671 674 passive_scalar, urban_surface 672 675 … … 761 764 !-- Use beta function for lad-profile construction 762 765 int_bpdf = 0.0_wp 763 canopy_height = pch_index * dz766 canopy_height = zw(pch_index) 764 767 765 768 DO k = 0, pch_index -
palm/trunk/SOURCE/pmc_interface_mod.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dz(1) 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 913 916 parent_grid_info_real(5) = lower_left_coord_x + ( nx + 1 ) * dx 914 917 parent_grid_info_real(6) = lower_left_coord_y + ( ny + 1 ) * dy 915 parent_grid_info_real(7) = dz 918 parent_grid_info_real(7) = dz(1) 916 919 917 920 parent_grid_info_int(1) = nx … … 1318 1321 fval(3) = dx 1319 1322 fval(4) = dy 1320 fval(5) = dz 1323 fval(5) = dz(1) 1321 1324 1322 1325 IF ( myid == 0 ) THEN -
palm/trunk/SOURCE/radiation_model_mod.f90
r3049 r3065 28 28 ! ----------------- 29 29 ! $Id$ 30 ! dz was replaced by dz(1), error message concerning vertical stretching was 31 ! added 32 ! 33 ! 3049 2018-05-29 13:52:36Z Giersch 30 34 ! Error messages revised 31 35 ! … … 330 334 ! ------------ 331 335 !> Radiation models and interfaces 336 !> @todo Replace dz(1) appropriatly to account for grid stretching 332 337 !> @todo move variable definitions used in radiation_init only to the subroutine 333 338 !> as they are no longer required after initialization. … … 4371 4376 IF ( plant_canopy ) THEN 4372 4377 pchf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp & 4373 / (cp * hyp(nzub:nzut) * dx*dy*dz ) !< equals to 1 / (rho * c_p * Vbox * T)4378 / (cp * hyp(nzub:nzut) * dx*dy*dz(1)) !< equals to 1 / (rho * c_p * Vbox * T) 4374 4379 pctf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp & 4375 4380 / (l_v * hyp(nzub:nzut) * dx*dy*dz) … … 4396 4401 !-- each dimension, so that this new direction vector will allow us 4397 4402 !-- to traverse the ray path within grid coordinates directly 4398 sunorig_grid = (/ sunorig(1)/dz , sunorig(2)/dy, sunorig(3)/dx /)4403 sunorig_grid = (/ sunorig(1)/dz(1), sunorig(2)/dy, sunorig(3)/dx /) 4399 4404 !-- sunorig_grid = sunorig_grid / norm2(sunorig_grid) 4400 4405 sunorig_grid = sunorig_grid / SQRT(SUM(sunorig_grid**2)) … … 4403 4408 !-- precompute effective box depth with prototype Leaf Area Density 4404 4409 pc_box_dimshift = MAXLOC(ABS(sunorig), 1) - 1 4405 CALL box_absorb(CSHIFT((/dz ,dy,dx/), pc_box_dimshift), &4410 CALL box_absorb(CSHIFT((/dz(1),dy,dx/), pc_box_dimshift), & 4406 4411 60, prototype_lad, & 4407 4412 CSHIFT(ABS(sunorig), pc_box_dimshift), & … … 4801 4806 IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx 4802 4807 IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy 4803 IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz 4808 IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1) 4804 4809 ENDDO 4805 4810 ! … … 5014 5019 SUBROUTINE radiation_interaction_init 5015 5020 5021 USE control_parameters, & 5022 ONLY: dz_stretch_level_start 5023 5016 5024 USE netcdf_data_input_mod, & 5017 5025 ONLY: leaf_area_density_f … … 5095 5103 nzpt = nzptl 5096 5104 #endif 5105 ! 5106 !-- Stretching (non-uniform grid spacing) is not considered in the radiation 5107 !-- model. Therefore, vertical stretching has to be applied above the area 5108 !-- where the parts of the radiation model which assume constant grid spacing 5109 !-- are active. ABS (...) is required because the default value of 5110 !-- dz_stretch_level_start is -9999999.9_wp (negative). 5111 IF ( ABS( dz_stretch_level_start(1) ) <= zw(nzut) ) THEN 5112 WRITE( message_string, * ) 'The lowest level where vertical ', & 5113 'stretching is applied have to be ', & 5114 'greater than ', zw(nzut) 5115 CALL message( 'init_grid', 'PA0496', 1, 2, 0, 6, 0 ) 5116 ENDIF 5097 5117 ! 5098 5118 !-- global number of urban and plant layers … … 5390 5410 IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx 5391 5411 IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy 5392 IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz 5412 IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1) 5393 5413 ENDDO 5394 5414 … … 5512 5532 5513 5533 !-- unit vector source -> target 5514 uv = (/ (ta(1)-sa(1))*dz , (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)5534 uv = (/ (ta(1)-sa(1))*dz(1), (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /) 5515 5535 sqdist = SUM(uv(:)**2) 5516 5536 uv = uv / SQRT(sqdist) … … 5677 5697 5678 5698 !-- unit vector source -> target 5679 uv = (/ (ta(1)-sa(1))*dz , (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /)5699 uv = (/ (ta(1)-sa(1))*dz(1), (ta(2)-sa(2))*dy, (ta(3)-sa(3))*dx /) 5680 5700 sqdist = SUM(uv(:)**2) 5681 5701 uv = uv / SQRT(sqdist) … … 6127 6147 ENDIF 6128 6148 uvect(:) = delta(:) / distance 6129 realdist = SQRT(SUM( (uvect(:)*(/dz ,dy,dx/))**2 ))6149 realdist = SQRT(SUM( (uvect(:)*(/dz(1),dy,dx/))**2 )) 6130 6150 6131 6151 lastdist = 0._wp … … 6509 6529 zb1 = CEILING(zexit * zsgn - .5_wp) - 1 ! because it must be smaller than exit 6510 6530 nz = MAX(zb1 - zb0 + 3, 2) 6511 rt2_dist(nz) = SQRT(((zexit-zorig)*dz )**2 + dxxyy)6531 rt2_dist(nz) = SQRT(((zexit-zorig)*dz(1))**2 + dxxyy) 6512 6532 qdist = rt2_dist(nz) / (zexit-zorig) 6513 6533 rt2_dist(2:nz-1) = (/( ((REAL(l, wp) + .5_wp) * zsgn - zorig) * qdist , l = zb0, zb1 )/) … … 6567 6587 zb1 = CEILING(zexit * zsgn - .5_wp) - 1 ! because it must be smaller than exit 6568 6588 nz = MAX(zb1 - zb0 + 3, 2) 6569 rt2_dist(nz) = SQRT(((zexit-zorig)*dz )**2 + dxxyy)6589 rt2_dist(nz) = SQRT(((zexit-zorig)*dz(1))**2 + dxxyy) 6570 6590 qdist = rt2_dist(nz) / (zexit-zorig) 6571 6591 rt2_dist(2:nz-1) = (/( ((REAL(l, wp) + .5_wp) * zsgn - zorig) * qdist , l = zb0, zb1 )/) -
palm/trunk/SOURCE/read_restart_data_mod.f90
r3056 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! New parameters concerning vertical grid stretching have been added 28 ! 29 ! 3056 2018-06-04 07:49:35Z Giersch 27 30 ! found variable has to be set to false inside overlap loop 28 31 ! … … 407 410 CASE ( 'dz_stretch_factor' ) 408 411 READ ( 13 ) dz_stretch_factor 412 CASE ( 'dz_stretch_factor_array' ) 413 READ ( 13 ) dz_stretch_factor_array 409 414 CASE ( 'dz_stretch_level' ) 410 415 READ ( 13 ) dz_stretch_level 416 CASE ( 'dz_stretch_level_end' ) 417 READ ( 13 ) dz_stretch_level_end 418 CASE ( 'dz_stretch_level_start' ) 419 READ ( 13 ) dz_stretch_level_start 411 420 CASE ( 'e_min' ) 412 421 READ ( 13 ) e_min … … 1251 1260 THEN 1252 1261 WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed ', & 1253 'in restart run PE', myid, ' will read from files ', &1262 'in restart run & PE', myid, ' will read from files ', & 1254 1263 file_list(1:files_to_be_opened) 1255 1264 CALL message( 'rrd_local', 'PA0285', 0, 0, 0, 6, 0 ) -
palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r3051 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Error message related to vertical stretching has been added, dz was replaced 28 ! by dz(1) 29 ! 30 ! 3051 2018-05-30 17:43:55Z suehring 27 31 ! Bugfix in calculation of initial Reynolds-stress tensor. 28 32 ! … … 313 317 314 318 USE control_parameters, & 315 ONLY: bc_lr, bc_ns, forcing, nest_domain, rans_mode, turbulent_inflow 319 ONLY: bc_lr, bc_ns, forcing, nest_domain, number_stretch_level_start, & 320 rans_mode, turbulent_inflow 316 321 317 322 USE pmc_interface, & … … 365 370 'is not allowed' 366 371 CALL message( 'stg_check_parameters', 'PA0039', 1, 2, 0, 6, 0 ) 372 ENDIF 373 374 IF ( number_stretch_level_start > 0 ) THEN 375 message_string = 'Using synthetic turbulence generator ' // & 376 'in combination with stretching is not allowed' 377 CALL message( 'stg_check_parameters', 'PA0420', 1, 2, 0, 6, 0 ) 367 378 ENDIF 368 379 … … 607 618 608 619 ! 609 !-- Convert length scales from meter to number of grid points 620 !-- Convert length scales from meter to number of grid points. Attention: 621 !-- Does not work if grid stretching is used 610 622 nuy(k) = INT( luy * ddy ) 611 nuz(k) = INT( luz / dz )623 nuz(k) = INT( luz / dz(1) ) 612 624 nvy(k) = INT( lvy * ddy ) 613 nvz(k) = INT( lvz / dz )625 nvz(k) = INT( lvz / dz(1) ) 614 626 nwy(k) = INT( lwy * ddy ) 615 nwz(k) = INT( lwz / dz )627 nwz(k) = INT( lwz / dz(1) ) 616 628 ! 617 629 !-- Workaround, assume isotropic turbulence -
palm/trunk/SOURCE/urban_surface_mod.f90
r3049 r3065 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Unused array dxdir was removed, dz was replaced by dzu to consider vertical 31 ! grid stretching 32 ! 33 ! 3049 2018-05-29 13:52:36Z Giersch 30 34 ! Error messages revised 31 35 ! … … 282 286 #if ! defined( __nopointer ) 283 287 USE arrays_3d, & 284 ONLY: hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend288 ONLY: dzu, hyp, zu, pt, pt_1, pt_2, p, u, v, w, hyp, tend 285 289 #endif 286 290 … … 292 296 293 297 USE control_parameters, & 294 ONLY: coupling_start_time, dz, topography, dt_3d,&298 ONLY: coupling_start_time, topography, dt_3d, & 295 299 intermediate_timestep_count, initializing_actions, & 296 300 intermediate_timestep_count_max, simulated_time, end_time, & … … 7207 7211 REAL(wp), DIMENSION(nzb:nzt) :: exn !< value of the Exner function in layers 7208 7212 7209 REAL(wp), DIMENSION(0:4) :: dxdir !< surface normal direction gridbox length7210 7213 REAL(wp) :: dtime !< simulated time of day (in UTC) 7211 7214 INTEGER(iwp) :: dhour !< simulated hour of day (in UTC) … … 7213 7216 7214 7217 7215 dxdir = (/dz,dy,dy,dx,dx/)7216 7218 #if ! defined( __nopointer ) 7217 7219 exn(nzb:nzt) = (hyp(nzb:nzt) / 100000.0_wp )**0.286_wp !< Exner function … … 7634 7636 (dtime/3600.0_wp-REAL(dhour,wp))*aheatprof(k,dhour+1) 7635 7637 IF ( aheat(k,j,i) > 0.0_wp ) THEN 7636 pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dz )7638 pt(k,j,i) = pt(k,j,i) + aheat(k,j,i)*acoef*dt_3d/(exn(k)*rho_cp*dzu(k)) 7637 7639 ENDIF 7638 7640 ENDIF -
palm/trunk/SOURCE/user_init_grid.f90
r2718 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! dz was replaced by dz(1) 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 108 111 !-- topography, bit is 1 for atmospheric grid point. 109 112 !-- The following example shows how to prescribe sine-like topography 110 !-- along x-direction with amplitude of 10 * dz and wavelength 10 * dy.113 !-- along x-direction with amplitude of 10 * dz(1) and wavelength 10 * dy. 111 114 ! DO i = nxlg, nxrg 112 ! h_topo = 10.0_wp * dz * (SIN(3.14_wp*0.5_wp)*i*dx / ( 5.0_wp * dy ) )**2115 ! h_topo = 10.0_wp * dz(1) * (SIN(3.14_wp*0.5_wp)*i*dx / ( 5.0_wp * dy ) )**2 113 116 ! 114 117 ! k_topo = MINLOC( ABS( zw - h_topo ), 1 ) - 1 -
palm/trunk/SOURCE/vertical_nesting_mod.f90
r3049 r3065 26 26 ! ----------------- 27 27 ! $Id$ 28 ! dz was replaced by dz(1), error messages related to vertical grid stretching 29 ! have been added 30 ! 31 ! 3049 2018-05-29 13:52:36Z Giersch 28 32 ! Error messages revised 29 33 ! … … 73 77 !> after spin-up of the CG 74 78 !> 79 !> @todo Replace dz(1) appropriatly to account for grid stretching 75 80 !> @todo Ensure that code can be compiled for serial and parallel mode. Please 76 81 !> check the placement of the directive "__parallel". … … 3639 3644 3640 3645 USE control_parameters, & 3641 ONLY: coupling_mode, coupling_mode_remote, coupling_topology, dz 3646 ONLY: coupling_mode, coupling_mode_remote, coupling_topology, dz, & 3647 dz_stretch_level_start, message_string 3642 3648 3643 3649 USE grid_variables, & … … 3678 3684 dxc = dx 3679 3685 dyc = dy 3680 dzc = dz 3686 dzc = dz(1) 3681 3687 cg_nprocs = numprocs 3682 3688 3683 3689 IF ( myid == 0 ) THEN 3684 3690 3685 CALL MPI_SEND( nxc, 1, MPI_INTEGER , numprocs, 1, comm_inter, 3691 CALL MPI_SEND( nxc, 1, MPI_INTEGER , numprocs, 1, comm_inter, & 3686 3692 ierr ) 3687 CALL MPI_SEND( nyc, 1, MPI_INTEGER , numprocs, 2, comm_inter, 3693 CALL MPI_SEND( nyc, 1, MPI_INTEGER , numprocs, 2, comm_inter, & 3688 3694 ierr ) 3689 CALL MPI_SEND( nzc, 1, MPI_INTEGER , numprocs, 3, comm_inter, 3695 CALL MPI_SEND( nzc, 1, MPI_INTEGER , numprocs, 3, comm_inter, & 3690 3696 ierr ) 3691 CALL MPI_SEND( dxc, 1, MPI_REAL , numprocs, 4, comm_inter, 3697 CALL MPI_SEND( dxc, 1, MPI_REAL , numprocs, 4, comm_inter, & 3692 3698 ierr ) 3693 CALL MPI_SEND( dyc, 1, MPI_REAL , numprocs, 5, comm_inter, 3699 CALL MPI_SEND( dyc, 1, MPI_REAL , numprocs, 5, comm_inter, & 3694 3700 ierr ) 3695 CALL MPI_SEND( dzc, 1, MPI_REAL , numprocs, 6, comm_inter, 3701 CALL MPI_SEND( dzc, 1, MPI_REAL , numprocs, 6, comm_inter, & 3696 3702 ierr ) 3697 CALL MPI_SEND( pdims, 2, MPI_INTEGER, numprocs, 7, comm_inter, 3703 CALL MPI_SEND( pdims, 2, MPI_INTEGER, numprocs, 7, comm_inter, & 3698 3704 ierr ) 3699 3705 CALL MPI_SEND( cg_nprocs, 1, MPI_INTEGER, numprocs, 8, comm_inter, & 3700 3706 ierr ) 3701 CALL MPI_RECV( nxf, 1, MPI_INTEGER, numprocs, 21, comm_inter, 3707 CALL MPI_RECV( nxf, 1, MPI_INTEGER, numprocs, 21, comm_inter, & 3702 3708 status, ierr ) 3703 CALL MPI_RECV( nyf, 1, MPI_INTEGER, numprocs, 22, comm_inter, 3709 CALL MPI_RECV( nyf, 1, MPI_INTEGER, numprocs, 22, comm_inter, & 3704 3710 status, ierr ) 3705 CALL MPI_RECV( nzf, 1, MPI_INTEGER, numprocs, 23, comm_inter, 3711 CALL MPI_RECV( nzf, 1, MPI_INTEGER, numprocs, 23, comm_inter, & 3706 3712 status, ierr ) 3707 CALL MPI_RECV( dxf, 1, MPI_REAL, numprocs, 24, comm_inter, 3713 CALL MPI_RECV( dxf, 1, MPI_REAL, numprocs, 24, comm_inter, & 3708 3714 status, ierr ) 3709 CALL MPI_RECV( dyf, 1, MPI_REAL, numprocs, 25, comm_inter, 3715 CALL MPI_RECV( dyf, 1, MPI_REAL, numprocs, 25, comm_inter, & 3710 3716 status, ierr ) 3711 CALL MPI_RECV( dzf, 1, MPI_REAL, numprocs, 26, comm_inter, 3717 CALL MPI_RECV( dzf, 1, MPI_REAL, numprocs, 26, comm_inter, & 3712 3718 status, ierr ) 3713 CALL MPI_RECV( pdims_partner, 2, MPI_INTEGER, 3719 CALL MPI_RECV( pdims_partner, 2, MPI_INTEGER, & 3714 3720 numprocs, 27, comm_inter, status, ierr ) 3715 CALL MPI_RECV( fg_nprocs, 1, MPI_INTEGER, &3721 CALL MPI_RECV( fg_nprocs, 1, MPI_INTEGER, & 3716 3722 numprocs, 28, comm_inter, status, ierr ) 3717 3723 ENDIF … … 3725 3731 CALL MPI_BCAST( pdims_partner, 2, MPI_INTEGER, 0, comm2d, ierr ) 3726 3732 CALL MPI_BCAST( fg_nprocs, 1, MPI_INTEGER, 0, comm2d, ierr ) 3733 3734 ! 3735 !-- Check if stretching is used within the nested domain. ABS(...) is 3736 !-- necessary because of the default value of -9999999.9_wp (negative) 3737 IF ( ABS( dz_stretch_level_start(1) ) <= (nzf+1)*dzf ) THEN 3738 message_string = 'Stretching in the parent domain is '// & 3739 'only allowed above the nested domain' 3740 CALL message( 'vertical_nesting_mod', 'PA0497', 1, 2, 0, 6, 0 ) 3741 ENDIF 3727 3742 3728 3743 ELSEIF ( coupling_mode == 'vnested_fine' ) THEN … … 3733 3748 dxf = dx 3734 3749 dyf = dy 3735 dzf = dz 3750 dzf = dz(1) 3736 3751 fg_nprocs = numprocs 3737 3752 … … 3774 3789 3775 3790 ENDIF 3776 3791 3777 3792 ngp_c = ( nxc+1 + 2 * nbgp ) * ( nyc+1 + 2 * nbgp ) 3778 3793 ngp_f = ( nxf+1 + 2 * nbgp ) * ( nyf+1 + 2 * nbgp ) … … 3910 3925 3911 3926 #if defined( __parallel ) 3912 USE arrays_3d, 3927 USE arrays_3d, & 3913 3928 ONLY: zu, zw 3914 3929 3915 USE control_parameters, 3916 ONLY: coupling_mode 3930 USE control_parameters, & 3931 ONLY: coupling_mode, message_string, number_stretch_level_start 3917 3932 3918 USE indices, 3933 USE indices, & 3919 3934 ONLY: nzt 3920 3935 … … 3925 3940 IMPLICIT NONE 3926 3941 3927 !-- Allocate and Exchange zuc and zuf, zwc and zwf 3942 ! 3943 !-- Allocate and Exchange zuc and zuf, zwc and zwf 3928 3944 IF ( coupling_mode(1:8) == 'vnested_' ) THEN 3929 3945 … … 3932 3948 3933 3949 IF ( coupling_mode == 'vnested_crse' ) THEN 3934 3935 zuc = zu 3936 zwc = zw 3950 3951 zuc = zu 3952 zwc = zw 3953 3937 3954 IF ( myid == 0 ) THEN 3938 3955 … … 3954 3971 ELSEIF ( coupling_mode == 'vnested_fine' ) THEN 3955 3972 3956 zuf = zu 3957 zwf = zw 3973 ! 3974 !-- Check if stretching is used within the nested domain 3975 IF ( number_stretch_level_start > 0 ) THEN 3976 message_string = 'Stretching in the nested domain is not '//& 3977 'allowed' 3978 CALL message( 'vertical_nesting_mod', 'PA0498', 1, 2, 0, 6, 0 ) 3979 ENDIF 3980 3981 zuf = zu 3982 zwf = zw 3983 3958 3984 IF ( myid == 0 ) THEN 3959 3985 -
palm/trunk/SOURCE/virtual_flight_mod.f90
r3049 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! IF statement revised to consider new vertical grid stretching 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 734 737 735 738 USE control_parameters, & 736 ONLY: dz, dz_stretch_level 739 ONLY: dz, dz_stretch_level_start 737 740 738 741 USE grid_variables, & … … 777 780 !-- if flight position is above the vertical grid-stretching level. 778 781 !-- fac=1 if variable is on scalar grid level, fac=0 for w-component. 779 IF ( z_pos(l) < dz_stretch_level ) THEN780 k = ( z_pos(l) + fac * 0.5_wp * dz ) / dz782 IF ( z_pos(l) < dz_stretch_level_start(1) ) THEN 783 k = ( z_pos(l) + fac * 0.5_wp * dz(1) ) / dz(1) 781 784 ELSE 782 785 ! -
palm/trunk/SOURCE/wind_turbine_model_mod.f90
r3049 r3065 26 26 ! ----------------- 27 27 ! $Id$ 28 ! dz was replaced by dz(1), error message concerning grid stretching has been 29 ! introduced 30 ! 31 ! 3049 2018-05-29 13:52:36Z Giersch 28 32 ! Error messages revised 29 33 ! … … 132 136 !> automatically). 133 137 !> 138 !> @todo Replace dz(1) appropriatly to account for grid stretching 134 139 !> @todo Revise code according to PALM Coding Standard 135 140 !> @todo Implement ADM and ALM turbine models … … 740 745 delta_r_factor = segment_width 741 746 delta_t_factor = segment_length 742 delta_r_init = delta_r_factor * MIN( dx, dy, dz )747 delta_r_init = delta_r_factor * MIN( dx, dy, dz(1)) 743 748 744 749 DO inot = 1, nturbines … … 937 942 938 943 944 USE control_parameters, & 945 ONLY: dz_stretch_level_start 946 939 947 IMPLICIT NONE 940 948 … … 992 1000 pi**( 1.0_wp / 6.0_wp ) * eps_kernel 993 1001 ! 1002 !-- Stretching (non-uniform grid spacing) is not considered in the wind 1003 !-- turbine model. Therefore, vertical stretching has to be applied above 1004 !-- the area where the wtm is active. ABS (...) is required because the 1005 !-- default value of dz_stretch_level_start is -9999999.9_wp (negative). 1006 IF ( ABS( dz_stretch_level_start(1) ) <= MAXVAL(rcz) + MAXVAL(rr) + & 1007 eps_min) THEN 1008 WRITE( message_string, * ) 'The lowest level where vertical ', & 1009 'stretching is applied &have to be ', & 1010 'greater than ', MAXVAL(rcz) + & 1011 MAXVAL(rr) + eps_min 1012 CALL message( 'init_grid', 'PA0484', 1, 2, 0, 6, 0 ) 1013 ENDIF 1014 ! 994 1015 !-- Square of eps_min: 995 1016 eps_min2 = eps_min**2 … … 1022 1043 i_hub(inot) = INT( rcx(inot) / dx ) 1023 1044 j_hub(inot) = INT( ( rcy(inot) + 0.5_wp * dy ) / dy ) 1024 k_hub(inot) = INT( ( rcz(inot) + 0.5_wp * dz ) / dz)1045 k_hub(inot) = INT( ( rcz(inot) + 0.5_wp * dz(1) ) / dz(1) ) 1025 1046 1026 1047 ! … … 1031 1052 i_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dx ) 1032 1053 j_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dy ) 1033 k_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dz )1054 k_smear(inot) = CEILING( ( rr(inot) + eps_min ) / dz(1) ) 1034 1055 1035 1056 ENDDO … … 1121 1142 !-- surface. All points between z=0 and z=dz/s would already be contained 1122 1143 !-- in grid box 1. 1123 index_nacb(inot) = INT( ( rcz(inot) - rnac(inot) ) / dz ) + 11124 index_nact(inot) = INT( ( rcz(inot) + rnac(inot) ) / dz ) + 11144 index_nacb(inot) = INT( ( rcz(inot) - rnac(inot) ) / dz(1) ) + 1 1145 index_nact(inot) = INT( ( rcz(inot) + rnac(inot) ) / dz(1) ) + 1 1125 1146 1126 1147 ! … … 1148 1169 IF ( j == tower_s ) THEN 1149 1170 tow_cd_surf(k,j,i) = ( rcz(inot) - & 1150 ( k_hub(inot) * dz - 0.5_wp * dz ) ) *& ! extension in z-direction1171 ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*& ! extension in z-direction 1151 1172 ( ( tower_s + 1.0_wp + 0.5_wp ) * dy - & 1152 1173 ( rcy(inot) - 0.5_wp * dtow(inot) ) ) * & ! extension in y-direction … … 1154 1175 ELSEIF ( j == tower_n ) THEN 1155 1176 tow_cd_surf(k,j,i) = ( rcz(inot) - & 1156 ( k_hub(inot) * dz - 0.5_wp * dz ) ) *& ! extension in z-direction1177 ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*& ! extension in z-direction 1157 1178 ( ( rcy(inot) + 0.5_wp * dtow(inot) ) - & 1158 1179 ( tower_n + 0.5_wp ) * dy ) * & ! extension in y-direction … … 1163 1184 ELSE 1164 1185 tow_cd_surf(k,j,i) = ( rcz(inot) - & 1165 ( k_hub(inot) * dz - 0.5_wp * dz ) ) *&1186 ( k_hub(inot) * dz(1) - 0.5_wp * dz(1) ) )*& 1166 1187 dy * turb_cd_tower(inot) 1167 1188 ENDIF … … 1169 1190 !-- tower lies completely within one grid box: 1170 1191 ELSE 1171 tow_cd_surf(k,j,i) = ( rcz(inot) -&1172 ( k_hub(inot) * dz - 0.5_wp * dz ) ) *&1192 tow_cd_surf(k,j,i) = ( rcz(inot) - ( k_hub(inot) * & 1193 dz(1) - 0.5_wp * dz(1) ) ) * & 1173 1194 dtow(inot) * turb_cd_tower(inot) 1174 1195 ENDIF … … 1182 1203 !-- leftmost and rightmost grid box: 1183 1204 IF ( j == tower_s ) THEN 1184 tow_cd_surf(k,j,i) = dz * (&1205 tow_cd_surf(k,j,i) = dz(1) * ( & 1185 1206 ( tower_s + 1 + 0.5_wp ) * dy - & 1186 1207 ( rcy(inot) - 0.5_wp * dtow(inot) ) & 1187 1208 ) * turb_cd_tower(inot) 1188 1209 ELSEIF ( j == tower_n ) THEN 1189 tow_cd_surf(k,j,i) = dz * (&1210 tow_cd_surf(k,j,i) = dz(1) * ( & 1190 1211 ( rcy(inot) + 0.5_wp * dtow(inot) ) - & 1191 1212 ( tower_n + 0.5_wp ) * dy & … … 1195 1216 !-- (where tow_cd_surf = grid box area): 1196 1217 ELSE 1197 tow_cd_surf(k,j,i) = dz * dy * turb_cd_tower(inot) 1218 tow_cd_surf(k,j,i) = dz(1) * dy * & 1219 turb_cd_tower(inot) 1198 1220 ENDIF 1199 1221 ! 1200 1222 !-- tower lies completely within one grid box: 1201 1223 ELSE 1202 tow_cd_surf(k,j,i) = dz * dtow(inot) *&1224 tow_cd_surf(k,j,i) = dz(1) * dtow(inot) * & 1203 1225 turb_cd_tower(inot) 1204 1226 ENDIF ! end if larger than grid box … … 1302 1324 ENDDO ! end of loop over turbines 1303 1325 1304 tow_cd_surf = tow_cd_surf / ( dx * dy * dz )! Normalize tower drag1305 nac_cd_surf = nac_cd_surf / ( dx * dy * dz ) ! Normalize nacelle drag1326 tow_cd_surf = tow_cd_surf / ( dx * dy * dz(1) ) ! Normalize tower drag 1327 nac_cd_surf = nac_cd_surf / ( dx * dy * dz(1) ) ! Normalize nacelle drag 1306 1328 1307 1329 CALL wtm_read_blade_tables … … 1746 1768 ii = rbx(ring,rseg) * ddx 1747 1769 jj = ( rby(ring,rseg) - 0.5_wp * dy ) * ddy 1748 kk = ( rbz(ring,rseg) + 0.5_wp * dz ) / dz1770 kk = ( rbz(ring,rseg) + 0.5_wp * dz(1) ) / dz(1) 1749 1771 ! 1750 1772 !-- Interpolate only if all required information is available on … … 1776 1798 1777 1799 u_int_1_l(inot,ring,rseg) = u_int_l + & 1778 ( rbz(ring,rseg) - zu(kk) ) / dz *&1800 ( rbz(ring,rseg) - zu(kk) ) / dz(1) * & 1779 1801 ( u_int_u - u_int_l ) 1780 1802 … … 1791 1813 ii = ( rbx(ring,rseg) - 0.5_wp * dx ) * ddx 1792 1814 jj = rby(ring,rseg) * ddy 1793 kk = ( rbz(ring,rseg) + 0.5_wp * dz ) / dz1815 kk = ( rbz(ring,rseg) + 0.5_wp * dz(1) ) / dz(1) 1794 1816 ! 1795 1817 !-- Interpolate only if all required information is available on … … 1821 1843 1822 1844 v_int_1_l(inot,ring,rseg) = v_int_l + & 1823 ( rbz(ring,rseg) - zu(kk) ) / dz *&1845 ( rbz(ring,rseg) - zu(kk) ) / dz(1) * & 1824 1846 ( v_int_u - v_int_l ) 1825 1847 … … 1836 1858 ii = ( rbx(ring,rseg) - 0.5_wp * dx ) * ddx 1837 1859 jj = ( rby(ring,rseg) - 0.5_wp * dy ) * ddy 1838 kk = rbz(ring,rseg) / dz 1860 kk = rbz(ring,rseg) / dz(1) 1839 1861 ! 1840 1862 !-- Interpolate only if all required information is available on … … 1866 1888 1867 1889 w_int_1_l(inot,ring,rseg) = w_int_l + & 1868 ( rbz(ring,rseg) - zw(kk) ) / dz *&1890 ( rbz(ring,rseg) - zw(kk) ) / dz(1) * & 1869 1891 ( w_int_u - w_int_l ) 1870 1892 ELSE … … 2273 2295 dist_u_3d = ( i * dx - rbx(ring,rseg) )**2 + & 2274 2296 ( j * dy + 0.5_wp * dy - rby(ring,rseg) )**2 + & 2275 ( k * dz - 0.5_wp * dz- rbz(ring,rseg) )**22297 ( k * dz(1) - 0.5_wp * dz(1) - rbz(ring,rseg) )**2 2276 2298 dist_v_3d = ( i * dx + 0.5_wp * dx - rbx(ring,rseg) )**2 + & 2277 2299 ( j * dy - rby(ring,rseg) )**2 + & 2278 ( k * dz - 0.5_wp * dz- rbz(ring,rseg) )**22300 ( k * dz(1) - 0.5_wp * dz(1) - rbz(ring,rseg) )**2 2279 2301 dist_w_3d = ( i * dx + 0.5_wp * dx - rbx(ring,rseg) )**2 + & 2280 2302 ( j * dy + 0.5_wp * dy - rby(ring,rseg) )**2 + & 2281 ( k * dz - rbz(ring,rseg) )**22303 ( k * dz(1) - rbz(ring,rseg) )**2 2282 2304 2283 2305 ! -
palm/trunk/SOURCE/write_restart_data_mod.f90
r3004 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! New parameters concerning vertical grid stretching have been added 28 ! 29 ! 3004 2018-04-27 12:33:25Z Giersch 27 30 ! precipitation_rate_av removed 28 31 ! … … 380 383 CALL wrd_write_string( 'dz' ) 381 384 WRITE ( 14 ) dz 382 385 383 386 CALL wrd_write_string( 'dz_max' ) 384 387 WRITE ( 14 ) dz_max … … 386 389 CALL wrd_write_string( 'dz_stretch_factor' ) 387 390 WRITE ( 14 ) dz_stretch_factor 388 391 392 CALL wrd_write_string( 'dz_stretch_factor_array' ) 393 WRITE ( 14 ) dz_stretch_factor_array 394 389 395 CALL wrd_write_string( 'dz_stretch_level' ) 390 396 WRITE ( 14 ) dz_stretch_level 391 397 398 CALL wrd_write_string( 'dz_stretch_level_end' ) 399 WRITE ( 14 ) dz_stretch_level_end 400 401 CALL wrd_write_string( 'dz_stretch_level_start' ) 402 WRITE ( 14 ) dz_stretch_level_start 403 392 404 CALL wrd_write_string( 'e_min' ) 393 405 WRITE ( 14 ) e_min
Note: See TracChangeset
for help on using the changeset viewer.
