Changeset 2508 for palm/trunk


Ignore:
Timestamp:
Oct 2, 2017 8:57:09 AM (7 years ago)
Author:
suehring
Message:

Bugfixes in SGS-TKE buoyancy production; revised initialization of vertical-gradient levels in case of ocean runs

Location:
palm/trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/INSTALL/example_cbl_rc

    r2331 r2508  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 4.0  Rev: 2330        *    atmosphere - run without 1D - prerun
     3 * PALM 4.0  Rev: 2507M       *    atmosphere - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:                 03-08-17    Run:       example_cbl         
    7  Time:                 17:10:03    Run-No.:   00
    8  Run on host:           lclocal
    9  Number of PEs:               8    Processor grid (x,y): (   2,   4) calculated
     6 Date:                 29-09-17    Run:       example_cbl         
     7 Time:                 15:33:18    Run-No.:   00
     8 Run on host:             lcmuk
     9 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
    1010 ------------------------------------------------------------------------------
    1111
     
    4141
    4242 Number of gridpoints (x,y,z):  (0:  39, 0:  39, 0:  41)
    43  Subdomain size (x,y,z):        (    20,     10,     42)
     43 Subdomain size (x,y,z):        (    20,     20,     42)
    4444
    4545
     
    123123
    124124
    125        Profile: pt, w"pt", w*pt*, wpt, w*2, pt*2,                           
     125       Profile: wpt, w"pt", w*pt*, w*2, pt, pt*2,                           
    126126       Output every               900.00 s
    127127       Time averaged over         600.00 s
     
    221221---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    222222  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
    223   0      1 00:00:20.00  20.0000A  -0.2125  -0.2273   0.1172  0.009  0.69  -0.131E+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
    224   0      2 00:00:40.00  20.0000A  -0.2117  -0.2252   0.1162  0.009  0.69  -0.132E+02    100.  0.102E-02  0.102E-02  0.939E-05  0.413E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    225   0      3 00:01:00.00  20.0000A  -0.2106  -0.2227   0.1152  0.009  0.69  -0.132E+02    100.  0.101E-02  0.101E-02  0.131E-04  0.417E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    226   0      4 00:01:20.00  20.0000A  -0.2093  -0.2197   0.1140  0.009  0.69  -0.133E+02    100.  0.100E-02  0.998E-03  0.168E-04  0.420E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    227   0      5 00:01:40.00  20.0000A  -0.2077  -0.2164   0.1128  0.009  0.69  -0.131E+02    100.  0.994E-03  0.990E-03  0.204E-04  0.424E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    228   0      6 00:02:00.00  20.0000A  -0.2058  -0.2126   0.1116  0.009  0.69  -0.128E+02    100.  0.986E-03  0.982E-03  0.241E-04  0.428E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    229   0      7 00:02:20.00  20.0000A  -0.2037  -0.2086   0.1102  0.009  0.69  -0.125E+02    100.  0.978E-03  0.974E-03  0.277E-04  0.432E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    230   0      8 00:02:40.00  20.0000A   0.3348D  0.3126D  0.1636  0.011  0.69  -0.102E+02    100.  0.202E-02  0.201E-02  0.602E-03  0.856E-13   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    231   0      9 00:03:00.00  20.0000A   0.3330   0.3039   0.1620  0.011  0.69  -0.997E+01    100.  0.200E-02  0.199E-02  0.378E-04  0.839E-15   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    232   0     10 00:03:20.00  20.0000A   0.3296  -0.2989   0.1602  0.012  0.69  -0.970E+01    100.  0.198E-02  0.196E-02  0.414E-04  0.845E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    233   0     11 00:03:40.00  20.0000A   0.3249  -0.2986   0.1582  0.012  0.69  -0.934E+01    100.  0.195E-02  0.194E-02  0.450E-04  0.857E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    234   0     12 00:04:00.00  20.0000A   0.3189  -0.2989   0.1561  0.013  0.69  -0.900E+01    100.  0.193E-02  0.192E-02  0.485E-04  0.876E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    235   0     13 00:04:20.00  20.0000D   0.3118  -0.3006   0.1539  0.013  0.93  -0.846E+01    250.  0.192E-02  0.190E-02  0.521E-04  0.901E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    236   0     14 00:04:40.00  20.0000D   0.3056  -0.3013   0.1516  0.014  0.93  -0.797E+01    250.  0.190E-02  0.189E-02  0.557E-04  0.933E-15   10   3   0    8   0   7    5  36  28    0.000    0.000       0
    237   0     15 00:05:00.00  20.0000D   0.4151D -0.3964D -0.1812  0.015  0.93  -0.761E+01    250.  0.296E-02  0.294E-02  0.629E-03  0.860E-13    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    238   0     16 00:05:20.00  20.0000D   0.4183  -0.4076   0.1999  0.016  0.93  -0.711E+01    250.  0.293E-02  0.291E-02  0.657E-04  0.144E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    239   0     17 00:05:40.00  20.0000D   0.4163  -0.4151   0.2287  0.018  0.93  -0.659E+01    250.  0.291E-02  0.289E-02  0.691E-04  0.150E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    240   0     18 00:06:00.00  20.0000D   0.4097  -0.4189   0.2650  0.019  1.22  -0.609E+01    550.  0.290E-02  0.288E-02  0.725E-04  0.157E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    241   0     19 00:06:20.00  20.0000D   0.3992  -0.4193   0.3104  0.021  1.22  -0.562E+01    550.  0.291E-02  0.289E-02  0.758E-04  0.166E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    242   0     20 00:06:40.00  20.0000D   0.3855  -0.4166   0.3665  0.023  1.22  -0.519E+01    550.  0.294E-02  0.291E-02  0.791E-04  0.178E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    243   0     21 00:07:00.00  20.0000A   0.3695  -0.4111   0.4352  0.025  1.22  -0.472E+01    550.  0.300E-02  0.297E-02  0.823E-04  0.192E-14    8  17   5    8   3   6    2   8  21    0.000    0.000       0
    244   0     22 00:07:20.00  20.0000A  -0.3702  -0.4033   0.5183  0.027  1.22  -0.433E+01    550.  0.310E-02  0.308E-02  0.853E-04  0.210E-14    8  13   9    8   3   6    2   8  21    0.000    0.000       0
    245   0     23 00:07:40.00  20.0000A  -0.4338D -0.4555D  0.6311  0.029  1.22  -0.394E+01    550.  0.431E-02  0.429E-02  0.659E-03  0.888E-13    4   4   4    8   1   8    2   8  21    0.000    0.000       0
    246   0     24 00:08:00.00  20.0000A  -0.4357  -0.4571   0.7474  0.032  1.22  -0.361E+01    550.  0.452E-02  0.450E-02  0.934E-04  0.295E-14    4   4   4    8   1   8    2   8  21    0.000    0.000       0
    247   0     25 00:08:20.00  20.0000A  -0.4388  -0.4525   0.8780  0.035  1.22  -0.332E+01    550.  0.484E-02  0.482E-02  0.956E-04  0.325E-14    4   4   4    8   1   8    2   8  21    0.000    0.000       0
    248   0     26 00:08:40.00  20.0000A  -0.4438  -0.4426   1.0186  0.038  1.22  -0.306E+01    550.  0.531E-02  0.529E-02  0.974E-04  0.361E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    249   0     27 00:09:00.00  20.0000A  -0.4667  -0.4423   1.1690  0.041  1.22  -0.283E+01    550.  0.598E-02  0.596E-02  0.987E-04  0.403E-14    1   3  23    1  29  18    2   6  12    0.000    0.000       0
    250   0     28 00:09:20.00  20.0000A  -0.5202  -0.5046   1.3508  0.044  1.18  -0.263E+01    500.  0.690E-02  0.688E-02  0.994E-04  0.453E-14    1   3  23    1  29  18    3   6  12    0.000    0.000       0
    251   0     29 00:09:40.00  20.0000A  -0.5757  -0.5716   1.6118  0.047  1.18  -0.245E+01    500.  0.814E-02  0.812E-02  0.997E-04  0.511E-14    1   4   7    1  29  18    3   6  12    0.000    0.000       0
    252   0     30 00:10:00.00  19.7000A  -0.6299D -0.6507D  1.8547  0.050  1.18  -0.228E+01    500.  0.109E-01  0.108E-01  0.651E-03  0.875E-13    1   4   7    1  29  17    3   6  12    0.000    0.000       0
    253   0     31 00:10:19.70  18.5000A  -0.6745  -0.7355   1.9910  0.054  1.18  -0.215E+01    500.  0.129E-01  0.129E-01  0.101E-03  0.680E-14    1  36  38    1  29  17    3   6  12    0.000    0.000       0
    254   0     32 00:10:38.20  17.6000A  -0.7292  -0.8133   2.2791  0.056  1.18  -0.203E+01    500.  0.153E-01  0.153E-01  0.959E-04  0.716E-14    1  36  38    1  29  17    4   6  12    0.000    0.000       0
    255   0     33 00:10:55.80  16.9000A  -0.7774  -0.8821   2.3959  0.059  1.18  -0.194E+01    500.  0.180E-01  0.180E-01  0.927E-04  0.759E-14    1  36  38    1  29  17    4   6  12    0.000    0.000       0
    256   0     34 00:11:12.70  16.5000A  -0.8188  -0.9395   2.6249  0.061  1.18  -0.186E+01    500.  0.212E-01  0.211E-01  0.911E-04  0.808E-14    1  36  38    1  29  17    5   6  12    0.000    0.000       0
    257   0     35 00:11:29.20  16.3000A  -0.8512  -0.9848   2.7031  0.064  1.18  -0.180E+01    500.  0.247E-01  0.247E-01  0.918E-04  0.870E-14    1  36  38    1  29  17    5   6  12    0.000    0.000       0
    258   0     36 00:11:45.50  15.7000A   0.8936  -1.0175   2.8625  0.066  1.18  -0.175E+01    500.  0.286E-01  0.286E-01  0.941E-04  0.945E-14    1  35  11    1  29  17    6   6  12    0.000    0.000       0
    259   0     37 00:12:01.20  15.6000A   0.9439  -1.0370   2.8835  0.067  1.28  -0.170E+01    650.  0.328E-01  0.328E-01  0.945E-04  0.995E-14    1  35  11    1  29  17    6   6  12    0.000    0.000       0
    260   0     38 00:12:16.80  15.1000A   0.9858  -1.0421   2.9773  0.069  1.28  -0.166E+01    650.  0.374E-01  0.374E-01  0.982E-04  0.108E-13    1  35  11    1  29  17    7   6  12    0.000    0.000       0
    261   0     39 00:12:31.90  14.8000A   1.0177  -1.0337   2.9369  0.071  1.28  -0.163E+01    650.  0.422E-01  0.422E-01  0.100E-03  0.113E-13    1  35  11    1  29  17    7   6  12    0.000    0.000       0
    262   0     40 00:12:46.70  14.6000A   1.0407  -1.0880   2.9904  0.072  1.29  -0.159E+01    650.  0.472E-01  0.471E-01  0.103E-03  0.119E-13    1  35  11    1  28  17    6  19  10    0.000    0.000       0
    263   0     41 00:13:01.30  14.3000A   1.0559  -1.1474   3.0119  0.073  1.29  -0.157E+01    650.  0.523E-01  0.523E-01  0.107E-03  0.127E-13    1  35  11    1  28  17    7   0  13    0.000    0.000       0
    264   0     42 00:13:15.60  14.1000A  -1.1136  -1.1912   3.1812  0.075  1.29  -0.156E+01    650.  0.576E-01  0.575E-01  0.110E-03  0.132E-13    1  24  19    1  28  17    7  19  10    0.000    0.000       0
    265   0     43 00:13:29.70  14.2000A  -1.1585  -1.2206   3.1558  0.076  1.29  -0.153E+01    650.  0.629E-01  0.629E-01  0.114E-03  0.138E-13    1  24  19    1  28  17    5  24   2    0.000    0.000       0
    266   0     44 00:13:43.90  13.6000A  -1.1885  -1.2385   3.3169  0.077  1.29  -0.151E+01    650.  0.684E-01  0.684E-01  0.120E-03  0.148E-13    1  24  19    1  28  17    8  19  10    0.000    0.000       0
    267   0     45 00:13:57.50  13.6000A  -1.2015  -1.2469   3.3181  0.077  1.29  -0.151E+01    650.  0.738E-01  0.738E-01  0.120E-03  0.150E-13    1  24  19    1  28  17    6  24   2    0.000    0.000       0
    268   0     46 00:14:11.10  13.2000A  -1.1966  -1.2485   3.3968  0.078  1.29  -0.148E+01    650.  0.792E-01  0.792E-01  0.125E-03  0.157E-13    1  24  19    1  28  17    6  24   2    0.000    0.000       0
    269   0     47 00:14:24.30  13.4000A  -1.1720  -1.2448   3.3599  0.079  1.29  -0.147E+01    650.  0.845E-01  0.845E-01  0.126E-03  0.160E-13    1  24  19    1  28  17    8  21  38    0.000    0.000       0
    270   0     48 00:14:37.70  12.5000A  -1.1307  -1.2622   3.5961  0.080  1.29  -0.147E+01    650.  0.900E-01  0.899E-01  0.133E-03  0.170E-13    1  24  19    1  28  18    7  34  12    0.000    0.000       0
    271   0     49 00:14:50.20  12.7000A  -1.0814  -1.2929   3.5398  0.080  1.38  -0.144E+01    800.  0.950E-01  0.950E-01  0.128E-03  0.165E-13    1  24  19    1  28  18    7  34  12    0.000    0.000       0
    272   0     50 00:15:02.90  12.1000A  -1.0800  -1.3177   3.7180  0.081  1.38  -0.143E+01    800.  0.100E+00  0.100E+00  0.133E-03  0.173E-13    5  34  31    1  28  18    8  34  12    0.000    0.000       0
    273   0     51 00:15:15.00  12.1000A  -1.1292  -1.3484   3.7037  0.082  1.38  -0.141E+01    800.  0.105E+00  0.105E+00  0.131E-03  0.171E-13    6  34  32    7   8  16    8  34  12    0.000    0.000       0
    274   0     52 00:15:27.10  11.9000A  -1.2222  -1.3571   3.7788  0.082  1.38  -0.140E+01    800.  0.110E+00  0.110E+00  0.134E-03  0.176E-13    6  34  31    7   7  16    9  34  12    0.000    0.000       0
    275   0     53 00:15:39.00  11.8000A  -1.3197  -1.3495   3.8225  0.083  1.38  -0.139E+01    800.  0.114E+00  0.114E+00  0.135E-03  0.178E-13    6  34  31    1  28  18    9  34  12    0.000    0.000       0
    276   0     54 00:15:50.80  11.8000A  -1.3314  -1.5600   3.8210  0.084  1.38  -0.138E+01    800.  0.119E+00  0.119E+00  0.137E-03  0.182E-13    6  34  31    8   7  16   10  34  12    0.000    0.000       0
    277   0     55 00:16:02.60  11.5000A  -1.3028  -1.5882   3.9183  0.084  1.38  -0.136E+01    800.  0.124E+00  0.124E+00  0.140E-03  0.186E-13    7  34  31    8   7  16   10  34  12    0.000    0.000       0
    278   0     56 00:16:14.10  11.5000A  -1.4626  -1.4300   3.9127  0.085  1.38  -0.135E+01    800.  0.128E+00  0.128E+00  0.139E-03  0.186E-13    7  34  31    9   7  16   10  34  12    0.000    0.000       0
    279   0     57 00:16:25.60  11.3000A  -1.5083  -1.6218   3.9839  0.086  1.38  -0.134E+01    800.  0.132E+00  0.132E+00  0.142E-03  0.190E-13    7  34  31    9   7  16   11  34  12    0.000    0.000       0
    280   0     58 00:16:36.90  11.1000A  -1.4260  -1.5809   4.0524  0.086  1.38  -0.132E+01    800.  0.137E+00  0.136E+00  0.142E-03  0.191E-13    7  34  31    9   7  16   10  34  12    0.000    0.000       0
    281   0     59 00:16:48.00  11.0000A  -1.3508  -1.3744   4.0932  0.087  1.38  -0.131E+01    800.  0.141E+00  0.141E+00  0.141E-03  0.191E-13    8  34  31    4  38  17   10  34  12    0.000    0.000       0
    282   0     60 00:16:59.00  10.8000A  -1.4557  -1.5001   4.1531  0.088  1.38  -0.130E+01    800.  0.145E+00  0.144E+00  0.143E-03  0.193E-13    8  34  31   10   7  16   11  34  12    0.000    0.000       0
    283   0     61 00:17:09.80  10.8000A  -1.4226  -1.4430   4.1810  0.088  1.38  -0.129E+01    800.  0.148E+00  0.148E+00  0.142E-03  0.193E-13    8  34  31   10   7  16   11  34  12    0.000    0.000       0
    284   0     62 00:17:20.60  10.7000A   1.3196  -1.3165   4.2106  0.089  1.38  -0.128E+01    800.  0.152E+00  0.152E+00  0.145E-03  0.197E-13    7  39   9    1  28  19   12  34  12    0.000    0.000       0
    285   0     63 00:17:31.30  10.7000A   1.3321  -1.3616   4.2218  0.090  1.38  -0.127E+01    800.  0.156E+00  0.156E+00  0.146E-03  0.199E-13    8   5  29   11   7  16   12  34  12    0.000    0.000       0
    286   0     64 00:17:42.00  10.6000A   1.4049   1.3774   4.2361  0.091  1.41  -0.126E+01    850.  0.159E+00  0.159E+00  0.149E-03  0.202E-13    7  14  30    9  14   7   13  34  12    0.000    0.000       0
    287   0     65 00:17:52.60  10.6000A   1.4742   1.3948   4.2412  0.091  1.41  -0.124E+01    850.  0.163E+00  0.162E+00  0.150E-03  0.204E-13    7  14  30    9  14   7   14  34  12    0.000    0.000       0
    288   0     66 00:18:03.20  10.6000A   1.5000  -1.4122   4.2394  0.092  1.41  -0.124E+01    850.  0.166E+00  0.166E+00  0.152E-03  0.207E-13    7  14  30    4  15  29   14  34  12    0.000    0.000       0
    289   0     67 00:18:13.79  10.6000A   1.4793  -1.5104   4.2445  0.093  1.41  -0.123E+01    850.  0.169E+00  0.169E+00  0.154E-03  0.210E-13    7  14  30    4  15  29   15  34  12    0.000    0.000       0
    290   0     68 00:18:24.39  10.7000A   1.5402  -1.5736   4.2092  0.094  1.41  -0.121E+01    850.  0.172E+00  0.172E+00  0.157E-03  0.213E-13    8  14  30    4  15  29   15  34  12    0.000    0.000       0
    291   0     69 00:18:35.09  10.8000A   1.5599  -1.5960   4.1702  0.095  1.41  -0.120E+01    850.  0.175E+00  0.174E+00  0.160E-03  0.218E-13    8  14  30    4  15  29   16  34  12    0.000    0.000       0
    292   0     70 00:18:45.89  11.0000A   1.5256  -1.5737   4.0998  0.096  1.41  -0.119E+01    850.  0.177E+00  0.177E+00  0.164E-03  0.222E-13    8  14  30    4  15  29   16  34  12    0.000    0.000       0
    293   0     71 00:18:56.89  11.3000A   1.4498  -1.6313   3.9764  0.096  1.41  -0.118E+01    850.  0.180E+00  0.179E+00  0.169E-03  0.229E-13    8  14  31    5  15  29   16  34  12    0.000    0.000       0
    294   0     72 00:19:08.19  11.7000A   1.4782  -1.7029   3.8426  0.097  1.41  -0.118E+01    850.  0.182E+00  0.182E+00  0.175E-03  0.238E-13    6  30   8    5  15  29   17  34  12    0.000    0.000       0
    295   0     73 00:19:19.89  12.3000A   1.4643  -1.7099   3.6654  0.098  1.41  -0.117E+01    850.  0.184E+00  0.184E+00  0.184E-03  0.249E-13    6  30   8    5  15  29   17  34  12    0.000    0.000       0
    296   0     74 00:19:32.19  13.0000A   1.3949   1.7055   3.4672  0.099  1.41  -0.116E+01    850.  0.186E+00  0.186E+00  0.195E-03  0.265E-13    1   2  35   12  15   7   16  34  12    0.000    0.000       0
    297   0     75 00:19:45.19  13.8000A   1.3741   1.8511   3.2689  0.100  1.41  -0.115E+01    850.  0.188E+00  0.188E+00  0.208E-03  0.284E-13    1   2  35   12  15   7   16  34  12    0.000    0.000       0
    298   0     76 00:19:58.99  14.0000A   1.3849   1.7958   3.2068  0.101  1.41  -0.113E+01    850.  0.190E+00  0.190E+00  0.224E-03  0.305E-13    1   3  35   12  14   7   16  17  29    0.000    0.000       0
    299   0     77 00:20:12.99  14.2000A   1.4296   1.9970   3.1645  0.102  1.41  -0.112E+01    850.  0.192E+00  0.191E+00  0.230E-03  0.313E-13    8   2   0   13  15   7   11  16  27    0.000    0.000       0
    300   0     78 00:20:27.19  14.3000A   1.4567   2.0398   3.1424  0.103  1.41  -0.111E+01    850.  0.193E+00  0.193E+00  0.236E-03  0.321E-13    8   2   0   13  15   7   11  16  27    0.000    0.000       0
    301   0     79 00:20:41.49  14.3000A   1.4272   1.8901   3.1549  0.104  1.38  -0.110E+01    800.  0.195E+00  0.195E+00  0.240E-03  0.328E-13    1   3  35   14  15   7   15   4  34    0.000    0.000       0
    302   0     80 00:20:55.79  14.8000A   1.4373   2.0496   3.0487  0.105  1.38  -0.108E+01    800.  0.197E+00  0.197E+00  0.244E-03  0.331E-13    1   3  35   14  15   7   11  23  15    0.000    0.000       0
    303   0     81 00:21:10.59  14.8000A   1.4433   2.0211   3.0399  0.105  1.38  -0.106E+01    800.  0.199E+00  0.199E+00  0.256E-03  0.346E-13    1   3  35   14  15   7   12  24  15    0.000    0.000       0
    304   0     82 00:21:25.39  14.5000A   1.4448   1.9147   3.1009  0.106  1.38  -0.105E+01    800.  0.202E+00  0.201E+00  0.259E-03  0.350E-13    1   3  35   15  16   7   13  24  15    0.000    0.000       0
    305   0     83 00:21:39.89  14.5000A   1.4802   2.0135   3.1045  0.107  1.38  -0.105E+01    800.  0.204E+00  0.203E+00  0.256E-03  0.347E-13    5  31  29   15  16   7   13  24  15    0.000    0.000       0
    306   0     84 00:21:54.39  14.4000A  -1.5293   1.9422   3.1286  0.107  1.38  -0.105E+01    800.  0.206E+00  0.205E+00  0.259E-03  0.351E-13    3  17  27   15  15   7   14  24  15    0.000    0.000       0
    307   0     85 00:22:08.79  14.9000A  -1.5794   1.9553   3.0169  0.108  1.38  -0.106E+01    800.  0.208E+00  0.207E+00  0.260E-03  0.351E-13    3  17  27   14  15   8   14  24  15    0.000    0.000       0
    308   0     86 00:22:23.69  15.1000A  -1.6061   1.8542   2.9895  0.109  1.38  -0.105E+01    800.  0.210E+00  0.209E+00  0.271E-03  0.366E-13    3  17  27   14  15   8    7   4   2    0.000    0.000       0
    309   0     87 00:22:38.79  14.7000A  -1.6133   1.7752   3.0537  0.109  1.38  -0.105E+01    800.  0.212E+00  0.211E+00  0.276E-03  0.374E-13    3  17  27   16  16   7   13  23  14    0.000    0.000       0
    310   0     88 00:22:53.49  14.5000A  -1.6394   1.8007   3.1072  0.110  1.38  -0.104E+01    800.  0.213E+00  0.213E+00  0.270E-03  0.366E-13    1  15  30   15  16   8    8  35  36    0.000    0.000       0
    311   0     89 00:23:07.99  14.3000A  -1.6763   1.8117   3.1424  0.110  1.38  -0.103E+01    800.  0.215E+00  0.215E+00  0.268E-03  0.363E-13    1  15  30   15  16   8    9  35  36    0.000    0.000       0
    312   0     90 00:23:22.29  14.4000A  -1.6778   1.7006   3.1271  0.111  1.38  -0.102E+01    800.  0.216E+00  0.216E+00  0.265E-03  0.359E-13    1  15  30   15  16   8    9   4   2    0.000    0.000       0
    313   0     91 00:23:36.69  14.0000A   1.6610  -1.6929   3.2213  0.112  1.38  -0.102E+01    800.  0.218E+00  0.217E+00  0.268E-03  0.362E-13    1  22   0    6  17  29   10  35  35    0.000    0.000       0
    314   0     92 00:23:50.69  13.8000A   1.7339  -1.6920   3.2528  0.113  1.41  -0.102E+01    850.  0.219E+00  0.219E+00  0.261E-03  0.353E-13    1  23  39    6  17  29   10  35  35    0.000    0.000       0
    315   0     93 00:24:04.49  13.7000A   1.7643  -1.6721   3.2950  0.113  1.41  -0.101E+01    850.  0.220E+00  0.220E+00  0.257E-03  0.349E-13    1  23  39    6  17  29   11  35  35    0.000    0.000       0
    316   0     94 00:24:18.19  13.8000A   1.7372  -1.6947   3.2679  0.114  1.41  -0.100E+01    850.  0.221E+00  0.221E+00  0.256E-03  0.347E-13    1  23  39    5  10   9   11  35  35    0.000    0.000       0
    317   0     95 00:24:31.99  13.8000A   1.7527  -1.7636   3.2723  0.115  1.41  -0.990E+00    850.  0.222E+00  0.222E+00  0.259E-03  0.350E-13    1  23   0    5  10   9   12  13  27    0.000    0.000       0
    318   0     96 00:24:45.79  13.7000A   1.8501  -1.7326   3.2878  0.115  1.41  -0.981E+00    850.  0.223E+00  0.222E+00  0.259E-03  0.350E-13    1  23   0    5  10   9   12  13  27    0.000    0.000       0
    319   0     97 00:24:59.49  13.8000A   1.8722   1.8059   3.2718  0.116  1.41  -0.977E+00    850.  0.223E+00  0.223E+00  0.258E-03  0.348E-13    1  23   0    8   2  38   13  13  27    0.000    0.000       0
    320   0     98 00:25:13.29  13.7000A   1.8364  -1.8480   3.2748  0.117  1.41  -0.965E+00    850.  0.224E+00  0.223E+00  0.261E-03  0.350E-13    1  23   0    6  10   9   13  13  27    0.000    0.000       0
    321   0     99 00:25:26.99  13.8000A   1.7995  -1.8245   3.2643  0.117  1.41  -0.960E+00    850.  0.224E+00  0.224E+00  0.260E-03  0.348E-13    1  23   1    6  10   9   12  36  35    0.000    0.000       0
    322   0    100 00:25:40.79  13.7000A   1.8721  -1.7917   3.2759  0.118  1.41  -0.957E+00    850.  0.224E+00  0.224E+00  0.262E-03  0.351E-13    1  23   1    6   9   9   12  36  35    0.000    0.000       0
    323   0    101 00:25:54.49  13.8000A   1.8790   1.7446   3.2551  0.119  1.41  -0.951E+00    850.  0.225E+00  0.224E+00  0.261E-03  0.349E-13    1  23   1    8   2  36   12  36  35    0.000    0.000       0
    324   0    102 00:26:08.29  14.1000A   1.8898   1.8065   3.1976  0.119  1.38  -0.946E+00    800.  0.225E+00  0.224E+00  0.263E-03  0.353E-13    4  24  28    9   3  39   12  36  35    0.000    0.000       0
    325   0    103 00:26:22.39  14.3000A   1.8885   1.8397   3.1462  0.120  1.38  -0.941E+00    800.  0.225E+00  0.224E+00  0.270E-03  0.361E-13    4  24  28    9   3  39   13  36  35    0.000    0.000       0
    326   0    104 00:26:36.69  14.5000A   1.9044   1.8226   3.1129  0.121  1.38  -0.932E+00    800.  0.225E+00  0.224E+00  0.274E-03  0.367E-13    4  24  29    9   3  38   13  36  35    0.000    0.000       0
    327   0    105 00:26:51.19  14.8000A   1.9263  -1.7545   3.0453  0.121  1.41  -0.928E+00    850.  0.225E+00  0.224E+00  0.278E-03  0.373E-13    1  21   1    5  10  23   13  36  35    0.000    0.000       0
    328   0    106 00:27:05.99  14.3000A   1.9735  -1.7593   3.1421  0.122  1.41  -0.930E+00    850.  0.225E+00  0.224E+00  0.284E-03  0.381E-13    1  34  25    5  10  23    9  23   1    0.000    0.000       0
    329   0    107 00:27:20.29  14.8000A   2.0489   1.7608   3.0483  0.123  1.41  -0.925E+00    850.  0.225E+00  0.224E+00  0.275E-03  0.370E-13    1  34  25    9   3  37   10  23   1    0.000    0.000       0
    330   0    108 00:27:35.09  14.4000A   2.0482  -1.7444   3.1179  0.123  1.41  -0.927E+00    850.  0.225E+00  0.224E+00  0.286E-03  0.384E-13    1  34  25    6  10  23    9  23   2    0.000    0.000       0
    331   0    109 00:27:49.49  14.4000A   2.0409  -1.7506   3.1273  0.124  1.41  -0.923E+00    850.  0.225E+00  0.225E+00  0.279E-03  0.375E-13    1  20  38    7  30   6    9  23   2    0.000    0.000       0
    332   0    110 00:28:03.89  14.7000A   2.0563  -1.8035   3.0671  0.125  1.41  -0.915E+00    850.  0.226E+00  0.225E+00  0.279E-03  0.376E-13    1  20  39    7  30   6   10  23   2    0.000    0.000       0
    333   0    111 00:28:18.59  14.4000A   2.1241  -1.8110   3.1324  0.125  1.41  -0.907E+00    850.  0.226E+00  0.225E+00  0.285E-03  0.384E-13    1  20  39    7  30   6    8  23   3    0.000    0.000       0
    334   0    112 00:28:32.99  14.1000A   2.1011  -1.7862   3.2014  0.126  1.38  -0.900E+00    800.  0.227E+00  0.226E+00  0.280E-03  0.377E-13    1  20  39    6   7  37    9  23   3    0.000    0.000       0
    335   0    113 00:28:47.09  13.8000A   2.1360  -1.8121   3.2494  0.126  1.38  -0.897E+00    800.  0.227E+00  0.226E+00  0.275E-03  0.369E-13    1  20   0    6   7  37    9  23   3    0.000    0.000       0
    336   0    114 00:29:00.89  13.9000A   2.1226  -1.8107   3.2363  0.127  1.38  -0.897E+00    800.  0.228E+00  0.227E+00  0.269E-03  0.361E-13    1  20   0    7  29   6    8  23   4    0.000    0.000       0
    337   0    115 00:29:14.79  13.7000A   2.1261  -1.8302   3.2966  0.127  1.38  -0.892E+00    800.  0.229E+00  0.228E+00  0.271E-03  0.363E-13    1  20   1    6   6  37    9  23   4    0.000    0.000       0
    338   0    116 00:29:28.49  13.6000A   2.1238  -1.8826   3.3087  0.127  1.38  -0.889E+00    800.  0.229E+00  0.229E+00  0.267E-03  0.357E-13    1  20   1    6   6  37    9  23   4    0.000    0.000       0
    339   0    117 00:29:42.09  13.4000A   2.1248  -1.8874   3.3662  0.128  1.38  -0.883E+00    800.  0.230E+00  0.230E+00  0.265E-03  0.353E-13    1  24   1    6   6  37   10  23   4    0.000    0.000       0
    340   0    118 00:29:55.49  13.6000A   2.1117  -1.8510   3.3099  0.128  1.38  -0.877E+00    800.  0.231E+00  0.230E+00  0.261E-03  0.348E-13    1  23   1    6   6  37   11  23   4    0.000    0.000       0
    341   0    119 00:30:09.09  13.5000A   2.1208  -1.8201   3.3295  0.128  1.38  -0.875E+00    800.  0.232E+00  0.231E+00  0.265E-03  0.352E-13    1  19   2    6   5  37   11  23   4    0.000    0.000       0
    342   0    120 00:30:22.59  13.7000A   2.1122  -1.8738   3.2737  0.128  1.38  -0.880E+00    800.  0.233E+00  0.232E+00  0.263E-03  0.349E-13    1  19   2    5   6  37   11  22   4    0.000    0.000       0
    343   0    121 00:30:36.29  13.6000A   2.1135  -1.9104   3.3098  0.127  1.38  -0.881E+00    800.  0.234E+00  0.233E+00  0.266E-03  0.353E-13    1  22   2    5   6  37   11  22   4    0.000    0.000       0
    344   0    122 00:30:49.89  13.7000A   2.1133  -1.9062   3.2907  0.127  1.38  -0.879E+00    800.  0.235E+00  0.234E+00  0.264E-03  0.349E-13    1  21   2    5   6  37   10  34  29    0.000    0.000       0
    345   0    123 00:31:03.59  13.5000A   2.0924  -1.8659   3.3411  0.127  1.41  -0.886E+00    850.  0.236E+00  0.235E+00  0.265E-03  0.350E-13    1  21   2    5   6  37   11  34  29    0.000    0.000       0
    346   0    124 00:31:17.09  13.7000A   2.0667  -1.9406   3.2893  0.127  1.41  -0.885E+00    850.  0.237E+00  0.236E+00  0.261E-03  0.344E-13    1  20   3   14   7   6   11  34  29    0.000    0.000       0
    347   0    125 00:31:30.79  14.2000A   2.0844  -1.9151   3.1657  0.127  1.43  -0.883E+00    900.  0.238E+00  0.237E+00  0.265E-03  0.347E-13    1  20   3   14   7   6   12  11   7    0.000    0.000       0
    348   0    126 00:31:44.99  14.0000A   2.0732   1.9057   3.2239  0.126  1.41  -0.886E+00    850.  0.239E+00  0.238E+00  0.274E-03  0.358E-13    1  20   3    5  32  22   11  35  30    0.000    0.000       0
    349   0    127 00:31:58.99  13.8000A   2.0331   1.9426   3.2629  0.126  1.41  -0.889E+00    850.  0.239E+00  0.239E+00  0.270E-03  0.352E-13    1  20   3    5  33  22   12  35  30    0.000    0.000       0
    350   0    128 00:32:12.79  13.7000A   1.9729   1.9580   3.2944  0.126  1.41  -0.892E+00    850.  0.240E+00  0.239E+00  0.266E-03  0.346E-13    1  20   3    5  33  22   12  35  30    0.000    0.000       0
    351   0    129 00:32:26.49  13.8000A   1.9539   1.9093   3.2598  0.125  1.41  -0.896E+00    850.  0.241E+00  0.240E+00  0.264E-03  0.342E-13    1  20   4    6  33  22   12  35  30    0.000    0.000       0
    352   0    130 00:32:40.29  14.1000A   1.9171   1.8826   3.1849  0.125  1.41  -0.901E+00    850.  0.242E+00  0.241E+00  0.267E-03  0.344E-13    1  20   4    6  34  22   12  35  30    0.000    0.000       0
    353   0    131 00:32:54.39  14.5000A   1.8873  -1.9132   3.1019  0.124  1.41  -0.906E+00    850.  0.243E+00  0.242E+00  0.273E-03  0.351E-13    1  19   4   11  29   6   13  35  30    0.000    0.000       0
    354   0    132 00:33:08.89  14.7000A   1.8593  -1.9441   3.0515  0.124  1.41  -0.911E+00    850.  0.244E+00  0.243E+00  0.281E-03  0.361E-13    1  19   4   11  29   6   15  35  31    0.000    0.000       0
    355   0    133 00:33:23.59  14.1000A   1.8374  -1.9285   3.1919  0.123  1.41  -0.919E+00    850.  0.244E+00  0.243E+00  0.285E-03  0.366E-13    1  22   5   11  29   6   13  12   6    0.000    0.000       0
    356   0    134 00:33:37.69  14.3000A   1.7985  -1.9136   3.1559  0.123  1.41  -0.924E+00    850.  0.245E+00  0.244E+00  0.275E-03  0.352E-13    5  32   1   11  28   6   14  12   6    0.000    0.000       0
    357   0    135 00:33:51.99  13.9000A   1.8252  -1.9198   3.2474  0.123  1.41  -0.923E+00    850.  0.245E+00  0.244E+00  0.279E-03  0.357E-13    5  32   2    9  32   5   14  12   6    0.000    0.000       0
    358   0    136 00:34:05.89  13.8000A   1.8535  -1.9075   3.2631  0.122  1.41  -0.927E+00    850.  0.245E+00  0.244E+00  0.272E-03  0.347E-13    5  32   2    9  32   5   15  12   6    0.000    0.000       0
    359   0    137 00:34:19.69  14.2000A   1.7996  -1.8627   3.1597  0.122  1.41  -0.929E+00    850.  0.245E+00  0.244E+00  0.272E-03  0.345E-13    5  32   2    9  31   5   15  12   6    0.000    0.000       0
    360   0    138 00:34:33.89  14.1000A   1.7124  -1.9516   3.1967  0.121  1.41  -0.931E+00    850.  0.245E+00  0.244E+00  0.281E-03  0.355E-13    1  21   6   11  28   5   16  12   6    0.000    0.000       0
    361   0    139 00:34:47.99  15.1000A  -1.7335  -1.9299   2.9803  0.121  1.43  -0.937E+00    900.  0.245E+00  0.244E+00  0.280E-03  0.352E-13    4  13  28   11  28   5   16  12   6    0.000    0.000       0
    362   0    140 00:35:03.09  15.4000A  -1.7852  -1.8970   2.9210  0.120  1.43  -0.938E+00    900.  0.245E+00  0.244E+00  0.300E-03  0.375E-13    4  13  28   10  31   5    6  27  23    0.000    0.000       0
    363   0    141 00:35:18.49  15.1000A  -1.8100  -1.9001   2.9761  0.120  1.41  -0.944E+00    850.  0.245E+00  0.244E+00  0.306E-03  0.383E-13    4  13  28   10  31   5    6  36   3    0.000    0.000       0
    364   0    142 00:35:33.59  14.6000A   1.8461  -1.8820   3.0852  0.120  1.41  -0.943E+00    850.  0.244E+00  0.243E+00  0.302E-03  0.375E-13   11  26  28   10  31   5    7  36   3    0.000    0.000       0
    365   0    143 00:35:48.19  14.6000A   1.8467  -1.8674   3.0885  0.119  1.41  -0.941E+00    850.  0.244E+00  0.243E+00  0.293E-03  0.363E-13   11  26  28   10  30   5    7  37   3    0.000    0.000       0
    366   0    144 00:36:02.79  14.4000A   1.8517  -1.8343   3.1235  0.119  1.41  -0.940E+00    850.  0.244E+00  0.243E+00  0.295E-03  0.364E-13   11  26  29   10  30   5    8  36   4    0.000    0.000       0
    367   0    145 00:36:17.19  14.6000A   1.8334   1.8697   3.0886  0.119  1.38  -0.936E+00    800.  0.244E+00  0.243E+00  0.292E-03  0.360E-13   11  26  29   11  15   2    8  36   4    0.000    0.000       0
    368   0    146 00:36:31.79  14.2000A   1.7882   1.8464   3.1627  0.119  1.38  -0.937E+00    800.  0.244E+00  0.243E+00  0.297E-03  0.365E-13   12  26  29   11  15   2    7  38   4    0.000    0.000       0
    369   0    147 00:36:45.99  14.3000A  -1.8382   1.8415   3.1561  0.118  1.38  -0.940E+00    800.  0.244E+00  0.243E+00  0.291E-03  0.356E-13    9  13   2   12  15   2    7  38   4    0.000    0.000       0
    370   0    148 00:37:00.29  14.1000A  -1.8092   1.8359   3.1965  0.118  1.38  -0.943E+00    800.  0.244E+00  0.243E+00  0.295E-03  0.359E-13    9  13   2   11  15   1    8  38   4    0.000    0.000       0
    371   0    149 00:37:14.39  14.1000A  -1.7676   1.8270   3.1950  0.118  1.38  -0.945E+00    800.  0.244E+00  0.243E+00  0.292E-03  0.355E-13   10  14   2   11  16   1    8  38   4    0.000    0.000       0
    372   0    150 00:37:28.49  13.8000A  -1.7841   1.7998   3.2689  0.118  1.38  -0.952E+00    800.  0.244E+00  0.243E+00  0.293E-03  0.355E-13   10  14   1   11  16   1    8  38   5    0.000    0.000       0
    373   0    151 00:37:42.29  13.5000A  -1.7949   1.7609   3.3327  0.118  1.35  -0.951E+00    750.  0.244E+00  0.243E+00  0.288E-03  0.348E-13   11  10   5   12  16   1    8  38   5    0.000    0.000       0
    374   0    152 00:37:55.79  13.2000A  -1.8109  -1.8123   3.4118  0.118  1.35  -0.951E+00    750.  0.245E+00  0.244E+00  0.283E-03  0.340E-13   11  10   5   11  18   6    9  38   5    0.000    0.000       0
    375   0    153 00:38:08.99  13.0000A  -1.8055  -1.8547   3.4589  0.117  1.35  -0.951E+00    750.  0.245E+00  0.244E+00  0.277E-03  0.332E-13   11   9   6   11  18   6    9  38   5    0.000    0.000       0
    376   0    154 00:38:21.99  12.9000A  -1.8636  -1.8821   3.4936  0.117  1.35  -0.948E+00    750.  0.245E+00  0.244E+00  0.274E-03  0.327E-13   11   9   6   11  18   6    9  37   5    0.000    0.000       0
    377   0    155 00:38:34.89  12.5000A  -1.8853  -1.8892   3.6025  0.117  1.38  -0.945E+00    800.  0.245E+00  0.244E+00  0.272E-03  0.324E-13   11   9   6   11  18   6    9  37   5    0.000    0.000       0
    378   0    156 00:38:47.39  12.2000A  -1.8821  -1.8719   3.6931  0.118  1.43  -0.944E+00    900.  0.245E+00  0.244E+00  0.265E-03  0.314E-13   15  36   1   11  18   6    9  37   5    0.000    0.000       0
    379   0    157 00:38:59.59  12.0000A  -1.8751  -1.8483   3.7411  0.118  1.43  -0.946E+00    900.  0.245E+00  0.244E+00  0.259E-03  0.307E-13   15  37   1    9  19   5    9  37   5    0.000    0.000       0
    380   0    158 00:39:11.59  11.8000A  -1.8629  -1.8273   3.8137  0.118  1.41  -0.947E+00    850.  0.245E+00  0.244E+00  0.255E-03  0.302E-13    7   9  26    9  19   5   10  37   5    0.000    0.000       0
    381   0    159 00:39:23.39  11.7000A  -1.8983  -1.7948   3.8512  0.118  1.41  -0.945E+00    850.  0.245E+00  0.244E+00  0.251E-03  0.297E-13    7   9  26   12  17   6   10  37   5    0.000    0.000       0
    382   0    160 00:39:35.09  11.7000A   1.9054  -1.7463   3.8374  0.118  1.41  -0.939E+00    850.  0.245E+00  0.244E+00  0.250E-03  0.294E-13   10  35  32   10  18   5   10  37   5    0.000    0.000       0
    383   0    161 00:39:46.79  11.6000A   1.9092  -1.7399   3.8678  0.118  1.41  -0.934E+00    850.  0.245E+00  0.244E+00  0.250E-03  0.294E-13   10  35  32   10  18   5   11  37   5    0.000    0.000       0
    384   0    162 00:39:58.39  11.7000A   1.9125  -1.6910   3.8576  0.118  1.41  -0.930E+00    850.  0.245E+00  0.244E+00  0.248E-03  0.292E-13    9  36  32   10  18   5   11  37   5    0.000    0.000       0
    385   0    163 00:40:10.09  11.8000A   1.9407  -1.6276   3.8074  0.119  1.41  -0.929E+00    850.  0.245E+00  0.244E+00  0.251E-03  0.294E-13    9  36  32   10  17   5   14  38   5    0.000    0.000       0
    386   0    164 00:40:21.89  11.8000A   1.9411  -1.6599   3.8090  0.119  1.41  -0.926E+00    850.  0.245E+00  0.244E+00  0.253E-03  0.297E-13    9  36  32   15  31   6   14  38   5    0.000    0.000       0
    387   0    165 00:40:33.69  12.0000A   1.9035  -1.7214   3.7627  0.119  1.41  -0.923E+00    850.  0.245E+00  0.243E+00  0.253E-03  0.297E-13    9  36  32   11  32  28   11  37   6    0.000    0.000       0
    388   0    166 00:40:45.69  11.7000A   1.8495  -1.7669   3.8369  0.119  1.41  -0.920E+00    850.  0.244E+00  0.243E+00  0.257E-03  0.302E-13   11  35  32   11  32  28   12  37   6    0.000    0.000       0
    389   0    167 00:40:57.39  11.6000A   1.8619  -1.7694   3.8693  0.120  1.41  -0.917E+00    850.  0.244E+00  0.243E+00  0.250E-03  0.294E-13   10  36  33   11  32  28   12  37   6    0.000    0.000       0
    390   0    168 00:41:08.99  11.5000A   1.8625  -1.7599   3.9156  0.120  1.43  -0.915E+00    900.  0.244E+00  0.243E+00  0.248E-03  0.291E-13   10  36  33   11  32  29   13  37   6    0.000    0.000       0
    391   0    169 00:41:20.49  11.5000A   1.8394  -1.7570   3.9212  0.120  1.43  -0.916E+00    900.  0.244E+00  0.242E+00  0.245E-03  0.288E-13   10  36  33   15  30   7   14  37   6    0.000    0.000       0
    392   0    170 00:41:31.99  11.5000A   1.8214  -1.7843   3.9015  0.120  1.43  -0.916E+00    900.  0.243E+00  0.242E+00  0.245E-03  0.288E-13   11  35  32   15  30   7   14  37   6    0.000    0.000       0
    393   0    171 00:41:43.49  11.6000A   1.8107  -1.7853   3.8732  0.121  1.43  -0.917E+00    900.  0.243E+00  0.241E+00  0.244E-03  0.287E-13   11  35  32   15  29   7   15  37   6    0.000    0.000       0
    394   0    172 00:41:55.09  11.8000A   1.7802  -1.8488   3.8285  0.121  1.43  -0.919E+00    900.  0.243E+00  0.241E+00  0.246E-03  0.289E-13    1  13   1   15  29   7   14  36   6    0.000    0.000       0
    395   0    173 00:42:06.89  11.8000A   1.7827  -1.8817   3.8277  0.121  1.43  -0.918E+00    900.  0.242E+00  0.241E+00  0.250E-03  0.293E-13    9  37  33   15  29   7   14  36   6    0.000    0.000       0
    396   0    174 00:42:18.69  11.9000A   1.7899  -1.8798   3.7684  0.121  1.43  -0.912E+00    900.  0.242E+00  0.240E+00  0.250E-03  0.292E-13    9  37  33   15  29   7   14  36   6    0.000    0.000       0
    397   0    175 00:42:30.59  12.2000A   1.7927  -1.8501   3.7036  0.121  1.46  -0.907E+00    950.  0.241E+00  0.240E+00  0.252E-03  0.294E-13    9  37  33   15  29   7   15  36   6    0.000    0.000       0
    398   0    176 00:42:42.79  12.4000A   1.8182  -1.8582   3.6160  0.122  1.43  -0.906E+00    900.  0.241E+00  0.239E+00  0.258E-03  0.301E-13   19  37  11   15  28   7   15  36   6    0.000    0.000       0
    399   0    177 00:42:55.19  12.7000A  -1.9094  -1.8648   3.5356  0.122  1.43  -0.905E+00    900.  0.241E+00  0.239E+00  0.262E-03  0.305E-13   12   7   7   15  28   7   15  36   7    0.000    0.000       0
    400   0    178 00:43:07.89  12.9000A  -1.9181  -1.8460   3.4780  0.122  1.43  -0.909E+00    900.  0.241E+00  0.239E+00  0.269E-03  0.312E-13   12   7   7   15  28   7   16  36   7    0.000    0.000       0
    401   0    179 00:43:20.79  13.2000A   1.8509  -1.8054   3.4190  0.122  1.41  -0.919E+00    850.  0.240E+00  0.239E+00  0.273E-03  0.317E-13   19  37  11   15  28   7   16  36   7    0.000    0.000       0
    402   0    180 00:43:33.99  13.3000A  -1.9409  -1.7742   3.3734  0.122  1.41  -0.922E+00    850.  0.240E+00  0.239E+00  0.280E-03  0.324E-13   13   7   6   12  34  33   10   3  19    0.000    0.000       0
    403   0    181 00:43:47.29  13.0000A  -1.9385  -1.8068   3.4501  0.122  1.38  -0.918E+00    800.  0.240E+00  0.239E+00  0.283E-03  0.328E-13   13   7   6   12  34  33   10   3  19    0.000    0.000       0
    404   0    182 00:44:00.29  12.7000A  -1.9523  -1.8310   3.5295  0.122  1.38  -0.912E+00    800.  0.241E+00  0.239E+00  0.277E-03  0.321E-13   14   7   6   12  34  33   11   3  19    0.000    0.000       0
    405   0    183 00:44:12.99  12.5000A   1.9783  -1.8452   3.5994  0.122  1.38  -0.910E+00    800.  0.241E+00  0.239E+00  0.272E-03  0.315E-13   14  34  32   12  34  33   11   3  19    0.000    0.000       0
    406   0    184 00:44:25.49  12.3000A   2.0174  -1.8479   3.6608  0.123  1.38  -0.908E+00    800.  0.241E+00  0.239E+00  0.269E-03  0.311E-13   14  34  32   12  34  33   12   3  19    0.000    0.000       0
    407   0    185 00:44:37.79  12.1000A   2.0369  -1.8378   3.7150  0.123  1.38  -0.907E+00    800.  0.241E+00  0.239E+00  0.266E-03  0.308E-13   14  34  32   12  34  33   12   3  19    0.000    0.000       0
    408   0    186 00:44:49.89  12.0000A   2.0359  -1.8412   3.7390  0.123  1.38  -0.907E+00    800.  0.241E+00  0.240E+00  0.262E-03  0.304E-13   14  34  32   12  33  33   13   3  19    0.000    0.000       0
    409   0    187 00:45:01.89  11.9000A   2.0462  -1.8555   3.7660  0.123  1.41  -0.907E+00    850.  0.242E+00  0.240E+00  0.261E-03  0.303E-13   14  35  32   12  33  33   13   3  19    0.000    0.000       0
    410   0    188 00:45:13.79  12.0000A   2.0171  -1.8590   3.7470  0.123  1.41  -0.907E+00    850.  0.242E+00  0.240E+00  0.260E-03  0.301E-13   14  35  32   12  33  33   14   3  19    0.000    0.000       0
    411   0    189 00:45:25.79  12.1000A   1.9931  -1.8505   3.7316  0.123  1.41  -0.907E+00    850.  0.242E+00  0.240E+00  0.263E-03  0.305E-13   14  34  33   12  33  33   14   3  19    0.000    0.000       0
    412   0    190 00:45:37.89  12.3000A   1.9619  -1.8306   3.6704  0.123  1.41  -0.904E+00    850.  0.242E+00  0.241E+00  0.265E-03  0.308E-13   14  34  33   12  33  33   15   3  19    0.000    0.000       0
    413   0    191 00:45:50.19  12.2000A   1.9310  -1.8022   3.6873  0.123  1.41  -0.901E+00    850.  0.243E+00  0.241E+00  0.270E-03  0.314E-13   15  34  33   12  33  33   14   2  18    0.000    0.000       0
    414   0    192 00:46:02.39  12.0000A   1.8724  -1.8250   3.7621  0.123  1.41  -0.900E+00    850.  0.243E+00  0.241E+00  0.269E-03  0.312E-13   15  34  33   16  29   5   14   2  18    0.000    0.000       0
    415   0    193 00:46:14.39  11.9000A   1.8731  -1.8366   3.7747  0.123  1.41  -0.900E+00    850.  0.243E+00  0.242E+00  0.265E-03  0.307E-13   15  35  33   16  29   5   15   2  18    0.000    0.000       0
    416   0    194 00:46:26.29  11.8000A   1.8315  -1.8382   3.8028  0.123  1.41  -0.901E+00    850.  0.244E+00  0.242E+00  0.263E-03  0.305E-13   15  35  33   16  28   5   15   2  18    0.000    0.000       0
    417   0    195 00:46:38.09  12.0000A   1.7918  -1.8795   3.7506  0.123  1.41  -0.900E+00    850.  0.244E+00  0.243E+00  0.261E-03  0.303E-13   19  34  10   16  28   5   15   2  18    0.000    0.000       0
    418   0    196 00:46:50.09  12.3000A   1.7749  -1.9053   3.6569  0.123  1.41  -0.900E+00    850.  0.245E+00  0.243E+00  0.266E-03  0.309E-13   19  34  10   16  28   5   15   2  18    0.000    0.000       0
    419   0    197 00:47:02.39  12.7000A  -1.7675  -1.9130   3.5520  0.123  1.41  -0.900E+00    850.  0.245E+00  0.244E+00  0.273E-03  0.317E-13    5  14  18   16  28   5   16   2  18    0.000    0.000       0
    420   0    198 00:47:15.09  12.9000A  -1.7553  -1.9008   3.4950  0.123  1.41  -0.899E+00    850.  0.246E+00  0.245E+00  0.282E-03  0.327E-13    5  14  18   16  28   5   15   1  18    0.000    0.000       0
    421   0    199 00:47:27.99  13.2000A   1.7921  -1.9003   3.4189  0.123  1.41  -0.899E+00    850.  0.247E+00  0.246E+00  0.286E-03  0.333E-13    3  19   7   16  27   5   15   1  18    0.000    0.000       0
    422   0    200 00:47:41.19  13.5000A   1.7846  -1.8990   3.3284  0.123  1.41  -0.900E+00    850.  0.247E+00  0.246E+00  0.293E-03  0.341E-13    3  19   7   16  27   5   16   1  18    0.000    0.000       0
    423   0    201 00:47:54.69  13.3000A   1.7871  -1.8790   3.3921  0.123  1.41  -0.902E+00    850.  0.248E+00  0.247E+00  0.301E-03  0.349E-13    7  18  26   13  32  33    9   1  11    0.000    0.000       0
    424   0    202 00:48:07.99  12.8000A   1.8224  -1.9230   3.5126  0.123  1.41  -0.904E+00    850.  0.249E+00  0.248E+00  0.297E-03  0.345E-13    7  18  26   13  31  33    9   1  11    0.000    0.000       0
    425   0    203 00:48:20.79  12.6000A   1.8448  -1.9635   3.5720  0.123  1.43  -0.907E+00    900.  0.250E+00  0.249E+00  0.287E-03  0.333E-13    7  18  26   13  31  33   10   1  11    0.000    0.000       0
    426   0    204 00:48:33.39  12.7000A   1.8596  -1.9610   3.5303  0.123  1.43  -0.917E+00    900.  0.250E+00  0.249E+00  0.283E-03  0.329E-13    7  18  26   13  31  33   10   1  11    0.000    0.000       0
    427   0    205 00:48:46.09  12.6000A   1.8777  -1.9223   3.5610  0.122  1.43  -0.926E+00    900.  0.251E+00  0.250E+00  0.286E-03  0.332E-13    1   7  10   13  31  33   11   1  12    0.000    0.000       0
    428   0    206 00:48:58.69  12.3000A   2.0056  -2.1292   3.6602  0.122  1.43  -0.911E+00    900.  0.252E+00  0.251E+00  0.285E-03  0.329E-13    1   7  10   16  39  11   11   1  12    0.000    0.000       0
    429   0    207 00:49:10.99  12.2000A   2.0492  -2.2028   3.6887  0.122  1.43  -0.911E+00    900.  0.252E+00  0.251E+00  0.279E-03  0.321E-13    1   7  10   16  39  11   12   1  12    0.000    0.000       0
    430   0    208 00:49:23.19  12.3000A   2.0151  -2.1455   3.6609  0.122  1.43  -0.912E+00    900.  0.253E+00  0.252E+00  0.277E-03  0.319E-13    1   7  10   11  16   8   12   1  12    0.000    0.000       0
    431   0    209 00:49:35.49  12.4000A   2.0876  -2.1920   3.6229  0.122  1.43  -0.914E+00    900.  0.253E+00  0.252E+00  0.280E-03  0.321E-13    1   7  11   16  39  12   12   0  11    0.000    0.000       0
    432   0    210 00:49:47.89  12.5000A   2.1261  -2.3248   3.5915  0.122  1.43  -0.917E+00    900.  0.254E+00  0.253E+00  0.283E-03  0.324E-13    1   7  11   17  38  11   13  37  28    0.000    0.000       0
    433   0    211 00:50:00.39  12.2000A   2.0887  -2.3431   3.6826  0.122  1.43  -0.919E+00    900.  0.254E+00  0.254E+00  0.287E-03  0.326E-13    1   7  11   17  38  11   13  37  28    0.000    0.000       0
    434   0    212 00:50:12.59  11.9000A   2.0070  -2.4010   3.7785  0.122  1.43  -0.925E+00    900.  0.255E+00  0.254E+00  0.281E-03  0.318E-13    1   7  11   17  37  11   14  37  28    0.000    0.000       0
    435   0    213 00:50:24.49  11.7000A   1.9207  -2.3776   3.8506  0.121  1.43  -0.923E+00    900.  0.256E+00  0.255E+00  0.275E-03  0.311E-13   15   1  14   17  37  11   14  37  28    0.000    0.000       0
    436   0    214 00:50:36.19  11.7000A  -1.8946  -2.3452   3.8577  0.121  1.43  -0.923E+00    900.  0.256E+00  0.255E+00  0.271E-03  0.306E-13   14  15   7   17  36  11   14  37  28    0.000    0.000       0
    437   0    215 00:50:47.89  11.7000A   1.8978  -2.3674   3.8622  0.121  1.41  -0.925E+00    850.  0.257E+00  0.256E+00  0.272E-03  0.306E-13   17   0  16   17  36  11   15  37  28    0.000    0.000       0
    438   0    216 00:50:59.59  11.7000A   1.9591  -2.3093   3.8353  0.121  1.41  -0.925E+00    850.  0.257E+00  0.257E+00  0.272E-03  0.307E-13   17   0  16   17  36  11   15  37  28    0.000    0.000       0
    439   0    217 00:51:11.29  11.9000A   1.9947  -2.2578   3.7742  0.121  1.41  -0.928E+00    850.  0.258E+00  0.257E+00  0.272E-03  0.308E-13   16   1  15   17  36  12   16  37  28    0.000    0.000       0
    440   0    218 00:51:23.19  12.1000A   1.9954  -2.2324   3.7174  0.121  1.41  -0.933E+00    850.  0.259E+00  0.258E+00  0.277E-03  0.313E-13   16   1  15   17  35  11   16  38  29    0.000    0.000       0
    441   0    219 00:51:35.29  12.2000A   1.9477  -2.2025   3.6739  0.121  1.41  -0.940E+00    850.  0.259E+00  0.259E+00  0.282E-03  0.319E-13   16   1  15   17  35  11   16  38  29    0.000    0.000       0
    442   0    220 00:51:47.49  12.7000A   1.8968  -2.1961   3.5448  0.121  1.41  -0.941E+00    850.  0.260E+00  0.259E+00  0.285E-03  0.322E-13   18  39  16   12  17  13   16  38  29    0.000    0.000       0
    443   0    221 00:52:00.19  12.6000A   1.9168  -2.2010   3.5678  0.120  1.41  -0.941E+00    850.  0.261E+00  0.260E+00  0.297E-03  0.335E-13   18  39  16   12  17  13    8  32  23    0.000    0.000       0
    444   0    222 00:52:12.79  12.4000A   1.9394  -2.1818   3.6289  0.120  1.43  -0.943E+00    900.  0.261E+00  0.260E+00  0.294E-03  0.332E-13   18  39  17   12  17  13   10  33  22    0.000    0.000       0
    445   0    223 00:52:25.19  11.9000A   1.9860  -2.2528   3.7676  0.120  1.43  -0.949E+00    900.  0.262E+00  0.261E+00  0.290E-03  0.327E-13   18  39  17   11  14   6   10  33  22    0.000    0.000       0
    446   0    224 00:52:37.09  11.9000A   2.0487  -2.2345   3.7951  0.120  1.43  -0.947E+00    900.  0.262E+00  0.261E+00  0.279E-03  0.314E-13   18  38  17   11  13   6   11  33  22    0.000    0.000       0
    447   0    225 00:52:48.99  11.4000A   2.0196  -2.3086   3.9566  0.120  1.43  -0.947E+00    900.  0.263E+00  0.262E+00  0.279E-03  0.314E-13   18  38  17   11  13   6   11  33  22    0.000    0.000       0
    448   0    226 00:53:00.39  11.2000A   1.9030  -2.2277   4.0313  0.120  1.43  -0.952E+00    900.  0.263E+00  0.262E+00  0.268E-03  0.301E-13   18  38  18   11  13   6   11  33  23    0.000    0.000       0
    449   0    227 00:53:11.59  10.6000A   1.9439  -2.2076   4.2396  0.120  1.43  -0.954E+00    900.  0.264E+00  0.263E+00  0.264E-03  0.297E-13   18  38  18   18  38   7   11  33  23    0.000    0.000       0
    450   0    228 00:53:22.19  10.6000A   1.8898  -2.2158   4.2257  0.120  1.43  -0.956E+00    900.  0.264E+00  0.263E+00  0.250E-03  0.282E-13   18  38  18   18  38   7   12  33  23    0.000    0.000       0
    451   0    229 00:53:32.79  10.3000A   1.8624  -2.1799   4.3623  0.120  1.43  -0.955E+00    900.  0.264E+00  0.263E+00  0.251E-03  0.283E-13    5  22   8   18  38   7   12  33  23    0.000    0.000       0
    452   0    230 00:53:43.09  10.5000A   1.8132  -2.1511   4.2981  0.120  1.43  -0.955E+00    900.  0.265E+00  0.264E+00  0.245E-03  0.276E-13   14   1  33   18  37   7   12  33  23    0.000    0.000       0
    453   0    231 00:53:53.59  10.4000A   1.8385  -2.1700   4.3442  0.120  1.43  -0.955E+00    900.  0.265E+00  0.264E+00  0.250E-03  0.282E-13   14   1  33   18  37   7   13  33  23    0.000    0.000       0
    454   0    232 00:54:03.99  10.5000A   1.8414  -2.1493   4.2731  0.120  1.40  -0.958E+00    850.  0.266E+00  0.265E+00  0.248E-03  0.280E-13   14   1  33   18  37   7   13  33  23    0.000    0.000       0
    455   0    233 00:54:14.49  10.7000A   1.8223  -2.1027   4.2171  0.120  1.40  -0.964E+00    850.  0.266E+00  0.265E+00  0.251E-03  0.283E-13   14   1  33   17  35   7   14  33  23    0.000    0.000       0
    456   0    234 00:54:25.19  10.8000A   1.8165  -2.1810   4.1756  0.120  1.40  -0.966E+00    850.  0.267E+00  0.266E+00  0.256E-03  0.289E-13    6  22   9   17  35   7   14  33  23    0.000    0.000       0
    457   0    235 00:54:35.99  10.9000A   1.8023  -2.2394   4.1198  0.120  1.40  -0.965E+00    850.  0.267E+00  0.266E+00  0.259E-03  0.292E-13    6  22  10   17  35   7   13  33  23    0.000    0.000       0
    458   0    236 00:54:46.89  11.0000A  -1.8409  -2.2713   4.1006  0.120  1.40  -0.966E+00    850.  0.268E+00  0.267E+00  0.261E-03  0.294E-13   17   8   1   17  35   7   14  33  23    0.000    0.000       0
    459   0    237 00:54:57.89  11.0000A  -1.8649  -2.2716   4.0996  0.120  1.40  -0.968E+00    850.  0.268E+00  0.268E+00  0.264E-03  0.297E-13   17   8   1   17  35   7   14  33  23    0.000    0.000       0
    460   0    238 00:55:08.89  11.1000A  -1.8709  -2.2383   4.0673  0.120  1.40  -0.972E+00    850.  0.269E+00  0.268E+00  0.264E-03  0.297E-13   17   8   1   17  35   7   14  33  23    0.000    0.000       0
    461   0    239 00:55:19.99  11.3000A  -1.8618  -2.2060   3.9913  0.120  1.40  -0.975E+00    850.  0.270E+00  0.269E+00  0.267E-03  0.300E-13   17   8   1   14  15  12   14  33  23    0.000    0.000       0
    462   0    240 00:55:31.29  11.3000A  -1.8770  -2.2286   3.9772  0.120  1.40  -0.972E+00    850.  0.270E+00  0.269E+00  0.272E-03  0.305E-13    1  14  24   14  15  12   15  33  23    0.000    0.000       0
    463   0    241 00:55:42.59  11.4000A  -1.8908  -2.2018   3.9534  0.120  1.40  -0.973E+00    850.  0.271E+00  0.270E+00  0.272E-03  0.305E-13    1  14  24   14  15  12   15  32  24    0.000    0.000       0
    464   0    242 00:55:53.99  11.3000A  -1.8989  -2.1232   3.9670  0.120  1.43  -0.981E+00    900.  0.271E+00  0.270E+00  0.275E-03  0.308E-13    1  14  24   16  34   5   15  32  24    0.000    0.000       0
    465   0    243 00:56:05.29  11.4000A  -1.9119  -2.1333   3.9427  0.120  1.43  -0.984E+00    900.  0.272E+00  0.271E+00  0.272E-03  0.305E-13    1  14  23   16  34   5   15  32  24    0.000    0.000       0
    466   0    244 00:56:16.69  11.4000A  -1.9295  -2.1183   3.9546  0.120  1.46  -0.987E+00    950.  0.272E+00  0.271E+00  0.275E-03  0.308E-13    1  14  23   15  33   5   16  32  24    0.000    0.000       0
    467   0    245 00:56:28.09  11.4000A  -1.9431  -2.1520   3.9533  0.120  1.46  -0.995E+00    950.  0.273E+00  0.272E+00  0.274E-03  0.308E-13    1  14  23   15  33   5   16  32  24    0.000    0.000       0
    468   0    246 00:56:39.49  11.5000A  -1.9630  -2.2336   3.9296  0.119  1.46  -0.992E+00    950.  0.273E+00  0.272E+00  0.274E-03  0.308E-13    1  13  23   15  32   5   16  32  24    0.000    0.000       0
    469   0    247 00:56:50.99  11.6000A   2.0059  -2.2697   3.8722  0.119  1.46  -0.992E+00    950.  0.274E+00  0.272E+00  0.276E-03  0.310E-13    4  21   4   15  32   5   16  32  24    0.000    0.000       0
    470   0    248 00:57:02.59  11.8000A   2.0182  -2.2392   3.8135  0.119  1.46  -0.994E+00    950.  0.274E+00  0.273E+00  0.279E-03  0.312E-13    4  21   4   15  32   5   17  32  24    0.000    0.000       0
    471   0    249 00:57:14.39  11.9000A   2.0614  -2.2080   3.7676  0.119  1.46  -0.995E+00    950.  0.274E+00  0.273E+00  0.284E-03  0.317E-13    4  21   5   15  31   5   17  32  24    0.000    0.000       0
    472   0    250 00:57:26.29  12.3000A   2.0752  -2.1482   3.6733  0.119  1.46  -0.999E+00    950.  0.275E+00  0.274E+00  0.286E-03  0.318E-13    4  21   5   15  31   5   17  32  24    0.000    0.000       0
    473   0    251 00:57:38.59  12.8000A   2.0545  -2.0109   3.5260  0.119  1.46  -0.100E+01    950.  0.275E+00  0.274E+00  0.295E-03  0.327E-13    4  21   5   14  31   5   17  32  24    0.000    0.000       0
    474   0    252 00:57:51.39  12.8000A   2.0072  -2.0460   3.5102  0.119  1.43  -0.101E+01    900.  0.275E+00  0.274E+00  0.307E-03  0.339E-13    4  21   5   14  30   5    9   0   5    0.000    0.000       0
    475   0    253 00:58:04.19  12.6000A  -2.0065  -2.1303   3.5677  0.119  1.43  -0.101E+01    900.  0.276E+00  0.275E+00  0.307E-03  0.337E-13    1  12  22   14  30   5   10  20  11    0.000    0.000       0
    476   0    254 00:58:16.79  12.5000A  -2.0028  -2.0665   3.6131  0.119  1.43  -0.102E+01    900.  0.277E+00  0.275E+00  0.302E-03  0.331E-13    1  12  22   14  30   5   11  20  11    0.000    0.000       0
    477   0    255 00:58:29.29  12.2000A   2.0277  -1.9795   3.6885  0.119  1.43  -0.101E+01    900.  0.277E+00  0.276E+00  0.300E-03  0.328E-13   12   8  25   14  29   5   11  20  11    0.000    0.000       0
    478   0    256 00:58:41.49  12.3000A   1.9976  -1.9238   3.6564  0.120  1.40  -0.103E+01    850.  0.278E+00  0.277E+00  0.293E-03  0.320E-13   12   8  25   14  30   4   12  20  11    0.000    0.000       0
    479   0    257 00:58:53.79  12.1000A   2.0332  -1.9483   3.7341  0.120  1.40  -0.100E+01    850.  0.279E+00  0.277E+00  0.296E-03  0.322E-13   13   8  26   14  30   4   12  20  11    0.000    0.000       0
    480   0    258 00:59:05.89  12.3000A   2.1001  -1.9169   3.6574  0.120  1.40  -0.100E+01    850.  0.279E+00  0.278E+00  0.292E-03  0.317E-13   13   8  26   14  12   6   12  20  11    0.000    0.000       0
    481   0    259 00:59:18.19  12.3000A   2.0587  -1.9158   3.6676  0.120  1.40  -0.991E+00    850.  0.280E+00  0.279E+00  0.297E-03  0.322E-13   13   8  26   14  29   4   13  20  11    0.000    0.000       0
    482   0    260 00:59:30.50  12.2000A   1.9849  -1.9554   3.6882  0.120  1.40  -0.989E+00    850.  0.281E+00  0.280E+00  0.298E-03  0.322E-13   14   8  26   16  18  13   12  21  11    0.000    0.000       0
    483   0    261 00:59:42.69  12.1000A   1.9925  -2.0021   3.7204  0.120  1.40  -0.988E+00    850.  0.282E+00  0.280E+00  0.296E-03  0.320E-13   14   8  27   16  18  13   12  21  11    0.000    0.000       0
    484   0    262 00:59:54.79  12.2000A   1.9466  -1.9601   3.7002  0.120  1.40  -0.984E+00    850.  0.283E+00  0.281E+00  0.294E-03  0.318E-13   14   8  27   16  18  13   12  21  11    0.000    0.000       0
    485   0    263 01:00:06.99  12.3000A   1.9426  -1.8653   3.6710  0.121  1.40  -0.986E+00    850.  0.284E+00  0.282E+00  0.297E-03  0.321E-13    1  14   9   17  17  13   13  21  11    0.000    0.000       0
     223  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
     224  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
     225  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
     226  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
     227  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
     228  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
     229  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
     230  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
     231  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
     232  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
     233  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
     234  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
     235  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
     236  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
     237  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
     238  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
     239  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
     240  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
     241  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
     242  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
     243  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
     244  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
     245  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
     246  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
     247  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
     248  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
     249  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
     250  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
     251  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
     252  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
     253  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
     254  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
     255  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
     256  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
     257  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
     258  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
     259  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
     260  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
     261  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
     262  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
     263  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
     264  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
     265  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
     266  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
     267  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
     268  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
     269  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
     270  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
     271  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
     272  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
     273  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
     274  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
     275  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
     276  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
     277  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
     278  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
     279  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
     280  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
     281  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
     282  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
     283  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
     284  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
     285  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
     286  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
     287  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
     288  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
     289  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
     290  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
     291  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
     292  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
     293  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
     294  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
     295  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
     296  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
     297  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
     298  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
     299  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
     300  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
     301  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
     302  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
     303  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
     304  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
     305  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
     306  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
     307  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
     308  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
     309  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
     310  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
     311  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
     312  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
     313  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
     314  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
     315  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
     316  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
     317  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
     318  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
     319  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
     320  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
     321  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
     322  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
     323  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
     324  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
     325  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
     326  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
     327  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
     328  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
     329  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
     330  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
     331  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
     332  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
     333  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
     334  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
     335  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
     336  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
     337  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
     338  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
     339  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
     340  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
     341  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
     342  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
     343  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
     344  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
     345  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
     346  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
     347  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
     348  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
     349  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
     350  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
     351  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
     352  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
     353  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
     354  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
     355  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
     356  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
     357  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
     358  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
     359  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
     360  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
     361  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
     362  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
     363  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
     364  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
     365  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
     366  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
     367  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
     368  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
     369  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
     370  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
     371  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
     372  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
     373  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
     374  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
     375  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
     376  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
     377  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
     378  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
     379  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
     380  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
     381  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
     382  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
     383  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
     384  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
     385  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
     386  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
     387  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
     388  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
     389  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
     390  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
     391  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
     392  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
     393  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
     394  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
     395  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
     396  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
     397  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
     398  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
     399  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
     400  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
     401  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
     402  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
     403  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
     404  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
     405  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
     406  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
     407  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
     408  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
     409  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
     410  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
     411  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
     412  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
     413  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
     414  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
     415  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
     416  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
     417  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
     418  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
     419  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
     420  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
     421  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
     422  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
     423  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
     424  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
     425  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
     426  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
     427  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
     428  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
     429  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
     430  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
     431  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
     432  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
     433  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
     434  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
     435  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
     436  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
     437  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
     438  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
     439  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
     440  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
     441  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
     442  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
     443  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
     444  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
     445  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
     446  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
     447  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
     448  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
     449  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
     450  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
     451  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
     452  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
     453  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
     454  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
     455  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
     456  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
     457  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
     458  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
     459  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
     460  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
     461  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
     462  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
     463  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
     464  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
     465  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
     466  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
     467  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
     468  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
     469  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
     470  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
     471  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
     472  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
     473  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
     474  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
     475  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
     476  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
  • palm/trunk/INSTALL/test_oceanml_rc

    r2492 r2508  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 4.0  Rev: 2486        *    ocean - run without 1D - prerun
     3 * PALM 4.0  Rev: 2507M       *    ocean - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:                 21-09-17    Run:       test_oceanml       
    7  Time:                 16:12:42    Run-No.:   00
    8  Run on host:             crayh
     6 Date:                 02-10-17    Run:       test_oceanml       
     7 Time:                 10:50:00    Run-No.:   00
     8 Run on host:             lcmuk
    99 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
    1010 ------------------------------------------------------------------------------
     
    9797    Characteristic levels of the initial temperature profile:
    9898
    99        Height:           0.0     0.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0  m
    100        Temperature:   300.00  300.00  300.00  300.00  300.00  300.00  300.00  300.00  300.00  300.00  300.00  K
    101        Gradient:      ------    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00  K/100m
    102        Gridpoint:          0      64      63      62      61      60      59      58      57      56      55
     99       Height:           0.0     0.0  m
     100       Temperature:   300.00  300.00  K
     101       Gradient:      ------    0.00  K/100m
     102       Gridpoint:          0      65
    103103
    104104    Characteristic levels of the initial salinity profile:
    105105
    106        Height:        0.0     0.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0    -1.0  m
    107        Salinity:    35.00   35.00   35.00   35.00   35.00   35.00   35.00   35.00   35.00   35.00   35.00  psu
    108        Gradient:   ------    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00  psu/100m
    109        Gridpoint:       0      64      63      62      61      60      59      58      57      56      55
     106       Height:        0.0     0.0  m
     107       Salinity:    35.00   35.00  psu
     108       Gradient:   ------    0.00  psu/100m
     109       Gridpoint:       0      65
    110110
    111111
     
    167167    Disturbance impulse (u,v) every :    60.00 s
    168168    Disturbance amplitude           :     0.03 m/s
    169     Lower disturbance level         :   -28.13 m (GP   42)
     169    Lower disturbance level         :   -28.12 m (GP   42)
    170170    Upper disturbance level         :    -4.38 m (GP   61)
    171171    Disturbances cease as soon as the disturbance energy exceeds 0.000 m**2/s**2
  • palm/trunk/SOURCE/check_parameters.f90

    r2422 r2508  
    2525! -----------------
    2626! $Id$
     27! Bugfix, change default value of vertical_gradient level in order to consider
     28! also ocean runs
     29!
     30! 2422 2017-09-08 08:25:41Z raasch
    2731! error message in case of missing "restart" file activation string
    2832!
     
    40314035!
    40324036!--    In case of no given gradients, choose zero gradient conditions
    4033        IF ( vertical_gradient_level(1) == -1.0_wp )  THEN
     4037       IF ( vertical_gradient_level(1) == -999999.9_wp )  THEN
    40344038          vertical_gradient_level(1) = 0.0_wp
    40354039       ENDIF
  • palm/trunk/SOURCE/modules.f90

    r2499 r2508  
    2525! -----------------
    2626! $Id$
     27! Change default value for pt/q/s/sa_vertical_gradient_level
     28!
     29! 2499 2017-09-22 16:47:58Z kanani
    2730! Default changed to fft_method = 'temperton-algorithm'
    2831!
     
    14101413    REAL(wp) ::  mask_scale(3)                                     !< collective array for mask_scale_x/y/z
    14111414    REAL(wp) ::  pt_vertical_gradient(10) = 0.0_wp                 !< namelist parameter
    1412     REAL(wp) ::  pt_vertical_gradient_level(10) = -1.0_wp          !< namelist parameter
     1415    REAL(wp) ::  pt_vertical_gradient_level(10) = -999999.9_wp   !< namelist parameter
    14131416    REAL(wp) ::  q_vertical_gradient(10) = 0.0_wp                  !< namelist parameter
    1414     REAL(wp) ::  q_vertical_gradient_level(10) = -1.0_wp           !< namelist parameter
     1417    REAL(wp) ::  q_vertical_gradient_level(10) = -999999.9_wp    !< namelist parameter
    14151418    REAL(wp) ::  s_vertical_gradient(10) = 0.0_wp                  !< namelist parameter
    1416     REAL(wp) ::  s_vertical_gradient_level(10) = -1.0_wp           !< namelist parameter
     1419    REAL(wp) ::  s_vertical_gradient_level(10) = -999999.9_wp    !< namelist parameter
    14171420    REAL(wp) ::  sa_vertical_gradient(10) = 0.0_wp                 !< namelist parameter
    1418     REAL(wp) ::  sa_vertical_gradient_level(10) = -1.0_wp          !< namelist parameter
     1421    REAL(wp) ::  sa_vertical_gradient_level(10) = -999999.9_wp   !< namelist parameter
    14191422    REAL(wp) ::  skip_time_domask(max_masks) = 9999999.9_wp        !< namelist parameter
    14201423    REAL(wp) ::  threshold(20) = 0.0_wp                            !< namelist parameter
  • palm/trunk/SOURCE/production_e.f90

    r2478 r2508  
    2525! -----------------
    2626! $Id$
     27! - Bugfix in buoyancy production term, wrong base state was set.
     28! - Consider use_single_reference_value case if humidity is used.
     29! - In case of use_top_fluxes, use correct inverse density at model top
     30! - Consider use_surface_fluxes and use_top_fluxes in ocean case 
     31!
     32! 2478 2017-09-18 13:37:24Z suehring
    2733! Bugfix, consider case where no constant-flux layer and no surfaces fluxes
    2834! are used
     
    152158                  humidity, kappa, neutral, ocean, pt_reference,               &
    153159                  rho_reference, use_single_reference_value,                   &
    154                   use_surface_fluxes, use_top_fluxes
     160                  use_surface_fluxes, use_top_fluxes, vpt_reference
    155161
    156162       USE grid_variables,                                                     &
     
    372378                DO  m = surf_s, surf_e
    373379                   k = surf_lsm_h%k(m)
    374 !
    375 !--                Please note, actually, an interpolation of u_0 and v_0
    376 !--                onto the grid center would be required. However, this
    377 !--                would require several data transfers between 2D-grid and
    378 !--                wall type. The effect of this missing interpolation is
    379 !--                negligible. (See also production_e_init).
     380
    380381                   dudz(k,j) = ( u(k+1,j,i) - surf_lsm_h%u_0(m) ) * dd2zu(k)   
    381382                   dvdz(k,j) = ( v(k+1,j,i) - surf_lsm_h%v_0(m) ) * dd2zu(k)
     
    388389                DO  m = surf_s, surf_e
    389390                   k = surf_usm_h%k(m)
    390 !
    391 !--                Please note, actually, an interpolation of u_0 and v_0
    392 !--                onto the grid center would be required. However, this
    393 !--                would require several data transfers between 2D-grid and
    394 !--                wall type. The effect of this missing interpolation is
    395 !--                negligible. (See also production_e_init).
     391
    396392                   dudz(k,j) = ( u(k+1,j,i) - surf_usm_h%u_0(m) ) * dd2zu(k)   
    397393                   dvdz(k,j) = ( v(k+1,j,i) - surf_usm_h%v_0(m) ) * dd2zu(k)
     
    405401                DO  m = surf_s, surf_e
    406402                   k = surf_def_h(1)%k(m)
    407 !
    408 !--                Please note, actually, an interpolation of u_0 and v_0
    409 !--                onto the grid center would be required. However, this
    410 !--                would require several data transfers between 2D-grid and
    411 !--                wall type. The effect of this missing interpolation is
    412 !--                negligible. (See also production_e_init).
     403
    413404                   dudz(k,j) = ( surf_def_h(1)%u_0(m) - u(k-1,j,i) ) * dd2zu(k)   
    414405                   dvdz(k,j) = ( surf_def_h(1)%v_0(m) - v(k-1,j,i) ) * dd2zu(k)
    415406
    416407                ENDDO
    417 
    418408             ENDDO
    419409
     
    493483             IF ( .NOT. humidity )  THEN
    494484
    495                 IF ( use_single_reference_value )  THEN
    496 
    497                    IF ( ocean )  THEN
    498 !
    499 !--                   So far in the ocean no special treatment of density flux
    500 !--                   in the bottom and top surface layer
    501                       DO  j = nys, nyn
    502                          DO  k = nzb+1, nzt
    503                             tend(k,j,i) = tend(k,j,i) +                        &
    504                                           kh(k,j,i) * g / rho_reference *      &
    505                                           ( prho(k+1,j,i) - prho(k-1,j,i) ) *  &
    506                                           dd2zu(k) *                           &
    507                                     MERGE( 1.0_wp, 0.0_wp,                     &
    508                                              BTEST( wall_flags_0(k,j,i), 0 )   &
    509                                           )
     485                IF ( ocean )  THEN
     486!
     487!--                So far in the ocean no special treatment of density flux
     488!--                in the bottom and top surface layer
     489                   DO  j = nys, nyn
     490                      DO  k = nzb+1, nzt
     491                         tend(k,j,i) = tend(k,j,i) +                           &
     492                                       kh(k,j,i) * g /                         &
     493                              MERGE( rho_reference, prho(k,j,i),               &
     494                                     use_single_reference_value ) *            &
     495                                       ( prho(k+1,j,i) - prho(k-1,j,i) ) *     &
     496                                       dd2zu(k) *                              &
     497                                   MERGE( 1.0_wp, 0.0_wp,                      &
     498                                          BTEST( wall_flags_0(k,j,i), 30 )     &
     499                                        )                            *         &
     500                                   MERGE( 1.0_wp, 0.0_wp,                      &
     501                                          BTEST( wall_flags_0(k,j,i), 9 )      &
     502                                        ) 
     503                      ENDDO
     504!
     505!--                   Treatment of near-surface grid points, at up- and down-
     506!--                   ward facing surfaces
     507                      IF ( use_surface_fluxes )  THEN
     508                         DO  l = 0, 1
     509                            surf_s = surf_def_h(l)%start_index(j,i)
     510                            surf_e = surf_def_h(l)%end_index(j,i)
     511                            DO  m = surf_s, surf_e
     512                               k = surf_def_h(l)%k(m)
     513                               tend(k,j,i) = tend(k,j,i) + g /                 &
     514                                         MERGE( rho_reference, prho(k,j,i),    &
     515                                                use_single_reference_value ) * &
     516                                         drho_air_zw(k-1) *                    &
     517                                         surf_def_h(l)%shf(m)
     518                            ENDDO
    510519                         ENDDO
     520
     521                      ENDIF
     522
     523                      IF ( use_top_fluxes )  THEN
     524                         surf_s = surf_def_h(2)%start_index(j,i)
     525                         surf_e = surf_def_h(2)%end_index(j,i)
     526                         DO  m = surf_s, surf_e
     527                            k = surf_def_h(2)%k(m)
     528                            tend(k,j,i) = tend(k,j,i) + g /                    &
     529                                         MERGE( rho_reference, prho(k,j,i),    &
     530                                                use_single_reference_value ) * &
     531                                         drho_air_zw(k) *                      &
     532                                         surf_def_h(2)%shf(m)
     533                         ENDDO
     534                      ENDIF
     535
     536                   ENDDO
     537
     538                ELSE
     539
     540                   DO  j = nys, nyn
     541                      DO  k = nzb+1, nzt
     542!
     543!--                      Flag 9 is used to mask top fluxes, flag 30 to mask
     544!--                      surface fluxes
     545                         tend(k,j,i) = tend(k,j,i) -                           &
     546                                       kh(k,j,i) * g /                         &
     547                                   MERGE( pt_reference, pt(k,j,i),             &
     548                                           use_single_reference_value ) *      &
     549                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) *         &
     550                                       dd2zu(k)                      *         &
     551                                   MERGE( 1.0_wp, 0.0_wp,                      &
     552                                          BTEST( wall_flags_0(k,j,i), 30 )     &
     553                                        )                            *         &
     554                                   MERGE( 1.0_wp, 0.0_wp,                      &
     555                                          BTEST( wall_flags_0(k,j,i), 9 )      &
     556                                        ) 
    511557                      ENDDO
    512558
    513                    ELSE
    514 
    515                       DO  j = nys, nyn
    516                          DO  k = nzb+1, nzt
    517 !
    518 !--                         Flag 9 is used to mask top fluxes, flag 30 to mask
    519 !--                         surface fluxes
    520                             tend(k,j,i) = tend(k,j,i) -                        &
    521                                           kh(k,j,i) * g / pt_reference  *      &
    522                                           ( pt(k+1,j,i) - pt(k-1,j,i) ) *      &
    523                                           dd2zu(k)                      *      &
    524                                       MERGE( 1.0_wp, 0.0_wp,                   &
    525                                              BTEST( wall_flags_0(k,j,i), 30 )  &
    526                                            )                            *      &
    527                                       MERGE( 1.0_wp, 0.0_wp,                   &
    528                                              BTEST( wall_flags_0(k,j,i), 9 )   &
    529                                            ) 
     559                      IF ( use_surface_fluxes )  THEN
     560!
     561!--                      Default surfaces, up- and downward-facing
     562                         DO  l = 0, 1
     563                            surf_s = surf_def_h(l)%start_index(j,i)
     564                            surf_e = surf_def_h(l)%end_index(j,i)
     565                            DO  m = surf_s, surf_e
     566                               k = surf_def_h(l)%k(m)
     567                               tend(k,j,i) = tend(k,j,i) + g /                 &
     568                                    MERGE( pt_reference, pt(k,j,i),            &
     569                                           use_single_reference_value )        &
     570                                                      * drho_air_zw(k-1)       &
     571                                                      * surf_def_h(l)%shf(m)   
     572                            ENDDO     
    530573                         ENDDO
    531 
    532                          IF ( use_surface_fluxes )  THEN
    533 !
    534 !--                         Default surfaces, up- and downward-facing
    535                             DO  l = 0, 1
    536                                surf_s = surf_def_h(l)%start_index(j,i)
    537                                surf_e = surf_def_h(l)%end_index(j,i)
    538                                DO  m = surf_s, surf_e
    539                                   k = surf_def_h(l)%k(m)
    540                                   tend(k,j,i) = tend(k,j,i) + g / pt_reference &
    541                                                          * drho_air_zw(k-1)    &
    542                                                          * surf_def_h(l)%shf(m)   
    543                                ENDDO   
    544                             ENDDO     
    545 !
    546 !--                         Natural surfaces
    547                             surf_s = surf_lsm_h%start_index(j,i)
    548                             surf_e = surf_lsm_h%end_index(j,i)
    549                             DO  m = surf_s, surf_e
    550                                k = surf_lsm_h%k(m)
    551                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    552                                                          * drho_air_zw(k-1)    &
    553                                                          * surf_lsm_h%shf(m)   
    554                             ENDDO
    555 !
    556 !--                         Urban surfaces
    557                             surf_s = surf_usm_h%start_index(j,i)
    558                             surf_e = surf_usm_h%end_index(j,i)
    559                             DO  m = surf_s, surf_e
    560                                k = surf_usm_h%k(m)
    561                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    562                                                          * drho_air_zw(k-1)    &
    563                                                          * surf_usm_h%shf(m)   
    564                             ENDDO                         
    565                          ENDIF
    566 
    567                          IF ( use_top_fluxes )  THEN
    568                             surf_s = surf_def_h(2)%start_index(j,i)
    569                             surf_e = surf_def_h(2)%end_index(j,i)
    570                             DO  m = surf_s, surf_e
    571                                k = surf_def_h(2)%k(m)
    572                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    573                                                          * drho_air_zw(k-1)    &
    574                                                          * surf_def_h(2)%shf(m)
    575                             ENDDO
    576                          ENDIF
    577                       ENDDO
    578 
    579                    ENDIF
    580 
    581                 ELSE
    582 
    583                    IF ( ocean )  THEN
    584 !
    585 !--                   So far in the ocean no special treatment of density flux
    586 !--                   in the bottom and top surface layer
    587                       DO  j = nys, nyn
    588                          DO  k = nzb+1, nzt
    589                             tend(k,j,i) = tend(k,j,i) +                        &
    590                                           kh(k,j,i) * g / prho(k,j,i) *        &
    591                                           ( prho(k+1,j,i) - prho(k-1,j,i) ) *  &
    592                                           dd2zu(k)                           * &
    593                                      MERGE( 1.0_wp, 0.0_wp,                    &
    594                                              BTEST( wall_flags_0(k,j,i), 0 )   &
    595                                           )
     574!
     575!--                      Natural surfaces
     576                         surf_s = surf_lsm_h%start_index(j,i)
     577                         surf_e = surf_lsm_h%end_index(j,i)
     578                         DO  m = surf_s, surf_e
     579                            k = surf_lsm_h%k(m)
     580                            tend(k,j,i) = tend(k,j,i) + g /                    &
     581                                    MERGE( pt_reference, pt(k,j,i),            &
     582                                           use_single_reference_value )        &
     583                                                      * drho_air_zw(k-1)       &
     584                                                      * surf_lsm_h%shf(m)   
    596585                         ENDDO
    597                       ENDDO
    598 
    599                    ELSE
    600 
    601                       DO  j = nys, nyn
    602                          DO  k = nzb+1, nzt
    603 !
    604 !--                         Flag 9 is used to mask top fluxes, flag 30 to mask
    605 !--                         surface fluxes
    606                             tend(k,j,i) = tend(k,j,i) -                        &
    607                                           kh(k,j,i) * g / pt(k,j,i) *          &
    608                                           ( pt(k+1,j,i) - pt(k-1,j,i) ) *      &
    609                                           dd2zu(k)                      *      &
    610                                       MERGE( 1.0_wp, 0.0_wp,                   &
    611                                              BTEST( wall_flags_0(k,j,i), 30 )  &
    612                                            )                            *      &
    613                                       MERGE( 1.0_wp, 0.0_wp,                   &
    614                                              BTEST( wall_flags_0(k,j,i), 9 )   &
    615                                            )
     586!
     587!--                      Urban surfaces
     588                         surf_s = surf_usm_h%start_index(j,i)
     589                         surf_e = surf_usm_h%end_index(j,i)
     590                         DO  m = surf_s, surf_e
     591                            k = surf_usm_h%k(m)
     592                            tend(k,j,i) = tend(k,j,i) + g /                    &
     593                                    MERGE( pt_reference, pt(k,j,i),            &
     594                                           use_single_reference_value )        &
     595                                                      * drho_air_zw(k-1)       &
     596                                                      * surf_usm_h%shf(m)   
     597                         ENDDO                         
     598                      ENDIF
     599
     600                      IF ( use_top_fluxes )  THEN
     601                         surf_s = surf_def_h(2)%start_index(j,i)
     602                         surf_e = surf_def_h(2)%end_index(j,i)
     603                         DO  m = surf_s, surf_e
     604                            k = surf_def_h(2)%k(m)
     605                            tend(k,j,i) = tend(k,j,i) + g /                    &
     606                                    MERGE( pt_reference, pt(k,j,i),            &
     607                                           use_single_reference_value )        &
     608                                                      * drho_air_zw(k)         &
     609                                                      * surf_def_h(2)%shf(m)
    616610                         ENDDO
    617 
    618                          IF ( use_surface_fluxes )  THEN
    619 !
    620 !--                         Default surfaces, up- and downwrd-facing
    621                             DO  l = 0, 1
    622                                surf_s = surf_def_h(l)%start_index(j,i)
    623                                surf_e = surf_def_h(l)%end_index(j,i)
    624                                DO  m = surf_s, surf_e
    625                                   k = surf_def_h(l)%k(m)
    626                                   tend(k,j,i) = tend(k,j,i) + g / pt_reference &
    627                                                          * drho_air_zw(k-1)    &
    628                                                          * surf_def_h(l)%shf(m)   
    629                                ENDDO 
    630                             ENDDO
    631 !
    632 !--                         Natural surfaces
    633                             surf_s = surf_lsm_h%start_index(j,i)
    634                             surf_e = surf_lsm_h%end_index(j,i)
    635                             DO  m = surf_s, surf_e
    636                                k = surf_lsm_h%k(m)
    637                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    638                                                          * drho_air_zw(k-1)    &
    639                                                          * surf_lsm_h%shf(m)   
    640                             ENDDO 
    641 !
    642 !--                         Urban surfaces
    643                             surf_s = surf_usm_h%start_index(j,i)
    644                             surf_e = surf_usm_h%end_index(j,i)
    645                             DO  m = surf_s, surf_e
    646                                k = surf_usm_h%k(m)
    647                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    648                                                          * drho_air_zw(k-1)    &
    649                                                          * surf_usm_h%shf(m)   
    650                             ENDDO 
    651                          ENDIF
    652 
    653                          IF ( use_top_fluxes )  THEN
    654                             surf_s = surf_def_h(2)%start_index(j,i)
    655                             surf_e = surf_def_h(2)%end_index(j,i)
    656                             DO  m = surf_s, surf_e
    657                                k = surf_def_h(2)%k(m)
    658                                tend(k,j,i) = tend(k,j,i) + g / pt_reference    &
    659                                                          * drho_air_zw(k-1)    &
    660                                                          * surf_def_h(2)%shf(m)
    661                             ENDDO
    662                          ENDIF
    663                       ENDDO
    664 
    665                    ENDIF
     611                      ENDIF
     612                   ENDDO
    666613
    667614                ENDIF
     
    679626                         k2 = 0.61_wp * pt(k,j,i)
    680627                         tend(k,j,i) = tend(k,j,i) - kh(k,j,i) *               &
    681                                          g / vpt(k,j,i) *                      &
     628                                         g /                                   &
     629                                    MERGE( vpt_reference, vpt(k,j,i),          &
     630                                           use_single_reference_value ) *      &
    682631                                         ( k1 * ( pt(k+1,j,i)-pt(k-1,j,i) ) +  &
    683632                                           k2 * ( q(k+1,j,i) - q(k-1,j,i) )    &
     
    704653                         ENDIF
    705654                         tend(k,j,i) = tend(k,j,i) - kh(k,j,i) *               &
    706                                          g / vpt(k,j,i) *                      &
     655                                         g /                                   &
     656                                    MERGE( vpt_reference, vpt(k,j,i),          &
     657                                           use_single_reference_value ) *      &
    707658                                         ( k1 * ( pt(k+1,j,i)-pt(k-1,j,i) ) +  &
    708659                                           k2 * ( q(k+1,j,i) - q(k-1,j,i) )    &
     
    718669                         k2 = 0.61_wp * pt(k,j,i)
    719670                         tend(k,j,i) = tend(k,j,i) -                           &
    720                                        kh(k,j,i) * g / vpt(k,j,i) *            &
     671                                       kh(k,j,i) * g /                         &
     672                                    MERGE( vpt_reference, vpt(k,j,i),          &
     673                                           use_single_reference_value ) *      &
    721674                                       ( k1 * ( pt(k+1,j,i)- pt(k-1,j,i) ) +   &
    722675                                         k2 * ( q(k+1,j,i) -  q(k-1,j,i) ) -   &
     
    739692                   DO  j = nys, nyn
    740693!
    741 !--                   Treat horizontal default surfaces, up- and downward-facing
     694!--                   Treat horizontal default surfaces
    742695                      DO  l = 0, 1
    743696                         surf_s = surf_def_h(l)%start_index(j,i)
     
    768721                            ENDIF
    769722
    770                             tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *       &
     723                            tend(k,j,i) = tend(k,j,i) + g /                    &
     724                                    MERGE( vpt_reference, vpt(k,j,i),          &
     725                                           use_single_reference_value ) *      &
    771726                                               ( k1 * surf_def_h(l)%shf(m) +   &
    772                                                  k2 * surf_def_h(l)%qsws(m)       &
     727                                                 k2 * surf_def_h(l)%qsws(m)    &
    773728                                               ) * drho_air_zw(k-1)
    774729                         ENDDO
     
    803758                         ENDIF
    804759
    805                          tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *          &
     760                         tend(k,j,i) = tend(k,j,i) + g /                       &
     761                                    MERGE( vpt_reference, vpt(k,j,i),          &
     762                                           use_single_reference_value ) *      &
    806763                                               ( k1 * surf_lsm_h%shf(m) +      &
    807764                                                 k2 * surf_lsm_h%qsws(m)       &
     
    837794                         ENDIF
    838795
    839                          tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *          &
     796                         tend(k,j,i) = tend(k,j,i) + g /                       &
     797                                    MERGE( vpt_reference, vpt(k,j,i),          &
     798                                           use_single_reference_value ) *      &
    840799                                               ( k1 * surf_usm_h%shf(m) +      &
    841800                                                 k2 * surf_usm_h%qsws(m)       &
     
    878837                         ENDIF
    879838
    880                          tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *          &
     839                         tend(k,j,i) = tend(k,j,i) + g /                       &
     840                                    MERGE( vpt_reference, vpt(k,j,i),          &
     841                                           use_single_reference_value ) *      &
    881842                                               ( k1 * surf_def_h(2)%shf(m) +   &
    882                                                  k2 * surf_def_h(2)%qsws(m)       &
    883                                                ) * drho_air_zw(k-1)
     843                                                 k2 * surf_def_h(2)%qsws(m)    &
     844                                               ) * drho_air_zw(k)
    884845
    885846                      ENDDO
     
    916877                  humidity, kappa, neutral, ocean, pt_reference,               &
    917878                  rho_reference, use_single_reference_value,                   &
    918                   use_surface_fluxes, use_top_fluxes
     879                  use_surface_fluxes, use_top_fluxes, vpt_reference
    919880
    920881       USE grid_variables,                                                     &
     
    10651026!--             -1.0 for right-facing wall, 1.0 for left-facing wall
    10661027                sign_dir = MERGE( 1.0_wp, -1.0_wp,                             &
    1067                                   BTEST( wall_flags_0(k,j,i-1), 0 ) ) 
     1028                                  BTEST( wall_flags_0(k,j,i-1), 0 ) )
    10681029                dvdx(k) = sign_dir * vsus / ( km_neutral + 1E-10_wp )
    10691030                dwdx(k) = sign_dir * wsus / ( km_neutral + 1E-10_wp )
     
    11291090          DO  m = surf_s, surf_e
    11301091             k = surf_lsm_h%k(m)
    1131 !
    1132 !--          Please note, actually, an interpolation of u_0 and v_0
    1133 !--          onto the grid center would be required. However, this
    1134 !--          would require several data transfers between 2D-grid and
    1135 !--          wall type. The effect of this missing interpolation is
    1136 !--          negligible. (See also production_e_init).
     1092
    11371093             dudz(k)     = ( u(k+1,j,i) - surf_lsm_h%u_0(m) ) * dd2zu(k)   
    11381094             dvdz(k)     = ( v(k+1,j,i) - surf_lsm_h%v_0(m) ) * dd2zu(k)
     
    11441100          DO  m = surf_s, surf_e
    11451101             k = surf_usm_h%k(m)
    1146 !
    1147 !--          Please note, actually, an interpolation of u_0 and v_0
    1148 !--          onto the grid center would be required. However, this
    1149 !--          would require several data transfers between 2D-grid and
    1150 !--          wall type. The effect of this missing interpolation is
    1151 !--          negligible. (See also production_e_init).
     1102
    11521103             dudz(k)     = ( u(k+1,j,i) - surf_usm_h%u_0(m) ) * dd2zu(k)   
    11531104             dvdz(k)     = ( v(k+1,j,i) - surf_usm_h%v_0(m) ) * dd2zu(k)
     
    11601111          DO  m = surf_s, surf_e
    11611112             k = surf_def_h(1)%k(m)
    1162 !
    1163 !--          Please note, actually, an interpolation of u_0 and v_0
    1164 !--          onto the grid center would be required. However, this
    1165 !--          would require several data transfers between 2D-grid and
    1166 !--          wall type. The effect of this missing interpolation is
    1167 !--          negligible. (See also production_e_init).
     1113
    11681114             dudz(k)     = ( surf_def_h(1)%u_0(m) - u(k-1,j,i) ) * dd2zu(k)   
    11691115             dvdz(k)     = ( surf_def_h(1)%v_0(m) - v(k-1,j,i) ) * dd2zu(k)
     
    12321178          IF ( .NOT. humidity )  THEN
    12331179
    1234              IF ( use_single_reference_value )  THEN
    1235 
    1236                 IF ( ocean )  THEN
    1237 !
    1238 !--                So far in the ocean no special treatment of density flux in
    1239 !--                the bottom and top surface layer
    1240                    DO  k = nzb+1, nzt
     1180             IF ( ocean )  THEN
     1181!
     1182!--             So far in the ocean no special treatment of density flux in
     1183!--             the bottom and top surface layer
     1184                DO  k = nzb+1, nzt
    12411185 
    1242                       tend(k,j,i) = tend(k,j,i) +                              &
    1243                                     kh(k,j,i) * g / rho_reference *            &
    1244                                     ( prho(k+1,j,i) - prho(k-1,j,i) ) *        &
    1245                                     dd2zu(k) *                                 &
    1246                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1247                                              BTEST( wall_flags_0(k,j,i), 0 )   &
    1248                                            )                   
     1186                   tend(k,j,i) = tend(k,j,i) +                                 &
     1187                                 kh(k,j,i) * g /                               &
     1188                                 MERGE( rho_reference, prho(k,j,i),            &
     1189                                        use_single_reference_value ) *         &
     1190                                 ( prho(k+1,j,i) - prho(k-1,j,i) ) *           &
     1191                                 dd2zu(k) *                                    &
     1192                                   MERGE( 1.0_wp, 0.0_wp,                      &
     1193                                          BTEST( wall_flags_0(k,j,i), 30 )     &
     1194                                        ) *                                    &
     1195                                   MERGE( 1.0_wp, 0.0_wp,                      &
     1196                                          BTEST( wall_flags_0(k,j,i), 9 )      &
     1197                                        )           
     1198                ENDDO
     1199
     1200                IF ( use_surface_fluxes )  THEN
     1201!
     1202!--                Default surfaces, up- and downward-facing
     1203                   DO  l = 0, 1
     1204                      surf_s = surf_def_h(l)%start_index(j,i)
     1205                      surf_e = surf_def_h(l)%end_index(j,i)
     1206                      DO  m = surf_s, surf_e
     1207                         k = surf_def_h(l)%k(m)
     1208                         tend(k,j,i) = tend(k,j,i) + g /                       &
     1209                                   MERGE( rho_reference, prho(k,j,i),          &
     1210                                          use_single_reference_value ) *       &
     1211                                   drho_air_zw(k-1) *                          &
     1212                                   surf_def_h(l)%shf(m)
     1213                      ENDDO
    12491214                   ENDDO
    12501215
    1251                 ELSE
    1252 
    1253                    DO  k = nzb+1, nzt
    1254 !
    1255 !--                   Flag 9 is used to mask top fluxes, flag 30 to mask
    1256 !--                   surface fluxes
    1257                       tend(k,j,i) = tend(k,j,i) -                              &
    1258                                     kh(k,j,i) * g / pt_reference *             &
    1259                                     ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k) * &
    1260                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1261                                              BTEST( wall_flags_0(k,j,i), 30 )  &
    1262                                            ) *                                 &
    1263                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1264                                              BTEST( wall_flags_0(k,j,i), 9 )   &
    1265                                            )
    1266 
     1216                ENDIF
     1217
     1218                IF ( use_top_fluxes )  THEN
     1219                   surf_s = surf_def_h(2)%start_index(j,i)
     1220                   surf_e = surf_def_h(2)%end_index(j,i)
     1221                   DO  m = surf_s, surf_e
     1222                      k = surf_def_h(2)%k(m)
     1223                      tend(k,j,i) = tend(k,j,i) + g /                          &
     1224                                   MERGE( rho_reference, prho(k,j,i),          &
     1225                                          use_single_reference_value ) *       &
     1226                                   drho_air_zw(k) *                            &
     1227                                   surf_def_h(2)%shf(m)
    12671228                   ENDDO
    1268 
    1269                    IF ( use_surface_fluxes )  THEN
    1270 !
    1271 !--                   Default surfaces, up- and downward-facing
    1272                       DO  l = 0, 1
    1273                          surf_s = surf_def_h(l)%start_index(j,i)
    1274                          surf_e = surf_def_h(l)%end_index(j,i)
    1275                          DO  m = surf_s, surf_e
    1276                             k = surf_def_h(l)%k(m)
    1277                             tend(k,j,i) = tend(k,j,i) + g / pt_reference *     &
    1278                                                            drho_air_zw(k-1) *  &
    1279                                                            surf_def_h(l)%shf(m)
    1280                          ENDDO
     1229                ENDIF
     1230
     1231             ELSE
     1232
     1233                DO  k = nzb+1, nzt
     1234!
     1235!--                Flag 9 is used to mask top fluxes, flag 30 to mask
     1236!--                surface fluxes
     1237                   tend(k,j,i) = tend(k,j,i) -                                 &
     1238                                 kh(k,j,i) * g /                               &
     1239                                   MERGE( pt_reference, pt(k,j,i),             &
     1240                                          use_single_reference_value ) *       &
     1241                                 ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k) *    &
     1242                                   MERGE( 1.0_wp, 0.0_wp,                      &
     1243                                          BTEST( wall_flags_0(k,j,i), 30 )     &
     1244                                        ) *                                    &
     1245                                   MERGE( 1.0_wp, 0.0_wp,                      &
     1246                                          BTEST( wall_flags_0(k,j,i), 9 )      &
     1247                                        )
     1248
     1249                ENDDO
     1250
     1251                IF ( use_surface_fluxes )  THEN
     1252!
     1253!--                Default surfaces, up- and downward-facing
     1254                   DO  l = 0, 1
     1255                      surf_s = surf_def_h(l)%start_index(j,i)
     1256                      surf_e = surf_def_h(l)%end_index(j,i)
     1257                      DO  m = surf_s, surf_e
     1258                         k = surf_def_h(l)%k(m)
     1259                         tend(k,j,i) = tend(k,j,i) + g /                       &
     1260                                   MERGE( pt_reference, pt(k,j,i),             &
     1261                                          use_single_reference_value ) *       &
     1262                                   drho_air_zw(k-1) *                          &
     1263                                   surf_def_h(l)%shf(m)
    12811264                      ENDDO
    1282 !
    1283 !--                   Natural surfaces
    1284                       surf_s = surf_lsm_h%start_index(j,i)
    1285                       surf_e = surf_lsm_h%end_index(j,i)
    1286                       DO  m = surf_s, surf_e
    1287                          k = surf_lsm_h%k(m)
    1288                          tend(k,j,i) = tend(k,j,i) + g / pt_reference *        &
    1289                                                            drho_air_zw(k-1) *  &
    1290                                                            surf_lsm_h%shf(m)
    1291                       ENDDO
    1292 !
    1293 !--                   Urban surfaces
    1294                       surf_s = surf_usm_h%start_index(j,i)
    1295                       surf_e = surf_usm_h%end_index(j,i)
    1296                       DO  m = surf_s, surf_e
    1297                          k = surf_usm_h%k(m)
    1298                          tend(k,j,i) = tend(k,j,i) + g / pt_reference *        &
    1299                                                            drho_air_zw(k-1) *  &
    1300                                                            surf_usm_h%shf(m)
    1301                       ENDDO
    1302                    ENDIF
    1303 
    1304                    IF ( use_top_fluxes )  THEN
    1305                       surf_s = surf_def_h(2)%start_index(j,i)
    1306                       surf_e = surf_def_h(2)%end_index(j,i)
    1307                       DO  m = surf_s, surf_e
    1308                          k = surf_def_h(2)%k(m)
    1309                          tend(k,j,i) = tend(k,j,i) + g / pt_reference *        &
    1310                                                            drho_air_zw(k-1) *  &
    1311                                                            surf_def_h(2)%shf(m)
    1312                       ENDDO
    1313                    ENDIF
    1314 
     1265                   ENDDO
     1266!
     1267!--                Natural surfaces
     1268                   surf_s = surf_lsm_h%start_index(j,i)
     1269                   surf_e = surf_lsm_h%end_index(j,i)
     1270                   DO  m = surf_s, surf_e
     1271                      k = surf_lsm_h%k(m)
     1272                      tend(k,j,i) = tend(k,j,i) + g /                          &
     1273                                   MERGE( pt_reference, pt(k,j,i),             &
     1274                                          use_single_reference_value ) *       &
     1275                                   drho_air_zw(k-1) *                          &
     1276                                   surf_lsm_h%shf(m)
     1277                   ENDDO
     1278!
     1279!--                Urban surfaces
     1280                   surf_s = surf_usm_h%start_index(j,i)
     1281                   surf_e = surf_usm_h%end_index(j,i)
     1282                   DO  m = surf_s, surf_e
     1283                      k = surf_usm_h%k(m)
     1284                      tend(k,j,i) = tend(k,j,i) + g /                          &
     1285                                   MERGE( pt_reference, pt(k,j,i),             &
     1286                                          use_single_reference_value ) *       &
     1287                                   drho_air_zw(k-1) *                          &
     1288                                   surf_usm_h%shf(m)
     1289                   ENDDO
    13151290                ENDIF
    13161291
    1317              ELSE
    1318 
    1319                 IF ( ocean )  THEN
    1320 !
    1321 !--                So far in the ocean no special treatment of density flux in
    1322 !--                the bottom and top surface layer
    1323                    DO  k = nzb+1, nzt
    1324                       tend(k,j,i) = tend(k,j,i) +                              &
    1325                                     kh(k,j,i) * g / prho(k,j,i)       *        &
    1326                                     ( prho(k+1,j,i) - prho(k-1,j,i) ) *        &
    1327                                     dd2zu(k) *                                 &
    1328                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1329                                              BTEST( wall_flags_0(k,j,i), 0 )   &
    1330                                            )           
     1292                IF ( use_top_fluxes )  THEN
     1293                   surf_s = surf_def_h(2)%start_index(j,i)
     1294                   surf_e = surf_def_h(2)%end_index(j,i)
     1295                   DO  m = surf_s, surf_e
     1296                      k = surf_def_h(2)%k(m)
     1297                      tend(k,j,i) = tend(k,j,i) + g /                          &
     1298                                   MERGE( pt_reference, pt(k,j,i),             &
     1299                                          use_single_reference_value ) *       &
     1300                                   drho_air_zw(k) *                            &
     1301                                   surf_def_h(2)%shf(m)
    13311302                   ENDDO
    1332 
    1333                 ELSE
    1334 
    1335                    DO  k = nzb+1, nzt
    1336 !
    1337 !--                   Flag 9 is used to mask top fluxes, flag 30 to mask
    1338 !--                   surface fluxes
    1339                       tend(k,j,i) = tend(k,j,i) -                              &
    1340                                     kh(k,j,i) * g / pt(k,j,i) *                &
    1341                                     ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k) * &
    1342                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1343                                              BTEST( wall_flags_0(k,j,i), 30 )  &
    1344                                            ) *                                 &
    1345                                       MERGE( 1.0_wp, 0.0_wp,                   &
    1346                                              BTEST( wall_flags_0(k,j,i), 9 )   &
    1347                                            )
    1348                    ENDDO
    1349 
    1350                    IF ( use_surface_fluxes )  THEN
    1351 !
    1352 !--                   Default surfaces, up- and downward-facing
    1353                       DO  l = 0, 1
    1354                          surf_s = surf_def_h(l)%start_index(j,i)
    1355                          surf_e = surf_def_h(l)%end_index(j,i)
    1356                          DO  m = surf_s, surf_e
    1357                             k = surf_def_h(l)%k(m)
    1358                             tend(k,j,i) = tend(k,j,i) + g / pt_reference       &
    1359                                                    * drho_air_zw(k-1)          &
    1360                                                    * surf_def_h(l)%shf(m)
    1361                          ENDDO 
    1362                       ENDDO
    1363 !
    1364 !--                   Natural surfaces
    1365                       surf_s = surf_lsm_h%start_index(j,i)
    1366                       surf_e = surf_lsm_h%end_index(j,i)
    1367                       DO  m = surf_s, surf_e
    1368                          k = surf_lsm_h%k(m)
    1369                          tend(k,j,i) = tend(k,j,i) + g / pt_reference          &
    1370                                                      * drho_air_zw(k-1)        &
    1371                                                      * surf_lsm_h%shf(m)
    1372                       ENDDO 
    1373 !
    1374 !--                   Urban surfaces
    1375                       surf_s = surf_usm_h%start_index(j,i)
    1376                       surf_e = surf_usm_h%end_index(j,i)
    1377                       DO  m = surf_s, surf_e
    1378                          k = surf_usm_h%k(m)
    1379                          tend(k,j,i) = tend(k,j,i) + g / pt_reference          &
    1380                                                      * drho_air_zw(k-1)        &
    1381                                                      * surf_usm_h%shf(m)
    1382                       ENDDO
    1383                    ENDIF
    1384 
    1385                    IF ( use_top_fluxes )  THEN
    1386                       surf_s = surf_def_h(2)%start_index(j,i)
    1387                       surf_e = surf_def_h(2)%end_index(j,i)
    1388                       DO  m = surf_s, surf_e
    1389                          k = surf_def_h(2)%k(m)
    1390                          tend(k,j,i) = tend(k,j,i) + g / pt_reference *        &
    1391                                                            drho_air_zw(k-1) *  &
    1392                                                            surf_def_h(2)%shf(m)
    1393                       ENDDO
    1394                    ENDIF
    1395 
    13961303                ENDIF
    13971304
     
    14061313                   k1 = 1.0_wp + 0.61_wp * q(k,j,i)
    14071314                   k2 = 0.61_wp * pt(k,j,i)
    1408                    tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / vpt(k,j,i) *    &
     1315                   tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g /                 &
     1316                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1317                                          use_single_reference_value ) *       &
    14091318                                         ( k1 * ( pt(k+1,j,i)-pt(k-1,j,i) ) +  &
    14101319                                           k2 * ( q(k+1,j,i) - q(k-1,j,i) )    &
     
    14311340                      k2 = theta * ( l_d_cp / temp * k1 - 1.0_wp )
    14321341                   ENDIF
    1433                    tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / vpt(k,j,i) *    &
     1342                   tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g /                 &
     1343                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1344                                          use_single_reference_value ) *       &
    14341345                                         ( k1 * ( pt(k+1,j,i)-pt(k-1,j,i) ) +  &
    14351346                                           k2 * ( q(k+1,j,i) - q(k-1,j,i) )    &
     
    14441355                   k1 = 1.0_wp + 0.61_wp * q(k,j,i) - ql(k,j,i)
    14451356                   k2 = 0.61_wp * pt(k,j,i)
    1446                    tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g / vpt(k,j,i) *    &
     1357                   tend(k,j,i) = tend(k,j,i) - kh(k,j,i) * g /                 &
     1358                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1359                                          use_single_reference_value ) *       &
    14471360                                     ( k1 * ( pt(k+1,j,i)-pt(k-1,j,i) ) +      &
    14481361                                       k2 * ( q(k+1,j,i) - q(k-1,j,i) ) -      &
     
    14891402                      ENDIF
    14901403
    1491                       tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *             &
     1404                      tend(k,j,i) = tend(k,j,i) + g /                          &
     1405                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1406                                          use_single_reference_value ) *       &
    14921407                                         ( k1 * surf_def_h(l)%shf(m) +         &
    14931408                                           k2 * surf_def_h(l)%qsws(m)          &
     
    15241439                   ENDIF
    15251440
    1526                    tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *                &
     1441                   tend(k,j,i) = tend(k,j,i) + g /                             &
     1442                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1443                                          use_single_reference_value ) *       &
    15271444                                            ( k1 * surf_lsm_h%shf(m) +         &
    15281445                                              k2 * surf_lsm_h%qsws(m)          &
     
    15581475                   ENDIF
    15591476
    1560                    tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *                &
     1477                   tend(k,j,i) = tend(k,j,i) + g /                             &
     1478                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1479                                          use_single_reference_value ) *       &
    15611480                                            ( k1 * surf_usm_h%shf(m) +         &
    15621481                                              k2 * surf_usm_h%qsws(m)          &
     
    15961515                   ENDIF
    15971516
    1598                    tend(k,j,i) = tend(k,j,i) + g / vpt(k,j,i) *                &
     1517                   tend(k,j,i) = tend(k,j,i) + g /                             &
     1518                                   MERGE( vpt_reference, vpt(k,j,i),           &
     1519                                          use_single_reference_value ) *       &
    15991520                               ( k1* surf_def_h(2)%shf(m) +                    &
    16001521                                 k2 * surf_def_h(2)%qsws(m)                    &
    1601                                ) * drho_air_zw(k-1)
     1522                               ) * drho_air_zw(k)
    16021523                ENDDO
    16031524
     
    16821603          ENDDO
    16831604!
    1684 !--       Default surfaces, downward-facing
     1605!--       Default surfaces, downward-facing surfaces
    16851606          !$OMP PARALLEL DO PRIVATE(i,j,k,m)
    16861607          DO  m = 1, surf_def_h(1)%ns
     
    16891610             j = surf_def_h(1)%j(m)
    16901611             k = surf_def_h(1)%k(m)
    1691 !
    1692 !--          Note, calculatione of u_0 and v_0 is not fully accurate, as u/v
    1693 !--          and km are not on the same grid. Actually, a further
    1694 !--          interpolation of km onto the u/v-grid is necessary. However, the
    1695 !--          effect of this error is negligible.
    1696 !--          In case of downward-facing surfaces, gradient is calculated
    1697 !--          between u_0 and u(k-1).
    1698              surf_def_h(1)%u_0(m) = u(k-1,j,i) - surf_def_h(1)%usws(m) *       &
    1699                                         drho_air_zw(k-1) *                     &
    1700                                         ( zu(k+1)    - zu(k-1)    )  /         &
     1612
     1613             surf_def_h(1)%u_0(m) = u(k-1,j,i) - surf_def_h(1)%usws(m) *    &
     1614                                        drho_air_zw(k-1) *                  &
     1615                                        ( zu(k+1)    - zu(k-1)    )  /      &
    17011616                                        ( km(k,j,i)  + 1.0E-20_wp ) 
    1702              surf_def_h(1)%v_0(m) = v(k-1,j,i) - surf_def_h(1)%vsws(m) *       &
    1703                                         drho_air_zw(k-1) *                     &
    1704                                         ( zu(k+1)    - zu(k-1)    )  /         &
     1617             surf_def_h(1)%v_0(m) = v(k-1,j,i) - surf_def_h(1)%vsws(m) *    &
     1618                                        drho_air_zw(k-1) *                  &
     1619                                        ( zu(k+1)    - zu(k-1)    )  /      &
    17051620                                        ( km(k,j,i)  + 1.0E-20_wp ) 
    17061621
    1707              IF ( ABS( surf_def_h(1)%u_0(m) - u(k-1,j,i) )  >                  &
    1708                   ABS( u(k+1,j,i)           - u(k-1,j,i) )                     &
     1622             IF ( ABS( surf_def_h(1)%u_0(m) - u(k-1,j,i) )  >               &
     1623                  ABS( u(k+1,j,i)           - u(k-1,j,i) )                  &
    17091624                )  surf_def_h(1)%u_0(m) = u(k+1,j,i)
    17101625
    1711              IF ( ABS( surf_def_h(1)%v_0(m) - v(k-1,j,i) )  >                  &
    1712                   ABS( v(k+1,j,i)           - v(k-1,j,i) )                     &
     1626             IF ( ABS( surf_def_h(1)%v_0(m) - v(k-1,j,i) )  >               &
     1627                  ABS( v(k+1,j,i)           - v(k-1,j,i) )                  &
    17131628                )  surf_def_h(1)%v_0(m) = v(k+1,j,i)
    17141629
  • palm/trunk/SOURCE/surface_mod.f90

    r2478 r2508  
    2525! -----------------
    2626! $Id$
     27! Minor formatting adjustment
     28!
     29! 2478 2017-09-18 13:37:24Z suehring
    2730! Bugfixes in initializing model top
    2831!
     
    795798       ALLOCATE ( surfaces%j(1:surfaces%ns)  )
    796799       ALLOCATE ( surfaces%k(1:surfaces%ns)  )
     800
     801       IF ( .NOT. constant_diffusion )  THEN   
     802          ALLOCATE ( surfaces%u_0(1:surfaces%ns) ) 
     803          ALLOCATE ( surfaces%v_0(1:surfaces%ns) )
     804       ENDIF
    797805!
    798806!--    Vertical momentum fluxes of u and v
     
    15041512!--          Prescribe latent heat flux at the top     
    15051513             IF ( humidity )  THEN
    1506              surf%qsws(num_h) = 0.0_wp
     1514                surf%qsws(num_h) = 0.0_wp
    15071515                IF ( cloud_physics  .AND.  microphysics_morrison ) THEN
    15081516                   surf%ncsws(num_h) = 0.0_wp
Note: See TracChangeset for help on using the changeset viewer.