Changeset 3083 for palm/trunk
- Timestamp:
- Jun 19, 2018 2:03:12 PM (6 years ago)
- Location:
- palm/trunk
- Files:
-
- 15 edited
-
. (modified) (1 prop)
-
INSTALL/example_cbl_rc (modified) (2 diffs)
-
SOURCE (modified) (1 prop)
-
SOURCE/Makefile (modified) (2 diffs, 1 prop)
-
SOURCE/check_parameters.f90 (modified) (2 diffs, 1 prop)
-
SOURCE/data_output_mask.f90 (modified) (1 diff)
-
SOURCE/header.f90 (modified) (4 diffs)
-
SOURCE/init_3d_model.f90 (modified) (1 prop)
-
SOURCE/model_1d_mod.f90 (modified) (28 diffs)
-
SOURCE/modules.f90 (modified) (2 diffs)
-
SOURCE/parin.f90 (modified) (4 diffs)
-
SOURCE/pmc_interface_mod.f90 (modified) (1 prop)
-
SOURCE/timestep.f90 (modified) (8 diffs)
-
SOURCE/turbulence_closure_mod.f90 (modified) (83 diffs)
-
SOURCE/vertical_nesting_mod.f90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk
- Property svn:mergeinfo changed
/palm/branches/rans merged: 2919,2922,2928-2929,2960,2962,2966,2976,2982,2987-2988,2991,3008-3009,3023,3047,3050,3059,3062,3071-3082
- Property svn:mergeinfo changed
-
palm/trunk/INSTALL/example_cbl_rc
r3068 r3083 1 1 2 2 ****************************** -------------------------------------------- 3 * PALM 5.0 Rev: 30 63:30* atmosphere - run without 1D - prerun3 * PALM 5.0 Rev: 3082 * atmosphere - run without 1D - prerun 4 4 ****************************** -------------------------------------------- 5 5 6 Date: 1 2-06-18 Run: example_cbl7 Time: 15: 39:03Run-No.: 008 Run on host: imuk_defau6 Date: 19-06-18 Run: example_cbl 7 Time: 15:53:32 Run-No.: 00 8 Run on host: default 9 9 Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated 10 10 ------------------------------------------------------------------------------ … … 230 230 RUN ITER. HH:MM:SS.SS DT(E) UMAX VMAX WMAX U* W* THETA* Z_I ENERG. DISTENERG DIVOLD DIVNEW UMAX(KJI) VMAX(KJI) WMAX(KJI) ADVECX ADVECY MGCYC 231 231 --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 232 0 0 00:00:00.00 20.0000A -0.2131D -0.2288D 0.1182 0.000 1.38 0.000E+00 800. 0.104E-02 0.103E-02 0.597E-03 0.867E-13 7 7 11 11 7 19 5 26 7 0.000 0.000 0 233 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 234 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 235 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.421E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 236 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 237 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 238 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 239 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.431E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 240 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 241 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.830E-15 10 3 39 10 8 27 5 36 28 0.000 0.000 0 242 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 243 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 244 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 245 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 246 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.842E-15 10 3 0 8 0 7 5 36 28 0.000 0.000 0 247 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 248 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 249 0 17 00:05:40.00 20.0000D 0.4163 -0.4151 -0.1796 0.015 1.22 -0.773E+01 550. 0.283E-02 0.280E-02 0.672E-04 0.129E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 250 0 18 00:06:00.00 20.0000D 0.4097 -0.4189 -0.1772 0.015 1.22 -0.735E+01 550. 0.279E-02 0.276E-02 0.703E-04 0.130E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 251 0 19 00:06:20.00 20.0000D 0.3992 -0.4192 0.1802 0.016 1.22 -0.698E+01 550. 0.275E-02 0.273E-02 0.734E-04 0.131E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 252 0 20 00:06:40.00 20.0000D 0.3855 -0.4165 0.1983 0.017 1.22 -0.656E+01 550. 0.272E-02 0.269E-02 0.765E-04 0.133E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 253 0 21 00:07:00.00 20.0000D 0.3695 -0.4109 0.2198 0.018 1.22 -0.616E+01 550. 0.269E-02 0.267E-02 0.795E-04 0.136E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 254 0 22 00:07:20.00 20.0000D -0.3706 -0.4030 0.2474 0.020 1.22 -0.577E+01 550. 0.267E-02 0.265E-02 0.824E-04 0.140E-14 8 13 9 8 3 6 2 6 12 0.000 0.000 0 255 0 23 00:07:40.00 20.0000D -0.4267D -0.4554D 0.2850 0.021 1.22 -0.539E+01 550. 0.371E-02 0.369E-02 0.653E-03 0.886E-13 4 4 4 8 1 8 2 6 12 0.000 0.000 0 256 0 24 00:08:00.00 20.0000D -0.4241 -0.4568 0.3233 0.023 1.22 -0.506E+01 550. 0.368E-02 0.366E-02 0.905E-04 0.190E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 257 0 25 00:08:20.00 20.0000D -0.4208 -0.4522 0.3672 0.024 1.22 -0.471E+01 550. 0.367E-02 0.365E-02 0.931E-04 0.196E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 258 0 26 00:08:40.00 20.0000D -0.4174 -0.4421 0.4169 0.026 1.22 -0.440E+01 550. 0.368E-02 0.366E-02 0.957E-04 0.204E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 259 0 27 00:09:00.00 20.0000D -0.4142 -0.4279 0.4727 0.028 1.22 -0.410E+01 550. 0.372E-02 0.370E-02 0.981E-04 0.214E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 260 0 28 00:09:20.00 20.0000D -0.4116 -0.4106 0.5348 0.030 1.18 -0.381E+01 500. 0.380E-02 0.378E-02 0.100E-03 0.227E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 261 0 29 00:09:40.00 20.0000D -0.4107 -0.4087 0.6032 0.032 1.18 -0.355E+01 500. 0.393E-02 0.391E-02 0.103E-03 0.242E-14 5 4 4 5 3 36 2 6 12 0.000 0.000 0 262 0 30 00:10:00.00 20.0000D 0.4798D 0.4695D 0.7094 0.035 1.18 -0.330E+01 500. 0.519E-02 0.516E-02 0.656E-03 0.867E-13 5 14 20 9 7 14 2 6 12 0.000 0.000 0 263 0 31 00:10:20.00 20.0000D 0.4589 -0.4757 0.8018 0.037 1.18 -0.308E+01 500. 0.543E-02 0.540E-02 0.108E-03 0.318E-14 5 14 20 8 2 6 3 6 12 0.000 0.000 0 264 0 32 00:10:40.00 20.0000D 0.4450 -0.4782 0.9236 0.040 1.18 -0.288E+01 500. 0.579E-02 0.576E-02 0.110E-03 0.344E-14 6 39 30 8 2 6 3 6 12 0.000 0.000 0 265 0 33 00:11:00.00 20.0000A 0.4349 -0.4767 1.0530 0.043 1.18 -0.266E+01 500. 0.629E-02 0.626E-02 0.111E-03 0.377E-14 6 39 30 8 2 6 3 6 12 0.000 0.000 0 266 0 34 00:11:20.00 20.0000A -0.4683 -0.5183 1.1853 0.046 1.18 -0.249E+01 500. 0.699E-02 0.696E-02 0.111E-03 0.416E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 267 0 35 00:11:40.00 20.0000A -0.5109 -0.5803 1.3145 0.049 1.18 -0.233E+01 500. 0.793E-02 0.791E-02 0.112E-03 0.464E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 268 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 269 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 270 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.673E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 271 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 272 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 273 0 41 00:13:40.00 20.0000A 0.8943 -0.9392 2.2218 0.069 1.18 -0.165E+01 500. 0.237E-01 0.236E-01 0.116E-03 0.100E-13 1 35 11 1 29 17 5 6 12 0.000 0.000 0 274 0 42 00:14:00.00 18.9000A 0.9547 -0.9776 2.3771 0.072 1.28 -0.157E+01 650. 0.290E-01 0.290E-01 0.120E-03 0.115E-13 1 35 11 1 29 17 6 6 12 0.000 0.000 0 275 0 43 00:14:18.89 18.4000A 1.0028 -1.0297 2.4487 0.075 1.28 -0.151E+01 650. 0.351E-01 0.351E-01 0.119E-03 0.123E-13 1 35 11 1 28 17 6 6 12 0.000 0.000 0 276 0 44 00:14:37.29 17.7000A 1.0384 -1.1110 2.5490 0.078 1.28 -0.147E+01 650. 0.421E-01 0.421E-01 0.122E-03 0.135E-13 1 35 11 1 28 17 7 6 12 0.000 0.000 0 277 0 45 00:14:55.00 16.2000A 1.0598 -1.1790 2.7722 0.081 1.28 -0.141E+01 650. 0.499E-01 0.498E-01 0.125E-03 0.146E-13 1 35 11 1 28 17 5 24 2 0.000 0.000 0 278 0 46 00:15:11.20 15.4000A 1.0660 -1.2306 2.9231 0.083 1.18 -0.138E+01 500. 0.578E-01 0.578E-01 0.122E-03 0.148E-13 1 35 11 1 28 17 6 24 2 0.000 0.000 0 279 0 47 00:15:26.60 14.3000A -1.0741 -1.2693 3.1400 0.084 1.28 -0.135E+01 650. 0.661E-01 0.661E-01 0.123E-03 0.155E-13 1 24 19 1 28 17 6 24 2 0.000 0.000 0 280 0 48 00:15:40.89 13.8000A 1.1031 -1.2960 3.2560 0.086 1.28 -0.132E+01 650. 0.745E-01 0.745E-01 0.122E-03 0.156E-13 1 17 38 1 28 17 6 24 2 0.000 0.000 0 281 0 49 00:15:54.69 13.2000A 1.1508 -1.3140 3.4217 0.087 1.28 -0.131E+01 650. 0.832E-01 0.832E-01 0.125E-03 0.163E-13 1 17 38 1 28 17 7 24 2 0.000 0.000 0 282 0 50 00:16:07.89 12.9000A 1.1879 -1.3247 3.4993 0.088 1.29 -0.130E+01 650. 0.919E-01 0.919E-01 0.126E-03 0.168E-13 1 17 38 1 28 17 7 24 2 0.000 0.000 0 283 0 51 00:16:20.79 12.4000A 1.2158 -1.3300 3.6230 0.089 1.29 -0.129E+01 650. 0.101E+00 0.101E+00 0.130E-03 0.175E-13 1 17 38 1 28 17 8 24 2 0.000 0.000 0 284 0 52 00:16:33.19 12.3000A 1.2346 -1.3449 3.6668 0.090 1.29 -0.128E+01 650. 0.110E+00 0.110E+00 0.131E-03 0.179E-13 1 17 38 1 28 18 8 34 13 0.000 0.000 0 285 0 53 00:16:45.49 11.9000A 1.2460 -1.3664 3.7680 0.091 1.29 -0.126E+01 650. 0.119E+00 0.119E+00 0.136E-03 0.188E-13 1 17 38 1 28 18 9 34 13 0.000 0.000 0 286 0 54 00:16:57.39 11.4000A 1.2899 -1.3815 3.9422 0.092 1.29 -0.126E+01 650. 0.128E+00 0.128E+00 0.138E-03 0.191E-13 5 31 10 1 28 18 9 34 13 0.000 0.000 0 287 0 55 00:17:08.79 11.2000A 1.3175 -1.3906 4.0211 0.092 1.29 -0.124E+01 650. 0.136E+00 0.136E+00 0.138E-03 0.192E-13 5 31 10 1 28 18 9 34 13 0.000 0.000 0 288 0 56 00:17:20.00 10.8000A -1.3135 -1.3943 4.1804 0.093 1.38 -0.123E+01 800. 0.145E+00 0.145E+00 0.141E-03 0.198E-13 5 29 17 1 28 18 10 34 13 0.000 0.000 0 289 0 57 00:17:30.79 10.6000A -1.3858 -1.3935 4.2483 0.093 1.38 -0.123E+01 800. 0.153E+00 0.153E+00 0.141E-03 0.199E-13 5 29 17 1 28 18 10 34 13 0.000 0.000 0 290 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 291 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 292 0 60 00:18:02.00 10.1000A -1.4413 1.4902 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 293 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 294 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 295 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 296 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 297 0 65 00:18:52.19 10.2000A -1.5591 -1.7454 4.4046 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 298 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 299 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 300 0 68 00:19:23.29 10.6000A -1.5827 1.6967 4.2637 0.100 1.41 -0.114E+01 850. 0.235E+00 0.234E+00 0.188E-03 0.273E-13 4 32 18 12 17 28 15 35 12 0.000 0.000 0 301 0 69 00:19:33.89 10.8000A -1.5957 -1.7194 4.1687 0.101 1.41 -0.113E+01 850. 0.241E+00 0.241E+00 0.194E-03 0.282E-13 7 12 6 9 7 16 16 35 13 0.000 0.000 0 302 0 70 00:19:44.69 10.7000A -1.5970 -1.7896 4.1867 0.102 1.41 -0.113E+01 850. 0.247E+00 0.247E+00 0.202E-03 0.293E-13 7 12 6 10 7 17 16 35 13 0.000 0.000 0 303 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 304 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 305 0 73 00:20:17.29 10.9000A 1.8396 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 306 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 307 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 308 0 76 00:20:49.60 10.7000A 1.9390 -1.9956 4.2224 0.107 1.41 -0.105E+01 850. 0.273E+00 0.273E+00 0.220E-03 0.314E-13 7 24 31 9 4 0 13 23 17 0.000 0.000 0 309 0 77 00:21:00.30 10.5000A 2.0200 -2.0694 4.2698 0.108 1.41 -0.104E+01 850. 0.276E+00 0.275E+00 0.223E-03 0.317E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 310 0 78 00:21:10.80 10.5000A 2.0530 -2.1072 4.3032 0.109 1.41 -0.103E+01 850. 0.278E+00 0.277E+00 0.221E-03 0.314E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 311 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 312 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 313 0 81 00:21:42.10 10.5000A 1.9801 -1.9716 4.2997 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 314 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 315 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 316 0 84 00:22:13.50 10.7000A 1.9413 -2.0334 4.2011 0.115 1.41 -0.961E+00 850. 0.283E+00 0.283E+00 0.229E-03 0.324E-13 14 4 0 10 4 0 15 23 17 0.000 0.000 0 317 0 85 00:22:24.20 10.8000A 1.9316 -2.0194 4.1554 0.116 1.41 -0.952E+00 850. 0.284E+00 0.283E+00 0.234E-03 0.332E-13 14 4 0 10 4 0 16 23 17 0.000 0.000 0 318 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 319 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 320 0 88 00:22:57.40 12.0000A 1.8992 -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 321 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 322 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 323 0 91 00:23:34.90 13.2000A 1.8498 -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 324 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 325 0 93 00:24:01.50 13.7000A 1.7141 -2.1684 3.2965 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 326 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.298E-03 0.426E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0 327 0 95 00:24:29.10 14.1000A 1.7197 -2.1962 3.1874 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 0 328 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 329 0 97 00:24:57.30 13.8000A 1.6954 -2.1302 3.2532 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 330 0 98 00:25:11.10 13.8000A -1.7060 -2.1049 3.2691 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 331 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.295E-03 0.416E-13 10 26 31 8 23 5 11 23 1 0.000 0.000 0 332 0 100 00:25:38.70 13.8000A 1.7646 -2.1973 3.2551 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 0 333 0 101 00:25:52.50 14.1000A 1.7850 -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 0 334 0 102 00:26:06.60 14.1000A 1.7815 -2.0922 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 335 0 103 00:26:20.70 14.4000A -1.8030 -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 336 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.303E-03 0.413E-13 1 36 26 8 22 5 12 23 19 0.000 0.000 0 337 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.276E+00 0.301E-03 0.410E-13 1 37 26 8 22 5 9 38 26 0.000 0.000 0 338 0 106 00:27:03.90 14.1000A 1.8658 1.8646 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 339 0 107 00:27:18.00 13.8000A 1.8273 1.8719 3.2544 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 0 340 0 108 00:27:31.80 13.8000A 1.8485 1.8998 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 0 341 0 109 00:27:45.60 13.9000A 1.9161 1.9299 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 0 342 0 110 00:27:59.50 13.9000A 1.9930 1.9488 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 0 343 0 111 00:28:13.40 13.9000A 2.0631 2.0386 3.2310 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 0 344 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.285E-03 0.386E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0 345 0 113 00:28:41.30 14.0000A 2.1239 2.1077 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 0 346 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.388E-13 1 16 0 9 36 23 13 37 37 0.000 0.000 0 347 0 115 00:29:09.30 13.9000A 2.1914 2.1430 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 0 348 0 116 00:29:23.20 13.9000A 2.1629 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 0 349 0 117 00:29:37.10 14.2000A 2.1088 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 0 350 0 118 00:29:51.30 14.3000A 2.1133 1.9030 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 0 351 0 119 00:30:05.60 14.4000A 2.0602 -1.9071 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 0 352 0 120 00:30:20.00 14.1000A 2.0168 -2.0168 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 0 353 0 121 00:30:34.10 13.8000A -2.0339 -2.0424 3.2554 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 0 354 0 122 00:30:47.90 13.4000A 1.9253 -2.1019 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 0 355 0 123 00:31:01.30 13.2000A -1.9151 -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 0 356 0 124 00:31:14.50 13.1000A -1.9606 -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 0 357 0 125 00:31:27.60 13.1000A -1.8503 -1.9343 3.4386 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 0 358 0 126 00:31:40.70 13.1000A -1.9153 -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 0 359 0 127 00:31:53.80 13.4000A -1.9319 -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 0 360 0 128 00:32:07.20 13.7000A -1.9610 -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 0 361 0 129 00:32:20.90 13.9000A -2.0531 -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 0 362 0 130 00:32:34.80 14.1000A -2.0925 -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 0 363 0 131 00:32:48.90 14.2000A -2.1181 -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 0 364 0 132 00:33:03.10 14.6000A -2.1477 -1.9040 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 0 365 0 133 00:33:17.70 14.7000A -2.2353 -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 0 366 0 134 00:33:32.40 14.7000A -2.3231 -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 0 367 0 135 00:33:47.10 14.6000A -2.4176 -1.8784 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 0 368 0 136 00:34:01.70 14.9000A -2.5034 -1.8777 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 0 369 0 137 00:34:16.60 15.2000A -2.5428 -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 0 370 0 138 00:34:31.80 15.3000A -2.4933 -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 0 371 0 139 00:34:47.10 15.4000A -2.3596 -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 0 372 0 140 00:35:02.50 15.5000A -2.3919 -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 0 373 0 141 00:35:18.00 15.8000A -2.4704 -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 0 374 0 142 00:35:33.80 15.8000A -2.4681 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 0 375 0 143 00:35:49.60 15.7000A -2.3624 -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 0 376 0 144 00:36:05.30 15.2000A -2.2471 -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 0 377 0 145 00:36:20.50 14.7000A -2.0989 -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 0 378 0 146 00:36:35.20 14.2000A -2.0337 -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 0 379 0 147 00:36:49.40 13.9000A -2.0495 -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 0 380 0 148 00:37:03.30 13.6000A -2.1261 -1.8874 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 0 381 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.258E-03 0.299E-13 11 12 13 1 5 34 8 38 32 0.000 0.000 0 382 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.255E-03 0.294E-13 11 12 13 1 6 35 8 38 32 0.000 0.000 0 383 0 151 00:37:43.70 13.5000A -2.1150 -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 0 384 0 152 00:37:57.20 13.6000A -2.0518 -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 0 385 0 153 00:38:10.80 13.9000A -1.9707 -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 0 386 0 154 00:38:24.70 14.2000A 1.9577 -2.0347 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 0 387 0 155 00:38:38.90 13.4000A 1.9885 -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 0 388 0 156 00:38:52.30 12.8000A 2.0049 -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 0 389 0 157 00:39:05.10 12.6000A 1.9495 -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 0 390 0 158 00:39:17.70 12.5000A -1.8994 -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 0 391 0 159 00:39:30.20 12.8000A -1.8343 -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 0 392 0 160 00:39:43.00 13.0000A -1.8423 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 0 393 0 161 00:39:56.00 13.1000A -1.8734 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 0 394 0 162 00:40:09.10 13.1000A -1.8857 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 0 395 0 163 00:40:22.20 12.9000A -1.8946 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 0 396 0 164 00:40:35.10 13.0000A -1.8988 2.1083 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 0 397 0 165 00:40:48.10 13.3000A -1.9469 2.1590 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 0 398 0 166 00:41:01.40 13.5000A -1.9768 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 0 399 0 167 00:41:14.90 13.7000A -1.9837 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 0 400 0 168 00:41:28.60 13.6000A -1.9913 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 0 401 0 169 00:41:42.20 13.4000A -1.9813 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 0 402 0 170 00:41:55.60 13.5000A -1.9944 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 0 403 0 171 00:42:09.10 13.5000A -2.0090 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 0 404 0 172 00:42:22.60 13.6000A -1.9963 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 0 405 0 173 00:42:36.20 13.9000A -1.9847 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 0 406 0 174 00:42:50.10 14.1000A -2.0020 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 0 407 0 175 00:43:04.20 14.3000A -1.9699 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 0 408 0 176 00:43:18.50 14.8000A -1.9215 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 0 409 0 177 00:43:33.30 15.2000A -1.9580 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 0 410 0 178 00:43:48.50 15.5000A -1.9076 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 0 411 0 179 00:44:04.00 15.0000A -1.8769 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 0 412 0 180 00:44:19.00 14.7000A -1.8140 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 0 413 0 181 00:44:33.70 14.6000A -1.7678 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 0 414 0 182 00:44:48.30 14.2000A 1.7560 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 0 415 0 183 00:45:02.50 13.9000A 1.7854 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 0 416 0 184 00:45:16.40 13.6000A 1.7977 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 0 417 0 185 00:45:30.00 13.3000A 1.8109 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 0 418 0 186 00:45:43.30 13.2000A 1.8581 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 0 419 0 187 00:45:56.50 13.1000A 1.8920 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 0 420 0 188 00:46:09.60 12.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 0 421 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 0 422 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 0 423 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 0 424 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 0 425 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 0 426 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 0 427 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 14 14 10 17 0.000 0.000 0 428 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 1 13 38 15 11 9 18 0.000 0.000 0 429 0 197 00:48:07.50 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 0 430 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 0 431 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 0 432 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 0 433 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 0 434 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 0 435 0 203 00:49:23.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 0 436 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 0 437 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 0 438 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 0 439 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 0 440 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 0 441 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 0 442 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 0 443 0 211 00:51:02.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 0 444 0 212 00:51:14.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 0 445 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 0 446 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 0 447 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 0 448 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 0 449 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 0 450 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 0 451 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 0 452 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 0 453 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 0 454 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 0 455 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 24 3 17 17 10 8 35 5 0.000 0.000 0 456 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 0 457 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 0 458 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 0 459 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 0 460 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 0 461 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 0 462 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 0 463 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 0 464 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 0 465 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 0 466 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 0 467 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 0 468 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 0 469 0 237 00:56:28.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 36 12 36 7 0.000 0.000 0 470 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 0 471 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 0 472 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 0 473 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 0 474 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 0 475 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 0 476 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 0 477 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 0 478 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 0 479 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 0 480 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 0 481 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 6 0.000 0.000 0 482 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 0 483 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 37 10 37 5 0.000 0.000 0 484 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 0 485 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 0 486 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 0 232 0 0 00:00:00.00 0.0200A -0.2131D -0.2288D 0.1182 0.000 1.38 0.000E+00 800. 0.104E-02 0.103E-02 0.597E-03 0.867E-13 7 7 11 11 7 19 5 26 7 0.000 0.000 0 233 0 1 00:00:00.02 0.0400A -0.2131 -0.2288 0.1182 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.326E-08 0.443E-18 7 7 11 11 7 19 5 26 7 0.000 0.000 0 234 0 2 00:00:00.06 0.0800A -0.2131 -0.2288 0.1182 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.609E-08 0.827E-18 7 7 11 11 7 19 5 26 7 0.000 0.000 0 235 0 3 00:00:00.14 0.1600A -0.2131 -0.2288 0.1182 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.121E-07 0.164E-17 7 7 11 11 7 19 5 26 7 0.000 0.000 0 236 0 4 00:00:00.30 0.3200A -0.2131 -0.2288 0.1181 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.243E-07 0.328E-17 7 7 11 11 7 19 5 26 7 0.000 0.000 0 237 0 5 00:00:00.62 0.6400A -0.2131 -0.2288 0.1181 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.490E-07 0.655E-17 7 7 11 11 7 19 5 26 7 0.000 0.000 0 238 0 6 00:00:01.26 1.2800A -0.2131 -0.2287 0.1181 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.101E-06 0.131E-16 7 7 11 11 7 19 5 26 7 0.000 0.000 0 239 0 7 00:00:02.54 2.5600A -0.2130 -0.2286 0.1180 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.214E-06 0.262E-16 7 7 11 11 7 19 5 26 7 0.000 0.000 0 240 0 8 00:00:05.09 5.1200D -0.2130 -0.2285 0.1179 0.009 0.69 -0.131E+02 100. 0.104E-02 0.103E-02 0.481E-06 0.525E-16 7 7 11 11 7 19 5 26 7 0.000 0.000 0 241 0 9 00:00:10.21 10.2000D -0.2129 -0.2281 0.1177 0.009 0.69 -0.131E+02 100. 0.103E-02 0.103E-02 0.117E-05 0.105E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 242 0 10 00:00:20.42 20.0000D -0.2125 -0.2272 0.1172 0.009 0.69 -0.132E+02 100. 0.103E-02 0.102E-02 0.317E-05 0.210E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 243 0 11 00:00:40.42 20.0000D -0.2117 -0.2252 0.1162 0.009 0.69 -0.133E+02 100. 0.102E-02 0.102E-02 0.942E-05 0.416E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 244 0 12 00:01:00.42 20.0000D -0.2106 -0.2227 0.1152 0.009 0.69 -0.134E+02 100. 0.101E-02 0.101E-02 0.131E-04 0.421E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 245 0 13 00:01:20.42 20.0000D -0.2093 -0.2197 0.1141 0.009 0.69 -0.134E+02 100. 0.100E-02 0.997E-03 0.167E-04 0.425E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 246 0 14 00:01:40.42 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.428E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 247 0 15 00:02:00.42 20.0000D -0.2059 -0.2126 0.1117 0.009 0.93 -0.133E+02 250. 0.983E-03 0.979E-03 0.240E-04 0.430E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 248 0 16 00:02:20.42 20.0000D -0.2038 -0.2086 0.1103 0.009 1.22 -0.131E+02 550. 0.975E-03 0.970E-03 0.276E-04 0.431E-15 7 7 11 11 7 19 5 26 7 0.000 0.000 0 249 0 17 00:02:40.42 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 250 0 18 00:03:00.42 20.0000D 0.3330 0.3039 0.1620 0.011 1.22 -0.105E+02 550. 0.199E-02 0.198E-02 0.374E-04 0.830E-15 10 3 39 10 8 27 5 36 28 0.000 0.000 0 251 0 19 00:03:20.42 20.0000D 0.3296 -0.2989 0.1601 0.011 1.22 -0.103E+02 550. 0.196E-02 0.195E-02 0.410E-04 0.828E-15 10 3 39 8 4 6 5 36 28 0.000 0.000 0 252 0 20 00:03:40.42 20.0000D 0.3249 -0.2987 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 253 0 21 00:04:00.42 20.0000D 0.3189 -0.2990 0.1559 0.012 1.22 -0.974E+01 550. 0.191E-02 0.190E-02 0.479E-04 0.831E-15 10 3 39 8 0 7 5 36 28 0.000 0.000 0 254 0 22 00:04:20.42 20.0000D 0.3119 -0.3006 0.1536 0.012 1.22 -0.940E+01 550. 0.189E-02 0.188E-02 0.513E-04 0.836E-15 10 3 39 8 0 7 5 36 28 0.000 0.000 0 255 0 23 00:04:40.42 20.0000D 0.3058 -0.3013 0.1512 0.013 1.22 -0.905E+01 550. 0.187E-02 0.186E-02 0.547E-04 0.842E-15 10 3 0 8 0 7 5 36 28 0.000 0.000 0 256 0 24 00:05:00.42 20.0000D 0.4152D -0.3965D -0.1813 0.013 1.22 -0.850E+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 257 0 25 00:05:20.42 20.0000D 0.4184 -0.4078 -0.1811 0.014 1.22 -0.811E+01 550. 0.287E-02 0.285E-02 0.641E-04 0.129E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 258 0 26 00:05:40.42 20.0000D 0.4164 -0.4153 -0.1797 0.015 1.22 -0.772E+01 550. 0.283E-02 0.280E-02 0.673E-04 0.129E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 259 0 27 00:06:00.42 20.0000D 0.4098 -0.4191 -0.1773 0.015 1.22 -0.735E+01 550. 0.279E-02 0.276E-02 0.704E-04 0.130E-14 8 17 5 8 3 6 8 19 3 0.000 0.000 0 260 0 28 00:06:20.42 20.0000D 0.3993 -0.4194 0.1804 0.016 1.22 -0.697E+01 550. 0.275E-02 0.273E-02 0.735E-04 0.131E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 261 0 29 00:06:40.42 20.0000D 0.3856 -0.4167 0.1985 0.017 1.22 -0.655E+01 550. 0.272E-02 0.270E-02 0.766E-04 0.133E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 262 0 30 00:07:00.42 20.0000D 0.3696 -0.4111 0.2201 0.018 1.22 -0.615E+01 550. 0.269E-02 0.267E-02 0.796E-04 0.136E-14 8 17 5 8 3 6 2 36 2 0.000 0.000 0 263 0 31 00:07:20.42 20.0000D -0.3708 -0.4031 0.2480 0.020 1.22 -0.576E+01 550. 0.267E-02 0.265E-02 0.825E-04 0.140E-14 8 13 9 8 3 6 2 6 12 0.000 0.000 0 264 0 32 00:07:40.42 20.0000D -0.4268D -0.4556D 0.2856 0.021 1.22 -0.538E+01 550. 0.372E-02 0.369E-02 0.653E-03 0.886E-13 4 4 4 8 1 8 2 6 12 0.000 0.000 0 265 0 33 00:08:00.42 20.0000D -0.4241 -0.4571 0.3241 0.023 1.22 -0.505E+01 550. 0.369E-02 0.366E-02 0.906E-04 0.190E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 266 0 34 00:08:20.42 20.0000D -0.4209 -0.4524 0.3681 0.024 1.22 -0.471E+01 550. 0.368E-02 0.365E-02 0.932E-04 0.196E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 267 0 35 00:08:40.42 20.0000D -0.4176 -0.4424 0.4179 0.026 1.22 -0.439E+01 550. 0.369E-02 0.366E-02 0.957E-04 0.205E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 268 0 36 00:09:00.42 20.0000D -0.4144 -0.4281 0.4739 0.028 1.22 -0.409E+01 550. 0.373E-02 0.371E-02 0.981E-04 0.215E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 269 0 37 00:09:20.42 20.0000D -0.4117 -0.4107 0.5361 0.030 1.18 -0.381E+01 500. 0.381E-02 0.378E-02 0.100E-03 0.227E-14 4 4 4 8 1 8 2 6 12 0.000 0.000 0 270 0 38 00:09:40.42 20.0000D -0.4112 -0.4088 0.6046 0.032 1.18 -0.354E+01 500. 0.393E-02 0.391E-02 0.103E-03 0.243E-14 5 4 4 5 3 36 2 6 12 0.000 0.000 0 271 0 39 00:10:00.42 20.0000D 0.4799D 0.4697D 0.7110 0.035 1.18 -0.329E+01 500. 0.520E-02 0.517E-02 0.656E-03 0.867E-13 5 14 20 9 7 14 2 6 12 0.000 0.000 0 272 0 40 00:10:20.42 20.0000D 0.4589 -0.4761 0.8042 0.037 1.18 -0.308E+01 500. 0.544E-02 0.541E-02 0.108E-03 0.318E-14 5 14 20 8 2 6 3 6 12 0.000 0.000 0 273 0 41 00:10:40.42 20.0000D 0.4452 -0.4785 0.9262 0.040 1.18 -0.287E+01 500. 0.580E-02 0.577E-02 0.110E-03 0.345E-14 6 39 30 8 2 6 3 6 12 0.000 0.000 0 274 0 42 00:11:00.42 20.0000A 0.4350 -0.4770 1.0557 0.043 1.18 -0.266E+01 500. 0.631E-02 0.628E-02 0.111E-03 0.377E-14 6 39 30 8 2 6 3 6 12 0.000 0.000 0 275 0 43 00:11:20.42 20.0000A -0.4690 -0.5194 1.1880 0.046 1.18 -0.248E+01 500. 0.701E-02 0.698E-02 0.111E-03 0.417E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 276 0 44 00:11:40.42 20.0000A -0.5115 -0.5815 1.3172 0.049 1.18 -0.232E+01 500. 0.796E-02 0.793E-02 0.112E-03 0.465E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 277 0 45 00:12:00.42 20.0000A -0.5527 -0.6472 1.4375 0.052 1.18 -0.218E+01 500. 0.922E-02 0.920E-02 0.112E-03 0.523E-14 1 3 23 1 29 17 3 6 12 0.000 0.000 0 278 0 46 00:12:20.42 20.0000A 0.6185 -0.7151 1.6172 0.055 1.18 -0.204E+01 500. 0.109E-01 0.109E-01 0.112E-03 0.593E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 279 0 47 00:12:40.42 20.0000A 0.6878 -0.7807 1.7894 0.059 1.18 -0.193E+01 500. 0.131E-01 0.130E-01 0.112E-03 0.674E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 280 0 48 00:13:00.42 20.0000A 0.7587 -0.8410 1.9157 0.062 1.18 -0.182E+01 500. 0.158E-01 0.158E-01 0.112E-03 0.769E-14 1 35 11 1 29 17 4 6 12 0.000 0.000 0 281 0 49 00:13:20.42 20.0000A 0.8287 -0.8944 2.1048 0.066 1.18 -0.173E+01 500. 0.194E-01 0.193E-01 0.114E-03 0.879E-14 1 35 11 1 29 17 5 6 12 0.000 0.000 0 282 0 50 00:13:40.42 20.0000A 0.8954 -0.9401 2.2237 0.069 1.18 -0.165E+01 500. 0.237E-01 0.237E-01 0.116E-03 0.101E-13 1 35 11 1 29 17 5 6 12 0.000 0.000 0 283 0 51 00:14:00.42 18.9000A 0.9557 -0.9784 2.3807 0.073 1.28 -0.157E+01 650. 0.291E-01 0.291E-01 0.120E-03 0.115E-13 1 35 11 1 29 17 6 6 12 0.000 0.000 0 284 0 52 00:14:19.32 18.4000A 1.0036 -1.0317 2.4493 0.076 1.28 -0.151E+01 650. 0.353E-01 0.352E-01 0.119E-03 0.123E-13 1 35 11 1 28 17 6 6 12 0.000 0.000 0 285 0 53 00:14:37.72 17.6000A 1.0389 -1.1128 2.5521 0.078 1.28 -0.147E+01 650. 0.423E-01 0.422E-01 0.122E-03 0.135E-13 1 35 11 1 28 17 7 6 12 0.000 0.000 0 286 0 54 00:14:55.32 16.2000A 1.0601 -1.1802 2.7746 0.081 1.28 -0.141E+01 650. 0.500E-01 0.499E-01 0.124E-03 0.145E-13 1 35 11 1 28 17 5 24 2 0.000 0.000 0 287 0 55 00:15:11.52 15.4000A 1.0660 -1.2317 2.9274 0.083 1.28 -0.138E+01 650. 0.579E-01 0.579E-01 0.122E-03 0.148E-13 1 35 11 1 28 17 6 24 2 0.000 0.000 0 288 0 56 00:15:26.92 14.3000A -1.0747 -1.2702 3.1422 0.084 1.28 -0.135E+01 650. 0.663E-01 0.663E-01 0.124E-03 0.155E-13 1 24 19 1 28 17 6 24 2 0.000 0.000 0 289 0 57 00:15:41.22 13.8000A 1.1042 -1.2968 3.2571 0.086 1.28 -0.132E+01 650. 0.747E-01 0.747E-01 0.122E-03 0.157E-13 1 17 38 1 28 17 6 24 2 0.000 0.000 0 290 0 58 00:15:55.01 13.1000A 1.1518 -1.3146 3.4236 0.087 1.28 -0.131E+01 650. 0.834E-01 0.833E-01 0.125E-03 0.164E-13 1 17 38 1 28 17 7 24 2 0.000 0.000 0 291 0 59 00:16:08.12 12.9000A 1.1885 -1.3252 3.4994 0.088 1.29 -0.130E+01 650. 0.920E-01 0.920E-01 0.126E-03 0.167E-13 1 17 38 1 28 17 7 24 2 0.000 0.000 0 292 0 60 00:16:21.01 12.4000A 1.2164 -1.3304 3.6240 0.089 1.29 -0.129E+01 650. 0.101E+00 0.101E+00 0.130E-03 0.175E-13 1 17 38 1 28 17 8 24 2 0.000 0.000 0 293 0 61 00:16:33.41 12.3000A 1.2351 -1.3457 3.6690 0.090 1.29 -0.127E+01 650. 0.110E+00 0.110E+00 0.131E-03 0.179E-13 1 17 38 1 28 18 8 34 13 0.000 0.000 0 294 0 62 00:16:45.71 11.9000A 1.2464 -1.3673 3.7728 0.091 1.29 -0.126E+01 650. 0.119E+00 0.119E+00 0.136E-03 0.188E-13 1 17 38 1 28 18 9 34 13 0.000 0.000 0 295 0 63 00:16:57.61 11.4000A 1.2904 -1.3825 3.9447 0.092 1.29 -0.126E+01 650. 0.128E+00 0.128E+00 0.138E-03 0.192E-13 5 31 10 1 28 18 9 34 13 0.000 0.000 0 296 0 64 00:17:09.01 11.2000A 1.3177 -1.3915 4.0232 0.092 1.29 -0.124E+01 650. 0.137E+00 0.137E+00 0.138E-03 0.193E-13 5 31 10 1 28 18 9 34 13 0.000 0.000 0 297 0 65 00:17:20.22 10.8000A -1.3137 -1.3953 4.1834 0.093 1.38 -0.123E+01 800. 0.145E+00 0.145E+00 0.141E-03 0.198E-13 5 29 17 1 28 18 10 34 13 0.000 0.000 0 298 0 66 00:17:31.01 10.6000A -1.3848 -1.3945 4.2505 0.093 1.38 -0.123E+01 800. 0.154E+00 0.153E+00 0.141E-03 0.200E-13 5 29 17 1 28 18 10 34 13 0.000 0.000 0 299 0 67 00:17:41.61 10.4000A -1.4187 -1.3899 4.3464 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 300 0 68 00:17:52.01 10.2000A -1.4319 1.3894 4.4010 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 301 0 69 00:18:02.22 10.1000A -1.4430 1.4896 4.4566 0.095 1.38 -0.120E+01 800. 0.178E+00 0.177E+00 0.148E-03 0.212E-13 4 32 19 7 37 19 12 34 13 0.000 0.000 0 302 0 70 00:18:12.31 10.0000A -1.4707 -1.5187 4.4834 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 303 0 71 00:18:22.31 10.0000A -1.4864 -1.6347 4.4882 0.096 1.38 -0.119E+01 800. 0.193E+00 0.193E+00 0.154E-03 0.222E-13 4 32 19 8 7 16 13 34 13 0.000 0.000 0 304 0 72 00:18:32.31 10.0000A -1.4923 -1.7223 4.4794 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 305 0 73 00:18:42.31 10.1000A -1.5096 -1.7607 4.4387 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 306 0 74 00:18:52.41 10.2000A -1.5638 -1.7362 4.4065 0.098 1.41 -0.116E+01 850. 0.214E+00 0.214E+00 0.169E-03 0.245E-13 4 32 18 8 7 16 13 34 13 0.000 0.000 0 307 0 75 00:19:02.61 10.4000A -1.5955 -1.6476 4.3472 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 308 0 76 00:19:13.01 10.5000A -1.6043 -1.6443 4.2951 0.099 1.41 -0.115E+01 850. 0.228E+00 0.228E+00 0.182E-03 0.265E-13 4 32 18 9 7 16 15 35 12 0.000 0.000 0 309 0 77 00:19:23.51 10.6000A -1.5910 -1.6991 4.2632 0.100 1.41 -0.114E+01 850. 0.235E+00 0.235E+00 0.188E-03 0.273E-13 4 32 18 9 7 16 15 35 12 0.000 0.000 0 310 0 78 00:19:34.11 10.8000A -1.5967 -1.7244 4.1654 0.101 1.41 -0.113E+01 850. 0.241E+00 0.241E+00 0.194E-03 0.282E-13 7 12 6 9 7 16 16 35 13 0.000 0.000 0 311 0 79 00:19:44.91 10.8000A -1.5986 -1.7950 4.1824 0.102 1.41 -0.113E+01 850. 0.247E+00 0.247E+00 0.202E-03 0.293E-13 7 12 6 10 7 17 16 35 13 0.000 0.000 0 312 0 80 00:19:55.71 11.0000A 1.6051 -1.8316 4.1076 0.102 1.41 -0.112E+01 850. 0.253E+00 0.253E+00 0.206E-03 0.298E-13 7 1 37 10 7 17 12 23 18 0.000 0.000 0 313 0 81 00:20:06.71 10.9000A 1.7144 -1.8050 4.1426 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 314 0 82 00:20:17.61 10.8000A 1.8454 1.7495 4.1494 0.104 1.41 -0.109E+01 850. 0.263E+00 0.262E+00 0.215E-03 0.310E-13 7 24 30 14 18 28 13 23 18 0.000 0.000 0 315 0 83 00:20:28.41 10.7000A 1.9170 1.8285 4.2009 0.105 1.41 -0.107E+01 850. 0.267E+00 0.266E+00 0.217E-03 0.310E-13 7 24 30 14 18 28 13 23 18 0.000 0.000 0 316 0 84 00:20:39.11 10.7000A 1.9339 -1.8947 4.1949 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 317 0 85 00:20:49.81 10.7000A 1.9419 -1.9959 4.2222 0.107 1.41 -0.105E+01 850. 0.273E+00 0.273E+00 0.220E-03 0.314E-13 7 24 31 9 4 0 13 23 17 0.000 0.000 0 318 0 86 00:21:00.51 10.5000A 2.0215 -2.0698 4.2704 0.108 1.41 -0.104E+01 850. 0.276E+00 0.275E+00 0.223E-03 0.317E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 319 0 87 00:21:11.01 10.5000A 2.0535 -2.1077 4.3029 0.109 1.41 -0.103E+01 850. 0.278E+00 0.277E+00 0.221E-03 0.314E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 320 0 88 00:21:21.51 10.4000A 2.0493 -2.1047 4.3239 0.110 1.41 -0.102E+01 850. 0.280E+00 0.279E+00 0.223E-03 0.316E-13 7 24 31 9 4 0 14 23 17 0.000 0.000 0 321 0 89 00:21:31.92 10.4000A 2.0209 -2.0588 4.3203 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 322 0 90 00:21:42.32 10.5000A 1.9795 -1.9725 4.2984 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 323 0 91 00:21:52.82 10.4000A 1.9323 -1.9477 4.3275 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 324 0 92 00:22:03.22 10.5000A 1.8929 -2.0067 4.3025 0.114 1.41 -0.971E+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 325 0 93 00:22:13.72 10.7000A 1.9405 -2.0350 4.2025 0.115 1.41 -0.961E+00 850. 0.283E+00 0.283E+00 0.229E-03 0.324E-13 14 4 0 10 4 0 15 23 17 0.000 0.000 0 326 0 94 00:22:24.42 10.8000A 1.9307 -2.0212 4.1558 0.116 1.41 -0.952E+00 850. 0.284E+00 0.283E+00 0.234E-03 0.332E-13 14 4 0 10 4 0 16 23 17 0.000 0.000 0 327 0 95 00:22:35.22 11.0000A 1.8885 -1.9605 4.0829 0.117 1.38 -0.945E+00 800. 0.284E+00 0.283E+00 0.237E-03 0.336E-13 14 4 0 10 4 0 16 23 17 0.000 0.000 0 328 0 96 00:22:46.22 11.4000A 1.8287 -1.8650 3.9342 0.118 1.38 -0.940E+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 329 0 97 00:22:57.62 12.0000A 1.8982 -1.7985 3.7406 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 330 0 98 00:23:09.62 12.6000A 1.9527 -1.8230 3.5839 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 331 0 99 00:23:22.22 12.9000A 1.9230 -1.9304 3.4804 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 332 0 100 00:23:35.12 13.2000A 1.8476 -2.0257 3.4144 0.121 1.38 -0.927E+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 333 0 101 00:23:48.32 13.4000A 1.7595 -2.1063 3.3486 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 334 0 102 00:24:01.72 13.7000A 1.7144 -2.1673 3.2912 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 335 0 103 00:24:15.42 13.9000A 1.7175 -2.2012 3.2333 0.122 1.35 -0.901E+00 750. 0.280E+00 0.279E+00 0.298E-03 0.426E-13 3 22 38 8 23 4 12 22 39 0.000 0.000 0 336 0 104 00:24:29.32 14.1000A 1.7200 -2.1940 3.1830 0.123 1.35 -0.895E+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 0 337 0 105 00:24:43.42 14.1000A 1.7151 -2.1666 3.1969 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 338 0 106 00:24:57.52 13.8000A 1.6954 -2.1250 3.2513 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 339 0 107 00:25:11.32 13.8000A -1.7019 -2.1113 3.2667 0.123 1.35 -0.890E+00 750. 0.279E+00 0.278E+00 0.295E-03 0.420E-13 8 25 1 8 23 5 11 23 1 0.000 0.000 0 340 0 108 00:25:25.12 13.8000A 1.7170 -2.1861 3.2552 0.123 1.35 -0.892E+00 750. 0.278E+00 0.278E+00 0.295E-03 0.417E-13 10 26 31 8 23 5 11 23 1 0.000 0.000 0 341 0 109 00:25:38.92 13.8000A 1.7572 -2.2020 3.2518 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 0 342 0 110 00:25:52.72 14.1000A 1.7772 -2.1520 3.1985 0.124 1.38 -0.902E+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 0 343 0 111 00:26:06.82 14.2000A 1.7739 -2.0928 3.1792 0.124 1.38 -0.896E+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 344 0 112 00:26:21.02 14.3000A -1.8071 -2.1438 3.1390 0.124 1.38 -0.894E+00 800. 0.277E+00 0.277E+00 0.300E-03 0.412E-13 16 22 38 8 22 5 12 23 19 0.000 0.000 0 345 0 113 00:26:35.32 14.4000A 1.7970 -2.1042 3.1245 0.124 1.38 -0.894E+00 800. 0.277E+00 0.276E+00 0.301E-03 0.411E-13 1 36 26 8 22 5 12 23 19 0.000 0.000 0 346 0 114 00:26:49.72 14.4000A 1.8521 -1.9823 3.1277 0.124 1.38 -0.894E+00 800. 0.276E+00 0.276E+00 0.301E-03 0.410E-13 1 37 26 8 22 5 9 38 26 0.000 0.000 0 347 0 115 00:27:04.12 14.1000A 1.8655 1.8641 3.2017 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 348 0 116 00:27:18.22 13.8000A 1.8259 1.8684 3.2623 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 0 349 0 117 00:27:32.02 13.7000A 1.8466 1.9007 3.2755 0.125 1.38 -0.899E+00 800. 0.274E+00 0.274E+00 0.285E-03 0.386E-13 1 16 38 8 35 23 10 38 26 0.000 0.000 0 350 0 118 00:27:45.72 13.8000A 1.9138 1.9282 3.2586 0.125 1.38 -0.900E+00 800. 0.274E+00 0.273E+00 0.282E-03 0.382E-13 1 16 39 8 35 23 10 38 26 0.000 0.000 0 351 0 119 00:27:59.52 13.8000A 1.9898 1.9444 3.2550 0.126 1.38 -0.899E+00 800. 0.273E+00 0.272E+00 0.283E-03 0.384E-13 1 16 39 8 36 23 11 38 26 0.000 0.000 0 352 0 120 00:28:13.32 13.9000A 2.0594 2.0384 3.2273 0.126 1.38 -0.897E+00 800. 0.272E+00 0.271E+00 0.283E-03 0.384E-13 1 16 39 9 36 24 11 38 26 0.000 0.000 0 353 0 121 00:28:27.22 14.0000A 2.1099 2.0924 3.2172 0.126 1.41 -0.895E+00 850. 0.271E+00 0.271E+00 0.285E-03 0.386E-13 1 16 39 9 36 24 13 37 36 0.000 0.000 0 354 0 122 00:28:41.22 14.1000A 2.1217 2.1131 3.1998 0.126 1.41 -0.895E+00 850. 0.270E+00 0.270E+00 0.287E-03 0.388E-13 1 16 39 9 36 23 13 37 37 0.000 0.000 0 355 0 123 00:28:55.32 14.0000A 2.1581 2.1139 3.2134 0.127 1.41 -0.892E+00 850. 0.269E+00 0.268E+00 0.288E-03 0.390E-13 1 16 0 9 36 23 13 37 37 0.000 0.000 0 356 0 124 00:29:09.32 13.9000A 2.1890 2.1442 3.2397 0.127 1.41 -0.890E+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 0 357 0 125 00:29:23.22 13.9000A 2.1613 2.1310 3.2296 0.127 1.41 -0.887E+00 850. 0.266E+00 0.266E+00 0.284E-03 0.383E-13 1 16 0 9 37 23 14 37 37 0.000 0.000 0 358 0 126 00:29:37.12 14.2000A 2.1070 2.0393 3.1733 0.128 1.38 -0.889E+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 0 359 0 127 00:29:51.32 14.3000A 2.1120 1.8897 3.1443 0.128 1.38 -0.900E+00 800. 0.263E+00 0.262E+00 0.291E-03 0.389E-13 1 16 1 9 37 23 13 6 16 0.000 0.000 0 360 0 128 00:30:05.62 14.4000A 2.0592 -1.9142 3.1322 0.128 1.38 -0.887E+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 361 0 129 00:30:20.02 14.0000A -2.0568 -2.0174 3.2119 0.128 1.38 -0.884E+00 800. 0.259E+00 0.258E+00 0.294E-03 0.391E-13 9 10 18 10 6 21 7 25 19 0.000 0.000 0 362 0 130 00:30:34.02 13.7000A -2.0451 -2.0359 3.2843 0.129 1.38 -0.882E+00 800. 0.257E+00 0.256E+00 0.285E-03 0.378E-13 9 10 18 10 6 21 8 25 19 0.000 0.000 0 363 0 131 00:30:47.72 13.3000A 1.9259 -2.1009 3.3730 0.129 1.41 -0.884E+00 850. 0.254E+00 0.254E+00 0.279E-03 0.368E-13 1 17 1 10 5 21 8 25 19 0.000 0.000 0 364 0 132 00:31:01.02 13.1000A -1.9186 -2.1189 3.4240 0.129 1.41 -0.878E+00 850. 0.252E+00 0.251E+00 0.270E-03 0.355E-13 10 9 18 10 5 21 9 25 19 0.000 0.000 0 365 0 133 00:31:14.12 13.0000A -1.9597 -2.0592 3.4698 0.129 1.41 -0.877E+00 850. 0.250E+00 0.249E+00 0.266E-03 0.348E-13 10 9 18 10 5 21 9 25 19 0.000 0.000 0 366 0 134 00:31:27.12 13.0000A -1.8415 -1.9297 3.4578 0.129 1.41 -0.880E+00 850. 0.247E+00 0.247E+00 0.263E-03 0.343E-13 10 9 18 10 5 21 10 25 19 0.000 0.000 0 367 0 135 00:31:40.12 13.0000A -1.9051 -1.9082 3.4511 0.128 1.41 -0.882E+00 850. 0.245E+00 0.244E+00 0.262E-03 0.340E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0 368 0 136 00:31:53.12 13.4000A -1.9185 -1.8703 3.3702 0.128 1.41 -0.889E+00 850. 0.242E+00 0.242E+00 0.261E-03 0.338E-13 13 26 20 6 12 21 10 25 19 0.000 0.000 0 369 0 137 00:32:06.52 13.7000A -1.9511 1.8787 3.2846 0.128 1.41 -0.895E+00 850. 0.240E+00 0.239E+00 0.268E-03 0.346E-13 13 26 19 10 10 12 11 25 19 0.000 0.000 0 370 0 138 00:32:20.22 13.9000A -2.0414 -1.8778 3.2356 0.127 1.41 -0.896E+00 850. 0.237E+00 0.236E+00 0.273E-03 0.351E-13 13 26 19 1 13 16 9 9 19 0.000 0.000 0 371 0 139 00:32:34.12 14.1000A -2.0833 -1.8864 3.2006 0.127 1.41 -0.897E+00 850. 0.234E+00 0.234E+00 0.275E-03 0.353E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0 372 0 140 00:32:48.22 14.2000A -2.1165 -1.8950 3.1741 0.127 1.41 -0.899E+00 850. 0.232E+00 0.231E+00 0.278E-03 0.355E-13 13 26 19 1 13 16 10 9 19 0.000 0.000 0 373 0 141 00:33:02.42 14.7000A -2.1634 -1.9011 3.0682 0.126 1.41 -0.901E+00 850. 0.229E+00 0.229E+00 0.278E-03 0.355E-13 14 26 19 1 13 16 10 8 19 0.000 0.000 0 374 0 142 00:33:17.12 14.7000A -2.2424 -1.9018 3.0542 0.126 1.41 -0.905E+00 850. 0.227E+00 0.226E+00 0.287E-03 0.363E-13 14 26 19 1 13 16 10 8 19 0.000 0.000 0 375 0 143 00:33:31.82 14.7000A -2.3446 -1.8939 3.0618 0.125 1.41 -0.909E+00 850. 0.225E+00 0.224E+00 0.286E-03 0.360E-13 14 26 19 1 13 16 11 8 19 0.000 0.000 0 376 0 144 00:33:46.52 14.7000A -2.4366 -1.8746 3.0702 0.125 1.41 -0.916E+00 850. 0.223E+00 0.222E+00 0.285E-03 0.356E-13 14 26 18 1 13 16 11 8 19 0.000 0.000 0 377 0 145 00:34:01.22 14.9000A -2.5232 -1.8732 3.0285 0.124 1.38 -0.914E+00 800. 0.221E+00 0.220E+00 0.284E-03 0.352E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0 378 0 146 00:34:16.12 15.2000A -2.5479 -1.8598 2.9686 0.124 1.38 -0.911E+00 800. 0.219E+00 0.219E+00 0.286E-03 0.354E-13 14 26 18 1 12 16 11 8 19 0.000 0.000 0 379 0 147 00:34:31.32 15.2000A -2.4724 -1.8309 2.9537 0.123 1.38 -0.912E+00 800. 0.218E+00 0.217E+00 0.291E-03 0.357E-13 14 26 18 1 12 16 12 8 19 0.000 0.000 0 380 0 148 00:34:46.52 15.3000A -2.3839 -1.7876 2.9361 0.123 1.38 -0.916E+00 800. 0.216E+00 0.216E+00 0.290E-03 0.354E-13 15 26 18 1 12 16 12 8 19 0.000 0.000 0 381 0 149 00:35:01.82 15.5000A -2.4539 -1.7567 2.9078 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 382 0 150 00:35:17.32 15.7000A -2.5193 1.7234 2.8580 0.121 1.38 -0.923E+00 800. 0.214E+00 0.214E+00 0.294E-03 0.354E-13 15 26 17 8 34 23 12 8 19 0.000 0.000 0 383 0 151 00:35:33.02 15.8000A -2.4903 1.7913 2.8419 0.121 1.38 -0.928E+00 800. 0.214E+00 0.213E+00 0.297E-03 0.355E-13 15 26 17 8 34 23 14 26 24 0.000 0.000 0 384 0 152 00:35:48.82 15.8000A -2.3635 -1.7756 2.8479 0.120 1.38 -0.929E+00 800. 0.213E+00 0.212E+00 0.299E-03 0.354E-13 15 26 17 1 7 34 6 38 32 0.000 0.000 0 385 0 153 00:36:04.62 15.3000A -2.2428 -1.7658 2.9398 0.119 1.38 -0.931E+00 800. 0.212E+00 0.212E+00 0.299E-03 0.352E-13 15 26 16 1 7 34 6 38 32 0.000 0.000 0 386 0 154 00:36:19.92 14.8000A -2.0822 -1.8215 3.0479 0.119 1.38 -0.939E+00 800. 0.212E+00 0.211E+00 0.290E-03 0.339E-13 15 26 16 1 6 34 7 38 32 0.000 0.000 0 387 0 155 00:36:34.72 14.3000A -2.0180 -1.8586 3.1546 0.118 1.38 -0.940E+00 800. 0.212E+00 0.211E+00 0.280E-03 0.327E-13 11 11 13 1 6 34 7 38 32 0.000 0.000 0 388 0 156 00:36:49.02 13.9000A -2.0312 -1.8315 3.2277 0.118 1.38 -0.945E+00 800. 0.212E+00 0.211E+00 0.271E-03 0.315E-13 9 12 14 1 6 34 7 38 32 0.000 0.000 0 389 0 157 00:37:02.92 13.7000A -2.1165 -1.9066 3.2855 0.117 1.38 -0.951E+00 800. 0.212E+00 0.211E+00 0.264E-03 0.306E-13 10 12 13 1 5 34 8 38 32 0.000 0.000 0 390 0 158 00:37:16.62 13.5000A -2.1128 -1.9114 3.3441 0.117 1.38 -0.953E+00 800. 0.211E+00 0.211E+00 0.260E-03 0.301E-13 11 12 13 1 5 34 8 38 32 0.000 0.000 0 391 0 159 00:37:30.12 13.4000A -2.1380 -1.8734 3.3540 0.116 1.38 -0.958E+00 800. 0.211E+00 0.210E+00 0.257E-03 0.296E-13 11 12 13 1 6 35 8 38 32 0.000 0.000 0 392 0 160 00:37:43.52 13.5000A -2.1108 -1.8757 3.3222 0.116 1.38 -0.963E+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 393 0 161 00:37:57.02 13.7000A -2.0627 -1.9699 3.2944 0.115 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 0 394 0 162 00:38:10.72 14.0000A -1.9774 -1.9542 3.2168 0.115 1.41 -0.970E+00 850. 0.211E+00 0.210E+00 0.261E-03 0.296E-13 13 12 12 1 5 35 9 38 32 0.000 0.000 0 395 0 163 00:38:24.72 14.1000A 1.9456 -2.0340 3.2004 0.115 1.41 -0.974E+00 850. 0.211E+00 0.210E+00 0.266E-03 0.301E-13 7 27 22 1 4 35 8 10 7 0.000 0.000 0 396 0 164 00:38:38.82 13.3000A -1.9424 -2.0740 3.3950 0.115 1.41 -0.973E+00 850. 0.211E+00 0.210E+00 0.269E-03 0.302E-13 12 13 12 1 4 35 8 9 7 0.000 0.000 0 397 0 165 00:38:52.12 12.6000A 1.9909 -2.0064 3.5667 0.114 1.41 -0.976E+00 850. 0.211E+00 0.210E+00 0.255E-03 0.284E-13 7 28 22 1 4 35 9 9 7 0.000 0.000 0 398 0 166 00:39:04.72 12.8000A 1.9909 -2.0271 3.5276 0.114 1.41 -0.980E+00 850. 0.211E+00 0.210E+00 0.242E-03 0.269E-13 7 28 22 1 3 35 9 9 7 0.000 0.000 0 399 0 167 00:39:17.52 12.8000A 1.9283 -2.0299 3.5283 0.114 1.41 -0.983E+00 850. 0.211E+00 0.210E+00 0.247E-03 0.273E-13 7 28 22 1 3 35 10 9 7 0.000 0.000 0 400 0 168 00:39:30.32 12.9000A -1.8906 -1.9463 3.4940 0.114 1.41 -0.984E+00 850. 0.211E+00 0.210E+00 0.247E-03 0.273E-13 8 13 15 1 3 35 10 8 6 0.000 0.000 0 401 0 169 00:39:43.22 13.4000A -1.9064 1.9409 3.3546 0.114 1.41 -0.985E+00 850. 0.212E+00 0.211E+00 0.250E-03 0.276E-13 8 13 14 17 15 11 10 8 6 0.000 0.000 0 402 0 170 00:39:56.62 13.5000A -1.9332 2.0092 3.3236 0.113 1.41 -0.989E+00 850. 0.212E+00 0.211E+00 0.261E-03 0.286E-13 8 13 14 17 15 11 9 8 7 0.000 0.000 0 403 0 171 00:40:10.12 13.0000A -1.9327 2.0173 3.4509 0.113 1.41 -0.989E+00 850. 0.212E+00 0.211E+00 0.263E-03 0.289E-13 8 13 14 17 15 11 10 9 8 0.000 0.000 0 404 0 172 00:40:23.12 13.1000A -1.9254 2.0677 3.4468 0.113 1.41 -0.989E+00 850. 0.212E+00 0.211E+00 0.255E-03 0.278E-13 9 13 14 17 16 11 10 9 8 0.000 0.000 0 405 0 173 00:40:36.22 13.2000A -1.9172 2.1591 3.4172 0.113 1.41 -0.991E+00 850. 0.212E+00 0.211E+00 0.258E-03 0.281E-13 9 13 14 17 16 11 11 8 7 0.000 0.000 0 406 0 174 00:40:49.42 13.6000A -1.9459 2.1733 3.3175 0.113 1.41 -0.993E+00 850. 0.213E+00 0.211E+00 0.261E-03 0.283E-13 18 13 9 17 16 11 11 8 7 0.000 0.000 0 407 0 175 00:41:03.02 13.8000A -1.9651 2.0985 3.2609 0.113 1.41 -0.995E+00 850. 0.213E+00 0.212E+00 0.270E-03 0.293E-13 18 13 9 17 17 11 11 8 8 0.000 0.000 0 408 0 176 00:41:16.82 13.6000A -1.9514 2.1932 3.3179 0.113 1.41 -0.999E+00 850. 0.213E+00 0.212E+00 0.275E-03 0.298E-13 18 13 9 17 17 11 11 8 8 0.000 0.000 0 409 0 177 00:41:30.42 13.6000A 1.9172 2.2178 3.3025 0.113 1.41 -0.100E+01 850. 0.213E+00 0.212E+00 0.272E-03 0.294E-13 12 5 12 17 17 11 12 8 8 0.000 0.000 0 410 0 178 00:41:44.02 13.5000A -1.9023 2.2163 3.3328 0.113 1.41 -0.101E+01 850. 0.213E+00 0.212E+00 0.273E-03 0.295E-13 18 14 5 18 18 9 12 8 8 0.000 0.000 0 411 0 179 00:41:57.52 13.7000A -1.9313 2.2157 3.2800 0.113 1.41 -0.101E+01 850. 0.214E+00 0.212E+00 0.272E-03 0.294E-13 18 15 5 17 18 10 12 8 8 0.000 0.000 0 412 0 180 00:42:11.22 13.9000A -1.9220 2.1850 3.2404 0.113 1.38 -0.101E+01 800. 0.214E+00 0.213E+00 0.277E-03 0.298E-13 18 15 5 18 18 9 13 8 8 0.000 0.000 0 413 0 181 00:42:25.12 14.1000A -1.9163 2.2124 3.1917 0.113 1.38 -0.102E+01 800. 0.215E+00 0.213E+00 0.281E-03 0.303E-13 18 14 7 17 19 10 13 8 8 0.000 0.000 0 414 0 182 00:42:39.22 14.6000A -1.9486 2.2088 3.0894 0.113 1.38 -0.101E+01 800. 0.216E+00 0.214E+00 0.286E-03 0.308E-13 18 14 7 18 19 8 13 8 8 0.000 0.000 0 415 0 183 00:42:53.82 15.0000A -1.8983 2.1797 3.0026 0.113 1.38 -0.101E+01 800. 0.216E+00 0.215E+00 0.296E-03 0.319E-13 18 14 7 17 20 9 14 8 8 0.000 0.000 0 416 0 184 00:43:08.82 15.3000A -1.8954 2.2107 2.9328 0.113 1.38 -0.100E+01 800. 0.217E+00 0.216E+00 0.304E-03 0.328E-13 18 15 7 17 20 9 14 8 8 0.000 0.000 0 417 0 185 00:43:24.12 15.8000A -1.8630 2.1607 2.8548 0.114 1.38 -0.100E+01 800. 0.218E+00 0.217E+00 0.310E-03 0.334E-13 18 15 6 17 20 9 14 8 8 0.000 0.000 0 418 0 186 00:43:39.92 15.6000A -1.8678 2.0965 2.8897 0.114 1.38 -0.999E+00 800. 0.220E+00 0.218E+00 0.320E-03 0.344E-13 17 14 4 18 18 12 9 9 18 0.000 0.000 0 419 0 187 00:43:55.52 15.0000A -1.8314 2.0690 2.9995 0.114 1.38 -0.994E+00 800. 0.221E+00 0.219E+00 0.316E-03 0.340E-13 17 14 4 18 18 12 10 9 18 0.000 0.000 0 420 0 188 00:44:10.52 14.8000A -1.7310 2.1297 3.0403 0.114 1.38 -0.994E+00 800. 0.222E+00 0.221E+00 0.304E-03 0.328E-13 17 14 4 18 19 11 11 9 18 0.000 0.000 0 421 0 189 00:44:25.32 14.6000A -1.7140 2.1094 3.0835 0.115 1.38 -0.990E+00 800. 0.224E+00 0.222E+00 0.301E-03 0.324E-13 16 15 3 18 17 13 13 11 16 0.000 0.000 0 422 0 190 00:44:39.92 14.2000A 1.7255 2.1319 3.1631 0.115 1.38 -0.991E+00 800. 0.225E+00 0.224E+00 0.298E-03 0.320E-13 9 27 23 18 17 13 13 11 16 0.000 0.000 0 423 0 191 00:44:54.12 13.9000A 1.7630 2.1271 3.2369 0.115 1.38 -0.993E+00 800. 0.227E+00 0.225E+00 0.291E-03 0.313E-13 9 27 23 18 17 13 13 11 16 0.000 0.000 0 424 0 192 00:45:08.02 13.6000A 1.7831 2.1534 3.2990 0.116 1.43 -0.991E+00 900. 0.228E+00 0.227E+00 0.286E-03 0.307E-13 9 27 23 18 18 13 14 11 16 0.000 0.000 0 425 0 193 00:45:21.62 13.4000A 1.7944 2.1802 3.3668 0.116 1.43 -0.100E+01 900. 0.230E+00 0.228E+00 0.281E-03 0.301E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0 426 0 194 00:45:35.02 13.2000A 1.8443 2.1889 3.4074 0.116 1.43 -0.996E+00 900. 0.231E+00 0.230E+00 0.278E-03 0.298E-13 9 27 24 18 18 13 14 11 16 0.000 0.000 0 427 0 195 00:45:48.22 13.2000A 1.8836 2.1814 3.4186 0.117 1.43 -0.991E+00 900. 0.233E+00 0.232E+00 0.274E-03 0.295E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 428 0 196 00:46:01.42 13.0000A 1.9089 2.1575 3.4497 0.117 1.43 -0.987E+00 900. 0.235E+00 0.233E+00 0.276E-03 0.296E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 429 0 197 00:46:14.42 13.1000A 1.9176 2.1154 3.4479 0.117 1.43 -0.985E+00 900. 0.236E+00 0.235E+00 0.273E-03 0.293E-13 9 27 24 18 18 13 15 11 16 0.000 0.000 0 430 0 198 00:46:27.52 13.2000A 1.9085 2.0760 3.4203 0.118 1.43 -0.986E+00 900. 0.238E+00 0.237E+00 0.276E-03 0.296E-13 9 27 24 18 19 12 15 11 16 0.000 0.000 0 431 0 199 00:46:40.72 13.2000A 1.8855 2.0566 3.4206 0.118 1.43 -0.987E+00 900. 0.240E+00 0.239E+00 0.279E-03 0.300E-13 9 26 25 18 19 12 16 11 16 0.000 0.000 0 432 0 200 00:46:53.92 13.2000A 1.9077 2.0374 3.4011 0.118 1.43 -0.983E+00 900. 0.242E+00 0.241E+00 0.280E-03 0.301E-13 9 27 25 8 3 0 16 11 16 0.000 0.000 0 433 0 201 00:47:07.12 13.4000A 1.9317 2.0179 3.3642 0.119 1.41 -0.978E+00 850. 0.244E+00 0.242E+00 0.281E-03 0.303E-13 9 27 25 8 4 0 16 11 16 0.000 0.000 0 434 0 202 00:47:20.52 13.4000A 1.9346 2.0385 3.3511 0.119 1.41 -0.974E+00 850. 0.245E+00 0.244E+00 0.287E-03 0.310E-13 9 27 25 19 18 14 15 10 16 0.000 0.000 0 435 0 203 00:47:33.92 13.3000A 1.9127 2.0929 3.3850 0.119 1.41 -0.977E+00 850. 0.247E+00 0.246E+00 0.288E-03 0.312E-13 9 27 25 19 18 14 15 9 17 0.000 0.000 0 436 0 204 00:47:47.22 13.2000A 1.9149 2.0377 3.4163 0.119 1.41 -0.974E+00 850. 0.249E+00 0.248E+00 0.288E-03 0.312E-13 9 27 26 19 18 14 15 9 17 0.000 0.000 0 437 0 205 00:48:00.42 13.2000A 1.9285 1.9118 3.4195 0.119 1.41 -0.972E+00 850. 0.251E+00 0.250E+00 0.287E-03 0.312E-13 9 27 26 19 19 13 15 9 17 0.000 0.000 0 438 0 206 00:48:13.62 13.1000A -2.0010 1.9771 3.4362 0.119 1.41 -0.973E+00 850. 0.253E+00 0.251E+00 0.288E-03 0.314E-13 6 3 1 17 13 10 16 9 17 0.000 0.000 0 439 0 207 00:48:26.72 13.1000A -2.0771 1.9674 3.4311 0.119 1.41 -0.982E+00 850. 0.255E+00 0.253E+00 0.287E-03 0.314E-13 6 3 0 17 13 10 16 9 17 0.000 0.000 0 440 0 208 00:48:39.82 13.2000A -2.1797 1.9858 3.4091 0.119 1.41 -0.977E+00 850. 0.257E+00 0.255E+00 0.289E-03 0.316E-13 6 3 0 17 14 10 13 9 18 0.000 0.000 0 441 0 209 00:48:53.01 13.2000A -2.1964 2.0063 3.4200 0.119 1.41 -0.973E+00 850. 0.259E+00 0.257E+00 0.293E-03 0.322E-13 6 3 0 17 14 10 8 39 1 0.000 0.000 0 442 0 210 00:49:06.21 13.1000A -2.2503 1.9578 3.4393 0.119 1.41 -0.973E+00 850. 0.261E+00 0.259E+00 0.295E-03 0.325E-13 7 2 1 17 15 10 9 39 1 0.000 0.000 0 443 0 211 00:49:19.31 13.3000A -2.2530 -1.9685 3.3805 0.119 1.41 -0.972E+00 850. 0.262E+00 0.260E+00 0.294E-03 0.325E-13 7 2 1 3 8 2 7 8 17 0.000 0.000 0 444 0 212 00:49:32.61 12.9000A -2.3251 -2.0228 3.4957 0.119 1.41 -0.974E+00 850. 0.264E+00 0.262E+00 0.300E-03 0.333E-13 7 3 0 3 8 2 9 0 0 0.000 0.000 0 445 0 213 00:49:45.51 12.7000A -2.3313 -2.0897 3.5383 0.119 1.41 -0.987E+00 850. 0.266E+00 0.264E+00 0.293E-03 0.326E-13 8 2 0 3 7 2 8 8 17 0.000 0.000 0 446 0 214 00:49:58.21 12.4000A -2.4812 -2.1653 3.6189 0.119 1.41 -0.995E+00 850. 0.267E+00 0.265E+00 0.290E-03 0.323E-13 8 1 0 3 7 2 10 7 17 0.000 0.000 0 447 0 215 00:50:10.61 12.0000A -2.5128 -2.1854 3.7471 0.119 1.41 -0.974E+00 850. 0.269E+00 0.267E+00 0.284E-03 0.318E-13 8 1 0 3 7 2 10 7 17 0.000 0.000 0 448 0 216 00:50:22.61 11.8000A -2.4764 -2.1261 3.7981 0.119 1.41 -0.973E+00 850. 0.270E+00 0.268E+00 0.276E-03 0.310E-13 9 1 39 3 7 2 10 7 17 0.000 0.000 0 449 0 217 00:50:34.42 11.9000A -2.6468 -2.0821 3.7849 0.118 1.43 -0.973E+00 900. 0.272E+00 0.269E+00 0.273E-03 0.307E-13 9 1 39 3 6 2 10 7 17 0.000 0.000 0 450 0 218 00:50:46.32 11.8000A -2.6822 -2.0757 3.8043 0.118 1.43 -0.970E+00 900. 0.273E+00 0.271E+00 0.276E-03 0.311E-13 9 1 39 3 6 2 11 7 17 0.000 0.000 0 451 0 219 00:50:58.12 11.9000A -2.6288 -2.0878 3.7871 0.118 1.43 -0.968E+00 900. 0.274E+00 0.272E+00 0.275E-03 0.310E-13 9 1 39 3 7 3 11 7 17 0.000 0.000 0 452 0 220 00:51:10.02 12.0000A -2.6844 -2.0617 3.7357 0.118 1.43 -0.969E+00 900. 0.275E+00 0.273E+00 0.278E-03 0.313E-13 9 1 38 3 7 3 14 6 18 0.000 0.000 0 453 0 221 00:51:22.02 11.9000A -2.6537 -1.9889 3.7962 0.118 1.43 -0.972E+00 900. 0.276E+00 0.274E+00 0.282E-03 0.317E-13 9 1 38 3 7 3 14 6 18 0.000 0.000 0 454 0 222 00:51:33.92 11.8000A -2.5385 -1.9568 3.7990 0.118 1.43 -0.977E+00 900. 0.278E+00 0.275E+00 0.280E-03 0.316E-13 9 1 38 3 6 3 14 6 18 0.000 0.000 0 455 0 223 00:51:45.72 11.8000A -2.4135 -1.9337 3.8177 0.118 1.43 -0.979E+00 900. 0.279E+00 0.276E+00 0.279E-03 0.314E-13 11 0 38 3 6 3 15 6 18 0.000 0.000 0 456 0 224 00:51:57.52 11.9000A -2.4027 -1.9037 3.7946 0.117 1.43 -0.982E+00 900. 0.280E+00 0.277E+00 0.279E-03 0.315E-13 11 0 37 15 32 38 15 6 18 0.000 0.000 0 457 0 225 00:52:09.42 12.1000A -2.4532 -1.8845 3.7186 0.117 1.43 -0.981E+00 900. 0.280E+00 0.278E+00 0.282E-03 0.319E-13 11 0 37 16 32 39 15 6 18 0.000 0.000 0 458 0 226 00:52:21.52 12.2000A -2.4193 -1.8701 3.6911 0.117 1.43 -0.982E+00 900. 0.281E+00 0.279E+00 0.288E-03 0.326E-13 11 0 37 16 32 39 16 5 18 0.000 0.000 0 459 0 227 00:52:33.72 12.4000A -2.3308 -1.8827 3.6176 0.117 1.43 -0.982E+00 900. 0.282E+00 0.279E+00 0.291E-03 0.329E-13 11 0 37 3 5 1 16 5 18 0.000 0.000 0 460 0 228 00:52:46.12 12.6000A -2.2282 1.9090 3.5577 0.117 1.43 -0.982E+00 900. 0.282E+00 0.280E+00 0.296E-03 0.335E-13 12 0 36 17 16 10 17 5 18 0.000 0.000 0 461 0 229 00:52:58.72 13.1000A -2.1808 1.9565 3.4457 0.117 1.43 -0.984E+00 900. 0.283E+00 0.280E+00 0.301E-03 0.341E-13 12 0 36 17 16 10 17 5 18 0.000 0.000 0 462 0 230 00:53:11.82 13.7000A -2.0728 1.9749 3.2822 0.116 1.43 -0.985E+00 900. 0.283E+00 0.280E+00 0.314E-03 0.354E-13 12 0 36 17 16 10 15 36 39 0.000 0.000 0 463 0 231 00:53:25.52 13.6000A -1.9391 1.9769 3.3165 0.116 1.43 -0.987E+00 900. 0.283E+00 0.280E+00 0.328E-03 0.369E-13 12 0 36 17 17 11 9 35 5 0.000 0.000 0 464 0 232 00:53:39.12 13.2000A -1.8656 1.9871 3.4090 0.116 1.43 -0.988E+00 900. 0.283E+00 0.281E+00 0.326E-03 0.365E-13 12 1 39 17 17 11 9 35 5 0.000 0.000 0 465 0 233 00:53:52.32 13.2000A -1.9476 1.9543 3.4211 0.116 1.43 -0.988E+00 900. 0.283E+00 0.281E+00 0.316E-03 0.354E-13 9 6 4 17 17 11 10 35 5 0.000 0.000 0 466 0 234 00:54:05.51 12.8000A -2.0044 1.9320 3.5149 0.116 1.43 -0.988E+00 900. 0.283E+00 0.280E+00 0.316E-03 0.353E-13 9 6 4 14 8 5 11 6 10 0.000 0.000 0 467 0 235 00:54:18.32 13.0000A -2.1071 1.9464 3.4637 0.115 1.43 -0.989E+00 900. 0.283E+00 0.280E+00 0.307E-03 0.342E-13 10 6 4 14 9 5 10 35 5 0.000 0.000 0 468 0 236 00:54:31.32 13.0000A -2.2080 2.0654 3.4627 0.115 1.43 -0.992E+00 900. 0.283E+00 0.280E+00 0.311E-03 0.346E-13 10 6 4 11 2 38 9 36 6 0.000 0.000 0 469 0 237 00:54:44.32 12.5000A -2.2464 2.1078 3.6015 0.115 1.43 -0.992E+00 900. 0.283E+00 0.280E+00 0.311E-03 0.346E-13 10 6 4 11 2 38 10 36 6 0.000 0.000 0 470 0 238 00:54:56.82 12.3000A -2.2361 2.0862 3.6655 0.115 1.40 -0.991E+00 850. 0.283E+00 0.280E+00 0.300E-03 0.332E-13 11 6 4 11 2 38 10 36 6 0.000 0.000 0 471 0 239 00:55:09.12 12.2000A -2.2916 2.0375 3.6744 0.115 1.40 -0.991E+00 850. 0.283E+00 0.280E+00 0.296E-03 0.327E-13 11 6 4 11 2 39 11 36 6 0.000 0.000 0 472 0 240 00:55:21.32 12.1000A -2.2994 2.0789 3.7339 0.115 1.40 -0.995E+00 850. 0.283E+00 0.280E+00 0.294E-03 0.324E-13 11 6 4 12 2 38 11 36 6 0.000 0.000 0 473 0 241 00:55:33.42 12.0000A -2.2669 2.1541 3.7608 0.115 1.40 -0.993E+00 850. 0.283E+00 0.280E+00 0.293E-03 0.321E-13 11 6 3 12 3 38 11 36 6 0.000 0.000 0 474 0 242 00:55:45.42 12.0000A -2.2433 2.2386 3.7470 0.115 1.40 -0.991E+00 850. 0.283E+00 0.280E+00 0.291E-03 0.319E-13 11 6 3 12 3 37 11 36 6 0.000 0.000 0 475 0 243 00:55:57.42 11.8000A -2.1748 2.2832 3.8004 0.115 1.40 -0.987E+00 850. 0.283E+00 0.280E+00 0.292E-03 0.319E-13 12 6 3 12 3 37 12 36 7 0.000 0.000 0 476 0 244 00:56:09.22 11.8000A -2.2169 2.2247 3.8031 0.115 1.40 -0.984E+00 850. 0.283E+00 0.280E+00 0.287E-03 0.314E-13 12 6 3 12 3 37 12 36 7 0.000 0.000 0 477 0 245 00:56:21.02 12.0000A -2.2126 2.2100 3.7549 0.115 1.40 -0.983E+00 850. 0.282E+00 0.280E+00 0.287E-03 0.314E-13 12 6 3 13 3 38 12 36 7 0.000 0.000 0 478 0 246 00:56:33.02 12.0000A -2.1607 2.2718 3.7616 0.115 1.43 -0.985E+00 900. 0.282E+00 0.279E+00 0.292E-03 0.319E-13 12 6 3 13 4 37 13 36 7 0.000 0.000 0 479 0 247 00:56:45.02 12.1000A -2.0640 2.3360 3.7201 0.115 1.43 -0.988E+00 900. 0.281E+00 0.279E+00 0.292E-03 0.320E-13 12 6 3 13 4 37 13 36 7 0.000 0.000 0 480 0 248 00:56:57.12 12.4000A -2.0537 2.3067 3.6170 0.115 1.43 -0.994E+00 900. 0.281E+00 0.278E+00 0.295E-03 0.323E-13 13 6 3 13 4 37 13 36 7 0.000 0.000 0 481 0 249 00:57:09.52 12.8000A -2.0386 2.2491 3.5073 0.115 1.43 -0.986E+00 900. 0.280E+00 0.278E+00 0.303E-03 0.331E-13 13 6 3 13 5 37 14 35 6 0.000 0.000 0 482 0 250 00:57:22.32 12.9000A 2.0738 2.1661 3.4923 0.115 1.43 -0.988E+00 900. 0.280E+00 0.277E+00 0.313E-03 0.343E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0 483 0 251 00:57:35.22 13.1000A 2.1173 2.2291 3.4375 0.115 1.43 -0.985E+00 900. 0.279E+00 0.277E+00 0.316E-03 0.347E-13 13 37 9 14 5 37 14 35 6 0.000 0.000 0 484 0 252 00:57:48.32 13.4000A 2.1033 2.2274 3.3545 0.115 1.43 -0.989E+00 900. 0.279E+00 0.276E+00 0.321E-03 0.353E-13 13 37 9 14 6 37 15 35 6 0.000 0.000 0 485 0 253 00:58:01.72 13.7000A 2.2294 2.2239 3.2838 0.115 1.46 -0.996E+00 950. 0.279E+00 0.276E+00 0.329E-03 0.362E-13 14 37 10 14 6 37 15 35 6 0.000 0.000 0 486 0 254 00:58:15.42 13.8000A 2.2842 2.1522 3.2636 0.115 1.40 -0.993E+00 850. 0.278E+00 0.276E+00 0.337E-03 0.371E-13 14 37 10 14 6 37 16 36 10 0.000 0.000 0 487 0 255 00:58:29.22 14.3000A 2.2566 2.1277 3.1445 0.115 1.40 -0.991E+00 850. 0.278E+00 0.276E+00 0.340E-03 0.375E-13 14 37 10 15 7 37 8 37 5 0.000 0.000 0 488 0 256 00:58:43.52 14.3000A 2.2013 2.1755 3.1545 0.115 1.38 -0.993E+00 800. 0.278E+00 0.275E+00 0.352E-03 0.389E-13 15 37 11 15 7 37 8 37 5 0.000 0.000 0 489 0 257 00:58:57.82 14.2000A 2.2568 2.1866 3.1704 0.115 1.38 -0.993E+00 800. 0.278E+00 0.275E+00 0.353E-03 0.389E-13 15 37 11 15 8 37 9 37 5 0.000 0.000 0 490 0 258 00:59:12.02 14.2000A 2.2312 2.1716 3.1684 0.115 1.38 -0.995E+00 800. 0.278E+00 0.275E+00 0.351E-03 0.387E-13 15 37 11 15 9 37 9 37 5 0.000 0.000 0 491 0 259 00:59:26.22 14.4000A 2.1733 2.1525 3.1350 0.115 1.38 -0.999E+00 800. 0.278E+00 0.276E+00 0.351E-03 0.387E-13 15 37 11 15 9 37 8 37 5 0.000 0.000 0 492 0 260 00:59:40.62 13.8000A 2.0889 2.1539 3.2687 0.115 1.38 -0.100E+01 800. 0.278E+00 0.276E+00 0.356E-03 0.392E-13 14 38 11 15 7 36 7 9 15 0.000 0.000 0 493 0 261 00:59:54.42 13.4000A 2.0555 2.1633 3.3538 0.115 1.38 -0.100E+01 800. 0.278E+00 0.276E+00 0.342E-03 0.375E-13 14 38 11 15 8 36 7 9 15 0.000 0.000 0 494 0 262 01:00:07.82 12.8000A 1.9438 2.1487 3.5253 0.115 1.38 -0.997E+00 800. 0.278E+00 0.276E+00 0.332E-03 0.363E-13 17 36 12 15 8 36 8 9 15 0.000 0.000 0 -
palm/trunk/SOURCE
- Property svn:mergeinfo changed
/palm/branches/rans/SOURCE merged: 2919,2922,2928-2929,2960,2962,2966,2976,2982,2987-2988,2991,3008-3009,3023,3047,3050,3059,3062,3071-3082
- Property svn:mergeinfo changed
-
palm/trunk/SOURCE/Makefile
- Property svn:mergeinfo changed
/palm/branches/rans/SOURCE/Makefile merged: 2919,2987,3081
r2963 r3083 25 25 # ----------------- 26 26 # $Id$ 27 # add turbulence_closure_mod to parin 28 # 29 # 2963 2018-04-12 14:47:44Z suehring 27 30 # Introduce index for vegetation/wall, pavement/green-wall and water/window 28 31 # surfaces, for clearer access of surface fraction, albedo, emissivity, etc. . … … 1230 1233 spectra_mod.o \ 1231 1234 synthetic_turbulence_generator_mod.o \ 1235 turbulence_closure_mod.o \ 1232 1236 uv_exposure_model_mod.o \ 1233 1237 vertical_nesting_mod.o \ - Property svn:mergeinfo changed
-
palm/trunk/SOURCE/check_parameters.f90
- Property svn:mergeinfo changed
/palm/branches/rans/SOURCE/check_parameters.f90 merged: 2919,2987,3081
r3065 r3083 24 24 ! Former revisions: 25 25 ! ----------------- 26 ! $Id$ 26 ! $Id: check_parameters.f90 2520 2017-10-05 13:50:26Z gronemeier & 27 ! Add inital profile output for e (TG) 28 ! 29 ! 3065 2018-06-12 07:03:02Z Giersch 27 30 ! dz was replaced by dz(1), error message revised 28 31 ! … … 2455 2458 ENDIF 2456 2459 2457 CASE ( 'e' )2460 CASE ( 'e', '#e' ) 2458 2461 dopr_index(i) = 8 2459 2462 dopr_unit(i) = 'm2/s2' 2460 2463 hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 ) 2461 2464 hom(nzb,2,8,:) = 0.0_wp 2465 IF ( data_output_pr(i)(1:1) == '#' ) THEN 2466 dopr_initial_index(i) = 8 2467 hom(:,2,8,:) = SPREAD( zu, 2, statistic_regions+1 ) 2468 data_output_pr(i) = data_output_pr(i)(2:) 2469 ENDIF 2462 2470 2463 2471 CASE ( 'km', '#km' ) - Property svn:mergeinfo changed
-
palm/trunk/SOURCE/data_output_mask.f90
r3045 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! 28 ! 29 ! 3045 2018-05-28 07:55:41Z Giersch 27 30 ! Error messages revised 28 31 ! -
palm/trunk/SOURCE/header.f90
r3065 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Print RANS-mode constants 28 ! 29 ! 3065 2018-06-12 07:03:02Z Giersch 27 30 ! Header output concerning stretching revised 28 31 ! … … 450 453 USE synthetic_turbulence_generator_mod, & 451 454 ONLY: stg_header 455 456 USE turbulence_closure_mod, & 457 ONLY: rans_const_c, rans_const_sigma 452 458 453 459 IMPLICIT NONE … … 1946 1952 IF ( wall_adjustment ) WRITE ( io, 453 ) wall_adjustment_factor 1947 1953 ENDIF 1948 1954 IF ( rans_mode ) THEN 1955 WRITE ( io, 457 ) rans_const_c, rans_const_sigma 1956 ENDIF 1949 1957 ! 1950 1958 !-- Special actions during the run … … 2388 2396 456 FORMAT (' Day of the year at model start : day_init = ',I3 & 2389 2397 /' UTC time at model start : time_utc_init = ',F7.1' s') 2398 457 FORMAT (' RANS-mode constants: c_0 = ',F9.5/ & 2399 ' c_1 = ',F9.5/ & 2400 ' c_2 = ',F9.5/ & 2401 ' c_3 = ',F9.5/ & 2402 ' c_4 = ',F9.5/ & 2403 ' sigma_e = ',F9.5/ & 2404 ' sigma_diss = ',F9.5) 2390 2405 470 FORMAT (//' Actions during the simulation:'/ & 2391 2406 ' -----------------------------'/) -
palm/trunk/SOURCE/init_3d_model.f90
- Property svn:mergeinfo changed
/palm/branches/rans/SOURCE/init_3d_model.f90 (added) merged: 2139,2244,2352-2353,2519-2521,2680,2760-2761,2842,2896,2917,3062,3081 /palm/trunk/SOURCE/init_3d_model.f90 removed
- Property svn:mergeinfo changed
-
palm/trunk/SOURCE/model_1d_mod.f90
r3049 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfixes: 28 ! - preset te_diss and te_e to avoid runtime errors 29 ! - implementation of buoyancy term to dissipation 30 ! according to Sogachev et al. (2012) 31 ! - where diss_p < 0 set diss_p = 0.1 diss 32 ! - calculate progn eq(diss) starting from nzb+1 33 ! Changes: 34 ! - add sig_e to TKE equation 35 ! - adjust prognostic equation of diss 36 ! - set model constants according to Koblitz (2013) 37 ! - renamed c_m to c_0 38 ! - rename l_black into l1d_init 39 ! - calculate l_grid within init_1d_model and save it as l1d_init 40 ! - calculate l1d according to DE85 if dissipation is a prognostic value 41 ! - made annotations doxygen-readable 42 ! 43 ! 3049 2018-05-29 13:52:36Z Giersch 27 44 ! Error messages revised 28 45 ! … … 164 181 USE kinds 165 182 166 USE pegrid 183 USE pegrid, & 184 ONLY: myid 167 185 168 186 … … 175 193 LOGICAL :: stop_dt_1d = .FALSE. !< termination flag, used in case of too small timestep (1d-model) 176 194 177 REAL(wp) :: c_1 = 1.44_wp !< model constant178 REAL(wp) :: c_ 2 = 1.92_wp !< model constant179 REAL(wp) :: c_ 3 = 1.44_wp !< model constant180 REAL(wp) :: c_ h = 0.0015_wp !< model constant according to Detering and Etling (1985)181 REAL(wp) :: c_ m = 0.4_wp !< model constant, 0.4 according to Detering and Etling (1985)182 REAL(wp) :: c_mu = 0.09_wp!< model constant195 REAL(wp) :: alpha_buoyancy !< model constant according to Koblitz (2013) 196 REAL(wp) :: c_0 = 0.03_wp**0.25_wp !< model constant according to Koblitz (2013) 197 REAL(wp) :: c_1 = 1.52_wp !< model constant according to Koblitz (2013) 198 REAL(wp) :: c_2 = 1.83_wp !< model constant according to Koblitz (2013) 199 REAL(wp) :: c_3 !< model constant 200 REAL(wp) :: c_mu !< model constant 183 201 REAL(wp) :: damp_level_1d = -1.0_wp !< namelist parameter 184 202 REAL(wp) :: dt_1d = 60.0_wp !< dynamic timestep (1d-model) … … 187 205 REAL(wp) :: dt_run_control_1d = 60.0_wp !< namelist parameter 188 206 REAL(wp) :: end_time_1d = 864000.0_wp !< namelist parameter 207 REAL(wp) :: lambda !< maximum mixing length 189 208 REAL(wp) :: qs1d !< characteristic humidity scale (1d-model) 190 209 REAL(wp) :: simulated_time_1d = 0.0_wp !< updated simulated time (1d-model) 191 REAL(wp) :: sig_diss = 1.3_wp !< model constant 210 REAL(wp) :: sig_diss = 2.95_wp !< model constant according to Koblitz (2013) 211 REAL(wp) :: sig_e = 2.95_wp !< model constant according to Koblitz (2013) 192 212 REAL(wp) :: time_pr_1d = 0.0_wp !< updated simulated time for profile output (1d-model) 193 213 REAL(wp) :: time_run_control_1d = 0.0_wp !< updated simulated time for run-control output (1d-model) … … 198 218 REAL(wp) :: z01d !< roughness length for momentum (1d-model) 199 219 REAL(wp) :: z0h1d !< roughness length for scalars (1d-model) 200 201 220 202 221 REAL(wp), DIMENSION(:), ALLOCATABLE :: diss1d !< tke dissipation rate (1d-model) … … 279 298 280 299 INTEGER(iwp) :: k !< loop index 281 282 REAL(wp) :: lambda !< maximum mixing length283 300 284 301 ! … … 324 341 ! 325 342 !-- Use the same mixing length as in 3D model (LES-mode) 326 ! @todo:rename (delete?) this option327 ! As the mixing length is different between RANS and LES mode, it328 ! must be distinguished here between these modes. For RANS mode,329 ! the mixing length is calculated accoding to Blackadar, which is330 ! the other option at this point.331 ! Maybe delete this option entirely (not appropriate in LES case)332 ! 2018-03-20, gronemeier343 !> @todo rename (delete?) this option 344 !> As the mixing length is different between RANS and LES mode, it 345 !> must be distinguished here between these modes. For RANS mode, 346 !> the mixing length is calculated accoding to Blackadar, which is 347 !> the other option at this point. 348 !> Maybe delete this option entirely (not appropriate in LES case) 349 !> 2018-03-20, gronemeier 333 350 DO k = nzb+1, nzt 334 351 l1d_init(k) = ( dx * dy * dzw(k) )**0.33333333333333_wp … … 359 376 us1d = 0.1_wp ! without initial friction the flow would not change 360 377 ELSE 361 diss1d(nzb+1) = 1.0_wp378 diss1d(nzb+1) = 0.001_wp 362 379 e1d(nzb+1) = 1.0_wp 363 380 km1d(nzb+1) = 1.0_wp … … 372 389 373 390 ! 374 !-- Tendencies must be preset in order to avoid runtime errors within the375 !-- first Runge-Kutta step391 !-- Tendencies must be preset in order to avoid runtime errors 392 te_diss = 0.0_wp 376 393 te_dissm = 0.0_wp 394 te_e = 0.0_wp 377 395 te_em = 0.0_wp 378 396 te_um = 0.0_wp … … 381 399 ! 382 400 !-- Set model constant 383 IF ( dissipation_1d == 'as_in_3d_model' ) c_m = 0.1_wp 401 IF ( dissipation_1d == 'as_in_3d_model' ) c_0 = 0.1_wp 402 c_mu = c_0**4 384 403 385 404 ! … … 423 442 !-- Determine the time step at the start of a 1D-simulation and 424 443 !-- determine and printout quantities used for run control 425 dt_1d = 1.0_wp444 dt_1d = 0.01_wp 426 445 CALL run_control_1d 427 446 … … 483 502 !-- dissipation rate 484 503 IF ( dissipation_1d == 'detering' ) THEN 485 diss1d(k) = c_ m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)504 diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k) 486 505 ELSEIF ( dissipation_1d == 'as_in_3d_model' ) THEN 487 506 diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) & … … 497 516 kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1) & 498 517 - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k) & 499 ) * ddzw(k) 518 ) * ddzw(k) / sig_e & 500 519 - diss1d(k) 501 520 … … 503 522 ! 504 523 !-- dissipation rate 505 te_diss(k) = km1d(k) * & 524 IF ( rif1d(k) >= 0.0_wp ) THEN 525 alpha_buoyancy = 1.0_wp - l1d(k) / lambda 526 ELSE 527 alpha_buoyancy = 1.0_wp - ( 1.0_wp + ( c_2 - 1.0_wp ) & 528 / ( c_2 - c_1 ) ) & 529 * l1d(k) / lambda 530 ENDIF 531 c_3 = ( c_1 - c_2 ) * alpha_buoyancy 532 te_diss(k) = ( km1d(k) * & 506 533 ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2 & 507 534 + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2 & 508 ) * c_1 * c_mu**0.75 / c_h * f / us1d & 509 * SQRT(e1d(k)) & 535 ) * ( c_1 + (c_2 - c_1) * l1d(k) / lambda ) & 510 536 - g / pt_0 * kh1d(k) * flux * c_3 & 511 * diss1d(k) / ( e1d(k) + 1.0E-20_wp ) & 512 + ( kmzp * ( diss1d(k+1) - diss1d(k) ) & 537 - c_2 * diss1d(k) & 538 ) * diss1d(k) / ( e1d(k) + 1.0E-20_wp ) & 539 + ( kmzp * ( diss1d(k+1) - diss1d(k) ) & 513 540 * ddzu(k+1) & 514 541 - kmzm * ( diss1d(k) - diss1d(k-1) ) & 515 542 * ddzu(k) & 516 ) * ddzw(k) / sig_diss & 517 - c_2 * diss1d(k)**2 / ( e1d(k) + 1.0E-20_wp ) 543 ) * ddzw(k) / sig_diss 518 544 519 545 ENDIF … … 546 572 !-- dissipation rate 547 573 IF ( dissipation_1d == 'detering' ) THEN 548 diss1d(k) = c_ m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)574 diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k) 549 575 ELSEIF ( dissipation_1d == 'as_in_3d_model' ) THEN 550 576 diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) ) & … … 565 591 !-- TKE 566 592 IF ( .NOT. dissipation_1d == 'prognostic' ) THEN 593 !> @query why integrate over 2dz 594 !> Why is it allowed to integrate over two delta-z for e 595 !> while for u and v it is not? 596 !> 2018-04-23, gronemeier 567 597 te_e(k) = km1d(k) * ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2& 568 598 + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2& … … 572 602 kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1) & 573 603 - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k) & 574 ) * ddzw(k) 604 ) * ddzw(k) / sig_e & 575 605 - diss1d(k) 576 606 ENDIF … … 595 625 ENDDO 596 626 597 IF ( dissipation_1d == 'prognostic' ) THEN598 DO k = nzb_diff, nzt599 diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + &600 tsc(3) * te_dissm(k) )601 ENDDO602 ENDIF603 627 ! 604 628 !-- Eliminate negative TKE values, which can result from the … … 606 630 !-- value is reduced to 10 percent of its old value. 607 631 WHERE ( e1d_p < 0.0_wp ) e1d_p = 0.1_wp * e1d 632 633 IF ( dissipation_1d == 'prognostic' ) THEN 634 DO k = nzb+1, nzt 635 diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + & 636 tsc(3) * te_dissm(k) ) 637 ENDDO 638 WHERE ( diss1d_p < 0.0_wp ) diss1d_p = 0.1_wp * diss1d 639 ENDIF 608 640 ENDIF 609 641 … … 643 675 IF ( dissipation_1d == 'prognostic' ) THEN 644 676 DO k = nzb+1, nzt 645 te_dissm(k) = -9.5625_wp * te_diss(k) + 5.3125_wp * te_dissm(k) 677 te_dissm(k) = -9.5625_wp * te_diss(k) & 678 + 5.3125_wp * te_dissm(k) 646 679 ENDDO 647 680 ENDIF … … 650 683 ENDIF 651 684 ENDIF 652 653 685 654 686 ! … … 701 733 702 734 ENDIF ! constant_flux_layer 703 735 !> @todo combine if clauses 736 !> The previous and following if clauses can be combined into a 737 !> single clause 738 !> 2018-04-23, gronemeier 704 739 ! 705 740 !-- Compute the Richardson-flux numbers, … … 789 824 !-- compatibility with the 3D model. 790 825 IF ( ibc_e_b == 2 ) THEN 791 e1d(nzb+1) = ( us1d / c_ m)**2826 e1d(nzb+1) = ( us1d / c_0 )**2 792 827 ENDIF 793 828 IF ( dissipation_1d == 'prognostic' ) THEN 794 e1d(nzb+1) = us1d**2 / SQRT( c_mu )829 e1d(nzb+1) = ( us1d / c_0 )**2 795 830 diss1d(nzb+1) = us1d**3 / ( kappa * zu(nzb+1) ) 796 831 diss1d(nzb) = diss1d(nzb+1) … … 829 864 !-- in the dissipation of TKE via l1d_diss. Otherwise, km1d would be 830 865 !-- too large. 831 IF ( mixing_length_1d == 'blackadar' ) THEN 866 IF ( dissipation_1d /= 'prognostic' ) THEN 867 IF ( mixing_length_1d == 'blackadar' ) THEN 868 DO k = nzb+1, nzt 869 IF ( rif1d(k) >= 0.0_wp ) THEN 870 l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) ) 871 l1d_diss(k) = l1d(k) 872 ELSE 873 l1d(k) = l1d_init(k) 874 l1d_diss(k) = l1d_init(k) * & 875 SQRT( 1.0_wp - 16.0_wp * rif1d(k) ) 876 ENDIF 877 ENDDO 878 ELSEIF ( mixing_length_1d == 'as_in_3d_model' ) THEN 879 DO k = nzb+1, nzt 880 dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k) 881 IF ( dpt_dz > 0.0_wp ) THEN 882 l_stable = 0.76_wp * SQRT( e1d(k) ) & 883 / SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp 884 ELSE 885 l_stable = l1d_init(k) 886 ENDIF 887 l1d(k) = MIN( l1d_init(k), l_stable ) 888 l1d_diss(k) = l1d(k) 889 ENDDO 890 ENDIF 891 ELSE 832 892 DO k = nzb+1, nzt 833 IF ( rif1d(k) >= 0.0_wp ) THEN 834 l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) ) 835 l1d_diss(k) = l1d(k) 836 ELSE 837 l1d(k) = l1d_init(k) 838 l1d_diss(k) = l1d_init(k) * & 839 SQRT( 1.0_wp - 16.0_wp * rif1d(k) ) 840 ENDIF 841 ENDDO 842 ELSEIF ( mixing_length_1d == 'as_in_3d_model' ) THEN 843 DO k = nzb+1, nzt 844 dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k) 845 IF ( dpt_dz > 0.0_wp ) THEN 846 l_stable = 0.76_wp * SQRT( e1d(k) ) / & 847 SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp 848 ELSE 849 l_stable = l1d_init(k) 850 ENDIF 851 l1d(k) = MIN( l1d_init(k), l_stable ) 852 l1d_diss(k) = l1d(k) 893 l1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) & 894 / ( diss1d(k) + 1.0E-30_wp ) 853 895 ENDDO 854 896 ENDIF … … 867 909 ENDIF 868 910 ENDIF 911 869 912 IF ( dissipation_1d == 'prognostic' ) THEN 870 913 DO k = nzb_diff, nzt … … 873 916 ELSE 874 917 DO k = nzb_diff, nzt 875 km1d(k) = c_ m* SQRT( e1d(k) ) * l1d(k)918 km1d(k) = c_0 * SQRT( e1d(k) ) * l1d(k) 876 919 ENDDO 877 920 ENDIF … … 1018 1061 1019 1062 REAL(wp) :: dt_diff !< time step accorind to diffusion criterion 1063 REAL(wp) :: dt_old !< previous time step 1020 1064 REAL(wp) :: fac !< factor of criterion 1021 1065 REAL(wp) :: value !< auxiliary variable 1022 1066 1023 1067 ! 1068 !-- Save previous time step 1069 dt_old = dt_1d 1070 1071 ! 1024 1072 !-- Compute the currently feasible time step according to the diffusion 1025 1073 !-- criterion. At nzb+1 the half grid length is used. 1026 fac = 0.125 !0.35_wp !### changed from 0.351074 fac = 0.125 1027 1075 dt_diff = dt_max_1d 1028 1076 DO k = nzb+2, nzt … … 1034 1082 1035 1083 ! 1084 !-- Limit the new time step to a maximum of 10 times the previous time step 1085 dt_1d = MIN( dt_old * 10.0_wp, dt_1d ) 1086 1087 ! 1036 1088 !-- Set flag when the time step becomes too small 1037 IF ( dt_1d < ( 0.00001_wp * dt_max_1d ) ) THEN1089 IF ( dt_1d < ( 1.0E-15_wp * dt_max_1d ) ) THEN 1038 1090 stop_dt_1d = .TRUE. 1039 1091 -
palm/trunk/SOURCE/modules.f90
r3065 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! set dt_3d = 0.01 28 ! 29 ! 3065 2018-06-12 07:03:02Z Giersch 27 30 ! Variables concerning stretching introduced or revised 28 31 ! … … 1422 1425 REAL(wp) :: dt_run_control = 60.0_wp !< namelist parameter 1423 1426 REAL(wp) :: dt_spinup = 60.0_wp !< namelist parameter 1424 REAL(wp) :: dt_3d = 1.0_wp!< time step1427 REAL(wp) :: dt_3d = 0.01_wp !< time step 1425 1428 REAL(wp) :: dz_max = 1000.0_wp !< namelist parameter 1426 1429 REAL(wp) :: dz_stretch_factor = 1.08_wp !< namelist parameter -
palm/trunk/SOURCE/parin.f90
r3065 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Added rans_const_c and rans_const_sigma as input parameters (TG) 28 ! 29 ! 3065 2018-06-12 07:03:02Z Giersch 27 30 ! New initialization parameters added 28 31 ! … … 467 470 USE synthetic_turbulence_generator_mod, & 468 471 ONLY: stg_parin 472 473 USE turbulence_closure_mod, & 474 ONLY: rans_const_c, rans_const_sigma 469 475 470 476 USE urban_surface_mod, & … … 531 537 pt_vertical_gradient_level, q_surface, q_surface_initial_change, & 532 538 q_vertical_gradient, q_vertical_gradient_level, & 533 random_generator, random_heatflux, rans_mode, & 539 random_generator, random_heatflux, rans_const_c, rans_const_sigma,& 540 rans_mode, & 534 541 rayleigh_damping_factor, rayleigh_damping_height, & 535 542 recycling_width, recycling_yshift, & … … 601 608 pt_vertical_gradient_level, q_surface, q_surface_initial_change, & 602 609 q_vertical_gradient, q_vertical_gradient_level, & 603 random_generator, random_heatflux, rans_mode, & 610 random_generator, random_heatflux, rans_const_c, rans_const_sigma,& 611 rans_mode, & 604 612 rayleigh_damping_factor, rayleigh_damping_height, & 605 613 recycling_width, recycling_yshift, & -
palm/trunk/SOURCE/pmc_interface_mod.f90
- Property svn:mergeinfo deleted
-
palm/trunk/SOURCE/timestep.f90
r3049 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! limit dt_3d to be at maximum 2*old_dt; define old_dt at beginning of routine 28 ! Add km/kh_max 29 ! 30 ! 3049 2018-05-29 13:52:36Z Giersch 27 31 ! Error messages revised 28 32 ! … … 161 165 INTEGER(iwp) :: j !< 162 166 INTEGER(iwp) :: k !< 167 INTEGER(iwp) :: km_max_ijk(3) = -1 !< index values (i,j,k) of location where km_max occurs 168 INTEGER(iwp) :: kh_max_ijk(3) = -1 !< index values (i,j,k) of location where kh_max occurs 163 169 164 170 LOGICAL :: stop_dt_local !< local switch for controlling the time stepping … … 173 179 REAL(wp) :: dt_w !< 174 180 REAL(wp) :: dt_w_l !< 181 REAL(wp) :: km_max !< maximum of Km in entire domain 182 REAL(wp) :: kh_max !< maximum of Kh in entire domain 175 183 REAL(wp) :: u_gtrans_l !< 176 184 REAL(wp) :: u_max_l !< … … 190 198 191 199 192 193 200 CALL cpu_log( log_point(12), 'calculate_timestep', 'start' ) 201 ! 202 !-- Save former time step as reference 203 old_dt = dt_3d 194 204 195 205 ! … … 318 328 !-- The time step is the minimum of the 3-4 components and the diffusion time 319 329 !-- step minus a reduction (cfl_factor) to be on the safe side. 320 !-- The time step must not exceed the maximum allowed value. 330 !-- The time step must not exceed the maximum allowed value and must not 331 !-- increase by more than a factor of 10. 321 332 dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w, dt_precipitation ) 322 dt_3d = MIN( dt_3d, dt_max )333 dt_3d = MIN( dt_3d, dt_max, 2.0_wp * old_dt ) 323 334 324 335 ! … … 334 345 IF ( dt_3d < ( 0.00001_wp * dt_max ) ) THEN 335 346 stop_dt = .TRUE. 347 348 ! 349 !-- Determine the maxima of the diffusion coefficients, including their 350 !-- grid index positions. 351 CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, km, 'abs', & 352 0.0_wp, km_max, km_max_ijk ) 353 CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, kh, 'abs', & 354 0.0_wp, kh_max, kh_max_ijk ) 336 355 337 356 WRITE( message_string, * ) 'Time step has reached minimum limit.', & … … 342 361 '&dt_w = ', dt_w, ' s', & 343 362 '&dt_diff = ', dt_diff, ' s', & 344 '&u_max = ', u_max, ' m/s k=', u_max_ijk(1),&363 '&u_max = ', u_max, ' m/s k=', u_max_ijk(1), & 345 364 ' j=', u_max_ijk(2), ' i=', u_max_ijk(3), & 346 '&v_max = ', v_max, ' m/s k=', v_max_ijk(1),&365 '&v_max = ', v_max, ' m/s k=', v_max_ijk(1), & 347 366 ' j=', v_max_ijk(2), ' i=', v_max_ijk(3), & 348 '&w_max = ', w_max, ' m/s k=', w_max_ijk(1), & 349 ' j=', w_max_ijk(2), ' i=', w_max_ijk(3) 367 '&w_max = ', w_max, ' m/s k=', w_max_ijk(1), & 368 ' j=', w_max_ijk(2), ' i=', w_max_ijk(3), & 369 '&km_max = ', km_max, ' m2/s2 k=', km_max_ijk(1), & 370 ' j=', km_max_ijk(2), ' i=', km_max_ijk(3), & 371 '&kh_max = ', kh_max, ' m2/s2 k=', kh_max_ijk(1), & 372 ' j=', kh_max_ijk(2), ' i=', kh_max_ijk(3) 350 373 CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 ) 351 374 ! … … 388 411 dt_3d = NINT( dt_3d * 100.0_wp / div ) * div / 100.0_wp 389 412 390 !391 !-- Adjust the time step392 old_dt = dt_3d393 394 413 ENDIF 395 414 -
palm/trunk/SOURCE/turbulence_closure_mod.f90
r3045 r3083 25 25 ! ----------------- 26 26 ! $Id$ 27 ! - set limits of diss at the end of prognostic equations 28 ! - call production_e to calculate production term of diss 29 ! - limit change of diss to -90% to +100% 30 ! - remove factor 0.5 from diffusion_diss_ij 31 ! - rename c_m into c_0, and c_h into c_4 32 ! - add rans_const_c and rans_const_sigma as namelist parameters 33 ! - add calculation of mixing length for profile output in case of rans_tke_e 34 ! - changed format of annotations to comply with doxygen standards 35 ! - calculate and save dissipation rate during rans_tke_l mode 36 ! - set bc at vertical walls for e, diss, km, kh 37 ! - bugfix: set l_wall = 0.0 within buildings 38 ! - set l_wall at bottom and top boundary (rans-mode) 39 ! - bugfix in production term for dissipation rate 40 ! - bugfix in diffusion of dissipation rate 41 ! - disable check for 1D model if rans_tke_e is used 42 ! - bugfixes for initialization (rans-mode): 43 ! - correction of dissipation-rate formula 44 ! - calculate km based on l_wall 45 ! - initialize diss if 1D model is not used 46 ! 47 ! 3045 2018-05-28 07:55:41Z Giersch 27 48 ! Error message revised 28 49 ! … … 36 57 ! Further todo's 37 58 ! 38 ! 2936 2018-03-27 14:49:27Z suehring59 ! 2936 2018-03-27 14:49:27Z gronemeier 39 60 ! - defined l_grid only within this module 40 61 ! - Moved l_wall definition from modules.f90 … … 81 102 !> add OpenMP directives whereever possible 82 103 !> remove debug output variables (dummy1, dummy2, dummy3) 83 !> @todo Move initialization of wall-mixing length from init_grid84 104 !> @todo Check for random disturbances 85 105 !> @note <Enter notes on the module> … … 146 166 147 167 148 REAL(wp) :: c_1 = 1.44_wp !< model constant for RANS mode 149 REAL(wp) :: c_2 = 1.92_wp !< model constant for RANS mode 150 REAL(wp) :: c_3 = 1.44_wp !< model constant for RANS mode 151 REAL(wp) :: c_h = 0.0015_wp !< model constant for RANS mode 152 REAL(wp) :: c_m !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES) 153 REAL(wp) :: c_mu = 0.09_wp !< model constant for RANS mode 154 REAL(wp) :: l_max !< maximum length scale for Blackadar mixing length 155 REAL(wp) :: sig_e = 1.0_wp !< factor to calculate Ke from Km 156 REAL(wp) :: sig_diss = 1.3_wp !< factor to calculate K_diss from Km 157 INTEGER(iwp) :: surf_e !< end index of surface elements at given i-j position 158 INTEGER(iwp) :: surf_s !< start index of surface elements at given i-j position 168 REAL(wp) :: c_0 !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES) 169 REAL(wp) :: c_1 !< model constant for RANS mode 170 REAL(wp) :: c_2 !< model constant for RANS mode 171 REAL(wp) :: c_3 !< model constant for RANS mode 172 REAL(wp) :: c_4 !< model constant for RANS mode 173 REAL(wp) :: l_max !< maximum length scale for Blackadar mixing length 174 REAL(wp) :: dsig_e = 1.0_wp !< factor to calculate Ke from Km (1/sigma_e) 175 REAL(wp) :: dsig_diss = 1.0_wp !< factor to calculate K_diss from Km (1/sigma_diss) 176 INTEGER(iwp) :: surf_e !< end index of surface elements at given i-j position 177 INTEGER(iwp) :: surf_s !< start index of surface elements at given i-j position 178 179 REAL(wp), DIMENSION(0:4) :: rans_const_c = & !< model constants for RANS mode (namelist param) 180 (/ 0.55_wp, 1.44_wp, 1.92_wp, 0.0_wp, 0.0_wp /) !> default values fit for standard-tke-e closure 181 182 REAL(wp), DIMENSION(2) :: rans_const_sigma = & !< model constants for RANS mode, sigma values (sigma_e, sigma_diss) (namelist param) 183 (/ 1.0_wp, 0.77_wp /) 159 184 160 185 REAL(wp), DIMENSION(:), ALLOCATABLE :: l_black !< mixing length according to Blackadar … … 163 188 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: l_wall !< near-wall mixing length 164 189 165 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy1 !< debug output variable 166 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy2 !< debug output variable 167 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy3 !< debug output variable 168 169 170 PUBLIC c_m, dummy1, dummy2, dummy3 190 !> @todo remove debug variables 191 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_prod1, diss_adve1, diss_diff1, & 192 diss_prod2, diss_adve2, diss_diff2, & 193 diss_prod3, diss_adve3, diss_diff3, & 194 dummy1, dummy2, dummy3 195 196 197 PUBLIC c_0, rans_const_c, rans_const_sigma 171 198 172 199 ! … … 294 321 ! ------------ 295 322 !> Check parameters routine for turbulence closure module. 323 !> @todo remove rans_mode from initialization namelist and rework checks 324 !> The way it is implemented at the moment, the user has to set two variables 325 !> so that the RANS mode is working. It would be better if only one parameter 326 !> has to be set. 327 !> 2018-06-18, gronemeier 296 328 !------------------------------------------------------------------------------! 297 329 SUBROUTINE tcm_check_parameters … … 307 339 IF ( rans_mode ) THEN 308 340 309 c_m = 0.4_wp !according to Detering and Etling (1985) 341 ! 342 !-- Assign values to constants for RANS mode 343 dsig_e = 1.0_wp / rans_const_sigma(1) 344 dsig_diss = 1.0_wp / rans_const_sigma(2) 345 346 c_0 = rans_const_c(0) 347 c_1 = rans_const_c(1) 348 c_2 = rans_const_c(2) 349 !c_3 = rans_const_c(3) !> @todo clarify how to switch between different models 350 c_4 = rans_const_c(4) 310 351 311 352 SELECT CASE ( TRIM( turbulence_closure ) ) … … 316 357 CASE ( 'TKE-e' ) 317 358 rans_tke_e = .TRUE. 318 319 IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) == 0 &320 .AND. .NOT. nest_domain ) THEN321 message_string = 'Initializing without 1D model while ' // &322 'using TKE-e closure&' // &323 'is not possible at the moment!'324 CALL message( 'tcm_check_parameters', 'TG0005', 1, 2, 0, 6, 0 )325 ENDIF326 359 327 360 CASE DEFAULT 328 361 message_string = 'Unknown turbulence closure: ' // & 329 362 TRIM( turbulence_closure ) 330 CALL message( 'tcm_check_parameters', ' TG0001', 1, 2, 0, 6, 0 )363 CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 ) 331 364 332 365 END SELECT 366 367 IF ( turbulent_inflow .OR. turbulent_outflow ) THEN 368 message_string = 'turbulent inflow/outflow is not yet '// & 369 'implemented for RANS mode' 370 CALL message( 'tcm_check_parameters', 'PA0501', 1, 2, 0, 6, 0 ) 371 ENDIF 333 372 334 373 message_string = 'RANS mode is still in development! ' // & 335 374 '&Not all features of PALM are yet compatible '// & 336 375 'with RANS mode. &Use at own risk!' 337 CALL message( 'tcm_check_parameters', ' TG0003', 0, 1, 0, 6, 0 )376 CALL message( 'tcm_check_parameters', 'PA0502', 0, 1, 0, 6, 0 ) 338 377 339 378 ELSE 340 379 341 c_m = 0.1_wp !according to Lilly (1967) and Deardorff (1980) 380 c_0 = 0.1_wp !according to Lilly (1967) and Deardorff (1980) 381 382 dsig_e = 1.0_wp !assure to use K_m to calculate TKE instead 383 !of K_e which is used in RANS mode 342 384 343 385 SELECT CASE ( TRIM( turbulence_closure ) ) … … 347 389 348 390 CASE DEFAULT 391 !> @todo rework this part so that only one call of this error exists 349 392 message_string = 'Unknown turbulence closure: ' // & 350 393 TRIM( turbulence_closure ) 351 CALL message( 'tcm_check_parameters', ' TG0001', 1, 2, 0, 6, 0 )394 CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 ) 352 395 353 396 END SELECT 354 355 ENDIF356 357 IF ( rans_tke_e ) THEN358 359 IF ( turbulent_inflow .OR. turbulent_outflow ) THEN360 message_string = 'turbulent inflow/outflow is not yet '// &361 'implemented for TKE-e closure'362 CALL message( 'tcm_check_parameters', 'TG0002', 1, 2, 0, 6, 0 )363 ENDIF364 397 365 398 ENDIF … … 379 412 IMPLICIT NONE 380 413 381 CHARACTER (LEN=*) :: unit !< 382 CHARACTER (LEN=*) :: var !< 383 384 INTEGER(iwp) :: i !< 385 INTEGER(iwp) :: ilen !< 386 INTEGER(iwp) :: k !< 414 CHARACTER (LEN=*) :: unit !< unit of output variable 415 CHARACTER (LEN=*) :: var !< name of output variable 416 417 INTEGER(iwp) :: i !< index of var in data_output 418 INTEGER(iwp) :: ilen !< length of var string 419 INTEGER(iwp) :: k !< flag if var contains one of '_xy', '_xz' or '_yz' 387 420 388 421 SELECT CASE ( TRIM( var ) ) 389 422 390 423 CASE ( 'diss' ) 391 IF ( .NOT. rans_tke_e ) THEN392 message_string = 'output of "' // TRIM( var ) // '" requi' // &393 'res TKE-e closure for RANS mode.'394 CALL message( 'tcm_check_data_output', 'TG0101', 1, 2, 0, 6, 0 )395 ENDIF396 424 unit = 'm2/s3' 397 425 398 CASE ( 'dummy2', 'dummy3', 'dummy1' ) 399 unit = 'mixing length' 426 CASE ( 'diss1', 'diss2', & !> @todo remove later 427 'diss_prod1', 'diss_adve1', 'diss_diff1', & 428 'diss_prod2', 'diss_adve2', 'diss_diff2', & 429 'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3' ) 430 unit = 'debug output' 400 431 401 432 CASE ( 'kh', 'km' ) … … 420 451 IMPLICIT NONE 421 452 422 CHARACTER (LEN=*), INTENT(OUT) :: grid_x !< 423 CHARACTER (LEN=*), INTENT(OUT) :: grid_y !< 424 CHARACTER (LEN=*), INTENT(OUT) :: grid_z !< 425 CHARACTER (LEN=*), INTENT(IN) :: var !< 426 427 LOGICAL, INTENT(OUT) :: found !< 428 453 CHARACTER (LEN=*), INTENT(OUT) :: grid_x !< x grid of output variable 454 CHARACTER (LEN=*), INTENT(OUT) :: grid_y !< y grid of output variable 455 CHARACTER (LEN=*), INTENT(OUT) :: grid_z !< z grid of output variable 456 CHARACTER (LEN=*), INTENT(IN) :: var !< name of output variable 457 458 LOGICAL, INTENT(OUT) :: found !< flag if output variable is found 459 429 460 found = .TRUE. 430 461 … … 438 469 grid_z = 'zu' 439 470 440 CASE ( 'dummy2', 'dummy3', 'dummy1' ) !### remove later 471 CASE ( 'diss1', 'diss2', & !> @todo remove later 472 'diss_prod1', 'diss_adve1', 'diss_diff1', & 473 'diss_prod2', 'diss_adve2', 'diss_diff2', & 474 'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3' ) 441 475 grid_x = 'x' 442 476 grid_y = 'y' … … 480 514 IMPLICIT NONE 481 515 482 CHARACTER (LEN=*) :: mode !< 483 CHARACTER (LEN=*) :: variable !< 484 485 INTEGER(iwp) :: i !< 486 INTEGER(iwp) :: j !< 487 INTEGER(iwp) :: k !< 516 CHARACTER (LEN=*) :: mode !< flag defining mode 'allocate', 'sum' or 'average' 517 CHARACTER (LEN=*) :: variable !< name of variable 518 519 INTEGER(iwp) :: i !< loop index 520 INTEGER(iwp) :: j !< loop index 521 INTEGER(iwp) :: k !< loop index 488 522 489 523 IF ( mode == 'allocate' ) THEN … … 616 650 IMPLICIT NONE 617 651 618 CHARACTER (LEN=*) :: grid !< 619 CHARACTER (LEN=*) :: mode !< 620 CHARACTER (LEN=*) :: variable !< 621 622 INTEGER(iwp) :: av !< 623 INTEGER(iwp) :: i !< 624 INTEGER(iwp) :: j !< 625 INTEGER(iwp) :: k !< 626 INTEGER(iwp) :: nzb_do !< 627 INTEGER(iwp) :: nzt_do !< 628 629 LOGICAL :: found !< 652 CHARACTER (LEN=*) :: grid !< name of vertical grid 653 CHARACTER (LEN=*) :: mode !< either 'xy', 'xz' or 'yz' 654 CHARACTER (LEN=*) :: variable !< name of variable 655 656 INTEGER(iwp) :: av !< flag for (non-)average output 657 INTEGER(iwp) :: i !< loop index 658 INTEGER(iwp) :: j !< loop index 659 INTEGER(iwp) :: k !< loop index 660 INTEGER(iwp) :: nzb_do !< vertical output index (bottom) 661 INTEGER(iwp) :: nzt_do !< vertical output index (top) 662 663 LOGICAL :: found !< flag if output variable is found 630 664 LOGICAL :: two_d !< flag parameter that indicates 2D variables (horizontal cross sections) 631 665 … … 737 771 IMPLICIT NONE 738 772 739 CHARACTER (LEN=*) :: variable !< 740 741 INTEGER(iwp) :: av !< 742 INTEGER(iwp) :: i !< 743 INTEGER(iwp) :: j !< 744 INTEGER(iwp) :: k !< 773 CHARACTER (LEN=*) :: variable !< name of variable 774 775 INTEGER(iwp) :: av !< flag for (non-)average output 776 INTEGER(iwp) :: i !< loop index 777 INTEGER(iwp) :: j !< loop index 778 INTEGER(iwp) :: k !< loop index 745 779 INTEGER(iwp) :: nzb_do !< lower limit of the data output (usually 0) 746 780 INTEGER(iwp) :: nzt_do !< vertical upper limit of the data output (usually nz_do3d) 747 781 748 LOGICAL :: found !< 782 LOGICAL :: found !< flag if output variable is found 749 783 750 784 REAL(wp) :: fill_value = -999.0_wp !< value for the _FillValue attribute … … 829 863 ENDIF 830 864 831 CASE ( 'dummy1' ) !### remove later 865 CASE ( 'dummy3' ) !> @todo remove later 866 IF ( av == 0 ) THEN 867 DO i = nxl, nxr 868 DO j = nys, nyn 869 DO k = nzb_do, nzt_do 870 local_pf(i,j,k) = dummy3(k,j,i) 871 ENDDO 872 ENDDO 873 ENDDO 874 ENDIF 875 876 CASE ( 'diss1' ) !> @todo remove later 832 877 IF ( av == 0 ) THEN 833 878 DO i = nxl, nxr … … 840 885 ENDIF 841 886 842 CASE ( 'd ummy2' ) !###remove later887 CASE ( 'diss2' ) !> @todo remove later 843 888 IF ( av == 0 ) THEN 844 889 DO i = nxl, nxr … … 851 896 ENDIF 852 897 853 CASE ( 'd ummy3' ) !###remove later898 CASE ( 'diss_prod1' ) !> @todo remove later 854 899 IF ( av == 0 ) THEN 855 900 DO i = nxl, nxr 856 901 DO j = nys, nyn 857 902 DO k = nzb_do, nzt_do 858 local_pf(i,j,k) = d ummy3(k,j,i)903 local_pf(i,j,k) = diss_prod1(k,j,i) 859 904 ENDDO 860 905 ENDDO … … 862 907 ENDIF 863 908 909 CASE ( 'diss_adve1' ) !> @todo remove later 910 IF ( av == 0 ) THEN 911 DO i = nxl, nxr 912 DO j = nys, nyn 913 DO k = nzb_do, nzt_do 914 local_pf(i,j,k) = diss_adve1(k,j,i) 915 ENDDO 916 ENDDO 917 ENDDO 918 ENDIF 919 920 CASE ( 'diss_diff1' ) !> @todo remove later 921 IF ( av == 0 ) THEN 922 DO i = nxl, nxr 923 DO j = nys, nyn 924 DO k = nzb_do, nzt_do 925 local_pf(i,j,k) = diss_diff1(k,j,i) 926 ENDDO 927 ENDDO 928 ENDDO 929 ENDIF 930 931 CASE ( 'diss_prod2' ) !> @todo remove later 932 IF ( av == 0 ) THEN 933 DO i = nxl, nxr 934 DO j = nys, nyn 935 DO k = nzb_do, nzt_do 936 local_pf(i,j,k) = diss_prod2(k,j,i) 937 ENDDO 938 ENDDO 939 ENDDO 940 ENDIF 941 942 CASE ( 'diss_adve2' ) !> @todo remove later 943 IF ( av == 0 ) THEN 944 DO i = nxl, nxr 945 DO j = nys, nyn 946 DO k = nzb_do, nzt_do 947 local_pf(i,j,k) = diss_adve2(k,j,i) 948 ENDDO 949 ENDDO 950 ENDDO 951 ENDIF 952 953 CASE ( 'diss_diff2' ) !> @todo remove later 954 IF ( av == 0 ) THEN 955 DO i = nxl, nxr 956 DO j = nys, nyn 957 DO k = nzb_do, nzt_do 958 local_pf(i,j,k) = diss_diff2(k,j,i) 959 ENDDO 960 ENDDO 961 ENDDO 962 ENDIF 963 964 CASE ( 'diss_prod3' ) !> @todo remove later 965 IF ( av == 0 ) THEN 966 DO i = nxl, nxr 967 DO j = nys, nyn 968 DO k = nzb_do, nzt_do 969 local_pf(i,j,k) = diss_prod3(k,j,i) 970 ENDDO 971 ENDDO 972 ENDDO 973 ENDIF 974 975 CASE ( 'diss_adve3' ) !> @todo remove later 976 IF ( av == 0 ) THEN 977 DO i = nxl, nxr 978 DO j = nys, nyn 979 DO k = nzb_do, nzt_do 980 local_pf(i,j,k) = diss_adve3(k,j,i) 981 ENDDO 982 ENDDO 983 ENDDO 984 ENDIF 985 986 CASE ( 'diss_diff3' ) !> @todo remove later 987 IF ( av == 0 ) THEN 988 DO i = nxl, nxr 989 DO j = nys, nyn 990 DO k = nzb_do, nzt_do 991 local_pf(i,j,k) = diss_diff3(k,j,i) 992 ENDDO 993 ENDDO 994 ENDDO 995 ENDIF 996 864 997 CASE DEFAULT 865 998 found = .FALSE. … … 891 1024 ALLOCATE( km(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 892 1025 893 ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) ! ###remove later1026 ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) !> @todo remove later 894 1027 ALLOCATE( dummy2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 895 1028 ALLOCATE( dummy3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 896 897 IF ( rans_mode ) ALLOCATE( l_black(nzb:nzt+1) ) 1029 ALLOCATE( diss_adve1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1030 ALLOCATE( diss_adve2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1031 ALLOCATE( diss_adve3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1032 ALLOCATE( diss_prod1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1033 ALLOCATE( diss_prod2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1034 ALLOCATE( diss_prod3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1035 ALLOCATE( diss_diff1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1036 ALLOCATE( diss_diff2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1037 ALLOCATE( diss_diff3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1038 dummy1 = 0.0_wp 1039 dummy2 = 0.0_wp 1040 dummy3 = 0.0_wp 1041 diss_adve1 = 0.0_wp 1042 diss_adve2 = 0.0_wp 1043 diss_adve3 = 0.0_wp 1044 diss_prod1 = 0.0_wp 1045 diss_prod2 = 0.0_wp 1046 diss_prod3 = 0.0_wp 1047 diss_diff1 = 0.0_wp 1048 diss_diff2 = 0.0_wp 1049 diss_diff3 = 0.0_wp 898 1050 899 1051 #if defined( __nopointer ) … … 912 1064 !-- they do not necessarily need to be transferred, which is attributed to 913 1065 !-- the design of the model coupler which allocates memory for each variable. 914 IF ( rans_ tke_e .OR. use_sgs_for_particles .OR. wang_kernel .OR.&1066 IF ( rans_mode .OR. use_sgs_for_particles .OR. wang_kernel .OR. & 915 1067 collision_turbulence .OR. nested_run ) THEN 916 1068 #if defined( __nopointer ) … … 934 1086 e => e_1; e_p => e_2; te_m => e_3 935 1087 936 IF ( rans_ tke_e .OR. use_sgs_for_particles .OR. &1088 IF ( rans_mode .OR. use_sgs_for_particles .OR. & 937 1089 wang_kernel .OR. collision_turbulence .OR. nested_run ) THEN 938 1090 diss => diss_1 … … 967 1119 INTEGER(iwp) :: j !< loop index 968 1120 INTEGER(iwp) :: k !< loop index 969 INTEGER(iwp) :: nz_s_shift !< 970 INTEGER(iwp) :: nz_s_shift_l !< 1121 INTEGER(iwp) :: nz_s_shift !< lower shift index for scalars 1122 INTEGER(iwp) :: nz_s_shift_l !< local lower shift index in case of turbulent inflow 971 1123 972 1124 ! 973 1125 !-- Initialize mixing length 974 1126 CALL tcm_init_mixing_length 1127 dummy3 = l_wall !> @todo remove later 975 1128 976 1129 ! … … 995 1148 996 1149 IF ( rans_tke_e ) THEN 997 IF ( dissipation_1d == 'prognostic' ) THEN ! ###Why must this be checked?998 DO i = nxlg, nxrg ! ###Should 'diss' not always999 DO j = nysg, nyng ! ###be prognostic in case rans_tke_e?1150 IF ( dissipation_1d == 'prognostic' ) THEN !> @query Why must this be checked? 1151 DO i = nxlg, nxrg !> Should 'diss' not always 1152 DO j = nysg, nyng !> be prognostic in case rans_tke_e? 1000 1153 diss(:,j,i) = diss1d 1001 1154 ENDDO … … 1005 1158 DO j = nysg, nyng 1006 1159 DO k = nzb+1, nzt 1007 diss(k,j,i) = e(k,j,i) * SQRT( e(k,j,i) ) / l1d(k)1160 diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km1d(k) 1008 1161 ENDDO 1009 1162 ENDDO … … 1016 1169 1017 1170 IF ( constant_diffusion ) THEN 1018 km 1019 kh 1020 e 1171 km = km_constant 1172 kh = km / prandtl_number 1173 e = 0.0_wp 1021 1174 ELSEIF ( e_init > 0.0_wp ) THEN 1022 DO k = nzb+1, nzt 1023 km(k,:,:) = c_m * l_grid(k) * SQRT( e_init ) 1175 DO i = nxlg, nxrg 1176 DO j = nysg, nyng 1177 DO k = nzb+1, nzt 1178 km(k,j,i) = c_0 * l_wall(k,j,i) * SQRT( e_init ) 1179 ENDDO 1180 ENDDO 1024 1181 ENDDO 1025 1182 km(nzb,:,:) = km(nzb+1,:,:) 1026 1183 km(nzt+1,:,:) = km(nzt,:,:) 1027 kh 1028 e 1184 kh = km / prandtl_number 1185 e = e_init 1029 1186 ELSE 1030 1187 IF ( .NOT. ocean ) THEN … … 1040 1197 ENDIF 1041 1198 1199 IF ( rans_tke_e ) THEN 1200 DO i = nxlg, nxrg 1201 DO j = nysg, nyng 1202 DO k = nzb+1, nzt 1203 diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km(k,j,i) 1204 ENDDO 1205 ENDDO 1206 ENDDO 1207 diss(nzb,:,:) = diss(nzb+1,:,:) 1208 diss(nzt+1,:,:) = diss(nzt,:,:) 1209 ENDIF 1210 1042 1211 ENDIF 1043 1212 ! … … 1073 1242 ENDDO 1074 1243 ENDDO 1244 IF ( rans_tke_e ) THEN 1245 DO i = nxlg, nxrg 1246 DO j = nysg, nyng 1247 nz_s_shift = get_topography_top_index_ji( j, i, 's' ) 1248 1249 diss(nz_s_shift:nzt+1,j,i) = diss(0:nzt+1-nz_s_shift,j,i) 1250 ENDDO 1251 ENDDO 1252 ENDIF 1075 1253 ENDIF 1076 1254 … … 1084 1262 !-- boundary and adjust mean inflow profiles 1085 1263 IF ( complex_terrain ) THEN 1086 IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 ) THEN 1264 IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. & 1265 nysg <= 0 .AND. nyng >= 0 ) THEN 1087 1266 nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' ) 1088 1267 ELSE … … 1095 1274 nz_s_shift = nz_s_shift_l 1096 1275 #endif 1097 mean_inflow_profiles(nz_s_shift:nzt+1,5) = hom_sum(0:nzt+1-nz_s_shift,8,0) ! e 1276 mean_inflow_profiles(nz_s_shift:nzt+1,5) = & 1277 hom_sum(0:nzt+1-nz_s_shift,8,0) ! e 1098 1278 ENDIF 1099 1279 ! … … 1114 1294 ! 1115 1295 !-- Inside buildings set TKE back to zero. 1116 !-- Other scalars (km, kh, diss,...) are ignored at present,1296 !-- Other scalars (km, kh,...) are ignored at present, 1117 1297 !-- maybe revise later. 1118 1298 DO i = nxlg, nxrg … … 1121 1301 e(k,j,i) = MERGE( e(k,j,i), 0.0_wp, & 1122 1302 BTEST( wall_flags_0(k,j,i), 0 ) ) 1123 te_m(k,j,i) = MERGE( te_m(k,j,i), 0.0_wp, &1124 BTEST( wall_flags_0(k,j,i), 0 ) )1125 1303 ENDDO 1126 1304 ENDDO 1127 1305 ENDDO 1128 1306 1307 IF ( rans_tke_e ) THEN 1308 DO i = nxlg, nxrg 1309 DO j = nysg, nyng 1310 DO k = nzb, nzt 1311 diss(k,j,i) = MERGE( diss(k,j,i), 0.0_wp, & 1312 BTEST( wall_flags_0(k,j,i), 0 ) ) 1313 ENDDO 1314 ENDDO 1315 ENDDO 1316 ENDIF 1129 1317 ENDIF 1130 1318 ! … … 1134 1322 ! 1135 1323 !-- Allthough tendency arrays are set in prognostic_equations, they have 1136 !-- have to be predefined here because they are used (but multiplied with 0)1137 !-- there before they are set.1324 !-- to be predefined here because there they are used (but multiplied with 0) 1325 !-- before they are set. 1138 1326 te_m = 0.0_wp 1327 1328 IF ( rans_tke_e ) THEN 1329 diss_p = diss 1330 tdiss_m = 0.0_wp 1331 ENDIF 1139 1332 1140 1333 ENDIF … … 1215 1408 1216 1409 DO k = 1, nzt 1217 IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR. &1410 IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR. & 1218 1411 l_grid(k) > 1.5_wp * dy * wall_adjustment_factor ) THEN 1219 WRITE( message_string, * ) 'grid anisotropy exceeds ', &1220 'threshold given by only local', &1221 ' horizontal reduction of near_wall ', &1222 'mixing length l_wall', &1223 ' starting from height level k = ', k, &1412 WRITE( message_string, * ) 'grid anisotropy exceeds ', & 1413 'threshold given by only local', & 1414 ' &horizontal reduction of near_wall ', & 1415 'mixing length l_wall', & 1416 ' &starting from height level k = ', k, & 1224 1417 '.' 1225 1418 CALL message( 'init_grid', 'PA0202', 0, 1, 0, 6, 0 ) … … 1320 1513 ! 1321 1514 !-- Initialize the mixing length in case of a RANS simulation 1515 ALLOCATE( l_black(nzb:nzt+1) ) 1322 1516 1323 1517 ! 1324 1518 !-- Calculate mixing length according to Blackadar (1962) 1325 1519 IF ( f /= 0.0_wp ) THEN 1326 l_max = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) /&1327 ABS( f ) + 1 E-10_wp1520 l_max = 2.7E-4_wp * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) / & 1521 ABS( f ) + 1.0E-10_wp 1328 1522 ELSE 1329 1523 l_max = 30.0_wp … … 1338 1532 ! 1339 1533 !-- Gather topography information of whole domain 1340 ! ## TODO:reduce amount of data sent by MPI call1341 ! ##By now, a whole global 3D-array is sent and received with1342 ! ##MPI_ALLREDUCE although most of the array is 0. This can be1343 ! ##drastically reduced if only the local subarray is sent and stored1344 ! ##in a global array. For that, an MPI data type or subarray must be1345 ! ##defined.1346 ! ##2018-03-19, gronemeier1534 !> @todo reduce amount of data sent by MPI call 1535 !> By now, a whole global 3D-array is sent and received with 1536 !> MPI_ALLREDUCE although most of the array is 0. This can be 1537 !> drastically reduced if only the local subarray is sent and stored 1538 !> in a global array. For that, an MPI data type or subarray must be 1539 !> defined. 1540 !> 2018-03-19, gronemeier 1347 1541 ALLOCATE( wall_flags_0_global(nzb:nzt+1,0:ny,0:nx) ) 1348 1542 … … 1373 1567 ENDDO 1374 1568 ENDDO 1569 1570 l_wall(nzb,:,:) = l_black(nzb) 1571 l_wall(nzt+1,:,:) = l_black(nzt+1) 1375 1572 ! 1376 1573 !-- Limit mixing length to either nearest wall or Blackadar mixing length. … … 1415 1612 IF ( rad_k_b /= 0 .OR. rad_k_t /= 0 ) THEN 1416 1613 1417 ! ## NOTE:shape of vicinity is larger in z direction1418 ! ##Shape of vicinity is two grid points larger than actual search1419 ! ##radius in vertical direction. The first and last grid point is1420 ! ##always set to 1 to asure correct detection of topography. See1421 ! ##function "shortest_distance" for details.1422 ! ##2018-03-16, gronemeier1614 !> @note shape of vicinity is larger in z direction 1615 !> Shape of vicinity is two grid points larger than actual search 1616 !> radius in vertical direction. The first and last grid point is 1617 !> always set to 1 to asure correct detection of topography. See 1618 !> function "shortest_distance" for details. 1619 !> 2018-03-16, gronemeier 1423 1620 ALLOCATE( vicinity(-rad_k-1:rad_k+1,-rad_j:rad_j,-rad_i:rad_i) ) 1424 1621 ALLOCATE( vic_yz(0:rad_k+1,0:rad_j) ) … … 1601 1798 ELSE !Check if (i,j,k) belongs to atmosphere 1602 1799 1603 l_wall(k,j,i) = -999.01800 l_wall(k,j,i) = l_black(k) 1604 1801 1605 1802 ENDIF … … 1630 1827 !> (pos_i/jj/kk), where (jj/kk) is the position of the maximum of 'array' 1631 1828 !> closest to the origin (0/0) of 'array'. 1829 !> @todo this part of PALM does not reproduce the same results for optimized 1830 !> and debug options for the compiler. This should be fixed 1632 1831 !------------------------------------------------------------------------------! 1633 1832 REAL FUNCTION shortest_distance( array, orientation, pos_i ) … … 1659 1858 IF ( orientation ) THEN !if array is oriented upwards 1660 1859 DO jj = 0, rad_j 1661 shortest_distance = MIN( shortest_distance, & 1662 SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2 & 1663 + MAX(jj*dy-0.5*dy,0.0)**2 & 1664 + MAX(zw(loc_k(jj)+k-1)-zu(k),0.0)**2 & 1665 ) & 1666 ) 1860 shortest_distance = & 1861 MIN( shortest_distance, & 1862 SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2 & 1863 + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2 & 1864 + MAX(zw(loc_k(jj)+k-1)-zu(k), 0.0_wp)**2 & 1865 ) & 1866 ) 1667 1867 ENDDO 1668 1868 ELSE !if array is oriented downwards 1669 ! ## NOTE:MAX within zw required to circumvent error at domain border1670 ! ##At the domain border, if non-cyclic boundary is present, the1671 ! ##index for zw could be -1, which will be errorneous (zw(-1) does1672 ! ##not exist). The MAX function limits the index to be at least 0.1869 !> @note MAX within zw required to circumvent error at domain border 1870 !> At the domain border, if non-cyclic boundary is present, the 1871 !> index for zw could be -1, which will be errorneous (zw(-1) does 1872 !> not exist). The MAX function limits the index to be at least 0. 1673 1873 DO jj = 0, rad_j 1674 shortest_distance = MIN( shortest_distance, & 1675 SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2 & 1676 + MAX(jj*dy-0.5*dy,0.0)**2 & 1677 + MAX(zu(k)-zw(MAX(k-loc_k(jj), & 1678 0_iwp)), & 1679 0.0)**2 & 1680 ) & 1681 ) 1874 shortest_distance = & 1875 MIN( shortest_distance, & 1876 SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2 & 1877 + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2 & 1878 + MAX(zu(k)-zw(MAX(k-loc_k(jj),0_iwp)), 0.0_wp)**2 & 1879 ) & 1880 ) 1682 1881 ENDDO 1683 1882 ENDIF 1684 1883 1685 1884 END FUNCTION 1686 1885 … … 2004 2203 ! 2005 2204 !-- Use special boundary condition in case of TKE-e closure 2205 !> @todo do the same for usm and lsm surfaces 2206 !> 2018-06-05, gronemeier 2006 2207 IF ( rans_tke_e ) THEN 2007 2208 DO i = nxl, nxr … … 2011 2212 DO m = surf_s, surf_e 2012 2213 k = surf_def_h(0)%k(m) 2013 e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_ m**22214 e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_0**2 2014 2215 ENDDO 2015 2216 ENDDO … … 2092 2293 DO k = nzb+1, nzt 2093 2294 ! tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k) 2094 tend(k,j,i) = tend(k,j,i) + c_1 * c_ mu * f / c_h & !###needs revision2295 tend(k,j,i) = tend(k,j,i) + c_1 * c_0**4 * f / c_4 & !> @todo needs revision 2095 2296 / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i)) & 2096 2297 * SQRT(e(k,j,i)) * produc(k,j,i) … … 2103 2304 ! 2104 2305 !-- Additional sink term for flows through plant canopies 2105 ! IF ( plant_canopy ) CALL pcm_tendency( ? ) ! ###what to do with this?2106 2107 ! CALL user_actions( 'diss-tendency' ) ! ###not yet implemented2306 ! IF ( plant_canopy ) CALL pcm_tendency( ? ) !> @query what to do with this? 2307 2308 ! CALL user_actions( 'diss-tendency' ) !> @todo not yet implemented 2108 2309 2109 2310 ! … … 2181 2382 USE arrays_3d, & 2182 2383 ONLY: ddzu, diss_l_diss, diss_l_e, diss_s_diss, diss_s_e, & 2183 flux_l_diss, flux_l_e, flux_s_diss, flux_s_e 2384 flux_l_diss, flux_l_e, flux_s_diss, flux_s_e,& 2385 u_p,v_p,w_p 2184 2386 2185 2387 USE control_parameters, & 2186 2388 ONLY: f, tsc 2389 2390 USE grid_variables, & 2391 ONLY: dx, dy 2187 2392 2188 2393 USE surface_mod, & … … 2190 2395 surf_usm_v 2191 2396 2397 use indices, only: nx, ny 2398 2192 2399 IMPLICIT NONE 2193 2400 2194 2401 INTEGER(iwp) :: i !< loop index x direction 2195 INTEGER(iwp) :: i_omp !< 2402 INTEGER(iwp) :: i_omp !< first loop index of i-loop in prognostic_equations 2196 2403 INTEGER(iwp) :: j !< loop index y direction 2197 2404 INTEGER(iwp) :: k !< loop index z direction 2405 INTEGER(iwp) :: l !< loop index 2198 2406 INTEGER(iwp) :: m !< loop index 2199 2407 INTEGER(iwp) :: surf_e !< end index of surface elements at given i-j position 2200 2408 INTEGER(iwp) :: surf_s !< start index of surface elements at given i-j position 2201 INTEGER(iwp) :: tn !< 2202 2203 REAL(wp), DIMENSION(nzb:nzt+1) :: advec !< advection term of TKE tendency 2204 REAL(wp), DIMENSION(nzb:nzt+1) :: produc !< production term of TKE tendency 2409 INTEGER(iwp) :: tn !< task number of openmp task 2410 2411 INTEGER(iwp) :: pis = 32 !< debug variable, print from i=pis !> @todo remove later 2412 INTEGER(iwp) :: pie = 32 !< debug variable, print until i=pie !> @todo remove later 2413 INTEGER(iwp) :: pjs = 26 !< debug variable, print from j=pjs !> @todo remove later 2414 INTEGER(iwp) :: pje = 26 !< debug variable, print until j=pje !> @todo remove later 2415 INTEGER(iwp) :: pkb = 1 !< debug variable, print from k=pkb !> @todo remove later 2416 INTEGER(iwp) :: pkt = 7 !< debug variable, print until k=pkt !> @todo remove later 2417 2418 REAL(wp), DIMENSION(nzb:nzt+1) :: dum_adv !< debug variable !> @todo remove later 2419 REAL(wp), DIMENSION(nzb:nzt+1) :: dum_pro !< debug variable !> @todo remove later 2420 REAL(wp), DIMENSION(nzb:nzt+1) :: dum_dif !< debug variable !> @todo remove later 2421 2422 5555 FORMAT(A,7(1X,E12.5)) !> @todo remove later 2205 2423 2206 2424 ! … … 2224 2442 ENDIF 2225 2443 2226 advec(:) = tend(:,j,i)2227 2228 CALL production_e( i, j )2229 2230 produc(:) = tend(:,j,i) - advec(:)2444 dum_adv = tend(:,j,i) !> @todo remove later 2445 2446 CALL production_e( i, j, .FALSE. ) 2447 2448 dum_pro = tend(:,j,i) - dum_adv !> @todo remove later 2231 2449 2232 2450 IF ( .NOT. humidity ) THEN … … 2239 2457 CALL diffusion_e( i, j, vpt, pt_reference ) 2240 2458 ENDIF 2459 2460 dum_dif = tend(:,j,i) - dum_adv - dum_pro !> @todo remove later 2241 2461 2242 2462 ! … … 2268 2488 DO m = surf_s, surf_e 2269 2489 k = surf_def_h(0)%k(m) 2270 e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_m**2 2490 e_p(k,j,i) = ( surf_def_h(0)%us(m) / c_0 )**2 2491 ENDDO 2492 2493 DO l = 0, 3 2494 surf_s = surf_def_v(l)%start_index(j,i) 2495 surf_e = surf_def_v(l)%end_index(j,i) 2496 DO m = surf_s, surf_e 2497 k = surf_def_v(l)%k(m) 2498 e_p(k,j,i) = ( surf_def_v(l)%us(m) / c_0 )**2 2499 ENDDO 2271 2500 ENDDO 2272 2501 ENDIF … … 2287 2516 ENDIF 2288 2517 ENDIF 2518 2519 ! if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then !> @todo remove later 2520 ! WRITE(9, *) '------' 2521 ! WRITE(9, '(A,F8.3,1X,F8.3,1X,I2)') 't, dt, int_ts:', simulated_time, dt_3d, intermediate_timestep_count 2522 ! WRITE(9, *) 'i:', i 2523 ! WRITE(9, *) 'j:', j 2524 ! WRITE(9, *) 'k:', pkb, ' - ', pkt 2525 ! WRITE(9, *) '---' 2526 ! WRITE(9, *) 'e:' 2527 ! WRITE(9, 5555) 'adv :', dum_adv(pkb:pkt) 2528 ! WRITE(9, 5555) 'pro :', dum_pro(pkb:pkt) 2529 ! WRITE(9, 5555) 'dif :', dum_dif(pkb:pkt) 2530 ! WRITE(9, 5555) 'tend:', tend(pkb:pkt,j,i) 2531 ! WRITE(9, 5555) 'e_p :', e_p(pkb:pkt,j,i) 2532 ! WRITE(9, 5555) 'e :', e(pkb:pkt,j,i) 2533 ! FLUSH(9) 2534 ! endif 2289 2535 2290 2536 ENDIF ! TKE equation … … 2309 2555 ENDIF 2310 2556 2557 IF ( intermediate_timestep_count == 1 ) diss_adve1(:,j,i) = tend(:,j,i) !> @todo remove later 2558 IF ( intermediate_timestep_count == 2 ) diss_adve2(:,j,i) = tend(:,j,i) 2559 IF ( intermediate_timestep_count == 3 ) diss_adve3(:,j,i) = tend(:,j,i) 2560 2311 2561 ! 2312 2562 !-- Production of TKE dissipation rate 2313 DO k = nzb+1, nzt 2314 ! tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k) 2315 tend(k,j,i) = tend(k,j,i) + c_1 * c_mu * f / c_h & !### needs revision 2316 / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i)) & 2317 * SQRT(e(k,j,i)) * produc(k) 2318 ENDDO 2319 2563 CALL production_e( i, j, .TRUE. ) 2564 2565 IF ( intermediate_timestep_count == 1 ) diss_prod1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) !> @todo remove later 2566 IF ( intermediate_timestep_count == 2 ) diss_prod2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i) 2567 IF ( intermediate_timestep_count == 3 ) diss_prod3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i) 2568 2569 dum_pro = tend(:,j,i) - dum_adv !> @todo remove later 2570 2571 ! 2572 !-- Diffusion term of TKE dissipation rate 2320 2573 CALL diffusion_diss( i, j ) 2321 2574 2575 IF ( intermediate_timestep_count == 1 ) diss_diff1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) - diss_prod1(:,j,i) !> @todo remove later 2576 IF ( intermediate_timestep_count == 2 ) diss_diff2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i) - diss_prod2(:,j,i) 2577 IF ( intermediate_timestep_count == 3 ) diss_diff3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i) - diss_prod3(:,j,i) 2578 IF ( intermediate_timestep_count == 3 ) dummy3(:,j,i) = km(:,j,i) 2579 2580 dum_dif = tend(:,j,i) - dum_adv - dum_pro !> @todo remove later 2581 2322 2582 ! 2323 2583 !-- Additional sink term for flows through plant canopies 2324 ! IF ( plant_canopy ) CALL pcm_tendency( i, j, ? ) ! ###not yet implemented2325 2326 ! CALL user_actions( i, j, 'diss-tendency' ) ! ###not yet implemented2584 ! IF ( plant_canopy ) CALL pcm_tendency( i, j, ? ) !> @todo not yet implemented 2585 2586 ! CALL user_actions( i, j, 'diss-tendency' ) !> @todo not yet implemented 2327 2587 2328 2588 ! … … 2338 2598 BTEST( wall_flags_0(k,j,i), 0 )& 2339 2599 ) 2340 IF ( diss_p(k,j,i) <= 0.0_wp ) diss_p(k,j,i) = 0.1_wp * diss(k,j,i)2341 2600 ENDDO 2342 2601 … … 2350 2609 ENDDO 2351 2610 2611 DO l = 0, 1 2612 surf_s = surf_def_v(l)%start_index(j,i) 2613 surf_e = surf_def_v(l)%end_index(j,i) 2614 DO m = surf_s, surf_e 2615 k = surf_def_v(l)%k(m) 2616 diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dy ) 2617 ENDDO 2618 ENDDO 2619 2620 DO l = 2, 3 2621 surf_s = surf_def_v(l)%start_index(j,i) 2622 surf_e = surf_def_v(l)%end_index(j,i) 2623 DO m = surf_s, surf_e 2624 k = surf_def_v(l)%k(m) 2625 diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dx ) 2626 ENDDO 2627 ENDDO 2352 2628 ! 2353 2629 !-- Calculate tendencies for the next Runge-Kutta step … … 2365 2641 ENDIF 2366 2642 ENDIF 2367 2368 ! IF ( intermediate_timestep_count == 1 ) dummy1(:,j,i) = e_p(:,j,i) 2369 ! IF ( intermediate_timestep_count == 2 ) dummy2(:,j,i) = e_p(:,j,i) 2370 ! IF ( intermediate_timestep_count == 3 ) dummy3(:,j,i) = e_p(:,j,i) 2643 ! 2644 !-- Limit change of diss to be between -90% and +100%. Also, set an absolute 2645 !-- minimum value 2646 DO k = nzb, nzt+1 2647 diss_p(k,j,i) = MIN( MAX( diss_p(k,j,i), & 2648 0.1_wp * diss(k,j,i), & 2649 0.0001_wp ), & 2650 2.0_wp * diss(k,j,i) ) 2651 ENDDO 2652 2653 IF ( intermediate_timestep_count == 1 ) dummy1(:,j,i) = diss_p(:,j,i) !> @todo remove later 2654 IF ( intermediate_timestep_count == 2 ) dummy2(:,j,i) = diss_p(:,j,i) 2655 2656 ! if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then !> @todo remove later 2657 ! WRITE(9, *) '---' 2658 ! WRITE(9, *) 'diss:' 2659 ! WRITE(9, 5555) 'adv :', dum_adv(pkb:pkt) 2660 ! WRITE(9, 5555) 'pro :', dum_pro(pkb:pkt) 2661 ! WRITE(9, 5555) 'dif :', dum_dif(pkb:pkt) 2662 ! WRITE(9, 5555) 'tend :', tend(pkb:pkt,j,i) 2663 ! WRITE(9, 5555) 'diss_p:', diss_p(pkb:pkt,j,i) 2664 ! WRITE(9, 5555) 'diss :', diss(pkb:pkt,j,i) 2665 ! WRITE(9, *) '---' 2666 ! WRITE(9, 5555) 'km :', km(pkb:pkt,j,i) 2667 ! flush(9) 2668 ! endif 2371 2669 2372 2670 ENDIF ! dissipation equation … … 2382 2680 !> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is 2383 2681 !> not considered well! 2682 !> @todo Adjust production term in case of rans_tke_e simulation 2384 2683 !------------------------------------------------------------------------------! 2385 2684 SUBROUTINE production_e … … 3099 3398 !> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is 3100 3399 !> not considered well! 3400 !> @todo non-neutral case is not yet considered for RANS mode 3101 3401 !------------------------------------------------------------------------------! 3102 SUBROUTINE production_e_ij( i, j )3402 SUBROUTINE production_e_ij( i, j, diss_production ) 3103 3403 3104 3404 USE arrays_3d, & … … 3121 3421 3122 3422 IMPLICIT NONE 3423 3424 LOGICAL :: diss_production 3123 3425 3124 3426 INTEGER(iwp) :: i !< running index x-direction … … 3153 3455 REAL(wp), DIMENSION(nzb+1:nzt) :: dwdy !< Gradient of w-component in y-direction 3154 3456 REAL(wp), DIMENSION(nzb+1:nzt) :: dwdz !< Gradient of w-component in z-direction 3155 3457 REAL(wp), DIMENSION(nzb+1:nzt) :: tend_temp !< temporal tendency 3156 3458 3157 3459 IF ( constant_flux_layer ) THEN … … 3348 3650 ENDDO 3349 3651 3652 ! IF ( .NOT. rans_tke_e ) THEN 3653 3350 3654 DO k = nzb+1, nzt 3351 3655 3352 def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) + & 3353 dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2 + & 3354 dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2 + & 3355 2.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) ) 3356 3357 ! 3358 !-- Production term according to Kato and Launder (1993) 3359 ! def = SQRT( ( dudx(k) + dudy(k) + dudz(k) + & 3360 ! dvdx(k) + dvdy(k) + dvdz(k) + & 3361 ! dwdx(k) + dwdy(k) + dwdz(k) & 3362 ! )**4 - & 3363 ! ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 + & 3364 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3365 ! dudz(k) * dwdx(k) + & 3366 ! dvdz(k) * dwdy(k) ) & 3367 ! )**2 & 3368 ! ) 3656 def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) + & 3657 dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2 + & 3658 dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2 + & 3659 2.0_wp * ( dvdx(k)*dudy(k) + & 3660 dwdx(k)*dudz(k) + & 3661 dwdy(k)*dvdz(k) ) 3369 3662 3370 3663 IF ( def < 0.0_wp ) def = 0.0_wp … … 3372 3665 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 3373 3666 3374 tend (k,j,i) = tend(k,j,i) +km(k,j,i) * def * flag3667 tend_temp(k) = km(k,j,i) * def * flag 3375 3668 3376 3669 ENDDO 3377 3670 3671 ! ELSE 3672 ! 3673 ! DO k = nzb+1, nzt 3674 ! ! 3675 ! !-- Production term according to Kato and Launder (1993) 3676 ! def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 + & 3677 ! dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 + & 3678 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3679 ! dvdz(k) * dwdy(k) + & 3680 ! dwdx(k) * dudz(k) ) ) & 3681 ! * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 + & 3682 ! dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 - & 3683 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3684 ! dvdz(k) * dwdy(k) + & 3685 ! dwdx(k) * dudz(k) ) ) & 3686 ! ) 3687 ! 3688 ! IF ( def < 0.0_wp ) def = 0.0_wp 3689 ! 3690 ! flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 3691 ! 3692 ! tend_temp(k) = km(k,j,i) * def * flag 3693 ! 3694 ! ENDDO 3695 ! 3696 ! ENDIF 3697 3698 ELSE ! not constant_flux_layer 3699 3700 ! IF ( .NOT. rans_tke_e ) THEN 3701 ! 3702 !-- Calculate TKE production by shear. Here, no additional 3703 !-- wall-bounded code is considered. 3704 !-- Why? 3705 DO k = nzb+1, nzt 3706 3707 dudx(k) = ( u(k,j,i+1) - u(k,j,i) ) * ddx 3708 dudy(k) = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) - & 3709 u(k,j-1,i) - u(k,j-1,i+1) ) * ddy 3710 dudz(k) = 0.5_wp * ( u(k+1,j,i) + u(k+1,j,i+1) - & 3711 u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k) 3712 3713 dvdx(k) = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) - & 3714 v(k,j,i-1) - v(k,j+1,i-1) ) * ddx 3715 dvdy(k) = ( v(k,j+1,i) - v(k,j,i) ) * ddy 3716 dvdz(k) = 0.5_wp * ( v(k+1,j,i) + v(k+1,j+1,i) - & 3717 v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k) 3718 3719 dwdx(k) = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) - & 3720 w(k,j,i-1) - w(k-1,j,i-1) ) * ddx 3721 dwdy(k) = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) - & 3722 w(k,j-1,i) - w(k-1,j-1,i) ) * ddy 3723 dwdz(k) = ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) 3724 3725 def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) + & 3726 dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2 + & 3727 dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2 + & 3728 2.0_wp * ( dvdx(k)*dudy(k) + & 3729 dwdx(k)*dudz(k) + & 3730 dwdy(k)*dvdz(k) ) 3731 3732 IF ( def < 0.0_wp ) def = 0.0_wp 3733 3734 flag = MERGE( 1.0_wp, 0.0_wp, & 3735 BTEST( wall_flags_0(k,j,i), 29 ) ) 3736 tend_temp(k) = km(k,j,i) * def * flag 3737 3738 ENDDO 3739 3740 ! ELSE 3741 ! 3742 ! DO k = nzb+1, nzt 3743 ! 3744 ! dudx(k) = ( u(k,j,i+1) - u(k,j,i) ) * ddx 3745 ! dudy(k) = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) - & 3746 ! u(k,j-1,i) - u(k,j-1,i+1) ) * ddy 3747 ! dudz(k) = 0.5_wp * ( u(k+1,j,i) + u(k+1,j,i+1) - & 3748 ! u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k) 3749 ! 3750 ! dvdx(k) = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) - & 3751 ! v(k,j,i-1) - v(k,j+1,i-1) ) * ddx 3752 ! dvdy(k) = ( v(k,j+1,i) - v(k,j,i) ) * ddy 3753 ! dvdz(k) = 0.5_wp * ( v(k+1,j,i) + v(k+1,j+1,i) - & 3754 ! v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k) 3755 ! 3756 ! dwdx(k) = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) - & 3757 ! w(k,j,i-1) - w(k-1,j,i-1) ) * ddx 3758 ! dwdy(k) = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) - & 3759 ! w(k,j-1,i) - w(k-1,j-1,i) ) * ddy 3760 ! dwdz(k) = ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) 3761 ! ! 3762 ! !-- Production term according to Kato and Launder (1993) 3763 ! def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 + & 3764 ! dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 + & 3765 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3766 ! dvdz(k) * dwdy(k) + & 3767 ! dwdx(k) * dudz(k) ) ) & 3768 ! * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 + & 3769 ! dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 - & 3770 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3771 ! dvdz(k) * dwdy(k) + & 3772 ! dwdx(k) * dudz(k) ) ) & 3773 ! ) 3774 ! 3775 ! IF ( def < 0.0_wp ) def = 0.0_wp 3776 ! 3777 ! flag = MERGE( 1.0_wp, 0.0_wp, & 3778 ! BTEST( wall_flags_0(k,j,i), 29 ) ) 3779 ! tend_temp(k) = km(k,j,i) * def * flag 3780 ! 3781 ! ENDDO 3782 ! 3783 ! ENDIF 3784 3785 ENDIF 3786 3787 IF ( .NOT. diss_production ) THEN 3788 ! 3789 !-- Production term in case of TKE production 3790 DO k = nzb+1, nzt 3791 tend(k,j,i) = tend(k,j,i) + tend_temp(k) 3792 ENDDO 3378 3793 ELSE 3379 3794 ! 3380 !-- Calculate TKE production by shear. Here, no additional 3381 !-- wall-bounded code is considered. 3382 !-- Why? 3795 !-- Production term in case of dissipation-rate production (rans_tke_e) 3383 3796 DO k = nzb+1, nzt 3384 3797 3385 dudx(k) = ( u(k,j,i+1) - u(k,j,i) ) * ddx 3386 dudy(k) = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) - & 3387 u(k,j-1,i) - u(k,j-1,i+1) ) * ddy 3388 dudz(k) = 0.5_wp * ( u(k+1,j,i) + u(k+1,j,i+1) - & 3389 u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k) 3390 3391 dvdx(k) = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) - & 3392 v(k,j,i-1) - v(k,j+1,i-1) ) * ddx 3393 dvdy(k) = ( v(k,j+1,i) - v(k,j,i) ) * ddy 3394 dvdz(k) = 0.5_wp * ( v(k+1,j,i) + v(k+1,j+1,i) - & 3395 v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k) 3396 3397 dwdx(k) = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) - & 3398 w(k,j,i-1) - w(k-1,j,i-1) ) * ddx 3399 dwdy(k) = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) - & 3400 w(k,j-1,i) - w(k-1,j-1,i) ) * ddy 3401 dwdz(k) = ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) 3402 3403 def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) + & 3404 dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2 + & 3405 dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2 + & 3406 2.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) ) 3407 3408 ! 3409 !-- Production term according to Kato and Launder (1993) 3410 ! def = SQRT( ( dudx(k) + dudy(k) + dudz(k) + & 3411 ! dvdx(k) + dvdy(k) + dvdz(k) + & 3412 ! dwdx(k) + dwdy(k) + dwdz(k) & 3413 ! )**4 - & 3414 ! ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 + & 3415 ! 2.0_wp * ( dudy(k) * dvdx(k) + & 3416 ! dudz(k) * dwdx(k) + & 3417 ! dvdz(k) * dwdy(k) ) & 3418 ! )**2 & 3419 ! ) 3420 3421 IF ( def < 0.0_wp ) def = 0.0_wp 3422 3423 flag = MERGE( 1.0_wp, 0.0_wp, & 3424 BTEST( wall_flags_0(k,j,i), 29 ) ) 3425 tend(k,j,i) = tend(k,j,i) + km(k,j,i) * def * flag 3426 3798 ! Standard TKE-e closure 3799 tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i) & 3800 /( e(k,j,i) + 1.0E-20_wp ) & 3801 * c_1 3802 ! ! Production according to Koblitz (2013) 3803 ! tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i) & 3804 ! /( e(k,j,i) + 1.0E-20_wp ) & 3805 ! * ( c_1 + ( c_2 - c_1 ) & 3806 ! * l_wall(k,j,i) / l_max ) 3807 ! ! Production according to Detering and Etling (1985) 3808 ! !> @todo us is not correct if there are vertical walls 3809 ! tend(k,j,i) = tend(k,j,i) + tend_temp(k) * SQRT(e(k,j,i)) & 3810 ! * c_1 * c_0**3 / c_4 * f & 3811 ! / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i)) 3427 3812 ENDDO 3428 3429 3813 ENDIF 3430 3814 … … 3822 4206 REAL(wp) :: l !< mixing length 3823 4207 REAL(wp) :: ll !< adjusted l 3824 REAL(wp) :: var_reference !< 4208 REAL(wp) :: var_reference !< reference temperature 3825 4209 3826 4210 #if defined( __nopointer ) 3827 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< 4211 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< temperature 3828 4212 #else 3829 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< 4213 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< temperature 3830 4214 #endif 3831 4215 REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) :: dissipation !< TKE dissipation … … 3837 4221 DO j = nys, nyn 3838 4222 DO k = nzb+1, nzt 4223 ! 4224 !-- Predetermine flag to mask topography 4225 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 3839 4226 3840 4227 ! … … 3851 4238 CALL mixing_length_rans( i, j, k, l, ll, var, var_reference ) 3852 4239 3853 dissipation(k,j) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll 4240 dissipation(k,j) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll 4241 4242 diss(k,j,i) = dissipation(k,j) * flag 3854 4243 3855 4244 ELSEIF ( rans_tke_e ) THEN … … 3859 4248 ENDIF 3860 4249 3861 ! 3862 !-- Predetermine flag to mask topography 3863 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 3864 3865 tend(k,j,i) = tend(k,j,i) & 3866 + ( & 4250 tend(k,j,i) = tend(k,j,i) + ( & 4251 ( & 3867 4252 ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) ) & 3868 4253 - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) ) & 3869 ) * ddx2 * flag&3870 + (&4254 ) * ddx2 * flag & 4255 + ( & 3871 4256 ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) ) & 3872 4257 - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) ) & 3873 ) * ddy2 * flag&3874 + (&4258 ) * ddy2 * flag & 4259 + ( & 3875 4260 ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) & 3876 4261 * rho_air_zw(k) & 3877 4262 - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k) & 3878 4263 * rho_air_zw(k-1) & 3879 ) * ddzw(k) * drho_air(k) * flag & 4264 ) * ddzw(k) * drho_air(k) & 4265 ) * flag * dsig_e & 3880 4266 - dissipation(k,j) * flag 3881 4267 … … 3959 4345 REAL(wp) :: l !< mixing length 3960 4346 REAL(wp) :: ll !< adjusted l 3961 REAL(wp) :: var_reference !< 4347 REAL(wp) :: var_reference !< reference temperature 3962 4348 3963 4349 #if defined( __nopointer ) 3964 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< 4350 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< temperature 3965 4351 #else 3966 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< 4352 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< temperature 3967 4353 #endif 3968 4354 REAL(wp), DIMENSION(nzb+1:nzt) :: dissipation !< dissipation of TKE … … 3991 4377 CALL mixing_length_rans( i, j, k, l, ll, var, var_reference ) 3992 4378 3993 dissipation(k) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll 4379 dissipation(k) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll 4380 4381 diss(k,j,i) = dissipation(k) * flag 3994 4382 3995 4383 ELSEIF ( rans_tke_e ) THEN … … 4001 4389 ! 4002 4390 !-- Calculate the tendency term 4003 tend(k,j,i) = tend(k,j,i) 4004 +( &4391 tend(k,j,i) = tend(k,j,i) + ( & 4392 ( & 4005 4393 ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) ) & 4006 4394 - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) ) & 4007 ) * ddx2 * flag / sig_e&4395 ) * ddx2 & 4008 4396 + ( & 4009 4397 ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) ) & 4010 4398 - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) ) & 4011 ) * ddy2 * flag / sig_e&4399 ) * ddy2 & 4012 4400 + ( & 4013 4401 ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) & … … 4015 4403 - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k) & 4016 4404 * rho_air_zw(k-1) & 4017 ) * ddzw(k) * drho_air(k) * flag / sig_e & 4018 - dissipation(k) * flag 4405 ) * ddzw(k) * drho_air(k) & 4406 ) * flag * dsig_e & 4407 - dissipation(k) * flag 4019 4408 4020 4409 ENDDO … … 4082 4471 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) ) 4083 4472 4084 tend(k,j,i) = tend(k,j,i) 4085 +( &4473 tend(k,j,i) = tend(k,j,i) + & 4474 ( ( & 4086 4475 ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) ) & 4087 4476 - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) ) & 4088 ) * ddx2 * flag&4477 ) * ddx2 & 4089 4478 + ( & 4090 4479 ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) ) & 4091 4480 - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) ) & 4092 ) * ddy2 * flag&4481 ) * ddy2 & 4093 4482 + ( & 4094 4483 ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1) & … … 4096 4485 - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k) & 4097 4486 * rho_air_zw(k-1) & 4098 ) * ddzw(k) * drho_air(k) * flag & 4099 - c_2 * diss(k,j,i)**2 & 4100 / ( e(k,j,i) + 1.0E-20_wp ) * flag 4487 ) * ddzw(k) * drho_air(k) & 4488 ) * flag * dsig_diss & 4489 - c_2 * diss(k,j,i)**2 & 4490 / ( e(k,j,i) + 1.0E-20_wp ) * flag 4101 4491 4102 4492 ENDDO … … 4123 4513 IMPLICIT NONE 4124 4514 4125 INTEGER(iwp) :: i !< running index x direction 4126 INTEGER(iwp) :: j !< running index y direction 4127 INTEGER(iwp) :: k !< running index z direction 4128 4129 REAL(wp) :: flag !< flag to mask topography 4130 4131 REAL(wp), DIMENSION(nzb+1:nzt) :: tend_temp 4515 INTEGER(iwp) :: i !< running index x direction 4516 INTEGER(iwp) :: j !< running index y direction 4517 INTEGER(iwp) :: k !< running index z direction 4518 4519 REAL(wp) :: flag !< flag to mask topography 4132 4520 4133 4521 ! … … 4141 4529 ! 4142 4530 !-- Calculate the tendency term 4143 tend_temp(k) = ( & 4531 tend(k,j,i) = tend(k,j,i) + & 4532 ( ( & 4144 4533 ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) ) & 4145 4534 - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) ) & 4146 ) * ddx2 * flag / sig_diss&4535 ) * ddx2 & 4147 4536 + ( & 4148 4537 ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) ) & 4149 4538 - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) ) & 4150 ) * ddy2 * flag / sig_diss&4539 ) * ddy2 & 4151 4540 + ( & 4152 4541 ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1) & … … 4154 4543 - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k) & 4155 4544 * rho_air_zw(k-1) & 4156 ) * ddzw(k) * drho_air(k) * flag / sig_diss & 4157 - c_2 * diss(k,j,i)**2 & 4158 / ( e(k,j,i) + 1.0E-20_wp ) * flag 4159 4160 tend(k,j,i) = tend(k,j,i) + tend_temp(k) 4545 ) * ddzw(k) * drho_air(k) & 4546 ) * flag * dsig_diss & 4547 - c_2 * diss(k,j,i)**2 / ( e(k,j,i) + 1.0E-20_wp ) * flag 4161 4548 4162 4549 ENDDO … … 4290 4677 !> Computation of the turbulent diffusion coefficients for momentum and heat 4291 4678 !> according to Prandtl-Kolmogorov. 4679 !> @todo consider non-default surfaces 4292 4680 !------------------------------------------------------------------------------! 4293 4681 SUBROUTINE tcm_diffusivities( var, var_reference ) … … 4297 4685 ONLY: e_min, outflow_l, outflow_n, outflow_r, outflow_s 4298 4686 4687 USE grid_variables, & 4688 ONLY: dx, dy 4689 4299 4690 USE statistics, & 4300 4691 ONLY : rmask, sums_l_l 4301 4692 4302 4693 USE surface_mod, & 4303 ONLY : bc_h, surf_def_h 4694 ONLY : bc_h, surf_def_h, surf_def_v 4304 4695 4305 4696 IMPLICIT NONE 4306 4697 4307 INTEGER(iwp) :: i !< 4308 INTEGER(iwp) :: j !< 4309 INTEGER(iwp) :: k !< 4310 INTEGER(iwp) :: m !< 4311 INTEGER(iwp) :: n !< 4312 INTEGER(iwp) :: omp_get_thread_num !< 4313 INTEGER(iwp) :: sr !< 4314 INTEGER(iwp) :: tn !< 4315 4316 REAL(wp) :: flag !< 4317 REAL(wp) :: l !< 4318 REAL(wp) :: ll !< 4319 REAL(wp) :: var_reference !< 4698 INTEGER(iwp) :: i !< loop index 4699 INTEGER(iwp) :: j !< loop index 4700 INTEGER(iwp) :: k !< loop index 4701 INTEGER(iwp) :: m !< loop index 4702 INTEGER(iwp) :: n !< loop index 4703 INTEGER(iwp) :: omp_get_thread_num !< opemmp function to get thread number 4704 INTEGER(iwp) :: sr !< statistic region 4705 INTEGER(iwp) :: tn !< thread number 4706 4707 REAL(wp) :: flag !< topography flag 4708 REAL(wp) :: l !< mixing length 4709 REAL(wp) :: ll !< adjusted mixing length 4710 REAL(wp) :: var_reference !< reference temperature 4320 4711 4321 4712 #if defined( __nopointer ) 4322 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< 4713 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: var !< temperature 4323 4714 #else 4324 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< 4715 REAL(wp), DIMENSION(:,:,:), POINTER :: var !< temperature 4325 4716 #endif 4326 4717 … … 4365 4756 ! 4366 4757 !-- Compute diffusion coefficients for momentum and heat 4367 km(k,j,i) = c_ m* l * SQRT( e(k,j,i) ) * flag4758 km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag 4368 4759 kh(k,j,i) = ( 1.0_wp + 2.0_wp * l / ll ) * km(k,j,i) * flag 4369 4760 ! … … 4391 4782 ! 4392 4783 !-- Compute diffusion coefficients for momentum and heat 4393 km(k,j,i) = c_ m* l * SQRT( e(k,j,i) ) * flag4784 km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag 4394 4785 kh(k,j,i) = km(k,j,i) / prandtl_number * flag 4395 4786 ! … … 4414 4805 ! 4415 4806 !-- Compute diffusion coefficients for momentum and heat 4416 km(k,j,i) = c_ mu * e(k,j,i)**2 / ( diss(k,j,i) + 1.E-10) * flag4807 km(k,j,i) = c_0**4 * e(k,j,i)**2 / ( diss(k,j,i) + 1.0E-30_wp ) * flag 4417 4808 kh(k,j,i) = km(k,j,i) / prandtl_number * flag 4809 ! 4810 !-- Summation for averaged profile of mixing length (cf. flow_statistics) 4811 DO sr = 0, statistic_regions 4812 sums_l_l(k,sr,tn) = sums_l_l(k,sr,tn) + & 4813 c_0**3 * e(k,j,i) * SQRT(e(k,j,i)) / & 4814 ( diss(k,j,i) + 1.0E-30_wp ) * rmask(j,i,sr) * flag 4815 ENDDO 4418 4816 4419 4817 ENDDO … … 4434 4832 !-- so far vertical surfaces require usage of a Prandtl-layer where the boundary 4435 4833 !-- values of the diffusivities are not needed. 4436 4437 4834 IF ( .NOT. rans_tke_e ) THEN 4438 4835 ! … … 4468 4865 ENDDO 4469 4866 ENDDO 4867 ! 4868 !-- North- and southward facing surfaces 4869 DO n = 0, 1 4870 DO m = 1, surf_def_v(n)%ns 4871 i = surf_def_v(n)%i(m) 4872 j = surf_def_v(n)%j(m) 4873 k = surf_def_v(n)%k(m) 4874 km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dy 4875 kh(k,j,i) = 1.35_wp * km(k,j,i) 4876 ENDDO 4877 ENDDO 4878 ! 4879 !-- West- and eastward facing surfaces 4880 DO n = 2, 3 4881 DO m = 1, surf_def_v(n)%ns 4882 i = surf_def_v(n)%i(m) 4883 j = surf_def_v(n)%j(m) 4884 k = surf_def_v(n)%k(m) 4885 km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dx 4886 kh(k,j,i) = 1.35_wp * km(k,j,i) 4887 ENDDO 4888 ENDDO 4470 4889 4471 4890 CALL exchange_horiz( km, nbgp ) … … 4473 4892 4474 4893 ENDIF 4894 4475 4895 ! 4476 4896 !-- Model top … … 4518 4938 INTEGER(iwp) :: j !< loop index y direction 4519 4939 INTEGER(iwp) :: k !< loop index z direction 4520 INTEGER, INTENT(IN) :: mod_count !< 4940 INTEGER, INTENT(IN) :: mod_count !< flag defining where pointers point to 4521 4941 4522 4942 #if defined( __nopointer ) -
palm/trunk/SOURCE/vertical_nesting_mod.f90
r3066 r3083 4058 4058 !-- Identical timestep for coarse and fine grids 4059 4059 dt_3d = MIN( dtc, dtf ) 4060 !> @fixme setting old_dt might be obsolete at this point 4061 !> Due to changes in timestep routine, setting of old_dt might be 4062 !> not necessary any more at this point. However, could not be 4063 !> tested so far. 4064 !> 2018-05-18, gronemeier 4060 4065 old_dt = dt_3d 4061 4066 #endif
Note: See TracChangeset
for help on using the changeset viewer.
