Changeset 3083 for palm/trunk


Ignore:
Timestamp:
Jun 19, 2018 2:03:12 PM (3 years ago)
Author:
gronemeier
Message:

merge with branch rans

Location:
palm/trunk
Files:
15 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk

  • palm/trunk/INSTALL/example_cbl_rc

    r3068 r3083  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 5.0  Rev: 3063:30     *    atmosphere - run without 1D - prerun
     3 * PALM 5.0  Rev: 3082        *    atmosphere - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:                 12-06-18    Run:       example_cbl         
    7  Time:                 15:39:03    Run-No.:   00
    8  Run on host:        imuk_defau
     6 Date:                 19-06-18    Run:       example_cbl         
     7 Time:                 15:53:32    Run-No.:   00
     8 Run on host:           default
    99 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
    1010 ------------------------------------------------------------------------------
     
    230230RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
    231231---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    232   0      0 00:00:00.00  20.0000A  -0.2131D -0.2288D  0.1182  0.000  1.38   0.000E+00    800.  0.104E-02  0.103E-02  0.597E-03  0.867E-13    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    233   0      1 00:00:20.00  20.0000D  -0.2125  -0.2273   0.1172  0.009  0.69  -0.132E+02    100.  0.103E-02  0.102E-02  0.570E-05  0.411E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    234   0      2 00:00:40.00  20.0000D  -0.2117  -0.2252   0.1162  0.009  0.69  -0.133E+02    100.  0.102E-02  0.102E-02  0.935E-05  0.416E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    235   0      3 00:01:00.00  20.0000D  -0.2106  -0.2227   0.1152  0.009  0.69  -0.134E+02    100.  0.101E-02  0.101E-02  0.130E-04  0.421E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    236   0      4 00:01:20.00  20.0000D  -0.2093  -0.2197   0.1141  0.009  0.69  -0.134E+02    100.  0.100E-02  0.997E-03  0.166E-04  0.424E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    237   0      5 00:01:40.00  20.0000D  -0.2077  -0.2164   0.1129  0.009  0.93  -0.134E+02    250.  0.992E-03  0.988E-03  0.203E-04  0.427E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    238   0      6 00:02:00.00  20.0000D  -0.2058  -0.2126   0.1117  0.009  0.93  -0.133E+02    250.  0.983E-03  0.979E-03  0.239E-04  0.429E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    239   0      7 00:02:20.00  20.0000D  -0.2037  -0.2086   0.1103  0.009  1.22  -0.131E+02    550.  0.974E-03  0.970E-03  0.275E-04  0.431E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
    240   0      8 00:02:40.00  20.0000D   0.3348D  0.3126D  0.1636  0.011  1.22  -0.105E+02    550.  0.202E-02  0.201E-02  0.601E-03  0.856E-13   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    241   0      9 00:03:00.00  20.0000D   0.3330   0.3039   0.1619  0.011  1.22  -0.105E+02    550.  0.199E-02  0.198E-02  0.374E-04  0.830E-15   10   3  39   10   8  27    5  36  28    0.000    0.000       0
    242   0     10 00:03:20.00  20.0000D   0.3296  -0.2989   0.1601  0.011  1.22  -0.103E+02    550.  0.196E-02  0.195E-02  0.409E-04  0.827E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    243   0     11 00:03:40.00  20.0000D   0.3249  -0.2986   0.1580  0.011  1.22  -0.101E+02    550.  0.194E-02  0.193E-02  0.444E-04  0.828E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
    244   0     12 00:04:00.00  20.0000D   0.3189  -0.2989   0.1558  0.012  1.22  -0.975E+01    550.  0.191E-02  0.190E-02  0.478E-04  0.831E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    245   0     13 00:04:20.00  20.0000D   0.3119  -0.3005   0.1535  0.012  1.22  -0.940E+01    550.  0.189E-02  0.188E-02  0.513E-04  0.835E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
    246   0     14 00:04:40.00  20.0000D   0.3056  -0.3012   0.1512  0.013  1.22  -0.905E+01    550.  0.187E-02  0.185E-02  0.546E-04  0.842E-15   10   3   0    8   0   7    5  36  28    0.000    0.000       0
    247   0     15 00:05:00.00  20.0000D   0.4151D -0.3964D -0.1812  0.013  1.22  -0.851E+01    550.  0.292E-02  0.289E-02  0.625E-03  0.860E-13    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    248   0     16 00:05:20.00  20.0000D   0.4183  -0.4076  -0.1810  0.014  1.22  -0.811E+01    550.  0.287E-02  0.285E-02  0.640E-04  0.129E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    249   0     17 00:05:40.00  20.0000D   0.4163  -0.4151  -0.1796  0.015  1.22  -0.773E+01    550.  0.283E-02  0.280E-02  0.672E-04  0.129E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    250   0     18 00:06:00.00  20.0000D   0.4097  -0.4189  -0.1772  0.015  1.22  -0.735E+01    550.  0.279E-02  0.276E-02  0.703E-04  0.130E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
    251   0     19 00:06:20.00  20.0000D   0.3992  -0.4192   0.1802  0.016  1.22  -0.698E+01    550.  0.275E-02  0.273E-02  0.734E-04  0.131E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
    252   0     20 00:06:40.00  20.0000D   0.3855  -0.4165   0.1983  0.017  1.22  -0.656E+01    550.  0.272E-02  0.269E-02  0.765E-04  0.133E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
    253   0     21 00:07:00.00  20.0000D   0.3695  -0.4109   0.2198  0.018  1.22  -0.616E+01    550.  0.269E-02  0.267E-02  0.795E-04  0.136E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
    254   0     22 00:07:20.00  20.0000D  -0.3706  -0.4030   0.2474  0.020  1.22  -0.577E+01    550.  0.267E-02  0.265E-02  0.824E-04  0.140E-14    8  13   9    8   3   6    2   6  12    0.000    0.000       0
    255   0     23 00:07:40.00  20.0000D  -0.4267D -0.4554D  0.2850  0.021  1.22  -0.539E+01    550.  0.371E-02  0.369E-02  0.653E-03  0.886E-13    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    256   0     24 00:08:00.00  20.0000D  -0.4241  -0.4568   0.3233  0.023  1.22  -0.506E+01    550.  0.368E-02  0.366E-02  0.905E-04  0.190E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    257   0     25 00:08:20.00  20.0000D  -0.4208  -0.4522   0.3672  0.024  1.22  -0.471E+01    550.  0.367E-02  0.365E-02  0.931E-04  0.196E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    258   0     26 00:08:40.00  20.0000D  -0.4174  -0.4421   0.4169  0.026  1.22  -0.440E+01    550.  0.368E-02  0.366E-02  0.957E-04  0.204E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    259   0     27 00:09:00.00  20.0000D  -0.4142  -0.4279   0.4727  0.028  1.22  -0.410E+01    550.  0.372E-02  0.370E-02  0.981E-04  0.214E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    260   0     28 00:09:20.00  20.0000D  -0.4116  -0.4106   0.5348  0.030  1.18  -0.381E+01    500.  0.380E-02  0.378E-02  0.100E-03  0.227E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
    261   0     29 00:09:40.00  20.0000D  -0.4107  -0.4087   0.6032  0.032  1.18  -0.355E+01    500.  0.393E-02  0.391E-02  0.103E-03  0.242E-14    5   4   4    5   3  36    2   6  12    0.000    0.000       0
    262   0     30 00:10:00.00  20.0000D   0.4798D  0.4695D  0.7094  0.035  1.18  -0.330E+01    500.  0.519E-02  0.516E-02  0.656E-03  0.867E-13    5  14  20    9   7  14    2   6  12    0.000    0.000       0
    263   0     31 00:10:20.00  20.0000D   0.4589  -0.4757   0.8018  0.037  1.18  -0.308E+01    500.  0.543E-02  0.540E-02  0.108E-03  0.318E-14    5  14  20    8   2   6    3   6  12    0.000    0.000       0
    264   0     32 00:10:40.00  20.0000D   0.4450  -0.4782   0.9236  0.040  1.18  -0.288E+01    500.  0.579E-02  0.576E-02  0.110E-03  0.344E-14    6  39  30    8   2   6    3   6  12    0.000    0.000       0
    265   0     33 00:11:00.00  20.0000A   0.4349  -0.4767   1.0530  0.043  1.18  -0.266E+01    500.  0.629E-02  0.626E-02  0.111E-03  0.377E-14    6  39  30    8   2   6    3   6  12    0.000    0.000       0
    266   0     34 00:11:20.00  20.0000A  -0.4683  -0.5183   1.1853  0.046  1.18  -0.249E+01    500.  0.699E-02  0.696E-02  0.111E-03  0.416E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
    267   0     35 00:11:40.00  20.0000A  -0.5109  -0.5803   1.3145  0.049  1.18  -0.233E+01    500.  0.793E-02  0.791E-02  0.112E-03  0.464E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
    268   0     36 00:12:00.00  20.0000A  -0.5522  -0.6460   1.4350  0.052  1.18  -0.218E+01    500.  0.920E-02  0.917E-02  0.112E-03  0.522E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
    269   0     37 00:12:20.00  20.0000A   0.6173  -0.7139   1.6131  0.055  1.18  -0.205E+01    500.  0.109E-01  0.108E-01  0.112E-03  0.591E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    270   0     38 00:12:40.00  20.0000A   0.6866  -0.7795   1.7862  0.059  1.18  -0.193E+01    500.  0.130E-01  0.130E-01  0.112E-03  0.673E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    271   0     39 00:13:00.00  20.0000A   0.7574  -0.8398   1.9133  0.062  1.18  -0.183E+01    500.  0.158E-01  0.158E-01  0.112E-03  0.767E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
    272   0     40 00:13:20.00  20.0000A   0.8276  -0.8933   2.1010  0.066  1.18  -0.173E+01    500.  0.193E-01  0.193E-01  0.114E-03  0.877E-14    1  35  11    1  29  17    5   6  12    0.000    0.000       0
    273   0     41 00:13:40.00  20.0000A   0.8943  -0.9392   2.2218  0.069  1.18  -0.165E+01    500.  0.237E-01  0.236E-01  0.116E-03  0.100E-13    1  35  11    1  29  17    5   6  12    0.000    0.000       0
    274   0     42 00:14:00.00  18.9000A   0.9547  -0.9776   2.3771  0.072  1.28  -0.157E+01    650.  0.290E-01  0.290E-01  0.120E-03  0.115E-13    1  35  11    1  29  17    6   6  12    0.000    0.000       0
    275   0     43 00:14:18.89  18.4000A   1.0028  -1.0297   2.4487  0.075  1.28  -0.151E+01    650.  0.351E-01  0.351E-01  0.119E-03  0.123E-13    1  35  11    1  28  17    6   6  12    0.000    0.000       0
    276   0     44 00:14:37.29  17.7000A   1.0384  -1.1110   2.5490  0.078  1.28  -0.147E+01    650.  0.421E-01  0.421E-01  0.122E-03  0.135E-13    1  35  11    1  28  17    7   6  12    0.000    0.000       0
    277   0     45 00:14:55.00  16.2000A   1.0598  -1.1790   2.7722  0.081  1.28  -0.141E+01    650.  0.499E-01  0.498E-01  0.125E-03  0.146E-13    1  35  11    1  28  17    5  24   2    0.000    0.000       0
    278   0     46 00:15:11.20  15.4000A   1.0660  -1.2306   2.9231  0.083  1.18  -0.138E+01    500.  0.578E-01  0.578E-01  0.122E-03  0.148E-13    1  35  11    1  28  17    6  24   2    0.000    0.000       0
    279   0     47 00:15:26.60  14.3000A  -1.0741  -1.2693   3.1400  0.084  1.28  -0.135E+01    650.  0.661E-01  0.661E-01  0.123E-03  0.155E-13    1  24  19    1  28  17    6  24   2    0.000    0.000       0
    280   0     48 00:15:40.89  13.8000A   1.1031  -1.2960   3.2560  0.086  1.28  -0.132E+01    650.  0.745E-01  0.745E-01  0.122E-03  0.156E-13    1  17  38    1  28  17    6  24   2    0.000    0.000       0
    281   0     49 00:15:54.69  13.2000A   1.1508  -1.3140   3.4217  0.087  1.28  -0.131E+01    650.  0.832E-01  0.832E-01  0.125E-03  0.163E-13    1  17  38    1  28  17    7  24   2    0.000    0.000       0
    282   0     50 00:16:07.89  12.9000A   1.1879  -1.3247   3.4993  0.088  1.29  -0.130E+01    650.  0.919E-01  0.919E-01  0.126E-03  0.168E-13    1  17  38    1  28  17    7  24   2    0.000    0.000       0
    283   0     51 00:16:20.79  12.4000A   1.2158  -1.3300   3.6230  0.089  1.29  -0.129E+01    650.  0.101E+00  0.101E+00  0.130E-03  0.175E-13    1  17  38    1  28  17    8  24   2    0.000    0.000       0
    284   0     52 00:16:33.19  12.3000A   1.2346  -1.3449   3.6668  0.090  1.29  -0.128E+01    650.  0.110E+00  0.110E+00  0.131E-03  0.179E-13    1  17  38    1  28  18    8  34  13    0.000    0.000       0
    285   0     53 00:16:45.49  11.9000A   1.2460  -1.3664   3.7680  0.091  1.29  -0.126E+01    650.  0.119E+00  0.119E+00  0.136E-03  0.188E-13    1  17  38    1  28  18    9  34  13    0.000    0.000       0
    286   0     54 00:16:57.39  11.4000A   1.2899  -1.3815   3.9422  0.092  1.29  -0.126E+01    650.  0.128E+00  0.128E+00  0.138E-03  0.191E-13    5  31  10    1  28  18    9  34  13    0.000    0.000       0
    287   0     55 00:17:08.79  11.2000A   1.3175  -1.3906   4.0211  0.092  1.29  -0.124E+01    650.  0.136E+00  0.136E+00  0.138E-03  0.192E-13    5  31  10    1  28  18    9  34  13    0.000    0.000       0
    288   0     56 00:17:20.00  10.8000A  -1.3135  -1.3943   4.1804  0.093  1.38  -0.123E+01    800.  0.145E+00  0.145E+00  0.141E-03  0.198E-13    5  29  17    1  28  18   10  34  13    0.000    0.000       0
    289   0     57 00:17:30.79  10.6000A  -1.3858  -1.3935   4.2483  0.093  1.38  -0.123E+01    800.  0.153E+00  0.153E+00  0.141E-03  0.199E-13    5  29  17    1  28  18   10  34  13    0.000    0.000       0
    290   0     58 00:17:41.39  10.4000A  -1.4206  -1.3889   4.3431  0.094  1.38  -0.122E+01    800.  0.162E+00  0.162E+00  0.144E-03  0.204E-13    5  29  17    1  28  18   11  34  13    0.000    0.000       0
    291   0     59 00:17:51.79  10.2000A  -1.4342   1.3908   4.3989  0.095  1.38  -0.121E+01    800.  0.170E+00  0.170E+00  0.146E-03  0.208E-13    5  29  17    7  37  19   11  34  13    0.000    0.000       0
    292   0     60 00:18:02.00  10.1000A  -1.4413   1.4902   4.4530  0.095  1.38  -0.120E+01    800.  0.177E+00  0.177E+00  0.148E-03  0.212E-13    4  32  19    7  37  19   12  34  13    0.000    0.000       0
    293   0     61 00:18:12.09  10.0000A  -1.4694  -1.5081   4.4813  0.096  1.38  -0.120E+01    800.  0.185E+00  0.185E+00  0.151E-03  0.217E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
    294   0     62 00:18:22.09  10.0000A  -1.4854  -1.6266   4.4842  0.096  1.38  -0.119E+01    800.  0.193E+00  0.192E+00  0.154E-03  0.222E-13    4  32  19    8   7  16   13  34  13    0.000    0.000       0
    295   0     63 00:18:32.09  10.0000A  -1.4915  -1.7186   4.4785  0.097  1.38  -0.118E+01    800.  0.200E+00  0.200E+00  0.158E-03  0.229E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
    296   0     64 00:18:42.09  10.1000A  -1.5053  -1.7631   4.4346  0.097  1.38  -0.116E+01    800.  0.207E+00  0.207E+00  0.163E-03  0.235E-13    4  32  18    8   7  16   14  34  13    0.000    0.000       0
    297   0     65 00:18:52.19  10.2000A  -1.5591  -1.7454   4.4046  0.098  1.41  -0.116E+01    850.  0.214E+00  0.214E+00  0.168E-03  0.244E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
    298   0     66 00:19:02.39  10.4000A  -1.5902  -1.6634   4.3471  0.099  1.41  -0.115E+01    850.  0.221E+00  0.221E+00  0.175E-03  0.253E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
    299   0     67 00:19:12.79  10.5000A  -1.5977  -1.6328   4.2947  0.099  1.41  -0.115E+01    850.  0.228E+00  0.228E+00  0.182E-03  0.264E-13    4  32  18    9   7  16   15  35  12    0.000    0.000       0
    300   0     68 00:19:23.29  10.6000A  -1.5827   1.6967   4.2637  0.100  1.41  -0.114E+01    850.  0.235E+00  0.234E+00  0.188E-03  0.273E-13    4  32  18   12  17  28   15  35  12    0.000    0.000       0
    301   0     69 00:19:33.89  10.8000A  -1.5957  -1.7194   4.1687  0.101  1.41  -0.113E+01    850.  0.241E+00  0.241E+00  0.194E-03  0.282E-13    7  12   6    9   7  16   16  35  13    0.000    0.000       0
    302   0     70 00:19:44.69  10.7000A  -1.5970  -1.7896   4.1867  0.102  1.41  -0.113E+01    850.  0.247E+00  0.247E+00  0.202E-03  0.293E-13    7  12   6   10   7  17   16  35  13    0.000    0.000       0
    303   0     71 00:19:55.39  11.0000A   1.6043  -1.8306   4.1022  0.102  1.41  -0.112E+01    850.  0.253E+00  0.252E+00  0.204E-03  0.295E-13    7   1  37   10   7  17   12  23  18    0.000    0.000       0
    304   0     72 00:20:06.39  10.9000A   1.7066  -1.8092   4.1388  0.103  1.41  -0.111E+01    850.  0.258E+00  0.258E+00  0.214E-03  0.308E-13    7  24  30   10   7  17   12  23  18    0.000    0.000       0
    305   0     73 00:20:17.29  10.9000A   1.8396   1.7475   4.1429  0.104  1.41  -0.109E+01    850.  0.262E+00  0.262E+00  0.215E-03  0.309E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
    306   0     74 00:20:28.20  10.7000A   1.9139   1.8281   4.1963  0.105  1.41  -0.107E+01    850.  0.267E+00  0.266E+00  0.219E-03  0.313E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
    307   0     75 00:20:38.90  10.7000A   1.9324  -1.8944   4.1945  0.106  1.41  -0.106E+01    850.  0.270E+00  0.270E+00  0.218E-03  0.311E-13    7  24  30    9   4   0   13  23  17    0.000    0.000       0
    308   0     76 00:20:49.60  10.7000A   1.9390  -1.9956   4.2224  0.107  1.41  -0.105E+01    850.  0.273E+00  0.273E+00  0.220E-03  0.314E-13    7  24  31    9   4   0   13  23  17    0.000    0.000       0
    309   0     77 00:21:00.30  10.5000A   2.0200  -2.0694   4.2698  0.108  1.41  -0.104E+01    850.  0.276E+00  0.275E+00  0.223E-03  0.317E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    310   0     78 00:21:10.80  10.5000A   2.0530  -2.1072   4.3032  0.109  1.41  -0.103E+01    850.  0.278E+00  0.277E+00  0.221E-03  0.314E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    311   0     79 00:21:21.30  10.4000A   2.0495  -2.1041   4.3252  0.110  1.41  -0.102E+01    850.  0.279E+00  0.279E+00  0.223E-03  0.316E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    312   0     80 00:21:31.70  10.4000A   2.0215  -2.0580   4.3215  0.111  1.41  -0.101E+01    850.  0.281E+00  0.280E+00  0.222E-03  0.315E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
    313   0     81 00:21:42.10  10.5000A   1.9801  -1.9716   4.2997  0.112  1.41  -0.100E+01    850.  0.282E+00  0.281E+00  0.224E-03  0.316E-13    7  24  31    9   4   0   15  23  17    0.000    0.000       0
    314   0     82 00:21:52.60  10.4000A   1.9327  -1.9463   4.3286  0.113  1.41  -0.982E+00    850.  0.283E+00  0.282E+00  0.227E-03  0.321E-13    7  24  31   10   4   0   15  23  17    0.000    0.000       0
    315   0     83 00:22:03.00  10.5000A   1.8937  -2.0053   4.3023  0.114  1.41  -0.972E+00    850.  0.283E+00  0.283E+00  0.226E-03  0.320E-13   14   4   0   10   4   0   15  23  17    0.000    0.000       0
    316   0     84 00:22:13.50  10.7000A   1.9413  -2.0334   4.2011  0.115  1.41  -0.961E+00    850.  0.283E+00  0.283E+00  0.229E-03  0.324E-13   14   4   0   10   4   0   15  23  17    0.000    0.000       0
    317   0     85 00:22:24.20  10.8000A   1.9316  -2.0194   4.1554  0.116  1.41  -0.952E+00    850.  0.284E+00  0.283E+00  0.234E-03  0.332E-13   14   4   0   10   4   0   16  23  17    0.000    0.000       0
    318   0     86 00:22:35.00  11.0000A   1.8894  -1.9584   4.0813  0.117  1.38  -0.945E+00    800.  0.283E+00  0.283E+00  0.237E-03  0.336E-13   14   4   0   10   4   0   16  23  17    0.000    0.000       0
    319   0     87 00:22:46.00  11.4000A   1.8295  -1.8632   3.9326  0.118  1.38  -0.941E+00    800.  0.283E+00  0.283E+00  0.242E-03  0.344E-13   14   4   0   10   3   0   16  23  17    0.000    0.000       0
    320   0     88 00:22:57.40  12.0000A   1.8992  -1.8015   3.7404  0.119  1.38  -0.929E+00    800.  0.283E+00  0.283E+00  0.252E-03  0.357E-13   15   4   1    7  24   4   16  23  17    0.000    0.000       0
    321   0     89 00:23:09.40  12.6000A   1.9542  -1.8235   3.5820  0.119  1.38  -0.923E+00    800.  0.283E+00  0.282E+00  0.265E-03  0.376E-13   15   4   1    8  23   4   17  23  17    0.000    0.000       0
    322   0     90 00:23:22.00  12.9000A   1.9250  -1.9310   3.4822  0.120  1.38  -0.920E+00    800.  0.282E+00  0.282E+00  0.278E-03  0.395E-13   15   4   1    8  23   4   16  12  29    0.000    0.000       0
    323   0     91 00:23:34.90  13.2000A   1.8498  -2.0264   3.4174  0.121  1.38  -0.928E+00    800.  0.282E+00  0.281E+00  0.284E-03  0.404E-13   15   4   1    8  23   4   11  22  39    0.000    0.000       0
    324   0     92 00:23:48.10  13.4000A   1.7617  -2.1071   3.3530  0.121  1.38  -0.909E+00    800.  0.281E+00  0.281E+00  0.290E-03  0.412E-13   15   4   1    8  23   4   12  22  39    0.000    0.000       0
    325   0     93 00:24:01.50  13.7000A   1.7141  -2.1684   3.2965  0.122  1.35  -0.903E+00    750.  0.280E+00  0.280E+00  0.293E-03  0.418E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    326   0     94 00:24:15.20  13.9000A   1.7172  -2.2027   3.2386  0.122  1.35  -0.901E+00    750.  0.280E+00  0.279E+00  0.298E-03  0.426E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    327   0     95 00:24:29.10  14.1000A   1.7197  -2.1962   3.1874  0.123  1.35  -0.896E+00    750.  0.279E+00  0.279E+00  0.300E-03  0.431E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
    328   0     96 00:24:43.20  14.1000A   1.7148  -2.1716   3.1978  0.123  1.35  -0.891E+00    750.  0.279E+00  0.278E+00  0.303E-03  0.435E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
    329   0     97 00:24:57.30  13.8000A   1.6954  -2.1302   3.2532  0.123  1.35  -0.890E+00    750.  0.279E+00  0.278E+00  0.302E-03  0.432E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
    330   0     98 00:25:11.10  13.8000A  -1.7060  -2.1049   3.2691  0.123  1.35  -0.891E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.420E-13    8  25   1    8  23   5   11  23   1    0.000    0.000       0
    331   0     99 00:25:24.90  13.8000A   1.7235  -2.1802   3.2581  0.123  1.35  -0.892E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.416E-13   10  26  31    8  23   5   11  23   1    0.000    0.000       0
    332   0    100 00:25:38.70  13.8000A   1.7646  -2.1973   3.2551  0.124  1.41  -0.893E+00    850.  0.278E+00  0.277E+00  0.294E-03  0.413E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
    333   0    101 00:25:52.50  14.1000A   1.7850  -2.1490   3.2018  0.124  1.38  -0.899E+00    800.  0.278E+00  0.277E+00  0.293E-03  0.409E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
    334   0    102 00:26:06.60  14.1000A   1.7815  -2.0922   3.1823  0.124  1.38  -0.897E+00    800.  0.278E+00  0.277E+00  0.299E-03  0.413E-13   10  26  31    8  22   5   13  23   1    0.000    0.000       0
    335   0    103 00:26:20.70  14.4000A  -1.8030  -2.1453   3.1344  0.124  1.38  -0.895E+00    800.  0.277E+00  0.277E+00  0.298E-03  0.409E-13   16  22  38    8  22   5   12  23  19    0.000    0.000       0
    336   0    104 00:26:35.10  14.4000A   1.7946  -2.1085   3.1234  0.124  1.38  -0.894E+00    800.  0.277E+00  0.276E+00  0.303E-03  0.413E-13    1  36  26    8  22   5   12  23  19    0.000    0.000       0
    337   0    105 00:26:49.50  14.4000A   1.8514  -1.9889   3.1209  0.124  1.38  -0.894E+00    800.  0.276E+00  0.276E+00  0.301E-03  0.410E-13    1  37  26    8  22   5    9  38  26    0.000    0.000       0
    338   0    106 00:27:03.90  14.1000A   1.8658   1.8646   3.1965  0.125  1.38  -0.898E+00    800.  0.276E+00  0.275E+00  0.299E-03  0.406E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
    339   0    107 00:27:18.00  13.8000A   1.8273   1.8719   3.2544  0.125  1.38  -0.900E+00    800.  0.275E+00  0.274E+00  0.292E-03  0.396E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
    340   0    108 00:27:31.80  13.8000A   1.8485   1.8998   3.2650  0.125  1.38  -0.899E+00    800.  0.274E+00  0.274E+00  0.285E-03  0.385E-13    1  16  38    8  35  23   10  38  26    0.000    0.000       0
    341   0    109 00:27:45.60  13.9000A   1.9161   1.9299   3.2455  0.125  1.38  -0.900E+00    800.  0.274E+00  0.273E+00  0.284E-03  0.384E-13    1  16  39    8  35  23   10  38  26    0.000    0.000       0
    342   0    110 00:27:59.50  13.9000A   1.9930   1.9488   3.2423  0.126  1.38  -0.902E+00    800.  0.273E+00  0.272E+00  0.285E-03  0.387E-13    1  16  39    8  36  23   11  38  26    0.000    0.000       0
    343   0    111 00:28:13.40  13.9000A   2.0631   2.0386   3.2310  0.126  1.38  -0.897E+00    800.  0.272E+00  0.271E+00  0.285E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
    344   0    112 00:28:27.30  14.0000A   2.1131   2.0896   3.2197  0.126  1.38  -0.896E+00    800.  0.271E+00  0.270E+00  0.285E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
    345   0    113 00:28:41.30  14.0000A   2.1239   2.1077   3.2062  0.126  1.41  -0.895E+00    850.  0.270E+00  0.269E+00  0.287E-03  0.388E-13    1  16  39    9  36  23   13  37  37    0.000    0.000       0
    346   0    114 00:28:55.30  14.0000A   2.1610   2.1093   3.2184  0.127  1.41  -0.893E+00    850.  0.269E+00  0.268E+00  0.286E-03  0.388E-13    1  16   0    9  36  23   13  37  37    0.000    0.000       0
    347   0    115 00:29:09.30  13.9000A   2.1914   2.1430   3.2457  0.127  1.41  -0.889E+00    850.  0.268E+00  0.267E+00  0.286E-03  0.387E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
    348   0    116 00:29:23.20  13.9000A   2.1629   2.1318   3.2336  0.127  1.41  -0.888E+00    850.  0.266E+00  0.266E+00  0.284E-03  0.384E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
    349   0    117 00:29:37.10  14.2000A   2.1088   2.0451   3.1754  0.127  1.38  -0.891E+00    800.  0.265E+00  0.264E+00  0.285E-03  0.382E-13    1  16   1    9  37  23   14  37  37    0.000    0.000       0
    350   0    118 00:29:51.30  14.3000A   2.1133   1.9030   3.1387  0.128  1.38  -0.904E+00    800.  0.263E+00  0.262E+00  0.291E-03  0.389E-13    1  16   1    9  37  23   14   9  10    0.000    0.000       0
    351   0    119 00:30:05.60  14.4000A   2.0602  -1.9071   3.1173  0.128  1.38  -0.888E+00    800.  0.261E+00  0.260E+00  0.293E-03  0.390E-13    1  16   1   10   6  21   13   6  16    0.000    0.000       0
    352   0    120 00:30:20.00  14.1000A   2.0168  -2.0168   3.1893  0.128  1.38  -0.885E+00    800.  0.259E+00  0.258E+00  0.294E-03  0.391E-13    1  17   1   10   6  21    7  25  19    0.000    0.000       0
    353   0    121 00:30:34.10  13.8000A  -2.0339  -2.0424   3.2554  0.128  1.38  -0.882E+00    800.  0.256E+00  0.256E+00  0.288E-03  0.381E-13    9  10  18   10   6  21    8  25  19    0.000    0.000       0
    354   0    122 00:30:47.90  13.4000A   1.9253  -2.1019   3.3500  0.129  1.41  -0.884E+00    850.  0.254E+00  0.253E+00  0.281E-03  0.371E-13    1  17   1   10   5  21    8  25  19    0.000    0.000       0
    355   0    123 00:31:01.30  13.2000A  -1.9151  -2.1226   3.3987  0.129  1.41  -0.879E+00    850.  0.252E+00  0.251E+00  0.272E-03  0.358E-13   10   8  19   10   5  21    9  25  19    0.000    0.000       0
    356   0    124 00:31:14.50  13.1000A  -1.9606  -2.0641   3.4475  0.129  1.41  -0.879E+00    850.  0.249E+00  0.249E+00  0.268E-03  0.350E-13   10   9  18   10   5  21    9  25  19    0.000    0.000       0
    357   0    125 00:31:27.60  13.1000A  -1.8503  -1.9343   3.4386  0.128  1.41  -0.885E+00    850.  0.247E+00  0.246E+00  0.265E-03  0.345E-13   10   9  18   10   5  21   10  25  19    0.000    0.000       0
    358   0    126 00:31:40.70  13.1000A  -1.9153  -1.8985   3.4354  0.128  1.41  -0.884E+00    850.  0.244E+00  0.244E+00  0.264E-03  0.343E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
    359   0    127 00:31:53.80  13.4000A  -1.9319  -1.8848   3.3614  0.128  1.41  -0.890E+00    850.  0.242E+00  0.241E+00  0.263E-03  0.340E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
    360   0    128 00:32:07.20  13.7000A  -1.9610  -1.8736   3.2883  0.128  1.41  -0.894E+00    850.  0.239E+00  0.239E+00  0.268E-03  0.346E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
    361   0    129 00:32:20.90  13.9000A  -2.0531  -1.8801   3.2398  0.127  1.41  -0.899E+00    850.  0.237E+00  0.236E+00  0.272E-03  0.351E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
    362   0    130 00:32:34.80  14.1000A  -2.0925  -1.8885   3.1975  0.127  1.41  -0.901E+00    850.  0.234E+00  0.233E+00  0.275E-03  0.353E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
    363   0    131 00:32:48.90  14.2000A  -2.1181  -1.8974   3.1694  0.127  1.41  -0.903E+00    850.  0.232E+00  0.231E+00  0.277E-03  0.355E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
    364   0    132 00:33:03.10  14.6000A  -2.1477  -1.9040   3.0733  0.126  1.41  -0.903E+00    850.  0.229E+00  0.228E+00  0.278E-03  0.354E-13   14  26  19    1  13  16   10   8  19    0.000    0.000       0
    365   0    133 00:33:17.70  14.7000A  -2.2353  -1.9052   3.0618  0.126  1.41  -0.908E+00    850.  0.227E+00  0.226E+00  0.285E-03  0.360E-13   14  26  18    1  13  16   10   8  19    0.000    0.000       0
    366   0    134 00:33:32.40  14.7000A  -2.3231  -1.8976   3.0628  0.125  1.41  -0.911E+00    850.  0.225E+00  0.224E+00  0.286E-03  0.359E-13   14  26  18    1  13  16   11   8  19    0.000    0.000       0
    367   0    135 00:33:47.10  14.6000A  -2.4176  -1.8784   3.0724  0.125  1.41  -0.913E+00    850.  0.223E+00  0.222E+00  0.284E-03  0.355E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    368   0    136 00:34:01.70  14.9000A  -2.5034  -1.8777   3.0298  0.124  1.38  -0.912E+00    800.  0.221E+00  0.220E+00  0.281E-03  0.349E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    369   0    137 00:34:16.60  15.2000A  -2.5428  -1.8641   2.9686  0.124  1.38  -0.911E+00    800.  0.219E+00  0.218E+00  0.286E-03  0.353E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
    370   0    138 00:34:31.80  15.3000A  -2.4933  -1.8350   2.9454  0.123  1.38  -0.912E+00    800.  0.218E+00  0.217E+00  0.290E-03  0.357E-13   14  26  18    1  12  16   12   8  19    0.000    0.000       0
    371   0    139 00:34:47.10  15.4000A  -2.3596  -1.7912   2.9241  0.122  1.38  -0.917E+00    800.  0.216E+00  0.215E+00  0.292E-03  0.356E-13   14  27  18    1  12  16   12   8  19    0.000    0.000       0
    372   0    140 00:35:02.50  15.5000A  -2.3919  -1.7717   2.8943  0.122  1.38  -0.917E+00    800.  0.215E+00  0.214E+00  0.293E-03  0.355E-13   15  26  18    1  11  15   12   8  19    0.000    0.000       0
    373   0    141 00:35:18.00  15.8000A  -2.4704  -1.7378   2.8456  0.121  1.38  -0.922E+00    800.  0.214E+00  0.213E+00  0.294E-03  0.354E-13   15  26  17    1  11  15   12   8  19    0.000    0.000       0
    374   0    142 00:35:33.80  15.8000A  -2.4681   1.7976   2.8489  0.121  1.38  -0.929E+00    800.  0.213E+00  0.212E+00  0.299E-03  0.357E-13   15  26  17    8  34  23   14  26  24    0.000    0.000       0
    375   0    143 00:35:49.60  15.7000A  -2.3624  -1.7673   2.8671  0.120  1.38  -0.927E+00    800.  0.213E+00  0.212E+00  0.299E-03  0.354E-13   15  26  17    1   8  35    6  38  32    0.000    0.000       0
    376   0    144 00:36:05.30  15.2000A  -2.2471  -1.7490   2.9525  0.119  1.38  -0.930E+00    800.  0.212E+00  0.211E+00  0.297E-03  0.349E-13   15  26  16    1   7  34    6  38  32    0.000    0.000       0
    377   0    145 00:36:20.50  14.7000A  -2.0989  -1.8137   3.0686  0.119  1.38  -0.939E+00    800.  0.212E+00  0.211E+00  0.288E-03  0.337E-13   15  26  16    1   6  34    7  38  32    0.000    0.000       0
    378   0    146 00:36:35.20  14.2000A  -2.0337  -1.8432   3.1685  0.118  1.38  -0.941E+00    800.  0.212E+00  0.211E+00  0.278E-03  0.325E-13   11  11  13    1   6  34    7  38  32    0.000    0.000       0
    379   0    147 00:36:49.40  13.9000A  -2.0495  -1.8095   3.2329  0.118  1.38  -0.945E+00    800.  0.211E+00  0.211E+00  0.269E-03  0.313E-13   10  12  14    1   6  34    7  38  32    0.000    0.000       0
    380   0    148 00:37:03.30  13.6000A  -2.1261  -1.8874   3.2988  0.117  1.38  -0.949E+00    800.  0.211E+00  0.210E+00  0.264E-03  0.306E-13   10  12  14    1   5  34    8  38  32    0.000    0.000       0
    381   0    149 00:37:16.90  13.4000A  -2.1436  -1.8857   3.3518  0.117  1.38  -0.952E+00    800.  0.211E+00  0.210E+00  0.258E-03  0.299E-13   11  12  13    1   5  34    8  38  32    0.000    0.000       0
    382   0    150 00:37:30.30  13.4000A  -2.1473  -1.8654   3.3558  0.116  1.38  -0.957E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.294E-13   11  12  13    1   6  35    8  38  32    0.000    0.000       0
    383   0    151 00:37:43.70  13.5000A  -2.1150  -1.8966   3.3332  0.116  1.38  -0.961E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.293E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
    384   0    152 00:37:57.20  13.6000A  -2.0518  -1.9679   3.3052  0.116  1.38  -0.964E+00    800.  0.211E+00  0.210E+00  0.257E-03  0.293E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
    385   0    153 00:38:10.80  13.9000A  -1.9707  -1.9302   3.2273  0.115  1.41  -0.968E+00    850.  0.211E+00  0.210E+00  0.259E-03  0.294E-13   13  12  12    1   5  35    9  38  32    0.000    0.000       0
    386   0    154 00:38:24.70  14.2000A   1.9577  -2.0347   3.1608  0.115  1.41  -0.973E+00    850.  0.211E+00  0.210E+00  0.265E-03  0.298E-13    7  28  21    1   4  35    8  10   7    0.000    0.000       0
    387   0    155 00:38:38.90  13.4000A   1.9885  -2.0495   3.3484  0.115  1.41  -0.974E+00    850.  0.211E+00  0.210E+00  0.271E-03  0.303E-13    7  28  22    1   4  35    8   9   7    0.000    0.000       0
    388   0    156 00:38:52.30  12.8000A   2.0049  -1.9578   3.5085  0.114  1.41  -0.976E+00    850.  0.211E+00  0.210E+00  0.256E-03  0.286E-13    7  28  22    1   4  35    9   9   7    0.000    0.000       0
    389   0    157 00:39:05.10  12.6000A   1.9495  -2.0192   3.5708  0.114  1.41  -0.980E+00    850.  0.211E+00  0.210E+00  0.246E-03  0.273E-13    7  28  22    1   3  35    9   9   7    0.000    0.000       0
    390   0    158 00:39:17.70  12.5000A  -1.8994  -1.9976   3.5936  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.243E-03  0.268E-13   14  13  11    1   3  35   10   9   7    0.000    0.000       0
    391   0    159 00:39:30.20  12.8000A  -1.8343  -1.8951   3.5263  0.114  1.41  -0.982E+00    850.  0.211E+00  0.210E+00  0.241E-03  0.266E-13   14  13  11    1   3  35   10   8   6    0.000    0.000       0
    392   0    160 00:39:43.00  13.0000A  -1.8423   1.8539   3.4660  0.114  1.41  -0.988E+00    850.  0.212E+00  0.211E+00  0.248E-03  0.272E-13    8  13  14   17  15  11   10   8   6    0.000    0.000       0
    393   0    161 00:39:56.00  13.1000A  -1.8734   1.9441   3.4252  0.113  1.41  -0.987E+00    850.  0.212E+00  0.211E+00  0.252E-03  0.276E-13    8  13  14   17  15  11   11   8   6    0.000    0.000       0
    394   0    162 00:40:09.10  13.1000A  -1.8857   1.9860   3.4418  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.255E-03  0.278E-13    8  13  14   17  15  11   10   9   8    0.000    0.000       0
    395   0    163 00:40:22.20  12.9000A  -1.8946   1.9867   3.4953  0.113  1.41  -0.992E+00    850.  0.212E+00  0.211E+00  0.256E-03  0.279E-13    9  13  14   17  16  11   10   9   8    0.000    0.000       0
    396   0    164 00:40:35.10  13.0000A  -1.8988   2.1083   3.4484  0.113  1.41  -0.994E+00    850.  0.213E+00  0.211E+00  0.253E-03  0.275E-13    9  13  14   17  16  11   11   9   8    0.000    0.000       0
    397   0    165 00:40:48.10  13.3000A  -1.9469   2.1590   3.3863  0.113  1.41  -0.994E+00    850.  0.213E+00  0.212E+00  0.256E-03  0.278E-13   18  13   9   17  16  11   11   8   7    0.000    0.000       0
    398   0    166 00:41:01.40  13.5000A  -1.9768   2.1184   3.3307  0.113  1.41  -0.996E+00    850.  0.213E+00  0.212E+00  0.263E-03  0.285E-13   18  13   9   17  16  11   12   8   7    0.000    0.000       0
    399   0    167 00:41:14.90  13.7000A  -1.9837   2.1241   3.2823  0.113  1.41  -0.100E+01    850.  0.213E+00  0.212E+00  0.268E-03  0.290E-13   18  14   9   17  17  11   11   8   8    0.000    0.000       0
    400   0    168 00:41:28.60  13.6000A  -1.9913   2.1833   3.3064  0.113  1.41  -0.100E+01    850.  0.214E+00  0.212E+00  0.274E-03  0.295E-13   18  14   9   17  17  11   12   8   8    0.000    0.000       0
    401   0    169 00:41:42.20  13.4000A  -1.9813   2.1737   3.3472  0.113  1.41  -0.101E+01    850.  0.214E+00  0.213E+00  0.273E-03  0.294E-13   18  14   9   18  18   9   12   8   8    0.000    0.000       0
    402   0    170 00:41:55.60  13.5000A  -1.9944   2.2409   3.3230  0.113  1.38  -0.101E+01    800.  0.214E+00  0.213E+00  0.270E-03  0.291E-13   18  14   8   17  18  10   12   8   8    0.000    0.000       0
    403   0    171 00:42:09.10  13.5000A  -2.0090   2.2335   3.3380  0.113  1.38  -0.101E+01    800.  0.215E+00  0.213E+00  0.272E-03  0.293E-13   18  14   8   17  18  10   13   8   8    0.000    0.000       0
    404   0    172 00:42:22.60  13.6000A  -1.9963   2.1836   3.3089  0.113  1.38  -0.101E+01    800.  0.215E+00  0.214E+00  0.273E-03  0.294E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
    405   0    173 00:42:36.20  13.9000A  -1.9847   2.2215   3.2261  0.113  1.38  -0.101E+01    800.  0.216E+00  0.214E+00  0.276E-03  0.297E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
    406   0    174 00:42:50.10  14.1000A  -2.0020   2.1920   3.1998  0.113  1.38  -0.102E+01    800.  0.217E+00  0.215E+00  0.282E-03  0.304E-13   18  15   7   17  19  10   14   8   8    0.000    0.000       0
    407   0    175 00:43:04.20  14.3000A  -1.9699   2.2250   3.1385  0.113  1.38  -0.101E+01    800.  0.218E+00  0.216E+00  0.286E-03  0.308E-13   18  15   7   17  20   9   14   8   8    0.000    0.000       0
    408   0    176 00:43:18.50  14.8000A  -1.9215   2.2241   3.0482  0.114  1.38  -0.100E+01    800.  0.219E+00  0.217E+00  0.290E-03  0.313E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
    409   0    177 00:43:33.30  15.2000A  -1.9580   2.1654   2.9561  0.114  1.38  -0.100E+01    800.  0.220E+00  0.218E+00  0.300E-03  0.324E-13   17  14   4   17  20   9   14   8   8    0.000    0.000       0
    410   0    178 00:43:48.50  15.5000A  -1.9076   2.1079   2.9114  0.114  1.38  -0.994E+00    800.  0.221E+00  0.219E+00  0.308E-03  0.332E-13   17  14   4   18  18  12   10   9  18    0.000    0.000       0
    411   0    179 00:44:04.00  15.0000A  -1.8769   2.1355   2.9998  0.114  1.38  -0.100E+01    800.  0.222E+00  0.221E+00  0.314E-03  0.339E-13   17  15   4   18  19  11   10   9  18    0.000    0.000       0
    412   0    180 00:44:19.00  14.7000A  -1.8140   2.1157   3.0702  0.115  1.38  -0.992E+00    800.  0.224E+00  0.222E+00  0.305E-03  0.328E-13   16  15   3   18  17  13   11   9  18    0.000    0.000       0
    413   0    181 00:44:33.70  14.6000A  -1.7678   2.1454   3.0915  0.115  1.43  -0.991E+00    900.  0.225E+00  0.223E+00  0.300E-03  0.323E-13   16  15   3   18  17  13   13  11  16    0.000    0.000       0
    414   0    182 00:44:48.30  14.2000A   1.7560   2.1502   3.1645  0.115  1.43  -0.995E+00    900.  0.226E+00  0.225E+00  0.299E-03  0.321E-13    9  27  23   18  17  13   13  11  16    0.000    0.000       0
    415   0    183 00:45:02.50  13.9000A   1.7854   2.1719   3.2326  0.116  1.43  -0.990E+00    900.  0.228E+00  0.226E+00  0.292E-03  0.314E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
    416   0    184 00:45:16.40  13.6000A   1.7977   2.1976   3.3144  0.116  1.43  -0.992E+00    900.  0.230E+00  0.228E+00  0.287E-03  0.308E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
    417   0    185 00:45:30.00  13.3000A   1.8109   2.2067   3.3797  0.116  1.43  -0.100E+01    900.  0.231E+00  0.230E+00  0.281E-03  0.302E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
    418   0    186 00:45:43.30  13.2000A   1.8581   2.2022   3.4134  0.117  1.43  -0.991E+00    900.  0.233E+00  0.231E+00  0.276E-03  0.297E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
    419   0    187 00:45:56.50  13.1000A   1.8920   2.1844   3.4439  0.117  1.43  -0.987E+00    900.  0.234E+00  0.233E+00  0.275E-03  0.295E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    420   0    188 00:46:09.60  12.9000A   1.9095   2.1525   3.4749  0.117  1.43  -0.985E+00    900.  0.236E+00  0.235E+00  0.274E-03  0.294E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    421   0    189 00:46:22.50  13.0000A   1.9088   2.1037   3.4713  0.118  1.43  -0.984E+00    900.  0.238E+00  0.237E+00  0.271E-03  0.291E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
    422   0    190 00:46:35.50  13.0000A   1.8892   2.0727   3.4501  0.118  1.43  -0.984E+00    900.  0.240E+00  0.238E+00  0.274E-03  0.294E-13    9  27  24   18  19  12   16  11  16    0.000    0.000       0
    423   0    191 00:46:48.50  13.0000A   1.8879   2.0544   3.4586  0.118  1.43  -0.986E+00    900.  0.242E+00  0.240E+00  0.275E-03  0.295E-13    9  26  25   18  19  12   16  11  16    0.000    0.000       0
    424   0    192 00:47:01.50  13.1000A   1.8897   2.0539   3.4431  0.119  1.41  -0.981E+00    850.  0.243E+00  0.242E+00  0.276E-03  0.296E-13    9  27  25    8   4   0   16  11  16    0.000    0.000       0
    425   0    193 00:47:14.60  13.2000A   1.9011   2.0297   3.4101  0.119  1.41  -0.978E+00    850.  0.245E+00  0.244E+00  0.279E-03  0.300E-13    9  27  25    8   4   0   13  10  17    0.000    0.000       0
    426   0    194 00:47:27.80  13.3000A   1.8908   2.0100   3.3868  0.119  1.41  -0.974E+00    850.  0.247E+00  0.246E+00  0.283E-03  0.305E-13    9  27  25   19  18  13   14  10  17    0.000    0.000       0
    427   0    195 00:47:41.10  13.2000A   1.8585   2.0275   3.4008  0.119  1.41  -0.970E+00    850.  0.249E+00  0.248E+00  0.286E-03  0.309E-13    9  27  25   19  18  14   14  10  17    0.000    0.000       0
    428   0    196 00:47:54.30  13.2000A  -1.8883  -1.9919   3.4148  0.119  1.41  -0.968E+00    850.  0.251E+00  0.250E+00  0.286E-03  0.310E-13    6   3   1   13  38  15   11   9  18    0.000    0.000       0
    429   0    197 00:48:07.50  12.9000A  -2.0109  -1.9287   3.4810  0.119  1.41  -0.973E+00    850.  0.253E+00  0.252E+00  0.288E-03  0.313E-13    6   3   1   13  38  15   12   9  18    0.000    0.000       0
    430   0    198 00:48:20.40  12.8000A  -2.0676   1.9262   3.5058  0.119  1.41  -0.971E+00    850.  0.255E+00  0.253E+00  0.283E-03  0.309E-13    6   3   1   19  19  13   12   9  18    0.000    0.000       0
    431   0    199 00:48:33.20  12.8000A  -2.0940   1.8888   3.5269  0.119  1.41  -0.972E+00    850.  0.257E+00  0.255E+00  0.283E-03  0.309E-13    6   4   0   19  19  13   13   9  18    0.000    0.000       0
    432   0    200 00:48:46.00  12.7000A   2.1500   1.8285   3.5411  0.119  1.41  -0.975E+00    850.  0.259E+00  0.257E+00  0.285E-03  0.312E-13    7   2  10   12   4  39   13   9  18    0.000    0.000       0
    433   0    201 00:48:58.70  12.7000A   2.1955   1.8433   3.5368  0.119  1.41  -0.974E+00    850.  0.261E+00  0.259E+00  0.284E-03  0.313E-13    7   2  10   12   4  39   14   9  18    0.000    0.000       0
    434   0    202 00:49:11.40  12.5000A  -2.2761  -1.8582   3.5983  0.119  1.41  -0.976E+00    850.  0.262E+00  0.261E+00  0.286E-03  0.316E-13    7   3   0    3   8   2    9  39   1    0.000    0.000       0
    435   0    203 00:49:23.90  12.9000A  -2.3204  -1.9224   3.4935  0.119  1.41  -0.973E+00    850.  0.264E+00  0.262E+00  0.284E-03  0.315E-13    7   3   0    3   8   2   15   9  18    0.000    0.000       0
    436   0    204 00:49:36.80  12.8000A  -2.2110  -1.9685   3.5230  0.119  1.41  -0.974E+00    850.  0.266E+00  0.264E+00  0.295E-03  0.327E-13    8   2   0    3   8   2    9   0   0    0.000    0.000       0
    437   0    205 00:49:49.60  12.6000A  -2.3818  -2.0373   3.5648  0.119  1.41  -0.975E+00    850.  0.268E+00  0.266E+00  0.294E-03  0.328E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
    438   0    206 00:50:02.20  12.4000A  -2.4161  -2.1012   3.6437  0.118  1.41  -0.973E+00    850.  0.269E+00  0.267E+00  0.291E-03  0.325E-13    8   1   0    3   7   2    9   7  17    0.000    0.000       0
    439   0    207 00:50:14.60  12.2000A  -2.3757  -2.1220   3.6851  0.118  1.41  -0.972E+00    850.  0.271E+00  0.269E+00  0.287E-03  0.322E-13    8   1   0    3   7   2   10   7  17    0.000    0.000       0
    440   0    208 00:50:26.80  12.0000A  -2.4802  -2.0799   3.7439  0.118  1.43  -0.972E+00    900.  0.272E+00  0.270E+00  0.284E-03  0.319E-13    9   1  39    3   7   2   10   7  17    0.000    0.000       0
    441   0    209 00:50:38.80  12.0000A  -2.5304  -2.0332   3.7519  0.118  1.43  -0.970E+00    900.  0.273E+00  0.271E+00  0.280E-03  0.315E-13    9   1  39    4   7   3   10   7  17    0.000    0.000       0
    442   0    210 00:50:50.80  12.1000A  -2.4939  -2.0329   3.7144  0.118  1.43  -0.968E+00    900.  0.275E+00  0.272E+00  0.281E-03  0.317E-13    9   1  39    3   6   2   10   7  17    0.000    0.000       0
    443   0    211 00:51:02.90  12.0000A  -2.4209  -1.9740   3.7581  0.118  1.43  -0.968E+00    900.  0.276E+00  0.274E+00  0.284E-03  0.321E-13    9   1  39    3   6   2   14   6  18    0.000    0.000       0
    444   0    212 00:51:14.90  11.7000A  -2.3846  -2.0341   3.8348  0.118  1.43  -0.970E+00    900.  0.277E+00  0.275E+00  0.283E-03  0.320E-13    9   1  38   15  33  38   14   6  18    0.000    0.000       0
    445   0    213 00:51:26.59  11.7000A  -2.3506  -2.1307   3.8576  0.117  1.43  -0.975E+00    900.  0.278E+00  0.276E+00  0.277E-03  0.314E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
    446   0    214 00:51:38.29  11.8000A  -2.3205  -2.1303   3.8242  0.117  1.43  -0.977E+00    900.  0.279E+00  0.277E+00  0.277E-03  0.315E-13   10   1  38   15  33  38   14   6  18    0.000    0.000       0
    447   0    215 00:51:50.09  12.0000A  -2.2789  -2.0978   3.7524  0.117  1.43  -0.981E+00    900.  0.280E+00  0.277E+00  0.281E-03  0.319E-13   11   0  38   15  32  38   15   5  18    0.000    0.000       0
    448   0    216 00:52:02.09  12.0000A  -2.2487  -2.0444   3.7616  0.117  1.43  -0.984E+00    900.  0.281E+00  0.278E+00  0.287E-03  0.326E-13   11   0  38   15  32  38   16   5  18    0.000    0.000       0
    449   0    217 00:52:14.09  12.0000A  -2.2217  -1.9955   3.7539  0.117  1.43  -0.980E+00    900.  0.281E+00  0.279E+00  0.288E-03  0.327E-13   11   0  37   15  32  37   16   5  18    0.000    0.000       0
    450   0    218 00:52:26.09  12.2000A  -2.1807  -2.0609   3.6997  0.117  1.43  -0.978E+00    900.  0.282E+00  0.279E+00  0.289E-03  0.329E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
    451   0    219 00:52:38.29  12.6000A  -2.0993  -2.1447   3.5845  0.117  1.43  -0.978E+00    900.  0.282E+00  0.280E+00  0.294E-03  0.335E-13   11   0  37   15  31  37   16   5  18    0.000    0.000       0
    452   0    220 00:52:50.89  13.2000A  -1.9941  -2.1586   3.4126  0.116  1.43  -0.979E+00    900.  0.283E+00  0.280E+00  0.304E-03  0.346E-13   13  23   4   15  31  37   16   5  18    0.000    0.000       0
    453   0    221 00:53:04.09  13.6000A  -1.9882  -2.0853   3.3005  0.116  1.43  -0.983E+00    900.  0.283E+00  0.280E+00  0.318E-03  0.362E-13   13  23   4   15  31  37   17   4  18    0.000    0.000       0
    454   0    222 00:53:17.69  13.6000A  -1.9653  -2.0588   3.2981  0.116  1.43  -0.993E+00    900.  0.283E+00  0.280E+00  0.328E-03  0.372E-13   13  23   4   15  30  37    8   4   4    0.000    0.000       0
    455   0    223 00:53:31.29  13.6000A  -1.9332   2.0510   3.3020  0.116  1.43  -0.994E+00    900.  0.283E+00  0.280E+00  0.328E-03  0.370E-13   14  24   3   17  17  10    8  35   5    0.000    0.000       0
    456   0    224 00:53:44.89  13.2000A  -1.9953   2.0070   3.4031  0.116  1.43  -0.987E+00    900.  0.283E+00  0.280E+00  0.327E-03  0.369E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
    457   0    225 00:53:58.09  13.1000A  -2.0501   1.9133   3.4417  0.116  1.43  -0.989E+00    900.  0.282E+00  0.280E+00  0.318E-03  0.357E-13    9   6   4   17  17  10    9  35   5    0.000    0.000       0
    458   0    226 00:54:11.19  13.1000A  -2.0938   1.8867   3.4369  0.115  1.43  -0.993E+00    900.  0.282E+00  0.279E+00  0.315E-03  0.352E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
    459   0    227 00:54:24.29  13.0000A  -2.1923   1.9127   3.4680  0.115  1.43  -0.100E+01    900.  0.282E+00  0.279E+00  0.315E-03  0.351E-13   10   6   4   13   8   4   10  35   5    0.000    0.000       0
    460   0    228 00:54:37.29  13.1000A  -2.2403   1.9712   3.4314  0.115  1.40  -0.995E+00    850.  0.282E+00  0.279E+00  0.312E-03  0.347E-13   10   6   4   14   9   5   10  35   5    0.000    0.000       0
    461   0    229 00:54:50.39  12.8000A  -2.2360   2.0200   3.5122  0.115  1.40  -0.994E+00    850.  0.282E+00  0.279E+00  0.315E-03  0.349E-13   10   6   4   14   9   5   10  36   6    0.000    0.000       0
    462   0    230 00:55:03.19  12.5000A  -2.2380   2.0059   3.5910  0.115  1.40  -0.996E+00    850.  0.282E+00  0.279E+00  0.308E-03  0.340E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
    463   0    231 00:55:15.69  12.4000A  -2.2664   1.9510   3.6299  0.115  1.40  -0.990E+00    850.  0.281E+00  0.279E+00  0.301E-03  0.332E-13   11   6   4   14   9   5   10  36   6    0.000    0.000       0
    464   0    232 00:55:28.09  12.3000A  -2.2644   1.9653   3.6637  0.115  1.40  -0.991E+00    850.  0.281E+00  0.279E+00  0.300E-03  0.329E-13   11   6   4   17  10   6   11  36   6    0.000    0.000       0
    465   0    233 00:55:40.39  12.1000A  -2.2635   1.9614   3.7037  0.115  1.40  -0.992E+00    850.  0.281E+00  0.279E+00  0.298E-03  0.326E-13   11   6   3   17  10   6   11  36   6    0.000    0.000       0
    466   0    234 00:55:52.49  12.1000A  -2.2125   1.9119   3.7274  0.115  1.40  -0.989E+00    850.  0.281E+00  0.278E+00  0.294E-03  0.321E-13   11   6   3   12   3  37   11  36   6    0.000    0.000       0
    467   0    235 00:56:04.59  12.1000A  -2.1711   1.8977   3.7116  0.115  1.40  -0.987E+00    850.  0.281E+00  0.278E+00  0.295E-03  0.321E-13   12   6   3   13   2  38   11  36   6    0.000    0.000       0
    468   0    236 00:56:16.69  12.3000A  -2.1905   1.9584   3.6493  0.115  1.40  -0.985E+00    850.  0.281E+00  0.278E+00  0.295E-03  0.321E-13   12   6   3   13   3  36   11  36   6    0.000    0.000       0
    469   0    237 00:56:28.99  12.3000A  -2.1734   1.9387   3.6577  0.115  1.40  -0.985E+00    850.  0.280E+00  0.278E+00  0.300E-03  0.327E-13   12   6   3   13   3  36   12  36   7    0.000    0.000       0
    470   0    238 00:56:41.29  12.4000A  -2.1181   1.9538   3.6195  0.115  1.40  -0.987E+00    850.  0.280E+00  0.278E+00  0.301E-03  0.327E-13   12   6   3   13   4  36   12  36   7    0.000    0.000       0
    471   0    239 00:56:53.69  12.5000A  -2.0250   1.9152   3.5962  0.115  1.38  -0.988E+00    800.  0.280E+00  0.277E+00  0.304E-03  0.330E-13   12   6   3   13   4  37   13  36   7    0.000    0.000       0
    472   0    240 00:57:06.19  12.7000A  -1.9906   1.9357   3.5442  0.115  1.38  -0.988E+00    800.  0.280E+00  0.277E+00  0.307E-03  0.334E-13   13   6   3   14   4  37   13  36   7    0.000    0.000       0
    473   0    241 00:57:18.89  12.8000A   2.0433   1.9496   3.5178  0.115  1.38  -0.990E+00    800.  0.279E+00  0.277E+00  0.313E-03  0.340E-13   12  37   9   14   5  37   14  35   6    0.000    0.000       0
    474   0    242 00:57:31.69  12.6000A   2.1147   1.9760   3.5706  0.115  1.40  -0.993E+00    850.  0.279E+00  0.277E+00  0.316E-03  0.343E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
    475   0    243 00:57:44.29  12.6000A   2.1795   1.9496   3.5723  0.115  1.40  -0.998E+00    850.  0.279E+00  0.276E+00  0.312E-03  0.339E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
    476   0    244 00:57:56.89  12.7000A   2.1966   1.9318   3.5502  0.114  1.40  -0.996E+00    850.  0.279E+00  0.276E+00  0.312E-03  0.340E-13   13  37   9   15   5  37   15  35   6    0.000    0.000       0
    477   0    245 00:58:09.59  12.7000A   2.1739   1.9579   3.5513  0.114  1.40  -0.996E+00    850.  0.279E+00  0.276E+00  0.315E-03  0.343E-13   13  37  10   15   5  37   15  35   6    0.000    0.000       0
    478   0    246 00:58:22.29  12.9000A   2.2250   1.9814   3.4946  0.114  1.38  -0.998E+00    800.  0.279E+00  0.276E+00  0.314E-03  0.344E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
    479   0    247 00:58:35.19  13.3000A   2.2578   1.9773   3.3900  0.114  1.38  -0.997E+00    800.  0.279E+00  0.276E+00  0.319E-03  0.350E-13   14  37  10   15   6  37   15  35   6    0.000    0.000       0
    480   0    248 00:58:48.49  13.6000A   2.2207   1.9381   3.3184  0.114  1.38  -0.998E+00    800.  0.279E+00  0.276E+00  0.329E-03  0.361E-13   14  37  10   15   6  37   15  34   3    0.000    0.000       0
    481   0    249 00:59:02.09  13.6000A   2.2069   1.9631   3.3056  0.114  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.336E-03  0.369E-13   14  37  11   16   6  37   12  37   6    0.000    0.000       0
    482   0    250 00:59:15.69  13.5000A   2.1185   1.9916   3.3423  0.114  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.336E-03  0.369E-13   14  37  11   16   7  37    9  37   5    0.000    0.000       0
    483   0    251 00:59:29.19  13.4000A   2.0242   2.0192   3.3527  0.113  1.43  -0.101E+01    900.  0.278E+00  0.276E+00  0.333E-03  0.365E-13   16  36  11   16   7  37   10  37   5    0.000    0.000       0
    484   0    252 00:59:42.59  13.2000A   1.9364   2.0137   3.4014  0.113  1.43  -0.101E+01    900.  0.278E+00  0.276E+00  0.331E-03  0.362E-13   15  37  12   16   7  37   10  37   5    0.000    0.000       0
    485   0    253 00:59:55.79  13.3000A   1.9199   2.0227   3.3908  0.113  1.43  -0.102E+01    900.  0.278E+00  0.276E+00  0.325E-03  0.355E-13   16  36  12   16   7  36   10  37   5    0.000    0.000       0
    486   0    254 01:00:09.09  13.2000A   1.9010   2.0017   3.4011  0.113  1.43  -0.102E+01    900.  0.278E+00  0.276E+00  0.328E-03  0.356E-13   16  36  12   16   7  36   11  37   5    0.000    0.000       0
     232  0      0 00:00:00.00   0.0200A  -0.2131D -0.2288D  0.1182  0.000  1.38   0.000E+00    800.  0.104E-02  0.103E-02  0.597E-03  0.867E-13    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     233  0      1 00:00:00.02   0.0400A  -0.2131  -0.2288   0.1182  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.326E-08  0.443E-18    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     234  0      2 00:00:00.06   0.0800A  -0.2131  -0.2288   0.1182  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.609E-08  0.827E-18    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     235  0      3 00:00:00.14   0.1600A  -0.2131  -0.2288   0.1182  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.121E-07  0.164E-17    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     236  0      4 00:00:00.30   0.3200A  -0.2131  -0.2288   0.1181  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.243E-07  0.328E-17    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     237  0      5 00:00:00.62   0.6400A  -0.2131  -0.2288   0.1181  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.490E-07  0.655E-17    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     238  0      6 00:00:01.26   1.2800A  -0.2131  -0.2287   0.1181  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.101E-06  0.131E-16    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     239  0      7 00:00:02.54   2.5600A  -0.2130  -0.2286   0.1180  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.214E-06  0.262E-16    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     240  0      8 00:00:05.09   5.1200D  -0.2130  -0.2285   0.1179  0.009  0.69  -0.131E+02    100.  0.104E-02  0.103E-02  0.481E-06  0.525E-16    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     241  0      9 00:00:10.21  10.2000D  -0.2129  -0.2281   0.1177  0.009  0.69  -0.131E+02    100.  0.103E-02  0.103E-02  0.117E-05  0.105E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     242  0     10 00:00:20.42  20.0000D  -0.2125  -0.2272   0.1172  0.009  0.69  -0.132E+02    100.  0.103E-02  0.102E-02  0.317E-05  0.210E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     243  0     11 00:00:40.42  20.0000D  -0.2117  -0.2252   0.1162  0.009  0.69  -0.133E+02    100.  0.102E-02  0.102E-02  0.942E-05  0.416E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     244  0     12 00:01:00.42  20.0000D  -0.2106  -0.2227   0.1152  0.009  0.69  -0.134E+02    100.  0.101E-02  0.101E-02  0.131E-04  0.421E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     245  0     13 00:01:20.42  20.0000D  -0.2093  -0.2197   0.1141  0.009  0.69  -0.134E+02    100.  0.100E-02  0.997E-03  0.167E-04  0.425E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     246  0     14 00:01:40.42  20.0000D  -0.2077  -0.2164   0.1129  0.009  0.93  -0.134E+02    250.  0.992E-03  0.988E-03  0.203E-04  0.428E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     247  0     15 00:02:00.42  20.0000D  -0.2059  -0.2126   0.1117  0.009  0.93  -0.133E+02    250.  0.983E-03  0.979E-03  0.240E-04  0.430E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     248  0     16 00:02:20.42  20.0000D  -0.2038  -0.2086   0.1103  0.009  1.22  -0.131E+02    550.  0.975E-03  0.970E-03  0.276E-04  0.431E-15    7   7  11   11   7  19    5  26   7    0.000    0.000       0
     249  0     17 00:02:40.42  20.0000D   0.3348D  0.3126D  0.1636  0.011  1.22  -0.105E+02    550.  0.202E-02  0.201E-02  0.601E-03  0.856E-13   10   3  39   10   8  27    5  36  28    0.000    0.000       0
     250  0     18 00:03:00.42  20.0000D   0.3330   0.3039   0.1620  0.011  1.22  -0.105E+02    550.  0.199E-02  0.198E-02  0.374E-04  0.830E-15   10   3  39   10   8  27    5  36  28    0.000    0.000       0
     251  0     19 00:03:20.42  20.0000D   0.3296  -0.2989   0.1601  0.011  1.22  -0.103E+02    550.  0.196E-02  0.195E-02  0.410E-04  0.828E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
     252  0     20 00:03:40.42  20.0000D   0.3249  -0.2987   0.1580  0.011  1.22  -0.101E+02    550.  0.194E-02  0.193E-02  0.444E-04  0.828E-15   10   3  39    8   4   6    5  36  28    0.000    0.000       0
     253  0     21 00:04:00.42  20.0000D   0.3189  -0.2990   0.1559  0.012  1.22  -0.974E+01    550.  0.191E-02  0.190E-02  0.479E-04  0.831E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
     254  0     22 00:04:20.42  20.0000D   0.3119  -0.3006   0.1536  0.012  1.22  -0.940E+01    550.  0.189E-02  0.188E-02  0.513E-04  0.836E-15   10   3  39    8   0   7    5  36  28    0.000    0.000       0
     255  0     23 00:04:40.42  20.0000D   0.3058  -0.3013   0.1512  0.013  1.22  -0.905E+01    550.  0.187E-02  0.186E-02  0.547E-04  0.842E-15   10   3   0    8   0   7    5  36  28    0.000    0.000       0
     256  0     24 00:05:00.42  20.0000D   0.4152D -0.3965D -0.1813  0.013  1.22  -0.850E+01    550.  0.292E-02  0.289E-02  0.625E-03  0.860E-13    8  17   5    8   3   6    8  19   3    0.000    0.000       0
     257  0     25 00:05:20.42  20.0000D   0.4184  -0.4078  -0.1811  0.014  1.22  -0.811E+01    550.  0.287E-02  0.285E-02  0.641E-04  0.129E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
     258  0     26 00:05:40.42  20.0000D   0.4164  -0.4153  -0.1797  0.015  1.22  -0.772E+01    550.  0.283E-02  0.280E-02  0.673E-04  0.129E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
     259  0     27 00:06:00.42  20.0000D   0.4098  -0.4191  -0.1773  0.015  1.22  -0.735E+01    550.  0.279E-02  0.276E-02  0.704E-04  0.130E-14    8  17   5    8   3   6    8  19   3    0.000    0.000       0
     260  0     28 00:06:20.42  20.0000D   0.3993  -0.4194   0.1804  0.016  1.22  -0.697E+01    550.  0.275E-02  0.273E-02  0.735E-04  0.131E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
     261  0     29 00:06:40.42  20.0000D   0.3856  -0.4167   0.1985  0.017  1.22  -0.655E+01    550.  0.272E-02  0.270E-02  0.766E-04  0.133E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
     262  0     30 00:07:00.42  20.0000D   0.3696  -0.4111   0.2201  0.018  1.22  -0.615E+01    550.  0.269E-02  0.267E-02  0.796E-04  0.136E-14    8  17   5    8   3   6    2  36   2    0.000    0.000       0
     263  0     31 00:07:20.42  20.0000D  -0.3708  -0.4031   0.2480  0.020  1.22  -0.576E+01    550.  0.267E-02  0.265E-02  0.825E-04  0.140E-14    8  13   9    8   3   6    2   6  12    0.000    0.000       0
     264  0     32 00:07:40.42  20.0000D  -0.4268D -0.4556D  0.2856  0.021  1.22  -0.538E+01    550.  0.372E-02  0.369E-02  0.653E-03  0.886E-13    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     265  0     33 00:08:00.42  20.0000D  -0.4241  -0.4571   0.3241  0.023  1.22  -0.505E+01    550.  0.369E-02  0.366E-02  0.906E-04  0.190E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     266  0     34 00:08:20.42  20.0000D  -0.4209  -0.4524   0.3681  0.024  1.22  -0.471E+01    550.  0.368E-02  0.365E-02  0.932E-04  0.196E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     267  0     35 00:08:40.42  20.0000D  -0.4176  -0.4424   0.4179  0.026  1.22  -0.439E+01    550.  0.369E-02  0.366E-02  0.957E-04  0.205E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     268  0     36 00:09:00.42  20.0000D  -0.4144  -0.4281   0.4739  0.028  1.22  -0.409E+01    550.  0.373E-02  0.371E-02  0.981E-04  0.215E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     269  0     37 00:09:20.42  20.0000D  -0.4117  -0.4107   0.5361  0.030  1.18  -0.381E+01    500.  0.381E-02  0.378E-02  0.100E-03  0.227E-14    4   4   4    8   1   8    2   6  12    0.000    0.000       0
     270  0     38 00:09:40.42  20.0000D  -0.4112  -0.4088   0.6046  0.032  1.18  -0.354E+01    500.  0.393E-02  0.391E-02  0.103E-03  0.243E-14    5   4   4    5   3  36    2   6  12    0.000    0.000       0
     271  0     39 00:10:00.42  20.0000D   0.4799D  0.4697D  0.7110  0.035  1.18  -0.329E+01    500.  0.520E-02  0.517E-02  0.656E-03  0.867E-13    5  14  20    9   7  14    2   6  12    0.000    0.000       0
     272  0     40 00:10:20.42  20.0000D   0.4589  -0.4761   0.8042  0.037  1.18  -0.308E+01    500.  0.544E-02  0.541E-02  0.108E-03  0.318E-14    5  14  20    8   2   6    3   6  12    0.000    0.000       0
     273  0     41 00:10:40.42  20.0000D   0.4452  -0.4785   0.9262  0.040  1.18  -0.287E+01    500.  0.580E-02  0.577E-02  0.110E-03  0.345E-14    6  39  30    8   2   6    3   6  12    0.000    0.000       0
     274  0     42 00:11:00.42  20.0000A   0.4350  -0.4770   1.0557  0.043  1.18  -0.266E+01    500.  0.631E-02  0.628E-02  0.111E-03  0.377E-14    6  39  30    8   2   6    3   6  12    0.000    0.000       0
     275  0     43 00:11:20.42  20.0000A  -0.4690  -0.5194   1.1880  0.046  1.18  -0.248E+01    500.  0.701E-02  0.698E-02  0.111E-03  0.417E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
     276  0     44 00:11:40.42  20.0000A  -0.5115  -0.5815   1.3172  0.049  1.18  -0.232E+01    500.  0.796E-02  0.793E-02  0.112E-03  0.465E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
     277  0     45 00:12:00.42  20.0000A  -0.5527  -0.6472   1.4375  0.052  1.18  -0.218E+01    500.  0.922E-02  0.920E-02  0.112E-03  0.523E-14    1   3  23    1  29  17    3   6  12    0.000    0.000       0
     278  0     46 00:12:20.42  20.0000A   0.6185  -0.7151   1.6172  0.055  1.18  -0.204E+01    500.  0.109E-01  0.109E-01  0.112E-03  0.593E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
     279  0     47 00:12:40.42  20.0000A   0.6878  -0.7807   1.7894  0.059  1.18  -0.193E+01    500.  0.131E-01  0.130E-01  0.112E-03  0.674E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
     280  0     48 00:13:00.42  20.0000A   0.7587  -0.8410   1.9157  0.062  1.18  -0.182E+01    500.  0.158E-01  0.158E-01  0.112E-03  0.769E-14    1  35  11    1  29  17    4   6  12    0.000    0.000       0
     281  0     49 00:13:20.42  20.0000A   0.8287  -0.8944   2.1048  0.066  1.18  -0.173E+01    500.  0.194E-01  0.193E-01  0.114E-03  0.879E-14    1  35  11    1  29  17    5   6  12    0.000    0.000       0
     282  0     50 00:13:40.42  20.0000A   0.8954  -0.9401   2.2237  0.069  1.18  -0.165E+01    500.  0.237E-01  0.237E-01  0.116E-03  0.101E-13    1  35  11    1  29  17    5   6  12    0.000    0.000       0
     283  0     51 00:14:00.42  18.9000A   0.9557  -0.9784   2.3807  0.073  1.28  -0.157E+01    650.  0.291E-01  0.291E-01  0.120E-03  0.115E-13    1  35  11    1  29  17    6   6  12    0.000    0.000       0
     284  0     52 00:14:19.32  18.4000A   1.0036  -1.0317   2.4493  0.076  1.28  -0.151E+01    650.  0.353E-01  0.352E-01  0.119E-03  0.123E-13    1  35  11    1  28  17    6   6  12    0.000    0.000       0
     285  0     53 00:14:37.72  17.6000A   1.0389  -1.1128   2.5521  0.078  1.28  -0.147E+01    650.  0.423E-01  0.422E-01  0.122E-03  0.135E-13    1  35  11    1  28  17    7   6  12    0.000    0.000       0
     286  0     54 00:14:55.32  16.2000A   1.0601  -1.1802   2.7746  0.081  1.28  -0.141E+01    650.  0.500E-01  0.499E-01  0.124E-03  0.145E-13    1  35  11    1  28  17    5  24   2    0.000    0.000       0
     287  0     55 00:15:11.52  15.4000A   1.0660  -1.2317   2.9274  0.083  1.28  -0.138E+01    650.  0.579E-01  0.579E-01  0.122E-03  0.148E-13    1  35  11    1  28  17    6  24   2    0.000    0.000       0
     288  0     56 00:15:26.92  14.3000A  -1.0747  -1.2702   3.1422  0.084  1.28  -0.135E+01    650.  0.663E-01  0.663E-01  0.124E-03  0.155E-13    1  24  19    1  28  17    6  24   2    0.000    0.000       0
     289  0     57 00:15:41.22  13.8000A   1.1042  -1.2968   3.2571  0.086  1.28  -0.132E+01    650.  0.747E-01  0.747E-01  0.122E-03  0.157E-13    1  17  38    1  28  17    6  24   2    0.000    0.000       0
     290  0     58 00:15:55.01  13.1000A   1.1518  -1.3146   3.4236  0.087  1.28  -0.131E+01    650.  0.834E-01  0.833E-01  0.125E-03  0.164E-13    1  17  38    1  28  17    7  24   2    0.000    0.000       0
     291  0     59 00:16:08.12  12.9000A   1.1885  -1.3252   3.4994  0.088  1.29  -0.130E+01    650.  0.920E-01  0.920E-01  0.126E-03  0.167E-13    1  17  38    1  28  17    7  24   2    0.000    0.000       0
     292  0     60 00:16:21.01  12.4000A   1.2164  -1.3304   3.6240  0.089  1.29  -0.129E+01    650.  0.101E+00  0.101E+00  0.130E-03  0.175E-13    1  17  38    1  28  17    8  24   2    0.000    0.000       0
     293  0     61 00:16:33.41  12.3000A   1.2351  -1.3457   3.6690  0.090  1.29  -0.127E+01    650.  0.110E+00  0.110E+00  0.131E-03  0.179E-13    1  17  38    1  28  18    8  34  13    0.000    0.000       0
     294  0     62 00:16:45.71  11.9000A   1.2464  -1.3673   3.7728  0.091  1.29  -0.126E+01    650.  0.119E+00  0.119E+00  0.136E-03  0.188E-13    1  17  38    1  28  18    9  34  13    0.000    0.000       0
     295  0     63 00:16:57.61  11.4000A   1.2904  -1.3825   3.9447  0.092  1.29  -0.126E+01    650.  0.128E+00  0.128E+00  0.138E-03  0.192E-13    5  31  10    1  28  18    9  34  13    0.000    0.000       0
     296  0     64 00:17:09.01  11.2000A   1.3177  -1.3915   4.0232  0.092  1.29  -0.124E+01    650.  0.137E+00  0.137E+00  0.138E-03  0.193E-13    5  31  10    1  28  18    9  34  13    0.000    0.000       0
     297  0     65 00:17:20.22  10.8000A  -1.3137  -1.3953   4.1834  0.093  1.38  -0.123E+01    800.  0.145E+00  0.145E+00  0.141E-03  0.198E-13    5  29  17    1  28  18   10  34  13    0.000    0.000       0
     298  0     66 00:17:31.01  10.6000A  -1.3848  -1.3945   4.2505  0.093  1.38  -0.123E+01    800.  0.154E+00  0.153E+00  0.141E-03  0.200E-13    5  29  17    1  28  18   10  34  13    0.000    0.000       0
     299  0     67 00:17:41.61  10.4000A  -1.4187  -1.3899   4.3464  0.094  1.38  -0.122E+01    800.  0.162E+00  0.162E+00  0.144E-03  0.204E-13    5  29  17    1  28  18   11  34  13    0.000    0.000       0
     300  0     68 00:17:52.01  10.2000A  -1.4319   1.3894   4.4010  0.095  1.38  -0.121E+01    800.  0.170E+00  0.170E+00  0.146E-03  0.208E-13    5  29  17    7  37  19   11  34  13    0.000    0.000       0
     301  0     69 00:18:02.22  10.1000A  -1.4430   1.4896   4.4566  0.095  1.38  -0.120E+01    800.  0.178E+00  0.177E+00  0.148E-03  0.212E-13    4  32  19    7  37  19   12  34  13    0.000    0.000       0
     302  0     70 00:18:12.31  10.0000A  -1.4707  -1.5187   4.4834  0.096  1.38  -0.120E+01    800.  0.185E+00  0.185E+00  0.151E-03  0.217E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
     303  0     71 00:18:22.31  10.0000A  -1.4864  -1.6347   4.4882  0.096  1.38  -0.119E+01    800.  0.193E+00  0.193E+00  0.154E-03  0.222E-13    4  32  19    8   7  16   13  34  13    0.000    0.000       0
     304  0     72 00:18:32.31  10.0000A  -1.4923  -1.7223   4.4794  0.097  1.38  -0.118E+01    800.  0.200E+00  0.200E+00  0.158E-03  0.229E-13    4  32  19    8   7  16   12  34  13    0.000    0.000       0
     305  0     73 00:18:42.31  10.1000A  -1.5096  -1.7607   4.4387  0.097  1.38  -0.116E+01    800.  0.207E+00  0.207E+00  0.163E-03  0.235E-13    4  32  18    8   7  16   14  34  13    0.000    0.000       0
     306  0     74 00:18:52.41  10.2000A  -1.5638  -1.7362   4.4065  0.098  1.41  -0.116E+01    850.  0.214E+00  0.214E+00  0.169E-03  0.245E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
     307  0     75 00:19:02.61  10.4000A  -1.5955  -1.6476   4.3472  0.099  1.41  -0.115E+01    850.  0.221E+00  0.221E+00  0.175E-03  0.253E-13    4  32  18    8   7  16   13  34  13    0.000    0.000       0
     308  0     76 00:19:13.01  10.5000A  -1.6043  -1.6443   4.2951  0.099  1.41  -0.115E+01    850.  0.228E+00  0.228E+00  0.182E-03  0.265E-13    4  32  18    9   7  16   15  35  12    0.000    0.000       0
     309  0     77 00:19:23.51  10.6000A  -1.5910  -1.6991   4.2632  0.100  1.41  -0.114E+01    850.  0.235E+00  0.235E+00  0.188E-03  0.273E-13    4  32  18    9   7  16   15  35  12    0.000    0.000       0
     310  0     78 00:19:34.11  10.8000A  -1.5967  -1.7244   4.1654  0.101  1.41  -0.113E+01    850.  0.241E+00  0.241E+00  0.194E-03  0.282E-13    7  12   6    9   7  16   16  35  13    0.000    0.000       0
     311  0     79 00:19:44.91  10.8000A  -1.5986  -1.7950   4.1824  0.102  1.41  -0.113E+01    850.  0.247E+00  0.247E+00  0.202E-03  0.293E-13    7  12   6   10   7  17   16  35  13    0.000    0.000       0
     312  0     80 00:19:55.71  11.0000A   1.6051  -1.8316   4.1076  0.102  1.41  -0.112E+01    850.  0.253E+00  0.253E+00  0.206E-03  0.298E-13    7   1  37   10   7  17   12  23  18    0.000    0.000       0
     313  0     81 00:20:06.71  10.9000A   1.7144  -1.8050   4.1426  0.103  1.41  -0.111E+01    850.  0.258E+00  0.258E+00  0.214E-03  0.308E-13    7  24  30   10   7  17   12  23  18    0.000    0.000       0
     314  0     82 00:20:17.61  10.8000A   1.8454   1.7495   4.1494  0.104  1.41  -0.109E+01    850.  0.263E+00  0.262E+00  0.215E-03  0.310E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
     315  0     83 00:20:28.41  10.7000A   1.9170   1.8285   4.2009  0.105  1.41  -0.107E+01    850.  0.267E+00  0.266E+00  0.217E-03  0.310E-13    7  24  30   14  18  28   13  23  18    0.000    0.000       0
     316  0     84 00:20:39.11  10.7000A   1.9339  -1.8947   4.1949  0.106  1.41  -0.106E+01    850.  0.270E+00  0.270E+00  0.218E-03  0.311E-13    7  24  30    9   4   0   13  23  17    0.000    0.000       0
     317  0     85 00:20:49.81  10.7000A   1.9419  -1.9959   4.2222  0.107  1.41  -0.105E+01    850.  0.273E+00  0.273E+00  0.220E-03  0.314E-13    7  24  31    9   4   0   13  23  17    0.000    0.000       0
     318  0     86 00:21:00.51  10.5000A   2.0215  -2.0698   4.2704  0.108  1.41  -0.104E+01    850.  0.276E+00  0.275E+00  0.223E-03  0.317E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
     319  0     87 00:21:11.01  10.5000A   2.0535  -2.1077   4.3029  0.109  1.41  -0.103E+01    850.  0.278E+00  0.277E+00  0.221E-03  0.314E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
     320  0     88 00:21:21.51  10.4000A   2.0493  -2.1047   4.3239  0.110  1.41  -0.102E+01    850.  0.280E+00  0.279E+00  0.223E-03  0.316E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
     321  0     89 00:21:31.92  10.4000A   2.0209  -2.0588   4.3203  0.111  1.41  -0.101E+01    850.  0.281E+00  0.280E+00  0.222E-03  0.315E-13    7  24  31    9   4   0   14  23  17    0.000    0.000       0
     322  0     90 00:21:42.32  10.5000A   1.9795  -1.9725   4.2984  0.112  1.41  -0.100E+01    850.  0.282E+00  0.281E+00  0.224E-03  0.316E-13    7  24  31    9   4   0   15  23  17    0.000    0.000       0
     323  0     91 00:21:52.82  10.4000A   1.9323  -1.9477   4.3275  0.113  1.41  -0.982E+00    850.  0.283E+00  0.282E+00  0.227E-03  0.321E-13    7  24  31   10   4   0   15  23  17    0.000    0.000       0
     324  0     92 00:22:03.22  10.5000A   1.8929  -2.0067   4.3025  0.114  1.41  -0.971E+00    850.  0.283E+00  0.283E+00  0.226E-03  0.320E-13   14   4   0   10   4   0   15  23  17    0.000    0.000       0
     325  0     93 00:22:13.72  10.7000A   1.9405  -2.0350   4.2025  0.115  1.41  -0.961E+00    850.  0.283E+00  0.283E+00  0.229E-03  0.324E-13   14   4   0   10   4   0   15  23  17    0.000    0.000       0
     326  0     94 00:22:24.42  10.8000A   1.9307  -2.0212   4.1558  0.116  1.41  -0.952E+00    850.  0.284E+00  0.283E+00  0.234E-03  0.332E-13   14   4   0   10   4   0   16  23  17    0.000    0.000       0
     327  0     95 00:22:35.22  11.0000A   1.8885  -1.9605   4.0829  0.117  1.38  -0.945E+00    800.  0.284E+00  0.283E+00  0.237E-03  0.336E-13   14   4   0   10   4   0   16  23  17    0.000    0.000       0
     328  0     96 00:22:46.22  11.4000A   1.8287  -1.8650   3.9342  0.118  1.38  -0.940E+00    800.  0.283E+00  0.283E+00  0.242E-03  0.344E-13   14   4   0   10   3   0   16  23  17    0.000    0.000       0
     329  0     97 00:22:57.62  12.0000A   1.8982  -1.7985   3.7406  0.119  1.38  -0.929E+00    800.  0.283E+00  0.283E+00  0.252E-03  0.357E-13   15   4   1    7  24   4   16  23  17    0.000    0.000       0
     330  0     98 00:23:09.62  12.6000A   1.9527  -1.8230   3.5839  0.119  1.38  -0.923E+00    800.  0.283E+00  0.282E+00  0.265E-03  0.376E-13   15   4   1    8  23   4   17  23  17    0.000    0.000       0
     331  0     99 00:23:22.22  12.9000A   1.9230  -1.9304   3.4804  0.120  1.38  -0.920E+00    800.  0.282E+00  0.282E+00  0.278E-03  0.395E-13   15   4   1    8  23   4   16  12  29    0.000    0.000       0
     332  0    100 00:23:35.12  13.2000A   1.8476  -2.0257   3.4144  0.121  1.38  -0.927E+00    800.  0.282E+00  0.281E+00  0.284E-03  0.404E-13   15   4   1    8  23   4   11  22  39    0.000    0.000       0
     333  0    101 00:23:48.32  13.4000A   1.7595  -2.1063   3.3486  0.121  1.38  -0.909E+00    800.  0.281E+00  0.281E+00  0.290E-03  0.412E-13   15   4   1    8  23   4   12  22  39    0.000    0.000       0
     334  0    102 00:24:01.72  13.7000A   1.7144  -2.1673   3.2912  0.122  1.35  -0.903E+00    750.  0.280E+00  0.280E+00  0.293E-03  0.418E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     335  0    103 00:24:15.42  13.9000A   1.7175  -2.2012   3.2333  0.122  1.35  -0.901E+00    750.  0.280E+00  0.279E+00  0.298E-03  0.426E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     336  0    104 00:24:29.32  14.1000A   1.7200  -2.1940   3.1830  0.123  1.35  -0.895E+00    750.  0.279E+00  0.279E+00  0.300E-03  0.431E-13    3  22  38    8  23   4   12  22  39    0.000    0.000       0
     337  0    105 00:24:43.42  14.1000A   1.7151  -2.1666   3.1969  0.123  1.35  -0.891E+00    750.  0.279E+00  0.278E+00  0.303E-03  0.435E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
     338  0    106 00:24:57.52  13.8000A   1.6954  -2.1250   3.2513  0.123  1.35  -0.890E+00    750.  0.279E+00  0.278E+00  0.302E-03  0.432E-13    3  22  38    7  23   5   10  23   1    0.000    0.000       0
     339  0    107 00:25:11.32  13.8000A  -1.7019  -2.1113   3.2667  0.123  1.35  -0.890E+00    750.  0.279E+00  0.278E+00  0.295E-03  0.420E-13    8  25   1    8  23   5   11  23   1    0.000    0.000       0
     340  0    108 00:25:25.12  13.8000A   1.7170  -2.1861   3.2552  0.123  1.35  -0.892E+00    750.  0.278E+00  0.278E+00  0.295E-03  0.417E-13   10  26  31    8  23   5   11  23   1    0.000    0.000       0
     341  0    109 00:25:38.92  13.8000A   1.7572  -2.2020   3.2518  0.124  1.41  -0.893E+00    850.  0.278E+00  0.277E+00  0.294E-03  0.413E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
     342  0    110 00:25:52.72  14.1000A   1.7772  -2.1520   3.1985  0.124  1.38  -0.902E+00    800.  0.278E+00  0.277E+00  0.293E-03  0.409E-13   10  26  31    8  23   5   12  23   1    0.000    0.000       0
     343  0    111 00:26:06.82  14.2000A   1.7739  -2.0928   3.1792  0.124  1.38  -0.896E+00    800.  0.278E+00  0.277E+00  0.299E-03  0.413E-13   10  26  31    8  22   5   13  23   1    0.000    0.000       0
     344  0    112 00:26:21.02  14.3000A  -1.8071  -2.1438   3.1390  0.124  1.38  -0.894E+00    800.  0.277E+00  0.277E+00  0.300E-03  0.412E-13   16  22  38    8  22   5   12  23  19    0.000    0.000       0
     345  0    113 00:26:35.32  14.4000A   1.7970  -2.1042   3.1245  0.124  1.38  -0.894E+00    800.  0.277E+00  0.276E+00  0.301E-03  0.411E-13    1  36  26    8  22   5   12  23  19    0.000    0.000       0
     346  0    114 00:26:49.72  14.4000A   1.8521  -1.9823   3.1277  0.124  1.38  -0.894E+00    800.  0.276E+00  0.276E+00  0.301E-03  0.410E-13    1  37  26    8  22   5    9  38  26    0.000    0.000       0
     347  0    115 00:27:04.12  14.1000A   1.8655   1.8641   3.2017  0.125  1.38  -0.898E+00    800.  0.276E+00  0.275E+00  0.299E-03  0.406E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
     348  0    116 00:27:18.22  13.8000A   1.8259   1.8684   3.2623  0.125  1.38  -0.900E+00    800.  0.275E+00  0.274E+00  0.292E-03  0.396E-13    1  37  26    8  35  24   10  38  26    0.000    0.000       0
     349  0    117 00:27:32.02  13.7000A   1.8466   1.9007   3.2755  0.125  1.38  -0.899E+00    800.  0.274E+00  0.274E+00  0.285E-03  0.386E-13    1  16  38    8  35  23   10  38  26    0.000    0.000       0
     350  0    118 00:27:45.72  13.8000A   1.9138   1.9282   3.2586  0.125  1.38  -0.900E+00    800.  0.274E+00  0.273E+00  0.282E-03  0.382E-13    1  16  39    8  35  23   10  38  26    0.000    0.000       0
     351  0    119 00:27:59.52  13.8000A   1.9898   1.9444   3.2550  0.126  1.38  -0.899E+00    800.  0.273E+00  0.272E+00  0.283E-03  0.384E-13    1  16  39    8  36  23   11  38  26    0.000    0.000       0
     352  0    120 00:28:13.32  13.9000A   2.0594   2.0384   3.2273  0.126  1.38  -0.897E+00    800.  0.272E+00  0.271E+00  0.283E-03  0.384E-13    1  16  39    9  36  24   11  38  26    0.000    0.000       0
     353  0    121 00:28:27.22  14.0000A   2.1099   2.0924   3.2172  0.126  1.41  -0.895E+00    850.  0.271E+00  0.271E+00  0.285E-03  0.386E-13    1  16  39    9  36  24   13  37  36    0.000    0.000       0
     354  0    122 00:28:41.22  14.1000A   2.1217   2.1131   3.1998  0.126  1.41  -0.895E+00    850.  0.270E+00  0.270E+00  0.287E-03  0.388E-13    1  16  39    9  36  23   13  37  37    0.000    0.000       0
     355  0    123 00:28:55.32  14.0000A   2.1581   2.1139   3.2134  0.127  1.41  -0.892E+00    850.  0.269E+00  0.268E+00  0.288E-03  0.390E-13    1  16   0    9  36  23   13  37  37    0.000    0.000       0
     356  0    124 00:29:09.32  13.9000A   2.1890   2.1442   3.2397  0.127  1.41  -0.890E+00    850.  0.268E+00  0.267E+00  0.286E-03  0.387E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
     357  0    125 00:29:23.22  13.9000A   2.1613   2.1310   3.2296  0.127  1.41  -0.887E+00    850.  0.266E+00  0.266E+00  0.284E-03  0.383E-13    1  16   0    9  37  23   14  37  37    0.000    0.000       0
     358  0    126 00:29:37.12  14.2000A   2.1070   2.0393   3.1733  0.128  1.38  -0.889E+00    800.  0.265E+00  0.264E+00  0.285E-03  0.382E-13    1  16   1    9  37  23   14  37  37    0.000    0.000       0
     359  0    127 00:29:51.32  14.3000A   2.1120   1.8897   3.1443  0.128  1.38  -0.900E+00    800.  0.263E+00  0.262E+00  0.291E-03  0.389E-13    1  16   1    9  37  23   13   6  16    0.000    0.000       0
     360  0    128 00:30:05.62  14.4000A   2.0592  -1.9142   3.1322  0.128  1.38  -0.887E+00    800.  0.261E+00  0.260E+00  0.292E-03  0.390E-13    1  16   1   10   6  21   13   6  16    0.000    0.000       0
     361  0    129 00:30:20.02  14.0000A  -2.0568  -2.0174   3.2119  0.128  1.38  -0.884E+00    800.  0.259E+00  0.258E+00  0.294E-03  0.391E-13    9  10  18   10   6  21    7  25  19    0.000    0.000       0
     362  0    130 00:30:34.02  13.7000A  -2.0451  -2.0359   3.2843  0.129  1.38  -0.882E+00    800.  0.257E+00  0.256E+00  0.285E-03  0.378E-13    9  10  18   10   6  21    8  25  19    0.000    0.000       0
     363  0    131 00:30:47.72  13.3000A   1.9259  -2.1009   3.3730  0.129  1.41  -0.884E+00    850.  0.254E+00  0.254E+00  0.279E-03  0.368E-13    1  17   1   10   5  21    8  25  19    0.000    0.000       0
     364  0    132 00:31:01.02  13.1000A  -1.9186  -2.1189   3.4240  0.129  1.41  -0.878E+00    850.  0.252E+00  0.251E+00  0.270E-03  0.355E-13   10   9  18   10   5  21    9  25  19    0.000    0.000       0
     365  0    133 00:31:14.12  13.0000A  -1.9597  -2.0592   3.4698  0.129  1.41  -0.877E+00    850.  0.250E+00  0.249E+00  0.266E-03  0.348E-13   10   9  18   10   5  21    9  25  19    0.000    0.000       0
     366  0    134 00:31:27.12  13.0000A  -1.8415  -1.9297   3.4578  0.129  1.41  -0.880E+00    850.  0.247E+00  0.247E+00  0.263E-03  0.343E-13   10   9  18   10   5  21   10  25  19    0.000    0.000       0
     367  0    135 00:31:40.12  13.0000A  -1.9051  -1.9082   3.4511  0.128  1.41  -0.882E+00    850.  0.245E+00  0.244E+00  0.262E-03  0.340E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
     368  0    136 00:31:53.12  13.4000A  -1.9185  -1.8703   3.3702  0.128  1.41  -0.889E+00    850.  0.242E+00  0.242E+00  0.261E-03  0.338E-13   13  26  20    6  12  21   10  25  19    0.000    0.000       0
     369  0    137 00:32:06.52  13.7000A  -1.9511   1.8787   3.2846  0.128  1.41  -0.895E+00    850.  0.240E+00  0.239E+00  0.268E-03  0.346E-13   13  26  19   10  10  12   11  25  19    0.000    0.000       0
     370  0    138 00:32:20.22  13.9000A  -2.0414  -1.8778   3.2356  0.127  1.41  -0.896E+00    850.  0.237E+00  0.236E+00  0.273E-03  0.351E-13   13  26  19    1  13  16    9   9  19    0.000    0.000       0
     371  0    139 00:32:34.12  14.1000A  -2.0833  -1.8864   3.2006  0.127  1.41  -0.897E+00    850.  0.234E+00  0.234E+00  0.275E-03  0.353E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
     372  0    140 00:32:48.22  14.2000A  -2.1165  -1.8950   3.1741  0.127  1.41  -0.899E+00    850.  0.232E+00  0.231E+00  0.278E-03  0.355E-13   13  26  19    1  13  16   10   9  19    0.000    0.000       0
     373  0    141 00:33:02.42  14.7000A  -2.1634  -1.9011   3.0682  0.126  1.41  -0.901E+00    850.  0.229E+00  0.229E+00  0.278E-03  0.355E-13   14  26  19    1  13  16   10   8  19    0.000    0.000       0
     374  0    142 00:33:17.12  14.7000A  -2.2424  -1.9018   3.0542  0.126  1.41  -0.905E+00    850.  0.227E+00  0.226E+00  0.287E-03  0.363E-13   14  26  19    1  13  16   10   8  19    0.000    0.000       0
     375  0    143 00:33:31.82  14.7000A  -2.3446  -1.8939   3.0618  0.125  1.41  -0.909E+00    850.  0.225E+00  0.224E+00  0.286E-03  0.360E-13   14  26  19    1  13  16   11   8  19    0.000    0.000       0
     376  0    144 00:33:46.52  14.7000A  -2.4366  -1.8746   3.0702  0.125  1.41  -0.916E+00    850.  0.223E+00  0.222E+00  0.285E-03  0.356E-13   14  26  18    1  13  16   11   8  19    0.000    0.000       0
     377  0    145 00:34:01.22  14.9000A  -2.5232  -1.8732   3.0285  0.124  1.38  -0.914E+00    800.  0.221E+00  0.220E+00  0.284E-03  0.352E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
     378  0    146 00:34:16.12  15.2000A  -2.5479  -1.8598   2.9686  0.124  1.38  -0.911E+00    800.  0.219E+00  0.219E+00  0.286E-03  0.354E-13   14  26  18    1  12  16   11   8  19    0.000    0.000       0
     379  0    147 00:34:31.32  15.2000A  -2.4724  -1.8309   2.9537  0.123  1.38  -0.912E+00    800.  0.218E+00  0.217E+00  0.291E-03  0.357E-13   14  26  18    1  12  16   12   8  19    0.000    0.000       0
     380  0    148 00:34:46.52  15.3000A  -2.3839  -1.7876   2.9361  0.123  1.38  -0.916E+00    800.  0.216E+00  0.216E+00  0.290E-03  0.354E-13   15  26  18    1  12  16   12   8  19    0.000    0.000       0
     381  0    149 00:35:01.82  15.5000A  -2.4539  -1.7567   2.9078  0.122  1.38  -0.917E+00    800.  0.215E+00  0.214E+00  0.291E-03  0.353E-13   15  26  18    1  11  15   12   8  19    0.000    0.000       0
     382  0    150 00:35:17.32  15.7000A  -2.5193   1.7234   2.8580  0.121  1.38  -0.923E+00    800.  0.214E+00  0.214E+00  0.294E-03  0.354E-13   15  26  17    8  34  23   12   8  19    0.000    0.000       0
     383  0    151 00:35:33.02  15.8000A  -2.4903   1.7913   2.8419  0.121  1.38  -0.928E+00    800.  0.214E+00  0.213E+00  0.297E-03  0.355E-13   15  26  17    8  34  23   14  26  24    0.000    0.000       0
     384  0    152 00:35:48.82  15.8000A  -2.3635  -1.7756   2.8479  0.120  1.38  -0.929E+00    800.  0.213E+00  0.212E+00  0.299E-03  0.354E-13   15  26  17    1   7  34    6  38  32    0.000    0.000       0
     385  0    153 00:36:04.62  15.3000A  -2.2428  -1.7658   2.9398  0.119  1.38  -0.931E+00    800.  0.212E+00  0.212E+00  0.299E-03  0.352E-13   15  26  16    1   7  34    6  38  32    0.000    0.000       0
     386  0    154 00:36:19.92  14.8000A  -2.0822  -1.8215   3.0479  0.119  1.38  -0.939E+00    800.  0.212E+00  0.211E+00  0.290E-03  0.339E-13   15  26  16    1   6  34    7  38  32    0.000    0.000       0
     387  0    155 00:36:34.72  14.3000A  -2.0180  -1.8586   3.1546  0.118  1.38  -0.940E+00    800.  0.212E+00  0.211E+00  0.280E-03  0.327E-13   11  11  13    1   6  34    7  38  32    0.000    0.000       0
     388  0    156 00:36:49.02  13.9000A  -2.0312  -1.8315   3.2277  0.118  1.38  -0.945E+00    800.  0.212E+00  0.211E+00  0.271E-03  0.315E-13    9  12  14    1   6  34    7  38  32    0.000    0.000       0
     389  0    157 00:37:02.92  13.7000A  -2.1165  -1.9066   3.2855  0.117  1.38  -0.951E+00    800.  0.212E+00  0.211E+00  0.264E-03  0.306E-13   10  12  13    1   5  34    8  38  32    0.000    0.000       0
     390  0    158 00:37:16.62  13.5000A  -2.1128  -1.9114   3.3441  0.117  1.38  -0.953E+00    800.  0.211E+00  0.211E+00  0.260E-03  0.301E-13   11  12  13    1   5  34    8  38  32    0.000    0.000       0
     391  0    159 00:37:30.12  13.4000A  -2.1380  -1.8734   3.3540  0.116  1.38  -0.958E+00    800.  0.211E+00  0.210E+00  0.257E-03  0.296E-13   11  12  13    1   6  35    8  38  32    0.000    0.000       0
     392  0    160 00:37:43.52  13.5000A  -2.1108  -1.8757   3.3222  0.116  1.38  -0.963E+00    800.  0.211E+00  0.210E+00  0.255E-03  0.292E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
     393  0    161 00:37:57.02  13.7000A  -2.0627  -1.9699   3.2944  0.115  1.38  -0.964E+00    800.  0.211E+00  0.210E+00  0.257E-03  0.293E-13   12  12  12    1   5  35    9  38  32    0.000    0.000       0
     394  0    162 00:38:10.72  14.0000A  -1.9774  -1.9542   3.2168  0.115  1.41  -0.970E+00    850.  0.211E+00  0.210E+00  0.261E-03  0.296E-13   13  12  12    1   5  35    9  38  32    0.000    0.000       0
     395  0    163 00:38:24.72  14.1000A   1.9456  -2.0340   3.2004  0.115  1.41  -0.974E+00    850.  0.211E+00  0.210E+00  0.266E-03  0.301E-13    7  27  22    1   4  35    8  10   7    0.000    0.000       0
     396  0    164 00:38:38.82  13.3000A  -1.9424  -2.0740   3.3950  0.115  1.41  -0.973E+00    850.  0.211E+00  0.210E+00  0.269E-03  0.302E-13   12  13  12    1   4  35    8   9   7    0.000    0.000       0
     397  0    165 00:38:52.12  12.6000A   1.9909  -2.0064   3.5667  0.114  1.41  -0.976E+00    850.  0.211E+00  0.210E+00  0.255E-03  0.284E-13    7  28  22    1   4  35    9   9   7    0.000    0.000       0
     398  0    166 00:39:04.72  12.8000A   1.9909  -2.0271   3.5276  0.114  1.41  -0.980E+00    850.  0.211E+00  0.210E+00  0.242E-03  0.269E-13    7  28  22    1   3  35    9   9   7    0.000    0.000       0
     399  0    167 00:39:17.52  12.8000A   1.9283  -2.0299   3.5283  0.114  1.41  -0.983E+00    850.  0.211E+00  0.210E+00  0.247E-03  0.273E-13    7  28  22    1   3  35   10   9   7    0.000    0.000       0
     400  0    168 00:39:30.32  12.9000A  -1.8906  -1.9463   3.4940  0.114  1.41  -0.984E+00    850.  0.211E+00  0.210E+00  0.247E-03  0.273E-13    8  13  15    1   3  35   10   8   6    0.000    0.000       0
     401  0    169 00:39:43.22  13.4000A  -1.9064   1.9409   3.3546  0.114  1.41  -0.985E+00    850.  0.212E+00  0.211E+00  0.250E-03  0.276E-13    8  13  14   17  15  11   10   8   6    0.000    0.000       0
     402  0    170 00:39:56.62  13.5000A  -1.9332   2.0092   3.3236  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.261E-03  0.286E-13    8  13  14   17  15  11    9   8   7    0.000    0.000       0
     403  0    171 00:40:10.12  13.0000A  -1.9327   2.0173   3.4509  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.263E-03  0.289E-13    8  13  14   17  15  11   10   9   8    0.000    0.000       0
     404  0    172 00:40:23.12  13.1000A  -1.9254   2.0677   3.4468  0.113  1.41  -0.989E+00    850.  0.212E+00  0.211E+00  0.255E-03  0.278E-13    9  13  14   17  16  11   10   9   8    0.000    0.000       0
     405  0    173 00:40:36.22  13.2000A  -1.9172   2.1591   3.4172  0.113  1.41  -0.991E+00    850.  0.212E+00  0.211E+00  0.258E-03  0.281E-13    9  13  14   17  16  11   11   8   7    0.000    0.000       0
     406  0    174 00:40:49.42  13.6000A  -1.9459   2.1733   3.3175  0.113  1.41  -0.993E+00    850.  0.213E+00  0.211E+00  0.261E-03  0.283E-13   18  13   9   17  16  11   11   8   7    0.000    0.000       0
     407  0    175 00:41:03.02  13.8000A  -1.9651   2.0985   3.2609  0.113  1.41  -0.995E+00    850.  0.213E+00  0.212E+00  0.270E-03  0.293E-13   18  13   9   17  17  11   11   8   8    0.000    0.000       0
     408  0    176 00:41:16.82  13.6000A  -1.9514   2.1932   3.3179  0.113  1.41  -0.999E+00    850.  0.213E+00  0.212E+00  0.275E-03  0.298E-13   18  13   9   17  17  11   11   8   8    0.000    0.000       0
     409  0    177 00:41:30.42  13.6000A   1.9172   2.2178   3.3025  0.113  1.41  -0.100E+01    850.  0.213E+00  0.212E+00  0.272E-03  0.294E-13   12   5  12   17  17  11   12   8   8    0.000    0.000       0
     410  0    178 00:41:44.02  13.5000A  -1.9023   2.2163   3.3328  0.113  1.41  -0.101E+01    850.  0.213E+00  0.212E+00  0.273E-03  0.295E-13   18  14   5   18  18   9   12   8   8    0.000    0.000       0
     411  0    179 00:41:57.52  13.7000A  -1.9313   2.2157   3.2800  0.113  1.41  -0.101E+01    850.  0.214E+00  0.212E+00  0.272E-03  0.294E-13   18  15   5   17  18  10   12   8   8    0.000    0.000       0
     412  0    180 00:42:11.22  13.9000A  -1.9220   2.1850   3.2404  0.113  1.38  -0.101E+01    800.  0.214E+00  0.213E+00  0.277E-03  0.298E-13   18  15   5   18  18   9   13   8   8    0.000    0.000       0
     413  0    181 00:42:25.12  14.1000A  -1.9163   2.2124   3.1917  0.113  1.38  -0.102E+01    800.  0.215E+00  0.213E+00  0.281E-03  0.303E-13   18  14   7   17  19  10   13   8   8    0.000    0.000       0
     414  0    182 00:42:39.22  14.6000A  -1.9486   2.2088   3.0894  0.113  1.38  -0.101E+01    800.  0.216E+00  0.214E+00  0.286E-03  0.308E-13   18  14   7   18  19   8   13   8   8    0.000    0.000       0
     415  0    183 00:42:53.82  15.0000A  -1.8983   2.1797   3.0026  0.113  1.38  -0.101E+01    800.  0.216E+00  0.215E+00  0.296E-03  0.319E-13   18  14   7   17  20   9   14   8   8    0.000    0.000       0
     416  0    184 00:43:08.82  15.3000A  -1.8954   2.2107   2.9328  0.113  1.38  -0.100E+01    800.  0.217E+00  0.216E+00  0.304E-03  0.328E-13   18  15   7   17  20   9   14   8   8    0.000    0.000       0
     417  0    185 00:43:24.12  15.8000A  -1.8630   2.1607   2.8548  0.114  1.38  -0.100E+01    800.  0.218E+00  0.217E+00  0.310E-03  0.334E-13   18  15   6   17  20   9   14   8   8    0.000    0.000       0
     418  0    186 00:43:39.92  15.6000A  -1.8678   2.0965   2.8897  0.114  1.38  -0.999E+00    800.  0.220E+00  0.218E+00  0.320E-03  0.344E-13   17  14   4   18  18  12    9   9  18    0.000    0.000       0
     419  0    187 00:43:55.52  15.0000A  -1.8314   2.0690   2.9995  0.114  1.38  -0.994E+00    800.  0.221E+00  0.219E+00  0.316E-03  0.340E-13   17  14   4   18  18  12   10   9  18    0.000    0.000       0
     420  0    188 00:44:10.52  14.8000A  -1.7310   2.1297   3.0403  0.114  1.38  -0.994E+00    800.  0.222E+00  0.221E+00  0.304E-03  0.328E-13   17  14   4   18  19  11   11   9  18    0.000    0.000       0
     421  0    189 00:44:25.32  14.6000A  -1.7140   2.1094   3.0835  0.115  1.38  -0.990E+00    800.  0.224E+00  0.222E+00  0.301E-03  0.324E-13   16  15   3   18  17  13   13  11  16    0.000    0.000       0
     422  0    190 00:44:39.92  14.2000A   1.7255   2.1319   3.1631  0.115  1.38  -0.991E+00    800.  0.225E+00  0.224E+00  0.298E-03  0.320E-13    9  27  23   18  17  13   13  11  16    0.000    0.000       0
     423  0    191 00:44:54.12  13.9000A   1.7630   2.1271   3.2369  0.115  1.38  -0.993E+00    800.  0.227E+00  0.225E+00  0.291E-03  0.313E-13    9  27  23   18  17  13   13  11  16    0.000    0.000       0
     424  0    192 00:45:08.02  13.6000A   1.7831   2.1534   3.2990  0.116  1.43  -0.991E+00    900.  0.228E+00  0.227E+00  0.286E-03  0.307E-13    9  27  23   18  18  13   14  11  16    0.000    0.000       0
     425  0    193 00:45:21.62  13.4000A   1.7944   2.1802   3.3668  0.116  1.43  -0.100E+01    900.  0.230E+00  0.228E+00  0.281E-03  0.301E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
     426  0    194 00:45:35.02  13.2000A   1.8443   2.1889   3.4074  0.116  1.43  -0.996E+00    900.  0.231E+00  0.230E+00  0.278E-03  0.298E-13    9  27  24   18  18  13   14  11  16    0.000    0.000       0
     427  0    195 00:45:48.22  13.2000A   1.8836   2.1814   3.4186  0.117  1.43  -0.991E+00    900.  0.233E+00  0.232E+00  0.274E-03  0.295E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     428  0    196 00:46:01.42  13.0000A   1.9089   2.1575   3.4497  0.117  1.43  -0.987E+00    900.  0.235E+00  0.233E+00  0.276E-03  0.296E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     429  0    197 00:46:14.42  13.1000A   1.9176   2.1154   3.4479  0.117  1.43  -0.985E+00    900.  0.236E+00  0.235E+00  0.273E-03  0.293E-13    9  27  24   18  18  13   15  11  16    0.000    0.000       0
     430  0    198 00:46:27.52  13.2000A   1.9085   2.0760   3.4203  0.118  1.43  -0.986E+00    900.  0.238E+00  0.237E+00  0.276E-03  0.296E-13    9  27  24   18  19  12   15  11  16    0.000    0.000       0
     431  0    199 00:46:40.72  13.2000A   1.8855   2.0566   3.4206  0.118  1.43  -0.987E+00    900.  0.240E+00  0.239E+00  0.279E-03  0.300E-13    9  26  25   18  19  12   16  11  16    0.000    0.000       0
     432  0    200 00:46:53.92  13.2000A   1.9077   2.0374   3.4011  0.118  1.43  -0.983E+00    900.  0.242E+00  0.241E+00  0.280E-03  0.301E-13    9  27  25    8   3   0   16  11  16    0.000    0.000       0
     433  0    201 00:47:07.12  13.4000A   1.9317   2.0179   3.3642  0.119  1.41  -0.978E+00    850.  0.244E+00  0.242E+00  0.281E-03  0.303E-13    9  27  25    8   4   0   16  11  16    0.000    0.000       0
     434  0    202 00:47:20.52  13.4000A   1.9346   2.0385   3.3511  0.119  1.41  -0.974E+00    850.  0.245E+00  0.244E+00  0.287E-03  0.310E-13    9  27  25   19  18  14   15  10  16    0.000    0.000       0
     435  0    203 00:47:33.92  13.3000A   1.9127   2.0929   3.3850  0.119  1.41  -0.977E+00    850.  0.247E+00  0.246E+00  0.288E-03  0.312E-13    9  27  25   19  18  14   15   9  17    0.000    0.000       0
     436  0    204 00:47:47.22  13.2000A   1.9149   2.0377   3.4163  0.119  1.41  -0.974E+00    850.  0.249E+00  0.248E+00  0.288E-03  0.312E-13    9  27  26   19  18  14   15   9  17    0.000    0.000       0
     437  0    205 00:48:00.42  13.2000A   1.9285   1.9118   3.4195  0.119  1.41  -0.972E+00    850.  0.251E+00  0.250E+00  0.287E-03  0.312E-13    9  27  26   19  19  13   15   9  17    0.000    0.000       0
     438  0    206 00:48:13.62  13.1000A  -2.0010   1.9771   3.4362  0.119  1.41  -0.973E+00    850.  0.253E+00  0.251E+00  0.288E-03  0.314E-13    6   3   1   17  13  10   16   9  17    0.000    0.000       0
     439  0    207 00:48:26.72  13.1000A  -2.0771   1.9674   3.4311  0.119  1.41  -0.982E+00    850.  0.255E+00  0.253E+00  0.287E-03  0.314E-13    6   3   0   17  13  10   16   9  17    0.000    0.000       0
     440  0    208 00:48:39.82  13.2000A  -2.1797   1.9858   3.4091  0.119  1.41  -0.977E+00    850.  0.257E+00  0.255E+00  0.289E-03  0.316E-13    6   3   0   17  14  10   13   9  18    0.000    0.000       0
     441  0    209 00:48:53.01  13.2000A  -2.1964   2.0063   3.4200  0.119  1.41  -0.973E+00    850.  0.259E+00  0.257E+00  0.293E-03  0.322E-13    6   3   0   17  14  10    8  39   1    0.000    0.000       0
     442  0    210 00:49:06.21  13.1000A  -2.2503   1.9578   3.4393  0.119  1.41  -0.973E+00    850.  0.261E+00  0.259E+00  0.295E-03  0.325E-13    7   2   1   17  15  10    9  39   1    0.000    0.000       0
     443  0    211 00:49:19.31  13.3000A  -2.2530  -1.9685   3.3805  0.119  1.41  -0.972E+00    850.  0.262E+00  0.260E+00  0.294E-03  0.325E-13    7   2   1    3   8   2    7   8  17    0.000    0.000       0
     444  0    212 00:49:32.61  12.9000A  -2.3251  -2.0228   3.4957  0.119  1.41  -0.974E+00    850.  0.264E+00  0.262E+00  0.300E-03  0.333E-13    7   3   0    3   8   2    9   0   0    0.000    0.000       0
     445  0    213 00:49:45.51  12.7000A  -2.3313  -2.0897   3.5383  0.119  1.41  -0.987E+00    850.  0.266E+00  0.264E+00  0.293E-03  0.326E-13    8   2   0    3   7   2    8   8  17    0.000    0.000       0
     446  0    214 00:49:58.21  12.4000A  -2.4812  -2.1653   3.6189  0.119  1.41  -0.995E+00    850.  0.267E+00  0.265E+00  0.290E-03  0.323E-13    8   1   0    3   7   2   10   7  17    0.000    0.000       0
     447  0    215 00:50:10.61  12.0000A  -2.5128  -2.1854   3.7471  0.119  1.41  -0.974E+00    850.  0.269E+00  0.267E+00  0.284E-03  0.318E-13    8   1   0    3   7   2   10   7  17    0.000    0.000       0
     448  0    216 00:50:22.61  11.8000A  -2.4764  -2.1261   3.7981  0.119  1.41  -0.973E+00    850.  0.270E+00  0.268E+00  0.276E-03  0.310E-13    9   1  39    3   7   2   10   7  17    0.000    0.000       0
     449  0    217 00:50:34.42  11.9000A  -2.6468  -2.0821   3.7849  0.118  1.43  -0.973E+00    900.  0.272E+00  0.269E+00  0.273E-03  0.307E-13    9   1  39    3   6   2   10   7  17    0.000    0.000       0
     450  0    218 00:50:46.32  11.8000A  -2.6822  -2.0757   3.8043  0.118  1.43  -0.970E+00    900.  0.273E+00  0.271E+00  0.276E-03  0.311E-13    9   1  39    3   6   2   11   7  17    0.000    0.000       0
     451  0    219 00:50:58.12  11.9000A  -2.6288  -2.0878   3.7871  0.118  1.43  -0.968E+00    900.  0.274E+00  0.272E+00  0.275E-03  0.310E-13    9   1  39    3   7   3   11   7  17    0.000    0.000       0
     452  0    220 00:51:10.02  12.0000A  -2.6844  -2.0617   3.7357  0.118  1.43  -0.969E+00    900.  0.275E+00  0.273E+00  0.278E-03  0.313E-13    9   1  38    3   7   3   14   6  18    0.000    0.000       0
     453  0    221 00:51:22.02  11.9000A  -2.6537  -1.9889   3.7962  0.118  1.43  -0.972E+00    900.  0.276E+00  0.274E+00  0.282E-03  0.317E-13    9   1  38    3   7   3   14   6  18    0.000    0.000       0
     454  0    222 00:51:33.92  11.8000A  -2.5385  -1.9568   3.7990  0.118  1.43  -0.977E+00    900.  0.278E+00  0.275E+00  0.280E-03  0.316E-13    9   1  38    3   6   3   14   6  18    0.000    0.000       0
     455  0    223 00:51:45.72  11.8000A  -2.4135  -1.9337   3.8177  0.118  1.43  -0.979E+00    900.  0.279E+00  0.276E+00  0.279E-03  0.314E-13   11   0  38    3   6   3   15   6  18    0.000    0.000       0
     456  0    224 00:51:57.52  11.9000A  -2.4027  -1.9037   3.7946  0.117  1.43  -0.982E+00    900.  0.280E+00  0.277E+00  0.279E-03  0.315E-13   11   0  37   15  32  38   15   6  18    0.000    0.000       0
     457  0    225 00:52:09.42  12.1000A  -2.4532  -1.8845   3.7186  0.117  1.43  -0.981E+00    900.  0.280E+00  0.278E+00  0.282E-03  0.319E-13   11   0  37   16  32  39   15   6  18    0.000    0.000       0
     458  0    226 00:52:21.52  12.2000A  -2.4193  -1.8701   3.6911  0.117  1.43  -0.982E+00    900.  0.281E+00  0.279E+00  0.288E-03  0.326E-13   11   0  37   16  32  39   16   5  18    0.000    0.000       0
     459  0    227 00:52:33.72  12.4000A  -2.3308  -1.8827   3.6176  0.117  1.43  -0.982E+00    900.  0.282E+00  0.279E+00  0.291E-03  0.329E-13   11   0  37    3   5   1   16   5  18    0.000    0.000       0
     460  0    228 00:52:46.12  12.6000A  -2.2282   1.9090   3.5577  0.117  1.43  -0.982E+00    900.  0.282E+00  0.280E+00  0.296E-03  0.335E-13   12   0  36   17  16  10   17   5  18    0.000    0.000       0
     461  0    229 00:52:58.72  13.1000A  -2.1808   1.9565   3.4457  0.117  1.43  -0.984E+00    900.  0.283E+00  0.280E+00  0.301E-03  0.341E-13   12   0  36   17  16  10   17   5  18    0.000    0.000       0
     462  0    230 00:53:11.82  13.7000A  -2.0728   1.9749   3.2822  0.116  1.43  -0.985E+00    900.  0.283E+00  0.280E+00  0.314E-03  0.354E-13   12   0  36   17  16  10   15  36  39    0.000    0.000       0
     463  0    231 00:53:25.52  13.6000A  -1.9391   1.9769   3.3165  0.116  1.43  -0.987E+00    900.  0.283E+00  0.280E+00  0.328E-03  0.369E-13   12   0  36   17  17  11    9  35   5    0.000    0.000       0
     464  0    232 00:53:39.12  13.2000A  -1.8656   1.9871   3.4090  0.116  1.43  -0.988E+00    900.  0.283E+00  0.281E+00  0.326E-03  0.365E-13   12   1  39   17  17  11    9  35   5    0.000    0.000       0
     465  0    233 00:53:52.32  13.2000A  -1.9476   1.9543   3.4211  0.116  1.43  -0.988E+00    900.  0.283E+00  0.281E+00  0.316E-03  0.354E-13    9   6   4   17  17  11   10  35   5    0.000    0.000       0
     466  0    234 00:54:05.51  12.8000A  -2.0044   1.9320   3.5149  0.116  1.43  -0.988E+00    900.  0.283E+00  0.280E+00  0.316E-03  0.353E-13    9   6   4   14   8   5   11   6  10    0.000    0.000       0
     467  0    235 00:54:18.32  13.0000A  -2.1071   1.9464   3.4637  0.115  1.43  -0.989E+00    900.  0.283E+00  0.280E+00  0.307E-03  0.342E-13   10   6   4   14   9   5   10  35   5    0.000    0.000       0
     468  0    236 00:54:31.32  13.0000A  -2.2080   2.0654   3.4627  0.115  1.43  -0.992E+00    900.  0.283E+00  0.280E+00  0.311E-03  0.346E-13   10   6   4   11   2  38    9  36   6    0.000    0.000       0
     469  0    237 00:54:44.32  12.5000A  -2.2464   2.1078   3.6015  0.115  1.43  -0.992E+00    900.  0.283E+00  0.280E+00  0.311E-03  0.346E-13   10   6   4   11   2  38   10  36   6    0.000    0.000       0
     470  0    238 00:54:56.82  12.3000A  -2.2361   2.0862   3.6655  0.115  1.40  -0.991E+00    850.  0.283E+00  0.280E+00  0.300E-03  0.332E-13   11   6   4   11   2  38   10  36   6    0.000    0.000       0
     471  0    239 00:55:09.12  12.2000A  -2.2916   2.0375   3.6744  0.115  1.40  -0.991E+00    850.  0.283E+00  0.280E+00  0.296E-03  0.327E-13   11   6   4   11   2  39   11  36   6    0.000    0.000       0
     472  0    240 00:55:21.32  12.1000A  -2.2994   2.0789   3.7339  0.115  1.40  -0.995E+00    850.  0.283E+00  0.280E+00  0.294E-03  0.324E-13   11   6   4   12   2  38   11  36   6    0.000    0.000       0
     473  0    241 00:55:33.42  12.0000A  -2.2669   2.1541   3.7608  0.115  1.40  -0.993E+00    850.  0.283E+00  0.280E+00  0.293E-03  0.321E-13   11   6   3   12   3  38   11  36   6    0.000    0.000       0
     474  0    242 00:55:45.42  12.0000A  -2.2433   2.2386   3.7470  0.115  1.40  -0.991E+00    850.  0.283E+00  0.280E+00  0.291E-03  0.319E-13   11   6   3   12   3  37   11  36   6    0.000    0.000       0
     475  0    243 00:55:57.42  11.8000A  -2.1748   2.2832   3.8004  0.115  1.40  -0.987E+00    850.  0.283E+00  0.280E+00  0.292E-03  0.319E-13   12   6   3   12   3  37   12  36   7    0.000    0.000       0
     476  0    244 00:56:09.22  11.8000A  -2.2169   2.2247   3.8031  0.115  1.40  -0.984E+00    850.  0.283E+00  0.280E+00  0.287E-03  0.314E-13   12   6   3   12   3  37   12  36   7    0.000    0.000       0
     477  0    245 00:56:21.02  12.0000A  -2.2126   2.2100   3.7549  0.115  1.40  -0.983E+00    850.  0.282E+00  0.280E+00  0.287E-03  0.314E-13   12   6   3   13   3  38   12  36   7    0.000    0.000       0
     478  0    246 00:56:33.02  12.0000A  -2.1607   2.2718   3.7616  0.115  1.43  -0.985E+00    900.  0.282E+00  0.279E+00  0.292E-03  0.319E-13   12   6   3   13   4  37   13  36   7    0.000    0.000       0
     479  0    247 00:56:45.02  12.1000A  -2.0640   2.3360   3.7201  0.115  1.43  -0.988E+00    900.  0.281E+00  0.279E+00  0.292E-03  0.320E-13   12   6   3   13   4  37   13  36   7    0.000    0.000       0
     480  0    248 00:56:57.12  12.4000A  -2.0537   2.3067   3.6170  0.115  1.43  -0.994E+00    900.  0.281E+00  0.278E+00  0.295E-03  0.323E-13   13   6   3   13   4  37   13  36   7    0.000    0.000       0
     481  0    249 00:57:09.52  12.8000A  -2.0386   2.2491   3.5073  0.115  1.43  -0.986E+00    900.  0.280E+00  0.278E+00  0.303E-03  0.331E-13   13   6   3   13   5  37   14  35   6    0.000    0.000       0
     482  0    250 00:57:22.32  12.9000A   2.0738   2.1661   3.4923  0.115  1.43  -0.988E+00    900.  0.280E+00  0.277E+00  0.313E-03  0.343E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
     483  0    251 00:57:35.22  13.1000A   2.1173   2.2291   3.4375  0.115  1.43  -0.985E+00    900.  0.279E+00  0.277E+00  0.316E-03  0.347E-13   13  37   9   14   5  37   14  35   6    0.000    0.000       0
     484  0    252 00:57:48.32  13.4000A   2.1033   2.2274   3.3545  0.115  1.43  -0.989E+00    900.  0.279E+00  0.276E+00  0.321E-03  0.353E-13   13  37   9   14   6  37   15  35   6    0.000    0.000       0
     485  0    253 00:58:01.72  13.7000A   2.2294   2.2239   3.2838  0.115  1.46  -0.996E+00    950.  0.279E+00  0.276E+00  0.329E-03  0.362E-13   14  37  10   14   6  37   15  35   6    0.000    0.000       0
     486  0    254 00:58:15.42  13.8000A   2.2842   2.1522   3.2636  0.115  1.40  -0.993E+00    850.  0.278E+00  0.276E+00  0.337E-03  0.371E-13   14  37  10   14   6  37   16  36  10    0.000    0.000       0
     487  0    255 00:58:29.22  14.3000A   2.2566   2.1277   3.1445  0.115  1.40  -0.991E+00    850.  0.278E+00  0.276E+00  0.340E-03  0.375E-13   14  37  10   15   7  37    8  37   5    0.000    0.000       0
     488  0    256 00:58:43.52  14.3000A   2.2013   2.1755   3.1545  0.115  1.38  -0.993E+00    800.  0.278E+00  0.275E+00  0.352E-03  0.389E-13   15  37  11   15   7  37    8  37   5    0.000    0.000       0
     489  0    257 00:58:57.82  14.2000A   2.2568   2.1866   3.1704  0.115  1.38  -0.993E+00    800.  0.278E+00  0.275E+00  0.353E-03  0.389E-13   15  37  11   15   8  37    9  37   5    0.000    0.000       0
     490  0    258 00:59:12.02  14.2000A   2.2312   2.1716   3.1684  0.115  1.38  -0.995E+00    800.  0.278E+00  0.275E+00  0.351E-03  0.387E-13   15  37  11   15   9  37    9  37   5    0.000    0.000       0
     491  0    259 00:59:26.22  14.4000A   2.1733   2.1525   3.1350  0.115  1.38  -0.999E+00    800.  0.278E+00  0.276E+00  0.351E-03  0.387E-13   15  37  11   15   9  37    8  37   5    0.000    0.000       0
     492  0    260 00:59:40.62  13.8000A   2.0889   2.1539   3.2687  0.115  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.356E-03  0.392E-13   14  38  11   15   7  36    7   9  15    0.000    0.000       0
     493  0    261 00:59:54.42  13.4000A   2.0555   2.1633   3.3538  0.115  1.38  -0.100E+01    800.  0.278E+00  0.276E+00  0.342E-03  0.375E-13   14  38  11   15   8  36    7   9  15    0.000    0.000       0
     494  0    262 01:00:07.82  12.8000A   1.9438   2.1487   3.5253  0.115  1.38  -0.997E+00    800.  0.278E+00  0.276E+00  0.332E-03  0.363E-13   17  36  12   15   8  36    8   9  15    0.000    0.000       0
  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/Makefile

    r2963 r3083  
    2525# -----------------
    2626# $Id$
     27# add turbulence_closure_mod to parin
     28#
     29# 2963 2018-04-12 14:47:44Z suehring
    2730# Introduce index for vegetation/wall, pavement/green-wall and water/window
    2831# surfaces, for clearer access of surface fraction, albedo, emissivity, etc. .
     
    12301233        spectra_mod.o \
    12311234        synthetic_turbulence_generator_mod.o \
     1235        turbulence_closure_mod.o \
    12321236        uv_exposure_model_mod.o \
    12331237        vertical_nesting_mod.o \
  • palm/trunk/SOURCE/check_parameters.f90

    r3065 r3083  
    2424! Former revisions:
    2525! -----------------
    26 ! $Id$
     26! $Id: check_parameters.f90 2520 2017-10-05 13:50:26Z gronemeier &
     27! Add inital profile output for e (TG)
     28!
     29! 3065 2018-06-12 07:03:02Z Giersch
    2730! dz was replaced by dz(1), error message revised
    2831!
     
    24552458             ENDIF
    24562459
    2457           CASE ( 'e' )
     2460          CASE ( 'e', '#e' )
    24582461             dopr_index(i)  = 8
    24592462             dopr_unit(i)   = 'm2/s2'
    24602463             hom(:,2,8,:)   = SPREAD( zu, 2, statistic_regions+1 )
    24612464             hom(nzb,2,8,:) = 0.0_wp
     2465             IF ( data_output_pr(i)(1:1) == '#' )  THEN
     2466                dopr_initial_index(i) = 8
     2467                hom(:,2,8,:)          = SPREAD( zu, 2, statistic_regions+1 )
     2468                data_output_pr(i)     = data_output_pr(i)(2:)
     2469             ENDIF
    24622470
    24632471          CASE ( 'km', '#km' )
  • palm/trunk/SOURCE/data_output_mask.f90

    r3045 r3083  
    2525! -----------------
    2626! $Id$
     27!
     28!
     29! 3045 2018-05-28 07:55:41Z Giersch
    2730! Error messages revised
    2831!
  • palm/trunk/SOURCE/header.f90

    r3065 r3083  
    2525! -----------------
    2626! $Id$
     27! Print RANS-mode constants
     28!
     29! 3065 2018-06-12 07:03:02Z Giersch
    2730! Header output concerning stretching revised
    2831!
     
    450453    USE synthetic_turbulence_generator_mod,                                    &
    451454        ONLY:  stg_header
     455
     456    USE turbulence_closure_mod,                                                &
     457        ONLY:  rans_const_c, rans_const_sigma
    452458
    453459    IMPLICIT NONE
     
    19461952       IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
    19471953    ENDIF
    1948 
     1954    IF ( rans_mode )  THEN
     1955       WRITE ( io, 457 )  rans_const_c, rans_const_sigma
     1956    ENDIF
    19491957!
    19501958!-- Special actions during the run
     
    23882396456 FORMAT ('    Day of the year at model start :   day_init = ',I3             &
    23892397            /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
     2398457 FORMAT ('    RANS-mode constants: c_0 = ',F9.5/         &
     2399            '                         c_1 = ',F9.5/         &
     2400            '                         c_2 = ',F9.5/         &
     2401            '                         c_3 = ',F9.5/         &
     2402            '                         c_4 = ',F9.5/         &
     2403            '                         sigma_e    = ',F9.5/  &
     2404            '                         sigma_diss = ',F9.5)
    23902405470 FORMAT (//' Actions during the simulation:'/ &
    23912406              ' -----------------------------'/)
  • palm/trunk/SOURCE/init_3d_model.f90

  • palm/trunk/SOURCE/model_1d_mod.f90

    r3049 r3083  
    2525! -----------------
    2626! $Id$
     27! Bugfixes:
     28!   - preset te_diss and te_e to avoid runtime errors
     29!   - implementation of buoyancy term to dissipation
     30!     according to Sogachev et al. (2012)
     31!   - where diss_p < 0 set diss_p = 0.1 diss
     32!   - calculate progn eq(diss) starting from nzb+1
     33! Changes:
     34!   - add sig_e to TKE equation
     35!   - adjust prognostic equation of diss
     36!   - set model constants according to Koblitz (2013)
     37!   - renamed c_m to c_0
     38!   - rename l_black into l1d_init
     39!   - calculate l_grid within init_1d_model and save it as l1d_init
     40!   - calculate l1d according to DE85 if dissipation is a prognostic value
     41!   - made annotations doxygen-readable
     42!
     43! 3049 2018-05-29 13:52:36Z Giersch
    2744! Error messages revised
    2845!
     
    164181    USE kinds
    165182
    166     USE pegrid
     183    USE pegrid,                                                                &
     184        ONLY:  myid
    167185       
    168186
     
    175193    LOGICAL ::  stop_dt_1d = .FALSE.             !< termination flag, used in case of too small timestep (1d-model)
    176194
    177     REAL(wp) ::  c_1 = 1.44_wp                 !< model constant
    178     REAL(wp) ::  c_2 = 1.92_wp                 !< model constant
    179     REAL(wp) ::  c_3 = 1.44_wp                 !< model constant
    180     REAL(wp) ::  c_h = 0.0015_wp               !< model constant according to Detering and Etling (1985)
    181     REAL(wp) ::  c_m = 0.4_wp                  !< model constant, 0.4 according to Detering and Etling (1985)
    182     REAL(wp) ::  c_mu = 0.09_wp                !< model constant
     195    REAL(wp) ::  alpha_buoyancy                !< model constant according to Koblitz (2013)
     196    REAL(wp) ::  c_0 = 0.03_wp**0.25_wp        !< model constant according to Koblitz (2013)
     197    REAL(wp) ::  c_1 = 1.52_wp                 !< model constant according to Koblitz (2013)
     198    REAL(wp) ::  c_2 = 1.83_wp                 !< model constant according to Koblitz (2013)
     199    REAL(wp) ::  c_3                           !< model constant
     200    REAL(wp) ::  c_mu                          !< model constant
    183201    REAL(wp) ::  damp_level_1d = -1.0_wp       !< namelist parameter
    184202    REAL(wp) ::  dt_1d = 60.0_wp               !< dynamic timestep (1d-model)
     
    187205    REAL(wp) ::  dt_run_control_1d = 60.0_wp   !< namelist parameter
    188206    REAL(wp) ::  end_time_1d = 864000.0_wp     !< namelist parameter
     207    REAL(wp) ::  lambda                        !< maximum mixing length
    189208    REAL(wp) ::  qs1d                          !< characteristic humidity scale (1d-model)
    190209    REAL(wp) ::  simulated_time_1d = 0.0_wp    !< updated simulated time (1d-model)
    191     REAL(wp) ::  sig_diss = 1.3_wp             !< model constant
     210    REAL(wp) ::  sig_diss = 2.95_wp            !< model constant according to Koblitz (2013)
     211    REAL(wp) ::  sig_e = 2.95_wp               !< model constant according to Koblitz (2013)
    192212    REAL(wp) ::  time_pr_1d = 0.0_wp           !< updated simulated time for profile output (1d-model)
    193213    REAL(wp) ::  time_run_control_1d = 0.0_wp  !< updated simulated time for run-control output (1d-model)
     
    198218    REAL(wp) ::  z01d                          !< roughness length for momentum (1d-model)
    199219    REAL(wp) ::  z0h1d                         !< roughness length for scalars (1d-model)
    200 
    201220
    202221    REAL(wp), DIMENSION(:), ALLOCATABLE ::  diss1d   !< tke dissipation rate (1d-model)
     
    279298
    280299    INTEGER(iwp) ::  k  !< loop index
    281    
    282     REAL(wp) ::  lambda !< maximum mixing length
    283300
    284301!
     
    324341!
    325342!--       Use the same mixing length as in 3D model (LES-mode)
    326           !@todo: rename (delete?) this option
    327           ! As the mixing length is different between RANS and LES mode, it
    328           ! must be distinguished here between these modes. For RANS mode,
    329           ! the mixing length is calculated accoding to Blackadar, which is
    330           ! the other option at this point.
    331           ! Maybe delete this option entirely (not appropriate in LES case)
    332           ! 2018-03-20, gronemeier
     343          !> @todo rename (delete?) this option
     344          !> As the mixing length is different between RANS and LES mode, it
     345          !> must be distinguished here between these modes. For RANS mode,
     346          !> the mixing length is calculated accoding to Blackadar, which is
     347          !> the other option at this point.
     348          !> Maybe delete this option entirely (not appropriate in LES case)
     349          !> 2018-03-20, gronemeier
    333350          DO  k = nzb+1, nzt
    334351             l1d_init(k)  = ( dx * dy * dzw(k) )**0.33333333333333_wp
     
    359376       us1d = 0.1_wp   ! without initial friction the flow would not change
    360377    ELSE
    361        diss1d(nzb+1) = 1.0_wp
     378       diss1d(nzb+1) = 0.001_wp
    362379       e1d(nzb+1)  = 1.0_wp
    363380       km1d(nzb+1) = 1.0_wp
     
    372389
    373390!
    374 !-- Tendencies must be preset in order to avoid runtime errors within the
    375 !-- first Runge-Kutta step
     391!-- Tendencies must be preset in order to avoid runtime errors
     392    te_diss  = 0.0_wp
    376393    te_dissm = 0.0_wp
     394    te_e  = 0.0_wp
    377395    te_em = 0.0_wp
    378396    te_um = 0.0_wp
     
    381399!
    382400!-- Set model constant
    383     IF ( dissipation_1d == 'as_in_3d_model' )  c_m = 0.1_wp
     401    IF ( dissipation_1d == 'as_in_3d_model' )  c_0 = 0.1_wp
     402    c_mu = c_0**4
    384403
    385404!
     
    423442!-- Determine the time step at the start of a 1D-simulation and
    424443!-- determine and printout quantities used for run control
    425     dt_1d = 1.0_wp
     444    dt_1d = 0.01_wp
    426445    CALL run_control_1d
    427446
     
    483502!--             dissipation rate
    484503                IF ( dissipation_1d == 'detering' )  THEN
    485                    diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
     504                   diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
    486505                ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
    487506                   diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) &
     
    497516                                     kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1)  &
    498517                                   - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k)    &
    499                                                  ) * ddzw(k)                   &
     518                                                 ) * ddzw(k) / sig_e           &
    500519                                   - diss1d(k)
    501520
     
    503522!
    504523!--                dissipation rate
    505                    te_diss(k) = km1d(k) *                                      &
     524                   IF ( rif1d(k) >= 0.0_wp )  THEN
     525                      alpha_buoyancy = 1.0_wp - l1d(k) / lambda
     526                   ELSE
     527                      alpha_buoyancy = 1.0_wp - ( 1.0_wp + ( c_2 - 1.0_wp )    &
     528                                                         / ( c_2 - c_1    ) )  &
     529                                              * l1d(k) / lambda
     530                   ENDIF
     531                   c_3 = ( c_1 - c_2 ) * alpha_buoyancy
     532                   te_diss(k) = ( km1d(k) *                                    &
    506533                                  ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2  &
    507534                                  + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2  &
    508                                   ) * c_1 * c_mu**0.75 / c_h * f / us1d        &
    509                                     * SQRT(e1d(k))                             &
     535                                  ) * ( c_1 + (c_2 - c_1) * l1d(k) / lambda )  &
    510536                                  - g / pt_0 * kh1d(k) * flux * c_3            &
    511                                     * diss1d(k) / ( e1d(k) + 1.0E-20_wp )      &
    512                                   + ( kmzp * ( diss1d(k+1) - diss1d(k) )       &
     537                                  - c_2 * diss1d(k)                            &
     538                                ) * diss1d(k) / ( e1d(k) + 1.0E-20_wp )        &
     539                                + (   kmzp * ( diss1d(k+1) - diss1d(k) )       &
    513540                                           * ddzu(k+1)                         &
    514541                                    - kmzm * ( diss1d(k) - diss1d(k-1) )       &
    515542                                           * ddzu(k)                           &
    516                                     ) * ddzw(k) / sig_diss                     &
    517                                   - c_2 * diss1d(k)**2 / ( e1d(k) + 1.0E-20_wp )
     543                                  ) * ddzw(k) / sig_diss
    518544
    519545                ENDIF
     
    546572!--          dissipation rate
    547573             IF ( dissipation_1d == 'detering' )  THEN
    548                 diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
     574                diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
    549575             ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
    550576                diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) )  &
     
    565591!--          TKE
    566592             IF ( .NOT. dissipation_1d == 'prognostic' )  THEN
     593                !> @query why integrate over 2dz
     594                !>   Why is it allowed to integrate over two delta-z for e
     595                !>   while for u and v it is not?
     596                !>   2018-04-23, gronemeier
    567597                te_e(k) = km1d(k) * ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2&
    568598                                    + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2&
     
    572602                                     kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1)  &
    573603                                   - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k)    &
    574                                                  ) * ddzw(k)                   &
     604                                                 ) * ddzw(k) / sig_e           &
    575605                                   - diss1d(k)
    576606             ENDIF
     
    595625             ENDDO
    596626
    597              IF ( dissipation_1d == 'prognostic' )  THEN
    598                 DO  k = nzb_diff, nzt
    599                    diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + &
    600                                                  tsc(3) * te_dissm(k) )
    601                 ENDDO
    602              ENDIF
    603627!
    604628!--          Eliminate negative TKE values, which can result from the
     
    606630!--          value is reduced to 10 percent of its old value.
    607631             WHERE ( e1d_p < 0.0_wp )  e1d_p = 0.1_wp * e1d
     632
     633             IF ( dissipation_1d == 'prognostic' )  THEN
     634                DO  k = nzb+1, nzt
     635                   diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + &
     636                                                       tsc(3) * te_dissm(k) )
     637                ENDDO
     638                WHERE ( diss1d_p < 0.0_wp )  diss1d_p = 0.1_wp * diss1d
     639             ENDIF
    608640          ENDIF
    609641
     
    643675                   IF ( dissipation_1d == 'prognostic' )  THEN
    644676                      DO k = nzb+1, nzt
    645                          te_dissm(k) = -9.5625_wp * te_diss(k) + 5.3125_wp * te_dissm(k)
     677                         te_dissm(k) = -9.5625_wp * te_diss(k)  &
     678                                     +  5.3125_wp * te_dissm(k)
    646679                      ENDDO
    647680                   ENDIF
     
    650683             ENDIF
    651684          ENDIF
    652 
    653685
    654686!
     
    701733
    702734             ENDIF    ! constant_flux_layer
    703 
     735             !> @todo combine if clauses
     736             !>   The previous and following if clauses can be combined into a
     737             !>   single clause
     738             !>   2018-04-23, gronemeier
    704739!
    705740!--          Compute the Richardson-flux numbers,
     
    789824!--             compatibility with the 3D model.
    790825                IF ( ibc_e_b == 2 )  THEN
    791                    e1d(nzb+1) = ( us1d / c_m )**2
     826                   e1d(nzb+1) = ( us1d / c_0 )**2
    792827                ENDIF
    793828                IF ( dissipation_1d == 'prognostic' )  THEN
    794                    e1d(nzb+1) = us1d**2 / SQRT( c_mu )
     829                   e1d(nzb+1) = ( us1d / c_0 )**2
    795830                   diss1d(nzb+1) = us1d**3 / ( kappa * zu(nzb+1) )
    796831                   diss1d(nzb) = diss1d(nzb+1)
     
    829864!--          in the dissipation of TKE via l1d_diss. Otherwise, km1d would be
    830865!--          too large.
    831              IF ( mixing_length_1d == 'blackadar' )  THEN
     866             IF ( dissipation_1d /= 'prognostic' )  THEN
     867                IF ( mixing_length_1d == 'blackadar' )  THEN
     868                   DO  k = nzb+1, nzt
     869                      IF ( rif1d(k) >= 0.0_wp )  THEN
     870                         l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
     871                         l1d_diss(k) = l1d(k)
     872                      ELSE
     873                         l1d(k) = l1d_init(k)
     874                         l1d_diss(k) = l1d_init(k) *                           &
     875                                       SQRT( 1.0_wp - 16.0_wp * rif1d(k) )
     876                      ENDIF
     877                   ENDDO
     878                ELSEIF ( mixing_length_1d == 'as_in_3d_model' )  THEN
     879                   DO  k = nzb+1, nzt
     880                      dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k)
     881                      IF ( dpt_dz > 0.0_wp )  THEN
     882                         l_stable = 0.76_wp * SQRT( e1d(k) )                   &
     883                                  / SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp
     884                      ELSE
     885                         l_stable = l1d_init(k)
     886                      ENDIF
     887                      l1d(k) = MIN( l1d_init(k), l_stable )
     888                      l1d_diss(k) = l1d(k)
     889                   ENDDO
     890                ENDIF
     891             ELSE
    832892                DO  k = nzb+1, nzt
    833                    IF ( rif1d(k) >= 0.0_wp )  THEN
    834                       l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
    835                       l1d_diss(k) = l1d(k)
    836                    ELSE
    837                       l1d(k) = l1d_init(k)
    838                       l1d_diss(k) = l1d_init(k) *                              &
    839                                     SQRT( 1.0_wp - 16.0_wp * rif1d(k) )
    840                    ENDIF
    841                 ENDDO
    842              ELSEIF ( mixing_length_1d == 'as_in_3d_model' )  THEN
    843                 DO  k = nzb+1, nzt
    844                    dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k)
    845                    IF ( dpt_dz > 0.0_wp )  THEN
    846                       l_stable = 0.76_wp * SQRT( e1d(k) ) /                    &
    847                                      SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp
    848                    ELSE
    849                       l_stable = l1d_init(k)
    850                    ENDIF
    851                    l1d(k) = MIN( l1d_init(k), l_stable )
    852                    l1d_diss(k) = l1d(k)
     893                   l1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) )                   &
     894                          / ( diss1d(k) + 1.0E-30_wp )
    853895                ENDDO
    854896             ENDIF
     
    867909                ENDIF
    868910             ENDIF
     911
    869912             IF ( dissipation_1d == 'prognostic' )  THEN
    870913                DO  k = nzb_diff, nzt
     
    873916             ELSE
    874917                DO  k = nzb_diff, nzt
    875                    km1d(k) = c_m * SQRT( e1d(k) ) * l1d(k)
     918                   km1d(k) = c_0 * SQRT( e1d(k) ) * l1d(k)
    876919                ENDDO
    877920             ENDIF
     
    10181061
    10191062    REAL(wp) ::  dt_diff  !< time step accorind to diffusion criterion
     1063    REAL(wp) ::  dt_old   !< previous time step
    10201064    REAL(wp) ::  fac      !< factor of criterion
    10211065    REAL(wp) ::  value    !< auxiliary variable
    10221066
    10231067!
     1068!-- Save previous time step
     1069    dt_old = dt_1d
     1070
     1071!
    10241072!-- Compute the currently feasible time step according to the diffusion
    10251073!-- criterion. At nzb+1 the half grid length is used.
    1026     fac = 0.125  !0.35_wp                                                       !### changed from 0.35
     1074    fac = 0.125
    10271075    dt_diff = dt_max_1d
    10281076    DO  k = nzb+2, nzt
     
    10341082
    10351083!
     1084!-- Limit the new time step to a maximum of 10 times the previous time step
     1085    dt_1d = MIN( dt_old * 10.0_wp, dt_1d )
     1086
     1087!
    10361088!-- Set flag when the time step becomes too small
    1037     IF ( dt_1d < ( 0.00001_wp * dt_max_1d ) )  THEN
     1089    IF ( dt_1d < ( 1.0E-15_wp * dt_max_1d ) )  THEN
    10381090       stop_dt_1d = .TRUE.
    10391091
  • palm/trunk/SOURCE/modules.f90

    r3065 r3083  
    2525! -----------------
    2626! $Id$
     27! set dt_3d = 0.01
     28!
     29! 3065 2018-06-12 07:03:02Z Giersch
    2730! Variables concerning stretching introduced or revised
    2831!
     
    14221425    REAL(wp) ::  dt_run_control = 60.0_wp                      !< namelist parameter
    14231426    REAL(wp) ::  dt_spinup = 60.0_wp                           !< namelist parameter
    1424     REAL(wp) ::  dt_3d = 1.0_wp                                !< time step
     1427    REAL(wp) ::  dt_3d = 0.01_wp                               !< time step
    14251428    REAL(wp) ::  dz_max = 1000.0_wp                            !< namelist parameter
    14261429    REAL(wp) ::  dz_stretch_factor = 1.08_wp                   !< namelist parameter
  • palm/trunk/SOURCE/parin.f90

    r3065 r3083  
    2525! -----------------
    2626! $Id$
     27! Added rans_const_c and rans_const_sigma as input parameters (TG)
     28!
     29! 3065 2018-06-12 07:03:02Z Giersch
    2730! New initialization parameters added
    2831!
     
    467470    USE synthetic_turbulence_generator_mod,                                    &
    468471        ONLY:  stg_parin
     472
     473    USE turbulence_closure_mod,                                                &
     474        ONLY:  rans_const_c, rans_const_sigma
    469475
    470476    USE urban_surface_mod,                                                     &
     
    531537             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
    532538             q_vertical_gradient, q_vertical_gradient_level,                   &
    533              random_generator, random_heatflux, rans_mode,                     &
     539             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
     540             rans_mode,                                                        &
    534541             rayleigh_damping_factor, rayleigh_damping_height,                 &
    535542             recycling_width, recycling_yshift,                                &
     
    601608             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
    602609             q_vertical_gradient, q_vertical_gradient_level,                   &
    603              random_generator, random_heatflux, rans_mode,                     &
     610             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
     611             rans_mode,                                                        &
    604612             rayleigh_damping_factor, rayleigh_damping_height,                 &
    605613             recycling_width, recycling_yshift,                                &
  • palm/trunk/SOURCE/pmc_interface_mod.f90

    • Property svn:mergeinfo deleted
  • palm/trunk/SOURCE/timestep.f90

    r3049 r3083  
    2525! -----------------
    2626! $Id$
     27! limit dt_3d to be at maximum 2*old_dt; define old_dt at beginning of routine
     28! Add km/kh_max
     29!
     30! 3049 2018-05-29 13:52:36Z Giersch
    2731! Error messages revised
    2832!
     
    161165    INTEGER(iwp) ::  j !<
    162166    INTEGER(iwp) ::  k !<
     167    INTEGER(iwp) ::  km_max_ijk(3) = -1  !< index values (i,j,k) of location where km_max occurs
     168    INTEGER(iwp) ::  kh_max_ijk(3) = -1  !< index values (i,j,k) of location where kh_max occurs
    163169
    164170    LOGICAL ::  stop_dt_local !< local switch for controlling the time stepping
     
    173179    REAL(wp) ::  dt_w              !<
    174180    REAL(wp) ::  dt_w_l            !<
     181    REAL(wp) ::  km_max            !< maximum of Km in entire domain
     182    REAL(wp) ::  kh_max            !< maximum of Kh in entire domain
    175183    REAL(wp) ::  u_gtrans_l        !<
    176184    REAL(wp) ::  u_max_l           !<
     
    190198
    191199
    192 
    193200    CALL cpu_log( log_point(12), 'calculate_timestep', 'start' )
     201!
     202!--    Save former time step as reference
     203       old_dt = dt_3d
    194204
    195205!
     
    318328!--    The time step is the minimum of the 3-4 components and the diffusion time
    319329!--    step minus a reduction (cfl_factor) to be on the safe side.
    320 !--    The time step must not exceed the maximum allowed value.
     330!--    The time step must not exceed the maximum allowed value and must not
     331!--    increase by more than a factor of 10.
    321332       dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w, dt_precipitation )
    322        dt_3d = MIN( dt_3d, dt_max )
     333       dt_3d = MIN( dt_3d, dt_max, 2.0_wp * old_dt )
    323334
    324335!
     
    334345       IF ( dt_3d < ( 0.00001_wp * dt_max ) )  THEN
    335346          stop_dt = .TRUE.
     347
     348!
     349!--       Determine the maxima of the diffusion coefficients, including their
     350!--       grid index positions.
     351          CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, km, 'abs',  &
     352                               0.0_wp, km_max, km_max_ijk )
     353          CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, kh, 'abs',  &
     354                               0.0_wp, kh_max, kh_max_ijk )
    336355
    337356          WRITE( message_string, * ) 'Time step has reached minimum limit.',   &
     
    342361               '&dt_w            = ', dt_w, ' s',                              &
    343362               '&dt_diff         = ', dt_diff, ' s',                           &
    344                '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),        &
     363               '&u_max           = ', u_max, ' m/s    k=', u_max_ijk(1),       &
    345364               '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                     &
    346                '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),        &
     365               '&v_max           = ', v_max, ' m/s    k=', v_max_ijk(1),       &
    347366               '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                     &
    348                '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),        &
    349                '  j=', w_max_ijk(2), '  i=', w_max_ijk(3)
     367               '&w_max           = ', w_max, ' m/s    k=', w_max_ijk(1),       &
     368               '  j=', w_max_ijk(2), '  i=', w_max_ijk(3),                     &
     369               '&km_max          = ', km_max, ' m2/s2  k=', km_max_ijk(1),     &
     370               '  j=', km_max_ijk(2), '  i=', km_max_ijk(3),                   &
     371               '&kh_max          = ', kh_max, ' m2/s2  k=', kh_max_ijk(1),     &
     372                '  j=', kh_max_ijk(2), '  i=', kh_max_ijk(3)
    350373          CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 )
    351374!
     
    388411       dt_3d = NINT( dt_3d * 100.0_wp / div ) * div / 100.0_wp
    389412
    390 !
    391 !--    Adjust the time step
    392        old_dt = dt_3d
    393 
    394413    ENDIF
    395414
  • palm/trunk/SOURCE/turbulence_closure_mod.f90

    r3045 r3083  
    2525! -----------------
    2626! $Id$
     27! - set limits of diss at the end of prognostic equations
     28! - call production_e to calculate production term of diss
     29! - limit change of diss to -90% to +100%
     30! - remove factor 0.5 from diffusion_diss_ij
     31! - rename c_m into c_0, and c_h into c_4
     32! - add rans_const_c and rans_const_sigma as namelist parameters
     33! - add calculation of mixing length for profile output in case of rans_tke_e
     34! - changed format of annotations to comply with doxygen standards
     35! - calculate and save dissipation rate during rans_tke_l mode
     36! - set bc at vertical walls for e, diss, km, kh
     37! - bugfix: set l_wall = 0.0 within buildings
     38! - set l_wall at bottom and top boundary (rans-mode)
     39! - bugfix in production term for dissipation rate
     40! - bugfix in diffusion of dissipation rate
     41! - disable check for 1D model if rans_tke_e is used
     42! - bugfixes for initialization (rans-mode):
     43!    - correction of dissipation-rate formula
     44!    - calculate km based on l_wall
     45!    - initialize diss if 1D model is not used
     46!
     47! 3045 2018-05-28 07:55:41Z Giersch
    2748! Error message revised
    2849!
     
    3657! Further todo's
    3758!
    38 ! 2936 2018-03-27 14:49:27Z suehring
     59! 2936 2018-03-27 14:49:27Z gronemeier
    3960! - defined l_grid only within this module
    4061! - Moved l_wall definition from modules.f90
     
    81102!>       add OpenMP directives whereever possible
    82103!>       remove debug output variables (dummy1, dummy2, dummy3)
    83 !> @todo Move initialization of wall-mixing length from init_grid
    84104!> @todo Check for random disturbances
    85105!> @note <Enter notes on the module>
     
    146166
    147167
    148     REAL(wp) ::  c_1  = 1.44_wp    !< model constant for RANS mode
    149     REAL(wp) ::  c_2  = 1.92_wp    !< model constant for RANS mode
    150     REAL(wp) ::  c_3  = 1.44_wp    !< model constant for RANS mode
    151     REAL(wp) ::  c_h  = 0.0015_wp  !< model constant for RANS mode
    152     REAL(wp) ::  c_m               !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES)
    153     REAL(wp) ::  c_mu = 0.09_wp    !< model constant for RANS mode
    154     REAL(wp) ::  l_max             !< maximum length scale for Blackadar mixing length
    155     REAL(wp) ::  sig_e = 1.0_wp    !< factor to calculate Ke from Km
    156     REAL(wp) ::  sig_diss = 1.3_wp !< factor to calculate K_diss from Km
    157     INTEGER(iwp) ::  surf_e        !< end index of surface elements at given i-j position
    158     INTEGER(iwp) ::  surf_s        !< start index of surface elements at given i-j position
     168    REAL(wp) ::  c_0                !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES)
     169    REAL(wp) ::  c_1                !< model constant for RANS mode
     170    REAL(wp) ::  c_2                !< model constant for RANS mode
     171    REAL(wp) ::  c_3                !< model constant for RANS mode
     172    REAL(wp) ::  c_4                !< model constant for RANS mode
     173    REAL(wp) ::  l_max              !< maximum length scale for Blackadar mixing length
     174    REAL(wp) ::  dsig_e = 1.0_wp    !< factor to calculate Ke from Km (1/sigma_e)
     175    REAL(wp) ::  dsig_diss = 1.0_wp !< factor to calculate K_diss from Km (1/sigma_diss)
     176    INTEGER(iwp) ::  surf_e         !< end index of surface elements at given i-j position
     177    INTEGER(iwp) ::  surf_s         !< start index of surface elements at given i-j position
     178
     179    REAL(wp), DIMENSION(0:4) :: rans_const_c = &       !< model constants for RANS mode (namelist param)
     180       (/ 0.55_wp, 1.44_wp, 1.92_wp, 0.0_wp, 0.0_wp /) !> default values fit for standard-tke-e closure
     181
     182    REAL(wp), DIMENSION(2) :: rans_const_sigma = &     !< model constants for RANS mode, sigma values (sigma_e, sigma_diss) (namelist param)
     183       (/ 1.0_wp, 0.77_wp /)
    159184
    160185    REAL(wp), DIMENSION(:), ALLOCATABLE ::  l_black    !< mixing length according to Blackadar
     
    163188    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  l_wall !< near-wall mixing length
    164189
    165     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy1 !< debug output variable
    166     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy2 !< debug output variable
    167     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy3 !< debug output variable
    168 
    169 
    170     PUBLIC c_m, dummy1, dummy2, dummy3
     190    !> @todo remove debug variables
     191    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_prod1, diss_adve1, diss_diff1, &
     192                                               diss_prod2, diss_adve2, diss_diff2, &
     193                                               diss_prod3, diss_adve3, diss_diff3, &
     194                                               dummy1, dummy2, dummy3
     195
     196
     197    PUBLIC c_0, rans_const_c, rans_const_sigma
    171198
    172199!
     
    294321! ------------
    295322!> Check parameters routine for turbulence closure module.
     323!> @todo remove rans_mode from initialization namelist and rework checks
     324!>   The way it is implemented at the moment, the user has to set two variables
     325!>   so that the RANS mode is working. It would be better if only one parameter
     326!>   has to be set.
     327!>   2018-06-18, gronemeier
    296328!------------------------------------------------------------------------------!
    297329 SUBROUTINE tcm_check_parameters
     
    307339    IF ( rans_mode )  THEN
    308340
    309        c_m = 0.4_wp  !according to Detering and Etling (1985)
     341!
     342!--    Assign values to constants for RANS mode
     343       dsig_e    = 1.0_wp / rans_const_sigma(1)
     344       dsig_diss = 1.0_wp / rans_const_sigma(2)
     345
     346       c_0 = rans_const_c(0)
     347       c_1 = rans_const_c(1)
     348       c_2 = rans_const_c(2)
     349       !c_3 = rans_const_c(3)   !> @todo clarify how to switch between different models
     350       c_4 = rans_const_c(4)
    310351
    311352       SELECT CASE ( TRIM( turbulence_closure ) )
     
    316357          CASE ( 'TKE-e' )
    317358             rans_tke_e = .TRUE.
    318 
    319              IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) == 0  &
    320                   .AND.  .NOT.  nest_domain )  THEN
    321                 message_string = 'Initializing without 1D model while ' //     &
    322                                  'using TKE-e closure&' //                     &
    323                                  'is not possible at the moment!'
    324                 CALL message( 'tcm_check_parameters', 'TG0005', 1, 2, 0, 6, 0 )
    325              ENDIF
    326359
    327360          CASE DEFAULT
    328361             message_string = 'Unknown turbulence closure: ' //                &
    329362                              TRIM( turbulence_closure )
    330              CALL message( 'tcm_check_parameters', 'TG0001', 1, 2, 0, 6, 0 )
     363             CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 )
    331364
    332365       END SELECT
     366
     367       IF ( turbulent_inflow .OR. turbulent_outflow )  THEN
     368          message_string = 'turbulent inflow/outflow is not yet '//            &
     369                           'implemented for RANS mode'
     370          CALL message( 'tcm_check_parameters', 'PA0501', 1, 2, 0, 6, 0 )
     371       ENDIF
    333372
    334373       message_string = 'RANS mode is still in development! ' //               &
    335374                        '&Not all features of PALM are yet compatible '//      &
    336375                        'with RANS mode. &Use at own risk!'
    337        CALL message( 'tcm_check_parameters', 'TG0003', 0, 1, 0, 6, 0 )
     376       CALL message( 'tcm_check_parameters', 'PA0502', 0, 1, 0, 6, 0 )
    338377
    339378    ELSE
    340379
    341        c_m = 0.1_wp !according to Lilly (1967) and Deardorff (1980)
     380       c_0 = 0.1_wp !according to Lilly (1967) and Deardorff (1980)
     381
     382       dsig_e = 1.0_wp !assure to use K_m to calculate TKE instead
     383                       !of K_e which is used in RANS mode
    342384
    343385       SELECT CASE ( TRIM( turbulence_closure ) )
     
    347389
    348390          CASE DEFAULT
     391             !> @todo rework this part so that only one call of this error exists
    349392             message_string = 'Unknown turbulence closure: ' //                &
    350393                              TRIM( turbulence_closure )
    351              CALL message( 'tcm_check_parameters', 'TG0001', 1, 2, 0, 6, 0 )
     394             CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 )
    352395
    353396       END SELECT
    354 
    355     ENDIF
    356 
    357     IF ( rans_tke_e )  THEN
    358 
    359        IF ( turbulent_inflow .OR. turbulent_outflow )  THEN
    360           message_string = 'turbulent inflow/outflow is not yet '//            &
    361                            'implemented for TKE-e closure'
    362           CALL message( 'tcm_check_parameters', 'TG0002', 1, 2, 0, 6, 0 )
    363        ENDIF
    364397
    365398    ENDIF
     
    379412    IMPLICIT NONE
    380413
    381     CHARACTER (LEN=*) ::  unit     !<
    382     CHARACTER (LEN=*) ::  var      !<
    383 
    384     INTEGER(iwp) ::  i      !<
    385     INTEGER(iwp) ::  ilen   !<
    386     INTEGER(iwp) ::  k      !<
     414    CHARACTER (LEN=*) ::  unit     !< unit of output variable
     415    CHARACTER (LEN=*) ::  var      !< name of output variable
     416
     417    INTEGER(iwp) ::  i      !< index of var in data_output
     418    INTEGER(iwp) ::  ilen   !< length of var string
     419    INTEGER(iwp) ::  k      !< flag if var contains one of '_xy', '_xz' or '_yz'
    387420
    388421    SELECT CASE ( TRIM( var ) )
    389422
    390423       CASE ( 'diss' )
    391           IF ( .NOT.  rans_tke_e )  THEN
    392              message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    393                       'res TKE-e closure for RANS mode.'
    394              CALL message( 'tcm_check_data_output', 'TG0101', 1, 2, 0, 6, 0 )
    395           ENDIF
    396424          unit = 'm2/s3'
    397425
    398        CASE ( 'dummy2', 'dummy3', 'dummy1' )
    399           unit = 'mixing length'
     426       CASE ( 'diss1', 'diss2',                         &                      !> @todo remove later
     427              'diss_prod1', 'diss_adve1', 'diss_diff1', &
     428              'diss_prod2', 'diss_adve2', 'diss_diff2', &
     429              'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3'  )
     430          unit = 'debug output'
    400431
    401432       CASE ( 'kh', 'km' )
     
    420451    IMPLICIT NONE
    421452
    422     CHARACTER (LEN=*), INTENT(OUT) ::  grid_x   !<
    423     CHARACTER (LEN=*), INTENT(OUT) ::  grid_y   !<
    424     CHARACTER (LEN=*), INTENT(OUT) ::  grid_z   !<
    425     CHARACTER (LEN=*), INTENT(IN)  ::  var      !<
    426    
    427     LOGICAL, INTENT(OUT) ::  found   !<
    428    
     453    CHARACTER (LEN=*), INTENT(OUT) ::  grid_x   !< x grid of output variable
     454    CHARACTER (LEN=*), INTENT(OUT) ::  grid_y   !< y grid of output variable
     455    CHARACTER (LEN=*), INTENT(OUT) ::  grid_z   !< z grid of output variable
     456    CHARACTER (LEN=*), INTENT(IN)  ::  var      !< name of output variable
     457
     458    LOGICAL, INTENT(OUT) ::  found   !< flag if output variable is found
     459
    429460    found  = .TRUE.
    430461
     
    438469          grid_z = 'zu'
    439470
    440        CASE ( 'dummy2', 'dummy3', 'dummy1' )                                    !### remove later
     471       CASE ( 'diss1', 'diss2',                         &                       !> @todo remove later
     472              'diss_prod1', 'diss_adve1', 'diss_diff1', &
     473              'diss_prod2', 'diss_adve2', 'diss_diff2', &
     474              'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3' )
    441475          grid_x = 'x'
    442476          grid_y = 'y'
     
    480514    IMPLICIT NONE
    481515
    482     CHARACTER (LEN=*) ::  mode       !<
    483     CHARACTER (LEN=*) ::  variable   !<
    484 
    485     INTEGER(iwp) ::  i   !<
    486     INTEGER(iwp) ::  j   !<
    487     INTEGER(iwp) ::  k   !<
     516    CHARACTER (LEN=*) ::  mode       !< flag defining mode 'allocate', 'sum' or 'average'
     517    CHARACTER (LEN=*) ::  variable   !< name of variable
     518
     519    INTEGER(iwp) ::  i   !< loop index
     520    INTEGER(iwp) ::  j   !< loop index
     521    INTEGER(iwp) ::  k   !< loop index
    488522
    489523    IF ( mode == 'allocate' )  THEN
     
    616650    IMPLICIT NONE
    617651
    618     CHARACTER (LEN=*) ::  grid       !<
    619     CHARACTER (LEN=*) ::  mode       !<
    620     CHARACTER (LEN=*) ::  variable   !<
    621 
    622     INTEGER(iwp) ::  av   !<
    623     INTEGER(iwp) ::  i    !<
    624     INTEGER(iwp) ::  j    !<
    625     INTEGER(iwp) ::  k    !<
    626     INTEGER(iwp) ::  nzb_do   !<
    627     INTEGER(iwp) ::  nzt_do   !<
    628 
    629     LOGICAL ::  found   !<
     652    CHARACTER (LEN=*) ::  grid       !< name of vertical grid
     653    CHARACTER (LEN=*) ::  mode       !< either 'xy', 'xz' or 'yz'
     654    CHARACTER (LEN=*) ::  variable   !< name of variable
     655
     656    INTEGER(iwp) ::  av   !< flag for (non-)average output
     657    INTEGER(iwp) ::  i    !< loop index
     658    INTEGER(iwp) ::  j    !< loop index
     659    INTEGER(iwp) ::  k    !< loop index
     660    INTEGER(iwp) ::  nzb_do   !< vertical output index (bottom)
     661    INTEGER(iwp) ::  nzt_do   !< vertical output index (top)
     662
     663    LOGICAL ::  found   !< flag if output variable is found
    630664    LOGICAL ::  two_d   !< flag parameter that indicates 2D variables (horizontal cross sections)
    631665
     
    737771    IMPLICIT NONE
    738772
    739     CHARACTER (LEN=*) ::  variable   !<
    740 
    741     INTEGER(iwp) ::  av     !<
    742     INTEGER(iwp) ::  i      !<
    743     INTEGER(iwp) ::  j      !<
    744     INTEGER(iwp) ::  k      !<
     773    CHARACTER (LEN=*) ::  variable   !< name of variable
     774
     775    INTEGER(iwp) ::  av     !< flag for (non-)average output
     776    INTEGER(iwp) ::  i      !< loop index
     777    INTEGER(iwp) ::  j      !< loop index
     778    INTEGER(iwp) ::  k      !< loop index
    745779    INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
    746780    INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)
    747781
    748     LOGICAL ::  found   !<
     782    LOGICAL ::  found   !< flag if output variable is found
    749783
    750784    REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
     
    829863          ENDIF
    830864
    831        CASE ( 'dummy1' )                                                        !### remove later
     865       CASE ( 'dummy3' )                                                        !> @todo remove later
     866          IF ( av == 0 )  THEN
     867             DO  i = nxl, nxr
     868                DO  j = nys, nyn
     869                   DO  k = nzb_do, nzt_do
     870                      local_pf(i,j,k) = dummy3(k,j,i)
     871                   ENDDO
     872                ENDDO
     873             ENDDO
     874          ENDIF
     875
     876       CASE ( 'diss1' )                                                         !> @todo remove later
    832877          IF ( av == 0 )  THEN
    833878             DO  i = nxl, nxr
     
    840885          ENDIF
    841886
    842        CASE ( 'dummy2' )                                                        !### remove later
     887       CASE ( 'diss2' )                                                         !> @todo remove later
    843888          IF ( av == 0 )  THEN
    844889             DO  i = nxl, nxr
     
    851896          ENDIF
    852897
    853        CASE ( 'dummy3' )                                                        !### remove later
     898       CASE ( 'diss_prod1' )                                                    !> @todo remove later
    854899          IF ( av == 0 )  THEN
    855900             DO  i = nxl, nxr
    856901                DO  j = nys, nyn
    857902                   DO  k = nzb_do, nzt_do
    858                       local_pf(i,j,k) = dummy3(k,j,i)
     903                      local_pf(i,j,k) = diss_prod1(k,j,i)
    859904                   ENDDO
    860905                ENDDO
     
    862907          ENDIF
    863908
     909       CASE ( 'diss_adve1' )                                                    !> @todo remove later
     910          IF ( av == 0 )  THEN
     911             DO  i = nxl, nxr
     912                DO  j = nys, nyn
     913                   DO  k = nzb_do, nzt_do
     914                      local_pf(i,j,k) = diss_adve1(k,j,i)
     915                   ENDDO
     916                ENDDO
     917             ENDDO
     918          ENDIF
     919
     920       CASE ( 'diss_diff1' )                                                    !> @todo remove later
     921          IF ( av == 0 )  THEN
     922             DO  i = nxl, nxr
     923                DO  j = nys, nyn
     924                   DO  k = nzb_do, nzt_do
     925                      local_pf(i,j,k) = diss_diff1(k,j,i)
     926                   ENDDO
     927                ENDDO
     928             ENDDO
     929          ENDIF
     930
     931       CASE ( 'diss_prod2' )                                                    !> @todo remove later
     932          IF ( av == 0 )  THEN
     933             DO  i = nxl, nxr
     934                DO  j = nys, nyn
     935                   DO  k = nzb_do, nzt_do
     936                      local_pf(i,j,k) = diss_prod2(k,j,i)
     937                   ENDDO
     938                ENDDO
     939             ENDDO
     940          ENDIF
     941
     942       CASE ( 'diss_adve2' )                                                    !> @todo remove later
     943          IF ( av == 0 )  THEN
     944             DO  i = nxl, nxr
     945                DO  j = nys, nyn
     946                   DO  k = nzb_do, nzt_do
     947                      local_pf(i,j,k) = diss_adve2(k,j,i)
     948                   ENDDO
     949                ENDDO
     950             ENDDO
     951          ENDIF
     952
     953       CASE ( 'diss_diff2' )                                                    !> @todo remove later
     954          IF ( av == 0 )  THEN
     955             DO  i = nxl, nxr
     956                DO  j = nys, nyn
     957                   DO  k = nzb_do, nzt_do
     958                      local_pf(i,j,k) = diss_diff2(k,j,i)
     959                   ENDDO
     960                ENDDO
     961             ENDDO
     962          ENDIF
     963
     964       CASE ( 'diss_prod3' )                                                    !> @todo remove later
     965          IF ( av == 0 )  THEN
     966             DO  i = nxl, nxr
     967                DO  j = nys, nyn
     968                   DO  k = nzb_do, nzt_do
     969                      local_pf(i,j,k) = diss_prod3(k,j,i)
     970                   ENDDO
     971                ENDDO
     972             ENDDO
     973          ENDIF
     974
     975       CASE ( 'diss_adve3' )                                                    !> @todo remove later
     976          IF ( av == 0 )  THEN
     977             DO  i = nxl, nxr
     978                DO  j = nys, nyn
     979                   DO  k = nzb_do, nzt_do
     980                      local_pf(i,j,k) = diss_adve3(k,j,i)
     981                   ENDDO
     982                ENDDO
     983             ENDDO
     984          ENDIF
     985
     986       CASE ( 'diss_diff3' )                                                    !> @todo remove later
     987          IF ( av == 0 )  THEN
     988             DO  i = nxl, nxr
     989                DO  j = nys, nyn
     990                   DO  k = nzb_do, nzt_do
     991                      local_pf(i,j,k) = diss_diff3(k,j,i)
     992                   ENDDO
     993                ENDDO
     994             ENDDO
     995          ENDIF
     996         
    864997       CASE DEFAULT
    865998          found = .FALSE.
     
    8911024    ALLOCATE( km(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    8921025
    893     ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )                           !### remove later
     1026    ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )                           !> @todo remove later
    8941027    ALLOCATE( dummy2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    8951028    ALLOCATE( dummy3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    896 
    897     IF ( rans_mode )  ALLOCATE( l_black(nzb:nzt+1) )
     1029    ALLOCATE( diss_adve1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1030    ALLOCATE( diss_adve2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1031    ALLOCATE( diss_adve3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1032    ALLOCATE( diss_prod1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1033    ALLOCATE( diss_prod2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1034    ALLOCATE( diss_prod3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1035    ALLOCATE( diss_diff1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1036    ALLOCATE( diss_diff2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1037    ALLOCATE( diss_diff3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1038    dummy1 = 0.0_wp
     1039    dummy2 = 0.0_wp
     1040    dummy3 = 0.0_wp
     1041    diss_adve1 = 0.0_wp
     1042    diss_adve2 = 0.0_wp
     1043    diss_adve3 = 0.0_wp
     1044    diss_prod1 = 0.0_wp
     1045    diss_prod2 = 0.0_wp
     1046    diss_prod3 = 0.0_wp
     1047    diss_diff1 = 0.0_wp
     1048    diss_diff2 = 0.0_wp
     1049    diss_diff3 = 0.0_wp
    8981050
    8991051#if defined( __nopointer )
     
    9121064!-- they do not necessarily need to be transferred, which is attributed to
    9131065!-- the design of the model coupler which allocates memory for each variable.
    914     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.  wang_kernel  .OR.      &
     1066    IF ( rans_mode  .OR.  use_sgs_for_particles  .OR.  wang_kernel  .OR.       &
    9151067         collision_turbulence  .OR.  nested_run )  THEN
    9161068#if defined( __nopointer )
     
    9341086    e  => e_1;   e_p  => e_2;   te_m  => e_3
    9351087
    936     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.     &
     1088    IF ( rans_mode  .OR.  use_sgs_for_particles  .OR.     &
    9371089         wang_kernel  .OR.  collision_turbulence  .OR.  nested_run )  THEN
    9381090       diss => diss_1
     
    9671119    INTEGER(iwp) :: j            !< loop index
    9681120    INTEGER(iwp) :: k            !< loop index
    969     INTEGER(iwp) :: nz_s_shift   !<
    970     INTEGER(iwp) :: nz_s_shift_l !<
     1121    INTEGER(iwp) :: nz_s_shift   !< lower shift index for scalars
     1122    INTEGER(iwp) :: nz_s_shift_l !< local lower shift index in case of turbulent inflow
    9711123
    9721124!
    9731125!-- Initialize mixing length
    9741126    CALL tcm_init_mixing_length
     1127    dummy3 = l_wall                 !> @todo remove later
    9751128
    9761129!
     
    9951148
    9961149          IF ( rans_tke_e )  THEN
    997              IF ( dissipation_1d == 'prognostic' )  THEN    !### Why must this be checked?
    998                 DO  i = nxlg, nxrg                          !### Should 'diss' not always
    999                    DO  j = nysg, nyng                       !### be prognostic in case rans_tke_e?
     1150             IF ( dissipation_1d == 'prognostic' )  THEN    !> @query Why must this be checked?
     1151                DO  i = nxlg, nxrg                          ! Should 'diss' not always
     1152                   DO  j = nysg, nyng                       ! be prognostic in case rans_tke_e?
    10001153                      diss(:,j,i) = diss1d
    10011154                   ENDDO
     
    10051158                   DO  j = nysg, nyng
    10061159                      DO  k = nzb+1, nzt
    1007                          diss(k,j,i) = e(k,j,i) * SQRT( e(k,j,i) ) / l1d(k)
     1160                         diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km1d(k)
    10081161                      ENDDO
    10091162                   ENDDO
     
    10161169
    10171170          IF ( constant_diffusion )  THEN
    1018              km   = km_constant
    1019              kh   = km / prandtl_number
    1020              e    = 0.0_wp
     1171             km = km_constant
     1172             kh = km / prandtl_number
     1173             e  = 0.0_wp
    10211174          ELSEIF ( e_init > 0.0_wp )  THEN
    1022              DO  k = nzb+1, nzt
    1023                 km(k,:,:) = c_m * l_grid(k) * SQRT( e_init )
     1175             DO  i = nxlg, nxrg
     1176                DO  j = nysg, nyng
     1177                   DO  k = nzb+1, nzt
     1178                      km(k,j,i) = c_0 * l_wall(k,j,i) * SQRT( e_init )
     1179                   ENDDO
     1180                ENDDO
    10241181             ENDDO
    10251182             km(nzb,:,:)   = km(nzb+1,:,:)
    10261183             km(nzt+1,:,:) = km(nzt,:,:)
    1027              kh   = km / prandtl_number
    1028              e    = e_init
     1184             kh = km / prandtl_number
     1185             e  = e_init
    10291186          ELSE
    10301187             IF ( .NOT. ocean )  THEN
     
    10401197          ENDIF
    10411198
     1199          IF ( rans_tke_e )  THEN
     1200             DO  i = nxlg, nxrg
     1201                DO  j = nysg, nyng
     1202                   DO  k = nzb+1, nzt
     1203                      diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km(k,j,i)
     1204                   ENDDO
     1205                ENDDO
     1206             ENDDO
     1207             diss(nzb,:,:) = diss(nzb+1,:,:)
     1208             diss(nzt+1,:,:) = diss(nzt,:,:)
     1209          ENDIF
     1210
    10421211       ENDIF
    10431212!
     
    10731242             ENDDO
    10741243          ENDDO
     1244          IF ( rans_tke_e )  THEN
     1245             DO  i = nxlg, nxrg
     1246                DO  j = nysg, nyng
     1247                   nz_s_shift = get_topography_top_index_ji( j, i, 's' )
     1248
     1249                   diss(nz_s_shift:nzt+1,j,i) = diss(0:nzt+1-nz_s_shift,j,i)
     1250                ENDDO
     1251             ENDDO
     1252          ENDIF
    10751253       ENDIF
    10761254
     
    10841262!--       boundary and adjust mean inflow profiles
    10851263          IF ( complex_terrain )  THEN
    1086              IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
     1264             IF ( nxlg <= 0 .AND. nxrg >= 0 .AND.  &
     1265                  nysg <= 0 .AND. nyng >= 0        )  THEN
    10871266                nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
    10881267             ELSE
     
    10951274             nz_s_shift = nz_s_shift_l
    10961275#endif
    1097              mean_inflow_profiles(nz_s_shift:nzt+1,5) = hom_sum(0:nzt+1-nz_s_shift,8,0)  ! e
     1276             mean_inflow_profiles(nz_s_shift:nzt+1,5) =  &
     1277                hom_sum(0:nzt+1-nz_s_shift,8,0)  ! e
    10981278          ENDIF
    10991279!
     
    11141294!
    11151295!--       Inside buildings set TKE back to zero.
    1116 !--       Other scalars (km, kh, diss, ...) are ignored at present,
     1296!--       Other scalars (km, kh,...) are ignored at present,
    11171297!--       maybe revise later.
    11181298          DO  i = nxlg, nxrg
     
    11211301                   e(k,j,i)     = MERGE( e(k,j,i), 0.0_wp,                     &
    11221302                                         BTEST( wall_flags_0(k,j,i), 0 ) )
    1123                    te_m(k,j,i)  = MERGE( te_m(k,j,i), 0.0_wp,                  &
    1124                                          BTEST( wall_flags_0(k,j,i), 0 ) )
    11251303                ENDDO
    11261304             ENDDO
    11271305          ENDDO
    11281306
     1307          IF ( rans_tke_e )  THEN
     1308             DO  i = nxlg, nxrg
     1309                DO  j = nysg, nyng
     1310                   DO  k = nzb, nzt
     1311                      diss(k,j,i)    = MERGE( diss(k,j,i), 0.0_wp,             &
     1312                                              BTEST( wall_flags_0(k,j,i), 0 ) )
     1313                   ENDDO
     1314                ENDDO
     1315             ENDDO
     1316          ENDIF
    11291317       ENDIF
    11301318!
     
    11341322!
    11351323!--    Allthough tendency arrays are set in prognostic_equations, they have
    1136 !--    have to be predefined here because they are used (but multiplied with 0)
    1137 !--    there before they are set.
     1324!--    to be predefined here because there they are used (but multiplied with 0)
     1325!--    before they are set.
    11381326       te_m = 0.0_wp
     1327
     1328       IF ( rans_tke_e )  THEN
     1329          diss_p = diss
     1330          tdiss_m = 0.0_wp
     1331       ENDIF
    11391332
    11401333    ENDIF
     
    12151408
    12161409       DO  k = 1, nzt
    1217           IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.           &
     1410          IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.            &
    12181411               l_grid(k) > 1.5_wp * dy * wall_adjustment_factor )  THEN
    1219              WRITE( message_string, * ) 'grid anisotropy exceeds ',            &
    1220                                         'threshold given by only local',       &
    1221                                         ' horizontal reduction of near_wall ', &
    1222                                         'mixing length l_wall',                &
    1223                                         ' starting from height level k = ', k, &
     1412             WRITE( message_string, * ) 'grid anisotropy exceeds ',             &
     1413                                        'threshold given by only local',        &
     1414                                        ' &horizontal reduction of near_wall ', &
     1415                                        'mixing length l_wall',                 &
     1416                                        ' &starting from height level k = ', k, &
    12241417                                        '.'
    12251418             CALL message( 'init_grid', 'PA0202', 0, 1, 0, 6, 0 )
     
    13201513!
    13211514!-- Initialize the mixing length in case of a RANS simulation
     1515       ALLOCATE( l_black(nzb:nzt+1) )
    13221516
    13231517!
    13241518!--    Calculate mixing length according to Blackadar (1962)
    13251519       IF ( f /= 0.0_wp )  THEN
    1326           l_max = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) /               &
    1327                   ABS( f ) + 1E-10_wp
     1520          l_max = 2.7E-4_wp * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) /            &
     1521                  ABS( f ) + 1.0E-10_wp
    13281522       ELSE
    13291523          l_max = 30.0_wp
     
    13381532!
    13391533!--    Gather topography information of whole domain
    1340        !## TODO: reduce amount of data sent by MPI call
    1341        !##  By now, a whole global 3D-array is sent and received with
    1342        !##  MPI_ALLREDUCE although most of the array is 0. This can be
    1343        !##  drastically reduced if only the local subarray is sent and stored
    1344        !##  in a global array. For that, an MPI data type or subarray must be
    1345        !##  defined.
    1346        !##  2018-03-19, gronemeier
     1534       !> @todo reduce amount of data sent by MPI call
     1535       !>   By now, a whole global 3D-array is sent and received with
     1536       !>   MPI_ALLREDUCE although most of the array is 0. This can be
     1537       !>   drastically reduced if only the local subarray is sent and stored
     1538       !>   in a global array. For that, an MPI data type or subarray must be
     1539       !>   defined.
     1540       !>   2018-03-19, gronemeier
    13471541       ALLOCATE( wall_flags_0_global(nzb:nzt+1,0:ny,0:nx) )
    13481542
     
    13731567          ENDDO
    13741568       ENDDO
     1569
     1570       l_wall(nzb,:,:) = l_black(nzb)
     1571       l_wall(nzt+1,:,:) = l_black(nzt+1)
    13751572!
    13761573!--    Limit mixing length to either nearest wall or Blackadar mixing length.
     
    14151612          IF ( rad_k_b /= 0 .OR. rad_k_t /= 0 )  THEN
    14161613
    1417              !## NOTE: shape of vicinity is larger in z direction
    1418              !##  Shape of vicinity is two grid points larger than actual search
    1419              !##  radius in vertical direction. The first and last grid point is
    1420              !##   always set to 1 to asure correct detection of topography. See
    1421              !##  function "shortest_distance" for details.
    1422              !##  2018-03-16, gronemeier
     1614             !> @note shape of vicinity is larger in z direction
     1615             !>   Shape of vicinity is two grid points larger than actual search
     1616             !>   radius in vertical direction. The first and last grid point is
     1617             !>   always set to 1 to asure correct detection of topography. See
     1618             !>   function "shortest_distance" for details.
     1619             !>   2018-03-16, gronemeier
    14231620             ALLOCATE( vicinity(-rad_k-1:rad_k+1,-rad_j:rad_j,-rad_i:rad_i) )
    14241621             ALLOCATE( vic_yz(0:rad_k+1,0:rad_j) )
     
    16011798                   ELSE  !Check if (i,j,k) belongs to atmosphere
    16021799
    1603                       l_wall(k,j,i) = -999.0
     1800                      l_wall(k,j,i) = l_black(k)
    16041801
    16051802                   ENDIF
     
    16301827!> (pos_i/jj/kk), where (jj/kk) is the position of the maximum of 'array'
    16311828!> closest to the origin (0/0) of 'array'.
     1829!> @todo this part of PALM does not reproduce the same results for optimized
     1830!>   and debug options for the compiler. This should be fixed
    16321831!------------------------------------------------------------------------------!
    16331832    REAL FUNCTION shortest_distance( array, orientation, pos_i )
     
    16591858       IF ( orientation ) THEN  !if array is oriented upwards
    16601859          DO  jj = 0, rad_j
    1661              shortest_distance = MIN( shortest_distance,                        &
    1662                                       SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2         &
    1663                                           + MAX(jj*dy-0.5*dy,0.0)**2            &
    1664                                           + MAX(zw(loc_k(jj)+k-1)-zu(k),0.0)**2 &
    1665                                           )                                     &
    1666                                     )
     1860             shortest_distance =                                               &
     1861                MIN( shortest_distance,                                        &
     1862                     SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2   &
     1863                         + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2      &
     1864                         + MAX(zw(loc_k(jj)+k-1)-zu(k), 0.0_wp)**2             &
     1865                         )                                                     &
     1866                   )
    16671867          ENDDO
    16681868       ELSE  !if array is oriented downwards
    1669           !## NOTE: MAX within zw required to circumvent error at domain border
    1670           !##  At the domain border, if non-cyclic boundary is present, the
    1671           !##  index for zw could be -1, which will be errorneous (zw(-1) does
    1672           !##  not exist). The MAX function limits the index to be at least 0.
     1869          !> @note MAX within zw required to circumvent error at domain border
     1870          !>   At the domain border, if non-cyclic boundary is present, the
     1871          !>   index for zw could be -1, which will be errorneous (zw(-1) does
     1872          !>   not exist). The MAX function limits the index to be at least 0.
    16731873          DO  jj = 0, rad_j
    1674              shortest_distance = MIN( shortest_distance,                       &
    1675                                       SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2        &
    1676                                           + MAX(jj*dy-0.5*dy,0.0)**2           &
    1677                                           + MAX(zu(k)-zw(MAX(k-loc_k(jj),      &
    1678                                                              0_iwp)),          &
    1679                                                 0.0)**2                        &
    1680                                           )                                    &
    1681                                     )
     1874             shortest_distance =                                               &
     1875                MIN( shortest_distance,                                        &
     1876                     SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2   &
     1877                         + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2      &
     1878                         + MAX(zu(k)-zw(MAX(k-loc_k(jj),0_iwp)), 0.0_wp)**2    &
     1879                         )                                                     &
     1880                   )
    16821881          ENDDO
    16831882       ENDIF
    1684 
     1883       
    16851884    END FUNCTION
    16861885
     
    20042203!
    20052204!--    Use special boundary condition in case of TKE-e closure
     2205       !> @todo do the same for usm and lsm surfaces
     2206       !>   2018-06-05, gronemeier
    20062207       IF ( rans_tke_e )  THEN
    20072208          DO  i = nxl, nxr
     
    20112212                DO  m = surf_s, surf_e
    20122213                   k = surf_def_h(0)%k(m)
    2013                    e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_m**2
     2214                   e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_0**2
    20142215                ENDDO
    20152216             ENDDO
     
    20922293             DO  k = nzb+1, nzt
    20932294!                tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k)
    2094                 tend(k,j,i) = tend(k,j,i) + c_1 * c_mu * f / c_h               &  !### needs revision
     2295                tend(k,j,i) = tend(k,j,i) + c_1 * c_0**4 * f / c_4               &  !> @todo needs revision
    20952296                      / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))       &
    20962297                      * SQRT(e(k,j,i)) * produc(k,j,i)
     
    21032304!
    21042305!--    Additional sink term for flows through plant canopies
    2105 !        IF ( plant_canopy )  CALL pcm_tendency( ? )                            !### what to do with this?
    2106 
    2107 !        CALL user_actions( 'diss-tendency' )                                   !### not yet implemented
     2306!        IF ( plant_canopy )  CALL pcm_tendency( ? )                            !> @query what to do with this?
     2307
     2308!        CALL user_actions( 'diss-tendency' )                                   !> @todo not yet implemented
    21082309
    21092310!
     
    21812382    USE arrays_3d,                                                             &
    21822383        ONLY:  ddzu, diss_l_diss, diss_l_e, diss_s_diss, diss_s_e,             &
    2183                flux_l_diss, flux_l_e, flux_s_diss, flux_s_e
     2384               flux_l_diss, flux_l_e, flux_s_diss, flux_s_e,&
     2385               u_p,v_p,w_p
    21842386
    21852387    USE control_parameters,                                                    &
    21862388        ONLY:  f, tsc
     2389
     2390    USE grid_variables,                                                        &
     2391        ONLY:  dx, dy
    21872392
    21882393    USE surface_mod,                                                           &
     
    21902395                surf_usm_v
    21912396
     2397    use indices, only: nx, ny
     2398
    21922399    IMPLICIT NONE
    21932400
    21942401    INTEGER(iwp) ::  i       !< loop index x direction
    2195     INTEGER(iwp) ::  i_omp   !<
     2402    INTEGER(iwp) ::  i_omp   !< first loop index of i-loop in prognostic_equations
    21962403    INTEGER(iwp) ::  j       !< loop index y direction
    21972404    INTEGER(iwp) ::  k       !< loop index z direction
     2405    INTEGER(iwp) ::  l       !< loop index
    21982406    INTEGER(iwp) ::  m       !< loop index
    21992407    INTEGER(iwp) ::  surf_e  !< end index of surface elements at given i-j position
    22002408    INTEGER(iwp) ::  surf_s  !< start index of surface elements at given i-j position
    2201     INTEGER(iwp) ::  tn      !<
    2202 
    2203     REAL(wp), DIMENSION(nzb:nzt+1) :: advec  !< advection term of TKE tendency
    2204     REAL(wp), DIMENSION(nzb:nzt+1) :: produc !< production term of TKE tendency
     2409    INTEGER(iwp) ::  tn      !< task number of openmp task
     2410
     2411    INTEGER(iwp) :: pis = 32 !< debug variable, print from i=pis                !> @todo remove later
     2412    INTEGER(iwp) :: pie = 32 !< debug variable, print until i=pie               !> @todo remove later
     2413    INTEGER(iwp) :: pjs = 26 !< debug variable, print from j=pjs                !> @todo remove later
     2414    INTEGER(iwp) :: pje = 26 !< debug variable, print until j=pje               !> @todo remove later
     2415    INTEGER(iwp) :: pkb = 1  !< debug variable, print from k=pkb                !> @todo remove later
     2416    INTEGER(iwp) :: pkt = 7  !< debug variable, print until k=pkt               !> @todo remove later
     2417
     2418    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_adv   !< debug variable               !> @todo remove later
     2419    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_pro   !< debug variable               !> @todo remove later
     2420    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_dif   !< debug variable               !> @todo remove later
     2421
     24225555 FORMAT(A,7(1X,E12.5))   !> @todo remove later
    22052423
    22062424!
     
    22242442       ENDIF
    22252443
    2226        advec(:) = tend(:,j,i)
    2227 
    2228        CALL production_e( i, j )
    2229 
    2230        produc(:) = tend(:,j,i) - advec(:)
     2444       dum_adv = tend(:,j,i)                                                    !> @todo remove later
     2445
     2446       CALL production_e( i, j, .FALSE. )
     2447
     2448       dum_pro = tend(:,j,i) - dum_adv                                          !> @todo remove later
    22312449
    22322450       IF ( .NOT. humidity )  THEN
     
    22392457          CALL diffusion_e( i, j, vpt, pt_reference )
    22402458       ENDIF
     2459
     2460       dum_dif = tend(:,j,i) - dum_adv - dum_pro                                !> @todo remove later
    22412461
    22422462!
     
    22682488          DO  m = surf_s, surf_e
    22692489             k = surf_def_h(0)%k(m)
    2270              e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_m**2
     2490             e_p(k,j,i) = ( surf_def_h(0)%us(m) / c_0 )**2
     2491          ENDDO
     2492
     2493          DO  l = 0, 3
     2494             surf_s = surf_def_v(l)%start_index(j,i)
     2495             surf_e = surf_def_v(l)%end_index(j,i)
     2496             DO  m = surf_s, surf_e
     2497                k = surf_def_v(l)%k(m)
     2498                e_p(k,j,i) = ( surf_def_v(l)%us(m) / c_0 )**2
     2499             ENDDO
    22712500          ENDDO
    22722501       ENDIF
     
    22872516          ENDIF
    22882517       ENDIF
     2518
     2519!        if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then        !> @todo remove later
     2520!           WRITE(9, *) '------'
     2521!           WRITE(9, '(A,F8.3,1X,F8.3,1X,I2)') 't, dt, int_ts:', simulated_time, dt_3d, intermediate_timestep_count
     2522!           WRITE(9, *) 'i:', i
     2523!           WRITE(9, *) 'j:', j
     2524!           WRITE(9, *) 'k:', pkb, ' - ', pkt
     2525!           WRITE(9, *) '---'
     2526!           WRITE(9, *) 'e:'
     2527!           WRITE(9, 5555) 'adv :', dum_adv(pkb:pkt)
     2528!           WRITE(9, 5555) 'pro :', dum_pro(pkb:pkt)
     2529!           WRITE(9, 5555) 'dif :', dum_dif(pkb:pkt)
     2530!           WRITE(9, 5555) 'tend:', tend(pkb:pkt,j,i)
     2531!           WRITE(9, 5555) 'e_p :', e_p(pkb:pkt,j,i)
     2532!           WRITE(9, 5555) 'e   :', e(pkb:pkt,j,i)
     2533!           FLUSH(9)
     2534!        endif
    22892535
    22902536    ENDIF   ! TKE equation
     
    23092555       ENDIF
    23102556
     2557       IF ( intermediate_timestep_count == 1 )  diss_adve1(:,j,i) = tend(:,j,i) !> @todo remove later
     2558       IF ( intermediate_timestep_count == 2 )  diss_adve2(:,j,i) = tend(:,j,i)
     2559       IF ( intermediate_timestep_count == 3 )  diss_adve3(:,j,i) = tend(:,j,i)
     2560
    23112561!
    23122562!--    Production of TKE dissipation rate
    2313        DO  k = nzb+1, nzt
    2314 !              tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k)
    2315           tend(k,j,i) = tend(k,j,i) + c_1 * c_mu * f / c_h                     &  !### needs revision
    2316                       / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))       &
    2317                       * SQRT(e(k,j,i)) * produc(k)
    2318        ENDDO
    2319 
     2563       CALL production_e( i, j, .TRUE. )
     2564
     2565       IF ( intermediate_timestep_count == 1 )  diss_prod1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) !> @todo remove later
     2566       IF ( intermediate_timestep_count == 2 )  diss_prod2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i)
     2567       IF ( intermediate_timestep_count == 3 )  diss_prod3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i)
     2568
     2569       dum_pro = tend(:,j,i) - dum_adv                                          !> @todo remove later
     2570
     2571!
     2572!--    Diffusion term of TKE dissipation rate
    23202573       CALL diffusion_diss( i, j )
    23212574
     2575       IF ( intermediate_timestep_count == 1 )  diss_diff1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) - diss_prod1(:,j,i) !> @todo remove later
     2576       IF ( intermediate_timestep_count == 2 )  diss_diff2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i) - diss_prod2(:,j,i)
     2577       IF ( intermediate_timestep_count == 3 )  diss_diff3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i) - diss_prod3(:,j,i)
     2578       IF ( intermediate_timestep_count == 3 )  dummy3(:,j,i) = km(:,j,i)
     2579
     2580       dum_dif = tend(:,j,i) - dum_adv - dum_pro                                !> @todo remove later
     2581
    23222582!
    23232583!--    Additional sink term for flows through plant canopies
    2324 !        IF ( plant_canopy )  CALL pcm_tendency( i, j, ? )                      !### not yet implemented
    2325 
    2326 !        CALL user_actions( i, j, 'diss-tendency' )                             !### not yet implemented
     2584!        IF ( plant_canopy )  CALL pcm_tendency( i, j, ? )                      !> @todo not yet implemented
     2585
     2586!        CALL user_actions( i, j, 'diss-tendency' )                             !> @todo not yet implemented
    23272587
    23282588!
     
    23382598                                                BTEST( wall_flags_0(k,j,i), 0 )&
    23392599                                               )
    2340           IF ( diss_p(k,j,i) <= 0.0_wp )  diss_p(k,j,i) = 0.1_wp * diss(k,j,i)
    23412600       ENDDO
    23422601
     
    23502609       ENDDO
    23512610
     2611       DO  l = 0, 1
     2612          surf_s = surf_def_v(l)%start_index(j,i)
     2613          surf_e = surf_def_v(l)%end_index(j,i)
     2614          DO  m = surf_s, surf_e
     2615             k = surf_def_v(l)%k(m)
     2616             diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dy )
     2617          ENDDO
     2618       ENDDO
     2619
     2620       DO  l = 2, 3
     2621          surf_s = surf_def_v(l)%start_index(j,i)
     2622          surf_e = surf_def_v(l)%end_index(j,i)
     2623          DO  m = surf_s, surf_e
     2624             k = surf_def_v(l)%k(m)
     2625             diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dx )
     2626          ENDDO
     2627       ENDDO
    23522628!
    23532629!--    Calculate tendencies for the next Runge-Kutta step
     
    23652641          ENDIF
    23662642       ENDIF
    2367 
    2368 !        IF ( intermediate_timestep_count == 1 )  dummy1(:,j,i) = e_p(:,j,i)
    2369 !        IF ( intermediate_timestep_count == 2 )  dummy2(:,j,i) = e_p(:,j,i)
    2370 !        IF ( intermediate_timestep_count == 3 )  dummy3(:,j,i) = e_p(:,j,i)
     2643!
     2644!--    Limit change of diss to be between -90% and +100%. Also, set an absolute
     2645!--    minimum value
     2646       DO  k = nzb, nzt+1
     2647          diss_p(k,j,i) = MIN( MAX( diss_p(k,j,i),         &
     2648                                    0.1_wp * diss(k,j,i),  &
     2649                                    0.0001_wp ),           &
     2650                               2.0_wp * diss(k,j,i) )
     2651       ENDDO
     2652
     2653       IF ( intermediate_timestep_count == 1 )  dummy1(:,j,i) = diss_p(:,j,i)   !> @todo remove later
     2654       IF ( intermediate_timestep_count == 2 )  dummy2(:,j,i) = diss_p(:,j,i)
     2655
     2656!        if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then        !> @todo remove later
     2657!           WRITE(9, *) '---'
     2658!           WRITE(9, *) 'diss:'
     2659!           WRITE(9, 5555) 'adv   :', dum_adv(pkb:pkt)
     2660!           WRITE(9, 5555) 'pro   :', dum_pro(pkb:pkt)
     2661!           WRITE(9, 5555) 'dif   :', dum_dif(pkb:pkt)
     2662!           WRITE(9, 5555) 'tend  :', tend(pkb:pkt,j,i)
     2663!           WRITE(9, 5555) 'diss_p:', diss_p(pkb:pkt,j,i)
     2664!           WRITE(9, 5555) 'diss  :', diss(pkb:pkt,j,i)
     2665!           WRITE(9, *) '---'
     2666!           WRITE(9, 5555) 'km    :', km(pkb:pkt,j,i)
     2667!           flush(9)
     2668!        endif
    23712669
    23722670    ENDIF   ! dissipation equation
     
    23822680!> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is
    23832681!>          not considered well!
     2682!> @todo Adjust production term in case of rans_tke_e simulation
    23842683!------------------------------------------------------------------------------!
    23852684 SUBROUTINE production_e
     
    30993398!> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is
    31003399!>          not considered well!
     3400!> @todo non-neutral case is not yet considered for RANS mode
    31013401!------------------------------------------------------------------------------!
    3102  SUBROUTINE production_e_ij( i, j )
     3402 SUBROUTINE production_e_ij( i, j, diss_production )
    31033403
    31043404    USE arrays_3d,                                                             &
     
    31213421
    31223422    IMPLICIT NONE
     3423
     3424    LOGICAL :: diss_production
    31233425
    31243426    INTEGER(iwp) ::  i       !< running index x-direction
     
    31533455    REAL(wp), DIMENSION(nzb+1:nzt)  ::  dwdy        !< Gradient of w-component in y-direction
    31543456    REAL(wp), DIMENSION(nzb+1:nzt)  ::  dwdz        !< Gradient of w-component in z-direction
    3155 
     3457    REAL(wp), DIMENSION(nzb+1:nzt)  ::  tend_temp   !< temporal tendency
    31563458
    31573459    IF ( constant_flux_layer )  THEN
     
    33483650       ENDDO
    33493651
     3652!        IF ( .NOT. rans_tke_e )  THEN
     3653
    33503654       DO  k = nzb+1, nzt
    33513655
    3352           def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +            &
    3353                            dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +            &
    3354                            dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +            &
    3355                 2.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) )
    3356 
    3357 !
    3358 !--       Production term according to Kato and Launder (1993)
    3359 !           def = SQRT( ( dudx(k) + dudy(k) + dudz(k) +                       &
    3360 !                         dvdx(k) + dvdy(k) + dvdz(k) +                       &
    3361 !                         dwdx(k) + dwdy(k) + dwdz(k)                         &
    3362 !                       )**4 -                                                &
    3363 !                       ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 +              &
    3364 !                         2.0_wp * ( dudy(k) * dvdx(k) +                      &
    3365 !                                    dudz(k) * dwdx(k) +                      &
    3366 !                                    dvdz(k) * dwdy(k) )                      &
    3367 !                       )**2                                                  &
    3368 !                     )
     3656          def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +         &
     3657                           dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +         &
     3658                           dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +         &
     3659                2.0_wp * ( dvdx(k)*dudy(k) +                                &
     3660                           dwdx(k)*dudz(k) +                                &
     3661                           dwdy(k)*dvdz(k) )
    33693662
    33703663          IF ( def < 0.0_wp )  def = 0.0_wp
     
    33723665          flag  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    33733666
    3374           tend(k,j,i) = tend(k,j,i) + km(k,j,i) * def * flag
     3667          tend_temp(k) = km(k,j,i) * def * flag
    33753668
    33763669       ENDDO
    33773670
     3671!        ELSE
     3672!
     3673!           DO  k = nzb+1, nzt
     3674! !
     3675! !--          Production term according to Kato and Launder (1993)
     3676!              def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
     3677!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 +              &
     3678!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
     3679!                                       dvdz(k) * dwdy(k) +                      &
     3680!                                       dwdx(k) * dudz(k) )       )              &
     3681!                        * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
     3682!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 -              &
     3683!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
     3684!                                       dvdz(k) * dwdy(k) +                      &
     3685!                                       dwdx(k) * dudz(k) )       )              &
     3686!                        )
     3687!
     3688!              IF ( def < 0.0_wp )  def = 0.0_wp
     3689!
     3690!              flag  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
     3691!
     3692!              tend_temp(k) = km(k,j,i) * def * flag
     3693!
     3694!           ENDDO
     3695!
     3696!        ENDIF
     3697
     3698    ELSE  ! not constant_flux_layer
     3699
     3700!        IF ( .NOT. rans_tke_e )  THEN
     3701!
     3702!--       Calculate TKE production by shear. Here, no additional
     3703!--       wall-bounded code is considered.
     3704!--       Why?
     3705          DO  k = nzb+1, nzt
     3706
     3707             dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
     3708             dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                &
     3709                                    u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
     3710             dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                &
     3711                                    u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
     3712
     3713             dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                &
     3714                                    v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
     3715             dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
     3716             dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                &
     3717                                    v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
     3718
     3719             dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                &
     3720                                    w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
     3721             dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                &
     3722                                    w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
     3723             dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
     3724
     3725             def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +         &
     3726                              dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +         &
     3727                              dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +         &
     3728                   2.0_wp * ( dvdx(k)*dudy(k) +                                &
     3729                              dwdx(k)*dudz(k) +                                &
     3730                              dwdy(k)*dvdz(k) )
     3731
     3732             IF ( def < 0.0_wp )  def = 0.0_wp
     3733
     3734             flag        = MERGE( 1.0_wp, 0.0_wp,                              &
     3735                                  BTEST( wall_flags_0(k,j,i), 29 ) )
     3736             tend_temp(k) = km(k,j,i) * def * flag
     3737
     3738          ENDDO
     3739
     3740!        ELSE
     3741!
     3742!           DO  k = nzb+1, nzt
     3743!
     3744!              dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
     3745!              dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                &
     3746!                                     u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
     3747!              dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                &
     3748!                                     u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
     3749!
     3750!              dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                &
     3751!                                     v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
     3752!              dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
     3753!              dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                &
     3754!                                     v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
     3755!
     3756!              dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                &
     3757!                                     w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
     3758!              dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                &
     3759!                                     w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
     3760!              dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
     3761! !
     3762! !--          Production term according to Kato and Launder (1993)
     3763!              def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
     3764!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 +              &
     3765!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
     3766!                                       dvdz(k) * dwdy(k) +                      &
     3767!                                       dwdx(k) * dudz(k) )       )              &
     3768!                        * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
     3769!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 -              &
     3770!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
     3771!                                       dvdz(k) * dwdy(k) +                      &
     3772!                                       dwdx(k) * dudz(k) )       )              &
     3773!                        )
     3774!
     3775!              IF ( def < 0.0_wp )  def = 0.0_wp
     3776!
     3777!              flag        = MERGE( 1.0_wp, 0.0_wp,                              &
     3778!                                   BTEST( wall_flags_0(k,j,i), 29 ) )
     3779!              tend_temp(k) = km(k,j,i) * def * flag
     3780!
     3781!           ENDDO
     3782!
     3783!        ENDIF
     3784
     3785    ENDIF
     3786
     3787    IF ( .NOT. diss_production )  THEN
     3788!
     3789!--    Production term in case of TKE production
     3790       DO  k = nzb+1, nzt
     3791          tend(k,j,i) = tend(k,j,i) + tend_temp(k)
     3792       ENDDO
    33783793    ELSE
    33793794!
    3380 !--    Calculate TKE production by shear. Here, no additional
    3381 !--    wall-bounded code is considered.
    3382 !--    Why?
     3795!--    Production term in case of dissipation-rate production (rans_tke_e)
    33833796       DO  k = nzb+1, nzt
    33843797
    3385           dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
    3386           dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                   &
    3387                                  u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
    3388           dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                   &
    3389                                  u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
    3390 
    3391           dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                   &
    3392                                  v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
    3393           dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
    3394           dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                   &
    3395                                  v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
    3396 
    3397           dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                   &
    3398                                  w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
    3399           dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                   &
    3400                                  w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
    3401           dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
    3402 
    3403           def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +            &
    3404                            dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +            &
    3405                            dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +            &
    3406                 2.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) )
    3407 
    3408 !
    3409 !--       Production term according to Kato and Launder (1993)
    3410 !           def = SQRT( ( dudx(k) + dudy(k) + dudz(k) +                        &
    3411 !                         dvdx(k) + dvdy(k) + dvdz(k) +                        &
    3412 !                         dwdx(k) + dwdy(k) + dwdz(k)                          &
    3413 !                       )**4 -                                                 &
    3414 !                       ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 +               &
    3415 !                         2.0_wp * ( dudy(k) * dvdx(k) +                       &
    3416 !                                    dudz(k) * dwdx(k) +                       &
    3417 !                                    dvdz(k) * dwdy(k) )                       &
    3418 !                       )**2                                                   &
    3419 !                     )
    3420 
    3421           IF ( def < 0.0_wp )  def = 0.0_wp
    3422 
    3423           flag        = MERGE( 1.0_wp, 0.0_wp,                                 &
    3424                                BTEST( wall_flags_0(k,j,i), 29 ) )
    3425           tend(k,j,i) = tend(k,j,i) + km(k,j,i) * def * flag
    3426 
     3798          ! Standard TKE-e closure
     3799          tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i)               &
     3800                                    /( e(k,j,i) + 1.0E-20_wp )                 &
     3801                                    * c_1
     3802!           ! Production according to Koblitz (2013)
     3803!           tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i)               &
     3804!                                     /( e(k,j,i) + 1.0E-20_wp )                 &
     3805!                                     * ( c_1 + ( c_2 - c_1 )                    &
     3806!                                             * l_wall(k,j,i) / l_max )
     3807!           ! Production according to Detering and Etling (1985)
     3808!           !> @todo us is not correct if there are vertical walls
     3809!           tend(k,j,i) = tend(k,j,i) + tend_temp(k) * SQRT(e(k,j,i))            &
     3810!                                     * c_1 * c_0**3 / c_4 * f                   &
     3811!              / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))
    34273812       ENDDO
    3428 
    34293813    ENDIF
    34303814
     
    38224206    REAL(wp)     ::  l              !< mixing length
    38234207    REAL(wp)     ::  ll             !< adjusted l
    3824     REAL(wp)     ::  var_reference  !<
     4208    REAL(wp)     ::  var_reference  !< reference temperature
    38254209
    38264210#if defined( __nopointer )
    3827     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
     4211    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
    38284212#else
    3829     REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !<
     4213    REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !< temperature
    38304214#endif
    38314215    REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  dissipation  !< TKE dissipation
     
    38374221       DO  j = nys, nyn
    38384222          DO  k = nzb+1, nzt
     4223!
     4224!--          Predetermine flag to mask topography
     4225             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    38394226
    38404227!
     
    38514238                CALL mixing_length_rans( i, j, k, l, ll, var, var_reference )
    38524239
    3853                 dissipation(k,j) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
     4240                dissipation(k,j) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
     4241
     4242                diss(k,j,i) = dissipation(k,j) * flag
    38544243
    38554244             ELSEIF ( rans_tke_e )  THEN
     
    38594248             ENDIF
    38604249
    3861 !
    3862 !--          Predetermine flag to mask topography
    3863              flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    3864 
    3865              tend(k,j,i) = tend(k,j,i)                                         &
    3866                                      + (                                       &
     4250             tend(k,j,i) = tend(k,j,i) + (                                     &
     4251                                           (                                   &
    38674252                       ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) )     &
    38684253                     - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) )     &
    3869                                        ) * ddx2  * flag                        &
    3870                                      + (                                       &
     4254                                           ) * ddx2  * flag                    &
     4255                                         + (                                   &
    38714256                       ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) )     &
    38724257                     - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) )     &
    3873                                        ) * ddy2  * flag                        &
    3874                                      + (                                       &
     4258                                           ) * ddy2  * flag                    &
     4259                                         + (                                   &
    38754260            ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1)    &
    38764261                                                          * rho_air_zw(k)      &
    38774262          - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)      &
    38784263                                                          * rho_air_zw(k-1)    &
    3879                                        ) * ddzw(k) * drho_air(k) * flag        &
     4264                                           ) * ddzw(k) * drho_air(k)           &
     4265                                         ) * flag * dsig_e                     &
    38804266                          - dissipation(k,j) * flag
    38814267
     
    39594345    REAL(wp)     ::  l              !< mixing length
    39604346    REAL(wp)     ::  ll             !< adjusted l
    3961     REAL(wp)     ::  var_reference  !<
     4347    REAL(wp)     ::  var_reference  !< reference temperature
    39624348
    39634349#if defined( __nopointer )
    3964     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
     4350    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
    39654351#else
    3966     REAL(wp), DIMENSION(:,:,:), POINTER ::  var     !<
     4352    REAL(wp), DIMENSION(:,:,:), POINTER ::  var     !< temperature
    39674353#endif
    39684354    REAL(wp), DIMENSION(nzb+1:nzt) ::  dissipation  !< dissipation of TKE
     
    39914377          CALL mixing_length_rans( i, j, k, l, ll, var, var_reference  )
    39924378
    3993           dissipation(k) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
     4379          dissipation(k) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
     4380
     4381          diss(k,j,i) = dissipation(k) * flag
    39944382
    39954383       ELSEIF ( rans_tke_e )  THEN
     
    40014389!
    40024390!--    Calculate the tendency term
    4003        tend(k,j,i) = tend(k,j,i)                                               &
    4004                                     + (                                        &
     4391       tend(k,j,i) = tend(k,j,i) + (                                           &
     4392                                      (                                        &
    40054393                      ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) )      &
    40064394                    - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) )      &
    4007                                       ) * ddx2  * flag / sig_e                 &
     4395                                      ) * ddx2                                 &
    40084396                                    + (                                        &
    40094397                      ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) )      &
    40104398                    - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) )      &
    4011                                       ) * ddy2  * flag / sig_e                 &
     4399                                      ) * ddy2                                 &
    40124400                                    + (                                        &
    40134401           ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1)     &
     
    40154403         - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)       &
    40164404                                                         * rho_air_zw(k-1)     &
    4017                                       ) * ddzw(k) * drho_air(k) * flag / sig_e &
    4018                                     - dissipation(k) * flag
     4405                                      ) * ddzw(k) * drho_air(k)                &
     4406                                   ) * flag * dsig_e                           &
     4407                                 - dissipation(k) * flag
    40194408
    40204409    ENDDO
     
    40824471             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    40834472
    4084              tend(k,j,i) = tend(k,j,i)                                         &
    4085                                + (                                             &
     4473             tend(k,j,i) = tend(k,j,i) +                                       &
     4474                         (      (                                             &
    40864475                 ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) )     &
    40874476               - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) )     &
    4088                                  ) * ddx2  * flag                              &
     4477                                 ) * ddx2                                      &
    40894478                               + (                                             &
    40904479                 ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) )     &
    40914480               - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) )     &
    4092                                  ) * ddy2  * flag                              &
     4481                                 ) * ddy2                                      &
    40934482                               + (                                             &
    40944483      ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1)    &
     
    40964485    - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k)      &
    40974486                                                    * rho_air_zw(k-1)          &
    4098                                  ) * ddzw(k) * drho_air(k) * flag              &
    4099                                - c_2 * diss(k,j,i)**2                          &
    4100                                      / ( e(k,j,i) + 1.0E-20_wp ) * flag
     4487                                 ) * ddzw(k) * drho_air(k)                     &
     4488                         ) * flag * dsig_diss                                  &
     4489                         - c_2 * diss(k,j,i)**2                                &
     4490                               / ( e(k,j,i) + 1.0E-20_wp ) * flag
    41014491
    41024492          ENDDO
     
    41234513    IMPLICIT NONE
    41244514
    4125     INTEGER(iwp) ::  i              !< running index x direction
    4126     INTEGER(iwp) ::  j              !< running index y direction
    4127     INTEGER(iwp) ::  k              !< running index z direction
    4128 
    4129     REAL(wp)     ::  flag           !< flag to mask topography
    4130 
    4131     REAL(wp), DIMENSION(nzb+1:nzt) :: tend_temp
     4515    INTEGER(iwp) ::  i         !< running index x direction
     4516    INTEGER(iwp) ::  j         !< running index y direction
     4517    INTEGER(iwp) ::  k         !< running index z direction
     4518
     4519    REAL(wp)     ::  flag      !< flag to mask topography
    41324520
    41334521!
     
    41414529!
    41424530!--    Calculate the tendency term
    4143        tend_temp(k) =           (                                              &
     4531       tend(k,j,i) =  tend(k,j,i) +                                            &
     4532                   (            (                                              &
    41444533                ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) )      &
    41454534              - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) )      &
    4146                                 ) * ddx2  * flag / sig_diss                    &
     4535                                ) * ddx2                                       &
    41474536                              + (                                              &
    41484537                ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) )      &
    41494538              - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) )      &
    4150                                 ) * ddy2  * flag / sig_diss                    &
     4539                                ) * ddy2                                       &
    41514540                              + (                                              &
    41524541     ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1)     &
     
    41544543   - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k)       &
    41554544                                                   * rho_air_zw(k-1)           &
    4156                                 ) * ddzw(k) * drho_air(k) * flag / sig_diss    &
    4157                               - c_2 * diss(k,j,i)**2                           &
    4158                                     / ( e(k,j,i) + 1.0E-20_wp ) * flag
    4159 
    4160        tend(k,j,i) = tend(k,j,i) + tend_temp(k)
     4545                                ) * ddzw(k) * drho_air(k)                      &
     4546                   ) * flag * dsig_diss                                        &
     4547                   - c_2 * diss(k,j,i)**2 / ( e(k,j,i) + 1.0E-20_wp ) * flag
    41614548
    41624549    ENDDO
     
    42904677!> Computation of the turbulent diffusion coefficients for momentum and heat
    42914678!> according to Prandtl-Kolmogorov.
     4679!> @todo consider non-default surfaces
    42924680!------------------------------------------------------------------------------!
    42934681 SUBROUTINE tcm_diffusivities( var, var_reference )
     
    42974685        ONLY:  e_min, outflow_l, outflow_n, outflow_r, outflow_s
    42984686
     4687    USE grid_variables,                                                        &
     4688        ONLY:  dx, dy
     4689
    42994690    USE statistics,                                                            &
    43004691        ONLY :  rmask, sums_l_l
    43014692
    43024693    USE surface_mod,                                                           &
    4303         ONLY :  bc_h, surf_def_h
     4694        ONLY :  bc_h, surf_def_h, surf_def_v
    43044695
    43054696    IMPLICIT NONE
    43064697
    4307     INTEGER(iwp) ::  i                   !<
    4308     INTEGER(iwp) ::  j                   !<
    4309     INTEGER(iwp) ::  k                   !<
    4310     INTEGER(iwp) ::  m                   !<
    4311     INTEGER(iwp) ::  n                   !<
    4312     INTEGER(iwp) ::  omp_get_thread_num  !<
    4313     INTEGER(iwp) ::  sr                  !<
    4314     INTEGER(iwp) ::  tn                  !<
    4315 
    4316     REAL(wp)     ::  flag                !<
    4317     REAL(wp)     ::  l                   !<
    4318     REAL(wp)     ::  ll                  !<
    4319     REAL(wp)     ::  var_reference       !<
     4698    INTEGER(iwp) ::  i                   !< loop index
     4699    INTEGER(iwp) ::  j                   !< loop index
     4700    INTEGER(iwp) ::  k                   !< loop index
     4701    INTEGER(iwp) ::  m                   !< loop index
     4702    INTEGER(iwp) ::  n                   !< loop index
     4703    INTEGER(iwp) ::  omp_get_thread_num  !< opemmp function to get thread number
     4704    INTEGER(iwp) ::  sr                  !< statistic region
     4705    INTEGER(iwp) ::  tn                  !< thread number
     4706
     4707    REAL(wp)     ::  flag                !< topography flag
     4708    REAL(wp)     ::  l                   !< mixing length
     4709    REAL(wp)     ::  ll                  !< adjusted mixing length
     4710    REAL(wp)     ::  var_reference       !< reference temperature
    43204711
    43214712#if defined( __nopointer )
    4322     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
     4713    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
    43234714#else
    4324     REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !<
     4715    REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !< temperature
    43254716#endif
    43264717
     
    43654756!
    43664757!--             Compute diffusion coefficients for momentum and heat
    4367                 km(k,j,i) = c_m * l * SQRT( e(k,j,i) ) * flag
     4758                km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag
    43684759                kh(k,j,i) = ( 1.0_wp + 2.0_wp * l / ll ) * km(k,j,i) * flag
    43694760!
     
    43914782!
    43924783!--             Compute diffusion coefficients for momentum and heat
    4393                 km(k,j,i) = c_m * l * SQRT( e(k,j,i) ) * flag
     4784                km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag
    43944785                kh(k,j,i) = km(k,j,i) / prandtl_number * flag
    43954786!
     
    44144805!
    44154806!--             Compute diffusion coefficients for momentum and heat
    4416                 km(k,j,i) = c_mu * e(k,j,i)**2 / ( diss(k,j,i) + 1.E-10 ) * flag
     4807                km(k,j,i) = c_0**4 * e(k,j,i)**2 / ( diss(k,j,i) + 1.0E-30_wp ) * flag
    44174808                kh(k,j,i) = km(k,j,i) / prandtl_number * flag
     4809!
     4810!--             Summation for averaged profile of mixing length (cf. flow_statistics)
     4811                DO  sr = 0, statistic_regions
     4812                   sums_l_l(k,sr,tn) = sums_l_l(k,sr,tn) +                     &
     4813                      c_0**3 * e(k,j,i) * SQRT(e(k,j,i)) /                     &
     4814                      ( diss(k,j,i) + 1.0E-30_wp ) * rmask(j,i,sr) * flag
     4815                ENDDO
    44184816
    44194817             ENDDO
     
    44344832!-- so far vertical surfaces require usage of a Prandtl-layer where the boundary
    44354833!-- values of the diffusivities are not needed.
    4436 
    44374834    IF ( .NOT. rans_tke_e )  THEN
    44384835!
     
    44684865          ENDDO
    44694866       ENDDO
     4867!
     4868!--    North- and southward facing surfaces
     4869       DO  n = 0, 1
     4870          DO  m = 1, surf_def_v(n)%ns
     4871             i = surf_def_v(n)%i(m)
     4872             j = surf_def_v(n)%j(m)
     4873             k = surf_def_v(n)%k(m)
     4874             km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dy
     4875             kh(k,j,i) = 1.35_wp * km(k,j,i)
     4876          ENDDO
     4877       ENDDO
     4878!
     4879!--    West- and eastward facing surfaces
     4880       DO  n = 2, 3
     4881          DO  m = 1, surf_def_v(n)%ns
     4882             i = surf_def_v(n)%i(m)
     4883             j = surf_def_v(n)%j(m)
     4884             k = surf_def_v(n)%k(m)
     4885             km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dx
     4886             kh(k,j,i) = 1.35_wp * km(k,j,i)
     4887          ENDDO
     4888       ENDDO
    44704889
    44714890       CALL exchange_horiz( km, nbgp )
     
    44734892
    44744893    ENDIF
     4894
    44754895!
    44764896!-- Model top
     
    45184938    INTEGER(iwp) ::  j      !< loop index y direction
    45194939    INTEGER(iwp) ::  k      !< loop index z direction
    4520     INTEGER, INTENT(IN) ::  mod_count  !<
     4940    INTEGER, INTENT(IN) ::  mod_count  !< flag defining where pointers point to
    45214941
    45224942#if defined( __nopointer )
  • palm/trunk/SOURCE/vertical_nesting_mod.f90

    r3066 r3083  
    40584058!-- Identical timestep for coarse and fine grids
    40594059          dt_3d = MIN( dtc, dtf )
     4060          !> @fixme setting old_dt might be obsolete at this point
     4061          !>   Due to changes in timestep routine, setting of old_dt might be
     4062          !>   not necessary any more at this point. However, could not be
     4063          !>   tested so far.
     4064          !>   2018-05-18, gronemeier
    40604065          old_dt = dt_3d
    40614066#endif
Note: See TracChangeset for help on using the changeset viewer.