Home

Context Navigation

← Previous Change
Next Change →

Changeset 3083 for palm

Timestamp:

Jun 19, 2018 2:03:12 PM (6 years ago)

Author:

gronemeier

Message:

merge with branch rans

Location:

palm/trunk

Files:

: 15 edited

. (modified) (1 prop)
INSTALL/example_cbl_rc (modified) (2 diffs)
SOURCE (modified) (1 prop)
SOURCE/Makefile (modified) (2 diffs, 1 prop)
SOURCE/check_parameters.f90 (modified) (2 diffs, 1 prop)
SOURCE/data_output_mask.f90 (modified) (1 diff)
SOURCE/header.f90 (modified) (4 diffs)
SOURCE/init_3d_model.f90 (modified) (1 prop)
SOURCE/model_1d_mod.f90 (modified) (28 diffs)
SOURCE/modules.f90 (modified) (2 diffs)
SOURCE/parin.f90 (modified) (4 diffs)
SOURCE/pmc_interface_mod.f90 (modified) (1 prop)
SOURCE/timestep.f90 (modified) (8 diffs)
SOURCE/turbulence_closure_mod.f90 (modified) (83 diffs)
SOURCE/vertical_nesting_mod.f90 (modified) (1 diff)

Legend:

: Unmodified
: Added
: Removed

palm/trunk
- Property svn:mergeinfo changed
  /palm/branches/rans merged: 2919,2922,2928-2929,2960,2962,2966,2976,2982,2987-2988,2991,3008-3009,3023,3047,3050,3059,3062,3071-3082

palm/trunk/INSTALL/example_cbl_rc

-                      r3068
+                      r3083
  ******************************    --------------------------------------------
  * PALM 5.0  Rev: 3063:30     *    atmosphere - run without 1D - prerun
+ * PALM 5.0  Rev: 3082        *    atmosphere - run without 1D - prerun
  ******************************    --------------------------------------------
  Date:                 12-06-18    Run:       example_cbl
  Time:                 15:39:03    Run-No.:   00
  Run on host:        imuk_defau
+ Date:                 19-06-18    Run:       example_cbl
+ Time:                 15:53:32    Run-No.:   00
+ Run on host:           default
  Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
  ------------------------------------------------------------------------------
 …
 RUN  ITER. HH:MM:SS.SS    DT(E)     UMAX     VMAX     WMAX     U*    W*      THETA*     Z_I     ENERG.   DISTENERG    DIVOLD     DIVNEW     UMAX(KJI)    VMAX(KJI)    WMAX(KJI)   ADVECX   ADVECY   MGCYC
 ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+      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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+     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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    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
+    262 01:00:07.82  12.8000A   1.9438   2.1487   3.5253  0.115  1.38  -0.997E+00    800.  0.278E+00  0.276E+00  0.332E-03  0.363E-13   17  36  12   15   8  36    8   9  15    0.000    0.000       0

palm/trunk/SOURCE
- Property svn:mergeinfo changed
  /palm/branches/rans/SOURCE merged: 2919,2922,2928-2929,2960,2962,2966,2976,2982,2987-2988,2991,3008-3009,3023,3047,3050,3059,3062,3071-3082

palm/trunk/SOURCE/Makefile

Property svn:mergeinfo changed
/palm/branches/rans/SOURCE/Makefile merged: 2919,2987,3081

-                      r2963
+                      r3083
 # -----------------
 # $Id$
+# add turbulence_closure_mod to parin
+#
+# 2963 2018-04-12 14:47:44Z suehring
 # Introduce index for vegetation/wall, pavement/green-wall and water/window
 # surfaces, for clearer access of surface fraction, albedo, emissivity, etc. .
 …
         spectra_mod.o \
         synthetic_turbulence_generator_mod.o \
+        turbulence_closure_mod.o \
         uv_exposure_model_mod.o \
         vertical_nesting_mod.o \

palm/trunk/SOURCE/check_parameters.f90

Property svn:mergeinfo changed
/palm/branches/rans/SOURCE/check_parameters.f90 merged: 2919,2987,3081

-                      r3065
+                      r3083
 ! Former revisions:
 ! -----------------
+! $Id$
+! $Id: check_parameters.f90 2520 2017-10-05 13:50:26Z gronemeier &
+! Add inital profile output for e (TG)
+!
+! 3065 2018-06-12 07:03:02Z Giersch
 ! dz was replaced by dz(1), error message revised
+!
 …
              ENDIF
           CASE ( 'e' )
+          CASE ( 'e', '#e' )
              dopr_index(i)  = 8
              dopr_unit(i)   = 'm2/s2'
              hom(:,2,8,:)   = SPREAD( zu, 2, statistic_regions+1 )
              hom(nzb,2,8,:) = 0.0_wp
+             IF ( data_output_pr(i)(1:1) == '#' )  THEN
+                dopr_initial_index(i) = 8
+                hom(:,2,8,:)          = SPREAD( zu, 2, statistic_regions+1 )
+                data_output_pr(i)     = data_output_pr(i)(2:)
+             ENDIF
           CASE ( 'km', '#km' )

palm/trunk/SOURCE/data_output_mask.f90

-                      r3045
+                      r3083
 ! -----------------
 ! $Id$
+!
+!
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Error messages revised
+!

palm/trunk/SOURCE/header.f90

-                      r3065
+                      r3083
 ! -----------------
 ! $Id$
+! Print RANS-mode constants
+!
+! 3065 2018-06-12 07:03:02Z Giersch
 ! Header output concerning stretching revised
+!
 …
     USE synthetic_turbulence_generator_mod,                                    &
         ONLY:  stg_header
+    USE turbulence_closure_mod,                                                &
+        ONLY:  rans_const_c, rans_const_sigma
     IMPLICIT NONE
 …
        IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
     ENDIF
+    IF ( rans_mode )  THEN
+       WRITE ( io, 457 )  rans_const_c, rans_const_sigma
+    ENDIF
+!
 !-- Special actions during the run
 …
 FORMAT ('    Day of the year at model start :   day_init = ',I3             &
             /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
+FORMAT ('    RANS-mode constants: c_0 = ',F9.5/         &
+            '                         c_1 = ',F9.5/         &
+            '                         c_2 = ',F9.5/         &
+            '                         c_3 = ',F9.5/         &
+            '                         c_4 = ',F9.5/         &
+            '                         sigma_e    = ',F9.5/  &
+            '                         sigma_diss = ',F9.5)
 FORMAT (//' Actions during the simulation:'/ &
               ' -----------------------------'/)

palm/trunk/SOURCE/init_3d_model.f90

Property svn:mergeinfo changed

/palm/branches/rans/SOURCE/init_3d_model.f90 (added)	merged: 2139,2244,2352-2353,2519-2521,2680,2760-2761,2842,2896,2917,3062,3081
/palm/trunk/SOURCE/init_3d_model.f90	removed

palm/trunk/SOURCE/model_1d_mod.f90

-                      r3049
+                      r3083
 ! -----------------
 ! $Id$
+! Bugfixes:
+!   - preset te_diss and te_e to avoid runtime errors
+!   - implementation of buoyancy term to dissipation
+!     according to Sogachev et al. (2012)
+!   - where diss_p < 0 set diss_p = 0.1 diss
+!   - calculate progn eq(diss) starting from nzb+1
+! Changes:
+!   - add sig_e to TKE equation
+!   - adjust prognostic equation of diss
+!   - set model constants according to Koblitz (2013)
+!   - renamed c_m to c_0
+!   - rename l_black into l1d_init
+!   - calculate l_grid within init_1d_model and save it as l1d_init
+!   - calculate l1d according to DE85 if dissipation is a prognostic value
+!   - made annotations doxygen-readable
+!
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
+!
 …
     USE kinds
+    USE pegrid
+    USE pegrid,                                                                &
+        ONLY:  myid
 …
     LOGICAL ::  stop_dt_1d = .FALSE.             !< termination flag, used in case of too small timestep (1d-model)
     REAL(wp) ::  c_1 = 1.44_wp                 !< model constant
     REAL(wp) ::  c_2 = 1.92_wp                 !< model constant
     REAL(wp) ::  c_3 = 1.44_wp                 !< model constant
     REAL(wp) ::  c_h = 0.0015_wp               !< model constant according to Detering and Etling (1985)
     REAL(wp) ::  c_m = 0.4_wp                  !< model constant, 0.4 according to Detering and Etling (1985)
     REAL(wp) ::  c_mu = 0.09_wp                !< model constant
+    REAL(wp) ::  alpha_buoyancy                !< model constant according to Koblitz (2013)
+    REAL(wp) ::  c_0 = 0.03_wp**0.25_wp        !< model constant according to Koblitz (2013)
+    REAL(wp) ::  c_1 = 1.52_wp                 !< model constant according to Koblitz (2013)
+    REAL(wp) ::  c_2 = 1.83_wp                 !< model constant according to Koblitz (2013)
+    REAL(wp) ::  c_3                           !< model constant
+    REAL(wp) ::  c_mu                          !< model constant
     REAL(wp) ::  damp_level_1d = -1.0_wp       !< namelist parameter
     REAL(wp) ::  dt_1d = 60.0_wp               !< dynamic timestep (1d-model)
 …
     REAL(wp) ::  dt_run_control_1d = 60.0_wp   !< namelist parameter
     REAL(wp) ::  end_time_1d = 864000.0_wp     !< namelist parameter
+    REAL(wp) ::  lambda                        !< maximum mixing length
     REAL(wp) ::  qs1d                          !< characteristic humidity scale (1d-model)
     REAL(wp) ::  simulated_time_1d = 0.0_wp    !< updated simulated time (1d-model)
+    REAL(wp) ::  sig_diss = 1.3_wp             !< model constant
+    REAL(wp) ::  sig_diss = 2.95_wp            !< model constant according to Koblitz (2013)
+    REAL(wp) ::  sig_e = 2.95_wp               !< model constant according to Koblitz (2013)
     REAL(wp) ::  time_pr_1d = 0.0_wp           !< updated simulated time for profile output (1d-model)
     REAL(wp) ::  time_run_control_1d = 0.0_wp  !< updated simulated time for run-control output (1d-model)
 …
     REAL(wp) ::  z01d                          !< roughness length for momentum (1d-model)
     REAL(wp) ::  z0h1d                         !< roughness length for scalars (1d-model)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  diss1d   !< tke dissipation rate (1d-model)
 …
     INTEGER(iwp) ::  k  !< loop index
-    REAL(wp) ::  lambda !< maximum mixing length
+!
 …
+!
 !--       Use the same mixing length as in 3D model (LES-mode)
           !@todo: rename (delete?) this option
           ! As the mixing length is different between RANS and LES mode, it
           ! must be distinguished here between these modes. For RANS mode,
           ! the mixing length is calculated accoding to Blackadar, which is
           ! the other option at this point.
           ! Maybe delete this option entirely (not appropriate in LES case)
           ! 2018-03-20, gronemeier
+          !> @todo rename (delete?) this option
+          !>  As the mixing length is different between RANS and LES mode, it
+          !>  must be distinguished here between these modes. For RANS mode,
+          !>  the mixing length is calculated accoding to Blackadar, which is
+          !>  the other option at this point.
+          !>  Maybe delete this option entirely (not appropriate in LES case)
+          !>  2018-03-20, gronemeier
           DO  k = nzb+1, nzt
              l1d_init(k)  = ( dx * dy * dzw(k) )**0.33333333333333_wp
 …
        us1d = 0.1_wp   ! without initial friction the flow would not change
     ELSE
        diss1d(nzb+1) = 1.0_wp
+       diss1d(nzb+1) = 0.001_wp
        e1d(nzb+1)  = 1.0_wp
        km1d(nzb+1) = 1.0_wp
 …
+!
 !-- Tendencies must be preset in order to avoid runtime errors within the
 !-- first Runge-Kutta step
+!-- Tendencies must be preset in order to avoid runtime errors
+    te_diss  = 0.0_wp
     te_dissm = 0.0_wp
+    te_e  = 0.0_wp
     te_em = 0.0_wp
     te_um = 0.0_wp
 …
+!
 !-- Set model constant
+    IF ( dissipation_1d == 'as_in_3d_model' )  c_m = 0.1_wp
+    IF ( dissipation_1d == 'as_in_3d_model' )  c_0 = 0.1_wp
+    c_mu = c_0**4
+!
 …
 !-- Determine the time step at the start of a 1D-simulation and
 !-- determine and printout quantities used for run control
     dt_1d = 1.0_wp
+    dt_1d = 0.01_wp
     CALL run_control_1d
 …
 !--             dissipation rate
                 IF ( dissipation_1d == 'detering' )  THEN
                    diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
+                   diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
                 ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
                    diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) &
 …
                                      kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1)  &
                                    - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k)    &
                                                  ) * ddzw(k)                   &
+                                                 ) * ddzw(k) / sig_e           &
                                    - diss1d(k)
 …
+!
 !--                dissipation rate
+                   te_diss(k) = km1d(k) *                                      &
+                   IF ( rif1d(k) >= 0.0_wp )  THEN
+                      alpha_buoyancy = 1.0_wp - l1d(k) / lambda
+                   ELSE
+                      alpha_buoyancy = 1.0_wp - ( 1.0_wp + ( c_2 - 1.0_wp )    &
+                                                         / ( c_2 - c_1    ) )  &
+                                              * l1d(k) / lambda
+                   ENDIF
+                   c_3 = ( c_1 - c_2 ) * alpha_buoyancy
+                   te_diss(k) = ( km1d(k) *                                    &
                                   ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2  &
                                   + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2  &
+                                  ) * c_1 * c_mu**0.75 / c_h * f / us1d        &
+                                    * SQRT(e1d(k))                             &
+                                  ) * ( c_1 + (c_2 - c_1) * l1d(k) / lambda )  &
                                   - g / pt_0 * kh1d(k) * flux * c_3            &
+                                    * diss1d(k) / ( e1d(k) + 1.0E-20_wp )      &
+                                  + ( kmzp * ( diss1d(k+1) - diss1d(k) )       &
+                                  - c_2 * diss1d(k)                            &
+                                ) * diss1d(k) / ( e1d(k) + 1.0E-20_wp )        &
+                                + (   kmzp * ( diss1d(k+1) - diss1d(k) )       &
                                            * ddzu(k+1)                         &
                                     - kmzm * ( diss1d(k) - diss1d(k-1) )       &
                                            * ddzu(k)                           &
+                                    ) * ddzw(k) / sig_diss                     &
+                                  - c_2 * diss1d(k)**2 / ( e1d(k) + 1.0E-20_wp )
+                                  ) * ddzw(k) / sig_diss
                 ENDIF
 …
 !--          dissipation rate
              IF ( dissipation_1d == 'detering' )  THEN
                 diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
+                diss1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
              ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
                 diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) )  &
 …
 !--          TKE
              IF ( .NOT. dissipation_1d == 'prognostic' )  THEN
+                !> @query why integrate over 2dz
+                !>   Why is it allowed to integrate over two delta-z for e
+                !>   while for u and v it is not?
+                !>   2018-04-23, gronemeier
                 te_e(k) = km1d(k) * ( ( ( u1d(k+1) - u1d(k-1) ) * dd2zu(k) )**2&
                                     + ( ( v1d(k+1) - v1d(k-1) ) * dd2zu(k) )**2&
 …
                                      kmzp * ( e1d(k+1) - e1d(k) ) * ddzu(k+1)  &
                                    - kmzm * ( e1d(k) - e1d(k-1) ) * ddzu(k)    &
                                                  ) * ddzw(k)                   &
+                                                 ) * ddzw(k) / sig_e           &
                                    - diss1d(k)
              ENDIF
 …
              ENDDO
-             IF ( dissipation_1d == 'prognostic' )  THEN
-                DO  k = nzb_diff, nzt
-                   diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + &
-                                                 tsc(3) * te_dissm(k) )
-                ENDDO
-             ENDIF
+!
 !--          Eliminate negative TKE values, which can result from the
 …
 !--          value is reduced to 10 percent of its old value.
              WHERE ( e1d_p < 0.0_wp )  e1d_p = 0.1_wp * e1d
+             IF ( dissipation_1d == 'prognostic' )  THEN
+                DO  k = nzb+1, nzt
+                   diss1d_p(k) = diss1d(k) + dt_1d * ( tsc(2) * te_diss(k) + &
+                                                       tsc(3) * te_dissm(k) )
+                ENDDO
+                WHERE ( diss1d_p < 0.0_wp )  diss1d_p = 0.1_wp * diss1d
+             ENDIF
           ENDIF
 …
                    IF ( dissipation_1d == 'prognostic' )  THEN
                       DO k = nzb+1, nzt
+                         te_dissm(k) = -9.5625_wp * te_diss(k) + 5.3125_wp * te_dissm(k)
+                         te_dissm(k) = -9.5625_wp * te_diss(k)  &
+                                     +  5.3125_wp * te_dissm(k)
                       ENDDO
                    ENDIF
 …
              ENDIF
           ENDIF
+!
 …
              ENDIF    ! constant_flux_layer
+             !> @todo combine if clauses
+             !>   The previous and following if clauses can be combined into a
+             !>   single clause
+             !>   2018-04-23, gronemeier
+!
 !--          Compute the Richardson-flux numbers,
 …
 !--             compatibility with the 3D model.
                 IF ( ibc_e_b == 2 )  THEN
                    e1d(nzb+1) = ( us1d / c_m )**2
+                   e1d(nzb+1) = ( us1d / c_0 )**2
                 ENDIF
                 IF ( dissipation_1d == 'prognostic' )  THEN
                    e1d(nzb+1) = us1d**2 / SQRT( c_mu )
+                   e1d(nzb+1) = ( us1d / c_0 )**2
                    diss1d(nzb+1) = us1d**3 / ( kappa * zu(nzb+1) )
                    diss1d(nzb) = diss1d(nzb+1)
 …
 !--          in the dissipation of TKE via l1d_diss. Otherwise, km1d would be
 !--          too large.
+             IF ( mixing_length_1d == 'blackadar' )  THEN
+             IF ( dissipation_1d /= 'prognostic' )  THEN
+                IF ( mixing_length_1d == 'blackadar' )  THEN
+                   DO  k = nzb+1, nzt
+                      IF ( rif1d(k) >= 0.0_wp )  THEN
+                         l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
+                         l1d_diss(k) = l1d(k)
+                      ELSE
+                         l1d(k) = l1d_init(k)
+                         l1d_diss(k) = l1d_init(k) *                           &
+                                       SQRT( 1.0_wp - 16.0_wp * rif1d(k) )
+                      ENDIF
+                   ENDDO
+                ELSEIF ( mixing_length_1d == 'as_in_3d_model' )  THEN
+                   DO  k = nzb+1, nzt
+                      dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k)
+                      IF ( dpt_dz > 0.0_wp )  THEN
+                         l_stable = 0.76_wp * SQRT( e1d(k) )                   &
+                                  / SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp
+                      ELSE
+                         l_stable = l1d_init(k)
+                      ENDIF
+                      l1d(k) = MIN( l1d_init(k), l_stable )
+                      l1d_diss(k) = l1d(k)
+                   ENDDO
+                ENDIF
+             ELSE
                 DO  k = nzb+1, nzt
+                   IF ( rif1d(k) >= 0.0_wp )  THEN
+                      l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
+                      l1d_diss(k) = l1d(k)
+                   ELSE
+                      l1d(k) = l1d_init(k)
+                      l1d_diss(k) = l1d_init(k) *                              &
+                                    SQRT( 1.0_wp - 16.0_wp * rif1d(k) )
+                   ENDIF
+                ENDDO
+             ELSEIF ( mixing_length_1d == 'as_in_3d_model' )  THEN
+                DO  k = nzb+1, nzt
+                   dpt_dz = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k)
+                   IF ( dpt_dz > 0.0_wp )  THEN
+                      l_stable = 0.76_wp * SQRT( e1d(k) ) /                    &
+                                     SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp
+                   ELSE
+                      l_stable = l1d_init(k)
+                   ENDIF
+                   l1d(k) = MIN( l1d_init(k), l_stable )
+                   l1d_diss(k) = l1d(k)
+                   l1d(k) = c_0**3 * e1d(k) * SQRT( e1d(k) )                   &
+                          / ( diss1d(k) + 1.0E-30_wp )
                 ENDDO
              ENDIF
 …
                 ENDIF
              ENDIF
              IF ( dissipation_1d == 'prognostic' )  THEN
                 DO  k = nzb_diff, nzt
 …
              ELSE
                 DO  k = nzb_diff, nzt
                    km1d(k) = c_m * SQRT( e1d(k) ) * l1d(k)
+                   km1d(k) = c_0 * SQRT( e1d(k) ) * l1d(k)
                 ENDDO
              ENDIF
 …
     REAL(wp) ::  dt_diff  !< time step accorind to diffusion criterion
+    REAL(wp) ::  dt_old   !< previous time step
     REAL(wp) ::  fac      !< factor of criterion
     REAL(wp) ::  value    !< auxiliary variable
+!
+!-- Save previous time step
+    dt_old = dt_1d
+!
 !-- Compute the currently feasible time step according to the diffusion
 !-- criterion. At nzb+1 the half grid length is used.
     fac = 0.125  !0.35_wp                                                       !### changed from 0.35
+    fac = 0.125
     dt_diff = dt_max_1d
     DO  k = nzb+2, nzt
 …
+!
+!-- Limit the new time step to a maximum of 10 times the previous time step
+    dt_1d = MIN( dt_old * 10.0_wp, dt_1d )
+!
 !-- Set flag when the time step becomes too small
     IF ( dt_1d < ( 0.00001_wp * dt_max_1d ) )  THEN
+    IF ( dt_1d < ( 1.0E-15_wp * dt_max_1d ) )  THEN
        stop_dt_1d = .TRUE.

palm/trunk/SOURCE/modules.f90

-                      r3065
+                      r3083
 ! -----------------
 ! $Id$
+! set dt_3d = 0.01
+!
+! 3065 2018-06-12 07:03:02Z Giersch
 ! Variables concerning stretching introduced or revised
+!
 …
     REAL(wp) ::  dt_run_control = 60.0_wp                      !< namelist parameter
     REAL(wp) ::  dt_spinup = 60.0_wp                           !< namelist parameter
     REAL(wp) ::  dt_3d = 1.0_wp                                !< time step
+    REAL(wp) ::  dt_3d = 0.01_wp                               !< time step
     REAL(wp) ::  dz_max = 1000.0_wp                            !< namelist parameter
     REAL(wp) ::  dz_stretch_factor = 1.08_wp                   !< namelist parameter

palm/trunk/SOURCE/parin.f90

-                      r3065
+                      r3083
 ! -----------------
 ! $Id$
+! Added rans_const_c and rans_const_sigma as input parameters (TG)
+!
+! 3065 2018-06-12 07:03:02Z Giersch
 ! New initialization parameters added
+!
 …
     USE synthetic_turbulence_generator_mod,                                    &
         ONLY:  stg_parin
+    USE turbulence_closure_mod,                                                &
+        ONLY:  rans_const_c, rans_const_sigma
     USE urban_surface_mod,                                                     &
 …
              pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
              q_vertical_gradient, q_vertical_gradient_level,                   &
+             random_generator, random_heatflux, rans_mode,                     &
+             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
+             rans_mode,                                                        &
              rayleigh_damping_factor, rayleigh_damping_height,                 &
              recycling_width, recycling_yshift,                                &
 …
              pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
              q_vertical_gradient, q_vertical_gradient_level,                   &
+             random_generator, random_heatflux, rans_mode,                     &
+             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
+             rans_mode,                                                        &
              rayleigh_damping_factor, rayleigh_damping_height,                 &
              recycling_width, recycling_yshift,                                &

palm/trunk/SOURCE/pmc_interface_mod.f90
- Property svn:mergeinfo deleted

palm/trunk/SOURCE/timestep.f90

-                      r3049
+                      r3083
 ! -----------------
 ! $Id$
+! limit dt_3d to be at maximum 2*old_dt; define old_dt at beginning of routine
+! Add km/kh_max
+!
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
+!
 …
     INTEGER(iwp) ::  j !<
     INTEGER(iwp) ::  k !<
+    INTEGER(iwp) ::  km_max_ijk(3) = -1  !< index values (i,j,k) of location where km_max occurs
+    INTEGER(iwp) ::  kh_max_ijk(3) = -1  !< index values (i,j,k) of location where kh_max occurs
     LOGICAL ::  stop_dt_local !< local switch for controlling the time stepping
 …
     REAL(wp) ::  dt_w              !<
     REAL(wp) ::  dt_w_l            !<
+    REAL(wp) ::  km_max            !< maximum of Km in entire domain
+    REAL(wp) ::  kh_max            !< maximum of Kh in entire domain
     REAL(wp) ::  u_gtrans_l        !<
     REAL(wp) ::  u_max_l           !<
 …
     CALL cpu_log( log_point(12), 'calculate_timestep', 'start' )
+!
+!--    Save former time step as reference
+       old_dt = dt_3d
+!
 …
 !--    The time step is the minimum of the 3-4 components and the diffusion time
 !--    step minus a reduction (cfl_factor) to be on the safe side.
+!--    The time step must not exceed the maximum allowed value.
+!--    The time step must not exceed the maximum allowed value and must not
+!--    increase by more than a factor of 10.
        dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w, dt_precipitation )
        dt_3d = MIN( dt_3d, dt_max )
+       dt_3d = MIN( dt_3d, dt_max, 2.0_wp * old_dt )
+!
 …
        IF ( dt_3d < ( 0.00001_wp * dt_max ) )  THEN
           stop_dt = .TRUE.
+!
+!--       Determine the maxima of the diffusion coefficients, including their
+!--       grid index positions.
+          CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, km, 'abs',  &
+.0_wp, km_max, km_max_ijk )
+          CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, kh, 'abs',  &
+.0_wp, kh_max, kh_max_ijk )
           WRITE( message_string, * ) 'Time step has reached minimum limit.',   &
 …
                '&dt_w            = ', dt_w, ' s',                              &
                '&dt_diff         = ', dt_diff, ' s',                           &
                '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),        &
+               '&u_max           = ', u_max, ' m/s    k=', u_max_ijk(1),       &
                '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                     &
                '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),        &
+               '&v_max           = ', v_max, ' m/s    k=', v_max_ijk(1),       &
                '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                     &
+               '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),        &
+               '  j=', w_max_ijk(2), '  i=', w_max_ijk(3)
+               '&w_max           = ', w_max, ' m/s    k=', w_max_ijk(1),       &
+               '  j=', w_max_ijk(2), '  i=', w_max_ijk(3),                     &
+               '&km_max          = ', km_max, ' m2/s2  k=', km_max_ijk(1),     &
+               '  j=', km_max_ijk(2), '  i=', km_max_ijk(3),                   &
+               '&kh_max          = ', kh_max, ' m2/s2  k=', kh_max_ijk(1),     &
+                '  j=', kh_max_ijk(2), '  i=', kh_max_ijk(3)
           CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 )
+!
 …
        dt_3d = NINT( dt_3d * 100.0_wp / div ) * div / 100.0_wp
+!
-!--    Adjust the time step
-       old_dt = dt_3d
     ENDIF

palm/trunk/SOURCE/turbulence_closure_mod.f90

-                      r3045
+                      r3083
 ! -----------------
 ! $Id$
+! - set limits of diss at the end of prognostic equations
+! - call production_e to calculate production term of diss
+! - limit change of diss to -90% to +100%
+! - remove factor 0.5 from diffusion_diss_ij
+! - rename c_m into c_0, and c_h into c_4
+! - add rans_const_c and rans_const_sigma as namelist parameters
+! - add calculation of mixing length for profile output in case of rans_tke_e
+! - changed format of annotations to comply with doxygen standards
+! - calculate and save dissipation rate during rans_tke_l mode
+! - set bc at vertical walls for e, diss, km, kh
+! - bugfix: set l_wall = 0.0 within buildings
+! - set l_wall at bottom and top boundary (rans-mode)
+! - bugfix in production term for dissipation rate
+! - bugfix in diffusion of dissipation rate
+! - disable check for 1D model if rans_tke_e is used
+! - bugfixes for initialization (rans-mode):
+!    - correction of dissipation-rate formula
+!    - calculate km based on l_wall
+!    - initialize diss if 1D model is not used
+!
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Error message revised
+!
 …
 ! Further todo's
+!
 ! 2936 2018-03-27 14:49:27Z suehring
+! 2936 2018-03-27 14:49:27Z gronemeier
 ! - defined l_grid only within this module
 ! - Moved l_wall definition from modules.f90
 …
 !>       add OpenMP directives whereever possible
 !>       remove debug output variables (dummy1, dummy2, dummy3)
-!> @todo Move initialization of wall-mixing length from init_grid
 !> @todo Check for random disturbances
 !> @note <Enter notes on the module>
 …
+    REAL(wp) ::  c_1  = 1.44_wp    !< model constant for RANS mode
+    REAL(wp) ::  c_2  = 1.92_wp    !< model constant for RANS mode
+    REAL(wp) ::  c_3  = 1.44_wp    !< model constant for RANS mode
+    REAL(wp) ::  c_h  = 0.0015_wp  !< model constant for RANS mode
+    REAL(wp) ::  c_m               !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES)
+    REAL(wp) ::  c_mu = 0.09_wp    !< model constant for RANS mode
+    REAL(wp) ::  l_max             !< maximum length scale for Blackadar mixing length
+    REAL(wp) ::  sig_e = 1.0_wp    !< factor to calculate Ke from Km
+    REAL(wp) ::  sig_diss = 1.3_wp !< factor to calculate K_diss from Km
+    INTEGER(iwp) ::  surf_e        !< end index of surface elements at given i-j position
+    INTEGER(iwp) ::  surf_s        !< start index of surface elements at given i-j position
+    REAL(wp) ::  c_0                !< constant used for diffusion coefficient and dissipation (dependent on mode RANS/LES)
+    REAL(wp) ::  c_1                !< model constant for RANS mode
+    REAL(wp) ::  c_2                !< model constant for RANS mode
+    REAL(wp) ::  c_3                !< model constant for RANS mode
+    REAL(wp) ::  c_4                !< model constant for RANS mode
+    REAL(wp) ::  l_max              !< maximum length scale for Blackadar mixing length
+    REAL(wp) ::  dsig_e = 1.0_wp    !< factor to calculate Ke from Km (1/sigma_e)
+    REAL(wp) ::  dsig_diss = 1.0_wp !< factor to calculate K_diss from Km (1/sigma_diss)
+    INTEGER(iwp) ::  surf_e         !< end index of surface elements at given i-j position
+    INTEGER(iwp) ::  surf_s         !< start index of surface elements at given i-j position
+    REAL(wp), DIMENSION(0:4) :: rans_const_c = &       !< model constants for RANS mode (namelist param)
+       (/ 0.55_wp, 1.44_wp, 1.92_wp, 0.0_wp, 0.0_wp /) !> default values fit for standard-tke-e closure
+    REAL(wp), DIMENSION(2) :: rans_const_sigma = &     !< model constants for RANS mode, sigma values (sigma_e, sigma_diss) (namelist param)
+       (/ 1.0_wp, 0.77_wp /)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  l_black    !< mixing length according to Blackadar
 …
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  l_wall !< near-wall mixing length
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy1 !< debug output variable
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy2 !< debug output variable
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dummy3 !< debug output variable
+    PUBLIC c_m, dummy1, dummy2, dummy3
+    !> @todo remove debug variables
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: diss_prod1, diss_adve1, diss_diff1, &
+                                               diss_prod2, diss_adve2, diss_diff2, &
+                                               diss_prod3, diss_adve3, diss_diff3, &
+                                               dummy1, dummy2, dummy3
+    PUBLIC c_0, rans_const_c, rans_const_sigma
+!
 …
 ! ------------
 !> Check parameters routine for turbulence closure module.
+!> @todo remove rans_mode from initialization namelist and rework checks
+!>   The way it is implemented at the moment, the user has to set two variables
+!>   so that the RANS mode is working. It would be better if only one parameter
+!>   has to be set.
+!>   2018-06-18, gronemeier
 !------------------------------------------------------------------------------!
  SUBROUTINE tcm_check_parameters
 …
     IF ( rans_mode )  THEN
+       c_m = 0.4_wp  !according to Detering and Etling (1985)
+!
+!--    Assign values to constants for RANS mode
+       dsig_e    = 1.0_wp / rans_const_sigma(1)
+       dsig_diss = 1.0_wp / rans_const_sigma(2)
+       c_0 = rans_const_c(0)
+       c_1 = rans_const_c(1)
+       c_2 = rans_const_c(2)
+       !c_3 = rans_const_c(3)   !> @todo clarify how to switch between different models
+       c_4 = rans_const_c(4)
        SELECT CASE ( TRIM( turbulence_closure ) )
 …
           CASE ( 'TKE-e' )
              rans_tke_e = .TRUE.
-             IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) == 0  &
-                  .AND.  .NOT.  nest_domain )  THEN
-                message_string = 'Initializing without 1D model while ' //     &
-                                 'using TKE-e closure&' //                     &
-                                 'is not possible at the moment!'
-                CALL message( 'tcm_check_parameters', 'TG0005', 1, 2, 0, 6, 0 )
-             ENDIF
           CASE DEFAULT
              message_string = 'Unknown turbulence closure: ' //                &
                               TRIM( turbulence_closure )
              CALL message( 'tcm_check_parameters', 'TG0001', 1, 2, 0, 6, 0 )
+             CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 )
        END SELECT
+       IF ( turbulent_inflow .OR. turbulent_outflow )  THEN
+          message_string = 'turbulent inflow/outflow is not yet '//            &
+                           'implemented for RANS mode'
+          CALL message( 'tcm_check_parameters', 'PA0501', 1, 2, 0, 6, 0 )
+       ENDIF
        message_string = 'RANS mode is still in development! ' //               &
                         '&Not all features of PALM are yet compatible '//      &
                         'with RANS mode. &Use at own risk!'
        CALL message( 'tcm_check_parameters', 'TG0003', 0, 1, 0, 6, 0 )
+       CALL message( 'tcm_check_parameters', 'PA0502', 0, 1, 0, 6, 0 )
     ELSE
+       c_m = 0.1_wp !according to Lilly (1967) and Deardorff (1980)
+       c_0 = 0.1_wp !according to Lilly (1967) and Deardorff (1980)
+       dsig_e = 1.0_wp !assure to use K_m to calculate TKE instead
+                       !of K_e which is used in RANS mode
        SELECT CASE ( TRIM( turbulence_closure ) )
 …
           CASE DEFAULT
+             !> @todo rework this part so that only one call of this error exists
              message_string = 'Unknown turbulence closure: ' //                &
                               TRIM( turbulence_closure )
              CALL message( 'tcm_check_parameters', 'TG0001', 1, 2, 0, 6, 0 )
+             CALL message( 'tcm_check_parameters', 'PA0500', 1, 2, 0, 6, 0 )
        END SELECT
-    ENDIF
-    IF ( rans_tke_e )  THEN
-       IF ( turbulent_inflow .OR. turbulent_outflow )  THEN
-          message_string = 'turbulent inflow/outflow is not yet '//            &
-                           'implemented for TKE-e closure'
-          CALL message( 'tcm_check_parameters', 'TG0002', 1, 2, 0, 6, 0 )
-       ENDIF
     ENDIF
 …
     IMPLICIT NONE
     CHARACTER (LEN=*) ::  unit     !<
     CHARACTER (LEN=*) ::  var      !<
     INTEGER(iwp) ::  i      !<
     INTEGER(iwp) ::  ilen   !<
     INTEGER(iwp) ::  k      !<
+    CHARACTER (LEN=*) ::  unit     !< unit of output variable
+    CHARACTER (LEN=*) ::  var      !< name of output variable
+    INTEGER(iwp) ::  i      !< index of var in data_output
+    INTEGER(iwp) ::  ilen   !< length of var string
+    INTEGER(iwp) ::  k      !< flag if var contains one of '_xy', '_xz' or '_yz'
     SELECT CASE ( TRIM( var ) )
        CASE ( 'diss' )
-          IF ( .NOT.  rans_tke_e )  THEN
-             message_string = 'output of "' // TRIM( var ) // '" requi' //  &
-                      'res TKE-e closure for RANS mode.'
-             CALL message( 'tcm_check_data_output', 'TG0101', 1, 2, 0, 6, 0 )
-          ENDIF
           unit = 'm2/s3'
+       CASE ( 'dummy2', 'dummy3', 'dummy1' )
+          unit = 'mixing length'
+       CASE ( 'diss1', 'diss2',                         &                      !> @todo remove later
+              'diss_prod1', 'diss_adve1', 'diss_diff1', &
+              'diss_prod2', 'diss_adve2', 'diss_diff2', &
+              'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3'  )
+          unit = 'debug output'
        CASE ( 'kh', 'km' )
 …
     IMPLICIT NONE
     CHARACTER (LEN=*), INTENT(OUT) ::  grid_x   !<
     CHARACTER (LEN=*), INTENT(OUT) ::  grid_y   !<
     CHARACTER (LEN=*), INTENT(OUT) ::  grid_z   !<
     CHARACTER (LEN=*), INTENT(IN)  ::  var      !<
     LOGICAL, INTENT(OUT) ::  found   !<
+    CHARACTER (LEN=*), INTENT(OUT) ::  grid_x   !< x grid of output variable
+    CHARACTER (LEN=*), INTENT(OUT) ::  grid_y   !< y grid of output variable
+    CHARACTER (LEN=*), INTENT(OUT) ::  grid_z   !< z grid of output variable
+    CHARACTER (LEN=*), INTENT(IN)  ::  var      !< name of output variable
+    LOGICAL, INTENT(OUT) ::  found   !< flag if output variable is found
     found  = .TRUE.
 …
           grid_z = 'zu'
+       CASE ( 'dummy2', 'dummy3', 'dummy1' )                                    !### remove later
+       CASE ( 'diss1', 'diss2',                         &                       !> @todo remove later
+              'diss_prod1', 'diss_adve1', 'diss_diff1', &
+              'diss_prod2', 'diss_adve2', 'diss_diff2', &
+              'diss_prod3', 'diss_adve3', 'diss_diff3', 'dummy3' )
           grid_x = 'x'
           grid_y = 'y'
 …
     IMPLICIT NONE
     CHARACTER (LEN=*) ::  mode       !<
     CHARACTER (LEN=*) ::  variable   !<
     INTEGER(iwp) ::  i   !<
     INTEGER(iwp) ::  j   !<
     INTEGER(iwp) ::  k   !<
+    CHARACTER (LEN=*) ::  mode       !< flag defining mode 'allocate', 'sum' or 'average'
+    CHARACTER (LEN=*) ::  variable   !< name of variable
+    INTEGER(iwp) ::  i   !< loop index
+    INTEGER(iwp) ::  j   !< loop index
+    INTEGER(iwp) ::  k   !< loop index
     IF ( mode == 'allocate' )  THEN
 …
     IMPLICIT NONE
     CHARACTER (LEN=*) ::  grid       !<
     CHARACTER (LEN=*) ::  mode       !<
     CHARACTER (LEN=*) ::  variable   !<
     INTEGER(iwp) ::  av   !<
     INTEGER(iwp) ::  i    !<
     INTEGER(iwp) ::  j    !<
     INTEGER(iwp) ::  k    !<
     INTEGER(iwp) ::  nzb_do   !<
     INTEGER(iwp) ::  nzt_do   !<
     LOGICAL ::  found   !<
+    CHARACTER (LEN=*) ::  grid       !< name of vertical grid
+    CHARACTER (LEN=*) ::  mode       !< either 'xy', 'xz' or 'yz'
+    CHARACTER (LEN=*) ::  variable   !< name of variable
+    INTEGER(iwp) ::  av   !< flag for (non-)average output
+    INTEGER(iwp) ::  i    !< loop index
+    INTEGER(iwp) ::  j    !< loop index
+    INTEGER(iwp) ::  k    !< loop index
+    INTEGER(iwp) ::  nzb_do   !< vertical output index (bottom)
+    INTEGER(iwp) ::  nzt_do   !< vertical output index (top)
+    LOGICAL ::  found   !< flag if output variable is found
     LOGICAL ::  two_d   !< flag parameter that indicates 2D variables (horizontal cross sections)
 …
     IMPLICIT NONE
     CHARACTER (LEN=*) ::  variable   !<
     INTEGER(iwp) ::  av     !<
     INTEGER(iwp) ::  i      !<
     INTEGER(iwp) ::  j      !<
     INTEGER(iwp) ::  k      !<
+    CHARACTER (LEN=*) ::  variable   !< name of variable
+    INTEGER(iwp) ::  av     !< flag for (non-)average output
+    INTEGER(iwp) ::  i      !< loop index
+    INTEGER(iwp) ::  j      !< loop index
+    INTEGER(iwp) ::  k      !< loop index
     INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
     INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)
     LOGICAL ::  found   !<
+    LOGICAL ::  found   !< flag if output variable is found
     REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
 …
           ENDIF
+       CASE ( 'dummy1' )                                                        !### remove later
+       CASE ( 'dummy3' )                                                        !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = dummy3(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss1' )                                                         !> @todo remove later
           IF ( av == 0 )  THEN
              DO  i = nxl, nxr
 …
           ENDIF
        CASE ( 'dummy2' )                                                        !### remove later
+       CASE ( 'diss2' )                                                         !> @todo remove later
           IF ( av == 0 )  THEN
              DO  i = nxl, nxr
 …
           ENDIF
        CASE ( 'dummy3' )                                                        !### remove later
+       CASE ( 'diss_prod1' )                                                    !> @todo remove later
           IF ( av == 0 )  THEN
              DO  i = nxl, nxr
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = dummy3(k,j,i)
+                      local_pf(i,j,k) = diss_prod1(k,j,i)
                    ENDDO
                 ENDDO
 …
           ENDIF
+       CASE ( 'diss_adve1' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_adve1(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_diff1' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_diff1(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_prod2' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_prod2(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_adve2' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_adve2(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_diff2' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_diff2(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_prod3' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_prod3(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_adve3' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_adve3(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
+       CASE ( 'diss_diff3' )                                                    !> @todo remove later
+          IF ( av == 0 )  THEN
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   DO  k = nzb_do, nzt_do
+                      local_pf(i,j,k) = diss_diff3(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
        CASE DEFAULT
           found = .FALSE.
 …
     ALLOCATE( km(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )                           !### remove later
+    ALLOCATE( dummy1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )                           !> @todo remove later
     ALLOCATE( dummy2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     ALLOCATE( dummy3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    IF ( rans_mode )  ALLOCATE( l_black(nzb:nzt+1) )
+    ALLOCATE( diss_adve1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_adve2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_adve3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_prod1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_prod2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_prod3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_diff1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_diff2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ALLOCATE( diss_diff3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    dummy1 = 0.0_wp
+    dummy2 = 0.0_wp
+    dummy3 = 0.0_wp
+    diss_adve1 = 0.0_wp
+    diss_adve2 = 0.0_wp
+    diss_adve3 = 0.0_wp
+    diss_prod1 = 0.0_wp
+    diss_prod2 = 0.0_wp
+    diss_prod3 = 0.0_wp
+    diss_diff1 = 0.0_wp
+    diss_diff2 = 0.0_wp
+    diss_diff3 = 0.0_wp
 #if defined( __nopointer )
 …
 !-- they do not necessarily need to be transferred, which is attributed to
 !-- the design of the model coupler which allocates memory for each variable.
     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.  wang_kernel  .OR.      &
+    IF ( rans_mode  .OR.  use_sgs_for_particles  .OR.  wang_kernel  .OR.       &
          collision_turbulence  .OR.  nested_run )  THEN
 #if defined( __nopointer )
 …
     e  => e_1;   e_p  => e_2;   te_m  => e_3
     IF ( rans_tke_e  .OR.  use_sgs_for_particles  .OR.     &
+    IF ( rans_mode  .OR.  use_sgs_for_particles  .OR.     &
          wang_kernel  .OR.  collision_turbulence  .OR.  nested_run )  THEN
        diss => diss_1
 …
     INTEGER(iwp) :: j            !< loop index
     INTEGER(iwp) :: k            !< loop index
     INTEGER(iwp) :: nz_s_shift   !<
     INTEGER(iwp) :: nz_s_shift_l !<
+    INTEGER(iwp) :: nz_s_shift   !< lower shift index for scalars
+    INTEGER(iwp) :: nz_s_shift_l !< local lower shift index in case of turbulent inflow
+!
 !-- Initialize mixing length
     CALL tcm_init_mixing_length
+    dummy3 = l_wall                 !> @todo remove later
+!
 …
           IF ( rans_tke_e )  THEN
              IF ( dissipation_1d == 'prognostic' )  THEN    !### Why must this be checked?
                 DO  i = nxlg, nxrg                          !### Should 'diss' not always
                    DO  j = nysg, nyng                       !### be prognostic in case rans_tke_e?
+             IF ( dissipation_1d == 'prognostic' )  THEN    !> @query Why must this be checked?
+                DO  i = nxlg, nxrg                          !>   Should 'diss' not always
+                   DO  j = nysg, nyng                       !>   be prognostic in case rans_tke_e?
                       diss(:,j,i) = diss1d
                    ENDDO
 …
                    DO  j = nysg, nyng
                       DO  k = nzb+1, nzt
                          diss(k,j,i) = e(k,j,i) * SQRT( e(k,j,i) ) / l1d(k)
+                         diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km1d(k)
                       ENDDO
                    ENDDO
 …
           IF ( constant_diffusion )  THEN
              km   = km_constant
              kh   = km / prandtl_number
              e    = 0.0_wp
+             km = km_constant
+             kh = km / prandtl_number
+             e  = 0.0_wp
           ELSEIF ( e_init > 0.0_wp )  THEN
+             DO  k = nzb+1, nzt
+                km(k,:,:) = c_m * l_grid(k) * SQRT( e_init )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb+1, nzt
+                      km(k,j,i) = c_0 * l_wall(k,j,i) * SQRT( e_init )
+                   ENDDO
+                ENDDO
              ENDDO
              km(nzb,:,:)   = km(nzb+1,:,:)
              km(nzt+1,:,:) = km(nzt,:,:)
              kh   = km / prandtl_number
              e    = e_init
+             kh = km / prandtl_number
+             e  = e_init
           ELSE
              IF ( .NOT. ocean )  THEN
 …
           ENDIF
+          IF ( rans_tke_e )  THEN
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb+1, nzt
+                      diss(k,j,i) = c_0**4 * e(k,j,i)**2 / km(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+             diss(nzb,:,:) = diss(nzb+1,:,:)
+             diss(nzt+1,:,:) = diss(nzt,:,:)
+          ENDIF
        ENDIF
+!
 …
              ENDDO
           ENDDO
+          IF ( rans_tke_e )  THEN
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   nz_s_shift = get_topography_top_index_ji( j, i, 's' )
+                   diss(nz_s_shift:nzt+1,j,i) = diss(0:nzt+1-nz_s_shift,j,i)
+                ENDDO
+             ENDDO
+          ENDIF
        ENDIF
 …
 !--       boundary and adjust mean inflow profiles
           IF ( complex_terrain )  THEN
+             IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
+             IF ( nxlg <= 0 .AND. nxrg >= 0 .AND.  &
+                  nysg <= 0 .AND. nyng >= 0        )  THEN
                 nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
              ELSE
 …
              nz_s_shift = nz_s_shift_l
 #endif
+             mean_inflow_profiles(nz_s_shift:nzt+1,5) = hom_sum(0:nzt+1-nz_s_shift,8,0)  ! e
+             mean_inflow_profiles(nz_s_shift:nzt+1,5) =  &
+                hom_sum(0:nzt+1-nz_s_shift,8,0)  ! e
           ENDIF
+!
 …
+!
 !--       Inside buildings set TKE back to zero.
 !--       Other scalars (km, kh, diss, ...) are ignored at present,
+!--       Other scalars (km, kh,...) are ignored at present,
 !--       maybe revise later.
           DO  i = nxlg, nxrg
 …
                    e(k,j,i)     = MERGE( e(k,j,i), 0.0_wp,                     &
                                          BTEST( wall_flags_0(k,j,i), 0 ) )
-                   te_m(k,j,i)  = MERGE( te_m(k,j,i), 0.0_wp,                  &
-                                         BTEST( wall_flags_0(k,j,i), 0 ) )
                 ENDDO
              ENDDO
           ENDDO
+          IF ( rans_tke_e )  THEN
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb, nzt
+                      diss(k,j,i)    = MERGE( diss(k,j,i), 0.0_wp,             &
+                                              BTEST( wall_flags_0(k,j,i), 0 ) )
+                   ENDDO
+                ENDDO
+             ENDDO
+          ENDIF
        ENDIF
+!
 …
+!
 !--    Allthough tendency arrays are set in prognostic_equations, they have
 !--    have to be predefined here because they are used (but multiplied with 0)
 !--    there before they are set.
+!--    to be predefined here because there they are used (but multiplied with 0)
+!--    before they are set.
        te_m = 0.0_wp
+       IF ( rans_tke_e )  THEN
+          diss_p = diss
+          tdiss_m = 0.0_wp
+       ENDIF
     ENDIF
 …
        DO  k = 1, nzt
           IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.           &
+          IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.            &
                l_grid(k) > 1.5_wp * dy * wall_adjustment_factor )  THEN
              WRITE( message_string, * ) 'grid anisotropy exceeds ',            &
                                         'threshold given by only local',       &
                                         ' horizontal reduction of near_wall ', &
                                         'mixing length l_wall',                &
                                         ' starting from height level k = ', k, &
+             WRITE( message_string, * ) 'grid anisotropy exceeds ',             &
+                                        'threshold given by only local',        &
+                                        ' &horizontal reduction of near_wall ', &
+                                        'mixing length l_wall',                 &
+                                        ' &starting from height level k = ', k, &
                                         '.'
              CALL message( 'init_grid', 'PA0202', 0, 1, 0, 6, 0 )
 …
+!
 !-- Initialize the mixing length in case of a RANS simulation
+       ALLOCATE( l_black(nzb:nzt+1) )
+!
 !--    Calculate mixing length according to Blackadar (1962)
        IF ( f /= 0.0_wp )  THEN
           l_max = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) /               &
                   ABS( f ) + 1E-10_wp
+          l_max = 2.7E-4_wp * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) /            &
+                  ABS( f ) + 1.0E-10_wp
        ELSE
           l_max = 30.0_wp
 …
+!
 !--    Gather topography information of whole domain
        !## TODO: reduce amount of data sent by MPI call
        !##  By now, a whole global 3D-array is sent and received with
        !##  MPI_ALLREDUCE although most of the array is 0. This can be
        !##  drastically reduced if only the local subarray is sent and stored
        !##  in a global array. For that, an MPI data type or subarray must be
        !##  defined.
        !##  2018-03-19, gronemeier
+       !> @todo reduce amount of data sent by MPI call
+       !>   By now, a whole global 3D-array is sent and received with
+       !>   MPI_ALLREDUCE although most of the array is 0. This can be
+       !>   drastically reduced if only the local subarray is sent and stored
+       !>   in a global array. For that, an MPI data type or subarray must be
+       !>   defined.
+       !>   2018-03-19, gronemeier
        ALLOCATE( wall_flags_0_global(nzb:nzt+1,0:ny,0:nx) )
 …
           ENDDO
        ENDDO
+       l_wall(nzb,:,:) = l_black(nzb)
+       l_wall(nzt+1,:,:) = l_black(nzt+1)
+!
 !--    Limit mixing length to either nearest wall or Blackadar mixing length.
 …
           IF ( rad_k_b /= 0 .OR. rad_k_t /= 0 )  THEN
              !## NOTE: shape of vicinity is larger in z direction
              !##  Shape of vicinity is two grid points larger than actual search
              !##  radius in vertical direction. The first and last grid point is
              !##   always set to 1 to asure correct detection of topography. See
              !##  function "shortest_distance" for details.
              !##  2018-03-16, gronemeier
+             !> @note shape of vicinity is larger in z direction
+             !>   Shape of vicinity is two grid points larger than actual search
+             !>   radius in vertical direction. The first and last grid point is
+             !>   always set to 1 to asure correct detection of topography. See
+             !>   function "shortest_distance" for details.
+             !>   2018-03-16, gronemeier
              ALLOCATE( vicinity(-rad_k-1:rad_k+1,-rad_j:rad_j,-rad_i:rad_i) )
              ALLOCATE( vic_yz(0:rad_k+1,0:rad_j) )
 …
                    ELSE  !Check if (i,j,k) belongs to atmosphere
                       l_wall(k,j,i) = -999.0
+                      l_wall(k,j,i) = l_black(k)
                    ENDIF
 …
 !> (pos_i/jj/kk), where (jj/kk) is the position of the maximum of 'array'
 !> closest to the origin (0/0) of 'array'.
+!> @todo this part of PALM does not reproduce the same results for optimized
+!>   and debug options for the compiler. This should be fixed
 !------------------------------------------------------------------------------!
     REAL FUNCTION shortest_distance( array, orientation, pos_i )
 …
        IF ( orientation ) THEN  !if array is oriented upwards
           DO  jj = 0, rad_j
+             shortest_distance = MIN( shortest_distance,                        &
+                                      SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2         &
+                                          + MAX(jj*dy-0.5*dy,0.0)**2            &
+                                          + MAX(zw(loc_k(jj)+k-1)-zu(k),0.0)**2 &
+                                          )                                     &
+                                    )
+             shortest_distance =                                               &
+                MIN( shortest_distance,                                        &
+                     SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2   &
+                         + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2      &
+                         + MAX(zw(loc_k(jj)+k-1)-zu(k), 0.0_wp)**2             &
+                         )                                                     &
+                   )
           ENDDO
        ELSE  !if array is oriented downwards
           !## NOTE: MAX within zw required to circumvent error at domain border
           !##  At the domain border, if non-cyclic boundary is present, the
           !##  index for zw could be -1, which will be errorneous (zw(-1) does
           !##  not exist). The MAX function limits the index to be at least 0.
+          !> @note MAX within zw required to circumvent error at domain border
+          !>   At the domain border, if non-cyclic boundary is present, the
+          !>   index for zw could be -1, which will be errorneous (zw(-1) does
+          !>   not exist). The MAX function limits the index to be at least 0.
           DO  jj = 0, rad_j
+             shortest_distance = MIN( shortest_distance,                       &
+                                      SQRT( MAX(pos_i*dx-0.5*dx,0.0)**2        &
+                                          + MAX(jj*dy-0.5*dy,0.0)**2           &
+                                          + MAX(zu(k)-zw(MAX(k-loc_k(jj),      &
+_iwp)),          &
+.0)**2                        &
+                                          )                                    &
+                                    )
+             shortest_distance =                                               &
+                MIN( shortest_distance,                                        &
+                     SQRT( MAX(REAL(pos_i, KIND=wp)*dx-0.5_wp*dx, 0.0_wp)**2   &
+                         + MAX(REAL(jj, KIND=wp)*dy-0.5_wp*dy, 0.0_wp)**2      &
+                         + MAX(zu(k)-zw(MAX(k-loc_k(jj),0_iwp)), 0.0_wp)**2    &
+                         )                                                     &
+                   )
           ENDDO
        ENDIF
     END FUNCTION
 …
+!
 !--    Use special boundary condition in case of TKE-e closure
+       !> @todo do the same for usm and lsm surfaces
+       !>   2018-06-05, gronemeier
        IF ( rans_tke_e )  THEN
           DO  i = nxl, nxr
 …
                 DO  m = surf_s, surf_e
                    k = surf_def_h(0)%k(m)
                    e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_m**2
+                   e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_0**2
                 ENDDO
              ENDDO
 …
              DO  k = nzb+1, nzt
 !                tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k)
                 tend(k,j,i) = tend(k,j,i) + c_1 * c_mu * f / c_h               &  !### needs revision
+                tend(k,j,i) = tend(k,j,i) + c_1 * c_0**4 * f / c_4               &  !> @todo needs revision
                       / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))       &
                       * SQRT(e(k,j,i)) * produc(k,j,i)
 …
+!
 !--    Additional sink term for flows through plant canopies
 !        IF ( plant_canopy )  CALL pcm_tendency( ? )                            !### what to do with this?
 !        CALL user_actions( 'diss-tendency' )                                   !### not yet implemented
+!        IF ( plant_canopy )  CALL pcm_tendency( ? )                            !> @query what to do with this?
+!        CALL user_actions( 'diss-tendency' )                                   !> @todo not yet implemented
+!
 …
     USE arrays_3d,                                                             &
         ONLY:  ddzu, diss_l_diss, diss_l_e, diss_s_diss, diss_s_e,             &
+               flux_l_diss, flux_l_e, flux_s_diss, flux_s_e
+               flux_l_diss, flux_l_e, flux_s_diss, flux_s_e,&
+               u_p,v_p,w_p
     USE control_parameters,                                                    &
         ONLY:  f, tsc
+    USE grid_variables,                                                        &
+        ONLY:  dx, dy
     USE surface_mod,                                                           &
 …
                 surf_usm_v
+    use indices, only: nx, ny
     IMPLICIT NONE
     INTEGER(iwp) ::  i       !< loop index x direction
     INTEGER(iwp) ::  i_omp   !<
+    INTEGER(iwp) ::  i_omp   !< first loop index of i-loop in prognostic_equations
     INTEGER(iwp) ::  j       !< loop index y direction
     INTEGER(iwp) ::  k       !< loop index z direction
+    INTEGER(iwp) ::  l       !< loop index
     INTEGER(iwp) ::  m       !< loop index
     INTEGER(iwp) ::  surf_e  !< end index of surface elements at given i-j position
     INTEGER(iwp) ::  surf_s  !< start index of surface elements at given i-j position
+    INTEGER(iwp) ::  tn      !<
+    REAL(wp), DIMENSION(nzb:nzt+1) :: advec  !< advection term of TKE tendency
+    REAL(wp), DIMENSION(nzb:nzt+1) :: produc !< production term of TKE tendency
+    INTEGER(iwp) ::  tn      !< task number of openmp task
+    INTEGER(iwp) :: pis = 32 !< debug variable, print from i=pis                !> @todo remove later
+    INTEGER(iwp) :: pie = 32 !< debug variable, print until i=pie               !> @todo remove later
+    INTEGER(iwp) :: pjs = 26 !< debug variable, print from j=pjs                !> @todo remove later
+    INTEGER(iwp) :: pje = 26 !< debug variable, print until j=pje               !> @todo remove later
+    INTEGER(iwp) :: pkb = 1  !< debug variable, print from k=pkb                !> @todo remove later
+    INTEGER(iwp) :: pkt = 7  !< debug variable, print until k=pkt               !> @todo remove later
+    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_adv   !< debug variable               !> @todo remove later
+    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_pro   !< debug variable               !> @todo remove later
+    REAL(wp), DIMENSION(nzb:nzt+1) :: dum_dif   !< debug variable               !> @todo remove later
+FORMAT(A,7(1X,E12.5))   !> @todo remove later
+!
 …
        ENDIF
        advec(:) = tend(:,j,i)
        CALL production_e( i, j )
        produc(:) = tend(:,j,i) - advec(:)
+       dum_adv = tend(:,j,i)                                                    !> @todo remove later
+       CALL production_e( i, j, .FALSE. )
+       dum_pro = tend(:,j,i) - dum_adv                                          !> @todo remove later
        IF ( .NOT. humidity )  THEN
 …
           CALL diffusion_e( i, j, vpt, pt_reference )
        ENDIF
+       dum_dif = tend(:,j,i) - dum_adv - dum_pro                                !> @todo remove later
+!
 …
           DO  m = surf_s, surf_e
              k = surf_def_h(0)%k(m)
+             e_p(k,j,i) = surf_def_h(0)%us(m)**2 / c_m**2
+             e_p(k,j,i) = ( surf_def_h(0)%us(m) / c_0 )**2
+          ENDDO
+          DO  l = 0, 3
+             surf_s = surf_def_v(l)%start_index(j,i)
+             surf_e = surf_def_v(l)%end_index(j,i)
+             DO  m = surf_s, surf_e
+                k = surf_def_v(l)%k(m)
+                e_p(k,j,i) = ( surf_def_v(l)%us(m) / c_0 )**2
+             ENDDO
           ENDDO
        ENDIF
 …
           ENDIF
        ENDIF
+!        if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then        !> @todo remove later
+!           WRITE(9, *) '------'
+!           WRITE(9, '(A,F8.3,1X,F8.3,1X,I2)') 't, dt, int_ts:', simulated_time, dt_3d, intermediate_timestep_count
+!           WRITE(9, *) 'i:', i
+!           WRITE(9, *) 'j:', j
+!           WRITE(9, *) 'k:', pkb, ' - ', pkt
+!           WRITE(9, *) '---'
+!           WRITE(9, *) 'e:'
+!           WRITE(9, 5555) 'adv :', dum_adv(pkb:pkt)
+!           WRITE(9, 5555) 'pro :', dum_pro(pkb:pkt)
+!           WRITE(9, 5555) 'dif :', dum_dif(pkb:pkt)
+!           WRITE(9, 5555) 'tend:', tend(pkb:pkt,j,i)
+!           WRITE(9, 5555) 'e_p :', e_p(pkb:pkt,j,i)
+!           WRITE(9, 5555) 'e   :', e(pkb:pkt,j,i)
+!           FLUSH(9)
+!        endif
     ENDIF   ! TKE equation
 …
        ENDIF
+       IF ( intermediate_timestep_count == 1 )  diss_adve1(:,j,i) = tend(:,j,i) !> @todo remove later
+       IF ( intermediate_timestep_count == 2 )  diss_adve2(:,j,i) = tend(:,j,i)
+       IF ( intermediate_timestep_count == 3 )  diss_adve3(:,j,i) = tend(:,j,i)
+!
 !--    Production of TKE dissipation rate
+       DO  k = nzb+1, nzt
+!              tend(k,j,i) = tend(k,j,i) + c_1 * diss(k,j,i) / ( e(k,j,i) + 1.0E-20_wp ) * produc(k)
+          tend(k,j,i) = tend(k,j,i) + c_1 * c_mu * f / c_h                     &  !### needs revision
+                      / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))       &
+                      * SQRT(e(k,j,i)) * produc(k)
+       ENDDO
+       CALL production_e( i, j, .TRUE. )
+       IF ( intermediate_timestep_count == 1 )  diss_prod1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) !> @todo remove later
+       IF ( intermediate_timestep_count == 2 )  diss_prod2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i)
+       IF ( intermediate_timestep_count == 3 )  diss_prod3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i)
+       dum_pro = tend(:,j,i) - dum_adv                                          !> @todo remove later
+!
+!--    Diffusion term of TKE dissipation rate
        CALL diffusion_diss( i, j )
+       IF ( intermediate_timestep_count == 1 )  diss_diff1(:,j,i) = tend(:,j,i) - diss_adve1(:,j,i) - diss_prod1(:,j,i) !> @todo remove later
+       IF ( intermediate_timestep_count == 2 )  diss_diff2(:,j,i) = tend(:,j,i) - diss_adve2(:,j,i) - diss_prod2(:,j,i)
+       IF ( intermediate_timestep_count == 3 )  diss_diff3(:,j,i) = tend(:,j,i) - diss_adve3(:,j,i) - diss_prod3(:,j,i)
+       IF ( intermediate_timestep_count == 3 )  dummy3(:,j,i) = km(:,j,i)
+       dum_dif = tend(:,j,i) - dum_adv - dum_pro                                !> @todo remove later
+!
 !--    Additional sink term for flows through plant canopies
 !        IF ( plant_canopy )  CALL pcm_tendency( i, j, ? )                      !### not yet implemented
 !        CALL user_actions( i, j, 'diss-tendency' )                             !### not yet implemented
+!        IF ( plant_canopy )  CALL pcm_tendency( i, j, ? )                      !> @todo not yet implemented
+!        CALL user_actions( i, j, 'diss-tendency' )                             !> @todo not yet implemented
+!
 …
                                                 BTEST( wall_flags_0(k,j,i), 0 )&
+                                               )
-          IF ( diss_p(k,j,i) <= 0.0_wp )  diss_p(k,j,i) = 0.1_wp * diss(k,j,i)
        ENDDO
 …
        ENDDO
+       DO  l = 0, 1
+          surf_s = surf_def_v(l)%start_index(j,i)
+          surf_e = surf_def_v(l)%end_index(j,i)
+          DO  m = surf_s, surf_e
+             k = surf_def_v(l)%k(m)
+             diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dy )
+          ENDDO
+       ENDDO
+       DO  l = 2, 3
+          surf_s = surf_def_v(l)%start_index(j,i)
+          surf_e = surf_def_v(l)%end_index(j,i)
+          DO  m = surf_s, surf_e
+             k = surf_def_v(l)%k(m)
+             diss_p(k,j,i) = surf_def_v(l)%us(m)**3 / ( kappa * 0.5_wp * dx )
+          ENDDO
+       ENDDO
+!
 !--    Calculate tendencies for the next Runge-Kutta step
 …
           ENDIF
        ENDIF
+!        IF ( intermediate_timestep_count == 1 )  dummy1(:,j,i) = e_p(:,j,i)
+!        IF ( intermediate_timestep_count == 2 )  dummy2(:,j,i) = e_p(:,j,i)
+!        IF ( intermediate_timestep_count == 3 )  dummy3(:,j,i) = e_p(:,j,i)
+!
+!--    Limit change of diss to be between -90% and +100%. Also, set an absolute
+!--    minimum value
+       DO  k = nzb, nzt+1
+          diss_p(k,j,i) = MIN( MAX( diss_p(k,j,i),         &
+.1_wp * diss(k,j,i),  &
+.0001_wp ),           &
+.0_wp * diss(k,j,i) )
+       ENDDO
+       IF ( intermediate_timestep_count == 1 )  dummy1(:,j,i) = diss_p(:,j,i)   !> @todo remove later
+       IF ( intermediate_timestep_count == 2 )  dummy2(:,j,i) = diss_p(:,j,i)
+!        if ( i >= pis .and. i <= pie .and. j >= pjs .and. j <= pje ) then        !> @todo remove later
+!           WRITE(9, *) '---'
+!           WRITE(9, *) 'diss:'
+!           WRITE(9, 5555) 'adv   :', dum_adv(pkb:pkt)
+!           WRITE(9, 5555) 'pro   :', dum_pro(pkb:pkt)
+!           WRITE(9, 5555) 'dif   :', dum_dif(pkb:pkt)
+!           WRITE(9, 5555) 'tend  :', tend(pkb:pkt,j,i)
+!           WRITE(9, 5555) 'diss_p:', diss_p(pkb:pkt,j,i)
+!           WRITE(9, 5555) 'diss  :', diss(pkb:pkt,j,i)
+!           WRITE(9, *) '---'
+!           WRITE(9, 5555) 'km    :', km(pkb:pkt,j,i)
+!           flush(9)
+!        endif
     ENDIF   ! dissipation equation
 …
 !> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is
 !>          not considered well!
+!> @todo Adjust production term in case of rans_tke_e simulation
 !------------------------------------------------------------------------------!
  SUBROUTINE production_e
 …
 !> @warning The case with constant_flux_layer = F and use_surface_fluxes = T is
 !>          not considered well!
+!> @todo non-neutral case is not yet considered for RANS mode
 !------------------------------------------------------------------------------!
  SUBROUTINE production_e_ij( i, j )
+ SUBROUTINE production_e_ij( i, j, diss_production )
     USE arrays_3d,                                                             &
 …
     IMPLICIT NONE
+    LOGICAL :: diss_production
     INTEGER(iwp) ::  i       !< running index x-direction
 …
     REAL(wp), DIMENSION(nzb+1:nzt)  ::  dwdy        !< Gradient of w-component in y-direction
     REAL(wp), DIMENSION(nzb+1:nzt)  ::  dwdz        !< Gradient of w-component in z-direction
+    REAL(wp), DIMENSION(nzb+1:nzt)  ::  tend_temp   !< temporal tendency
     IF ( constant_flux_layer )  THEN
 …
        ENDDO
+!        IF ( .NOT. rans_tke_e )  THEN
        DO  k = nzb+1, nzt
+          def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +            &
+                           dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +            &
+                           dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +            &
+.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) )
+!
+!--       Production term according to Kato and Launder (1993)
+!           def = SQRT( ( dudx(k) + dudy(k) + dudz(k) +                       &
+!                         dvdx(k) + dvdy(k) + dvdz(k) +                       &
+!                         dwdx(k) + dwdy(k) + dwdz(k)                         &
+!                       )**4 -                                                &
+!                       ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 +              &
+!                         2.0_wp * ( dudy(k) * dvdx(k) +                      &
+!                                    dudz(k) * dwdx(k) +                      &
+!                                    dvdz(k) * dwdy(k) )                      &
+!                       )**2                                                  &
+!                     )
+          def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +         &
+                           dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +         &
+                           dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +         &
+.0_wp * ( dvdx(k)*dudy(k) +                                &
+                           dwdx(k)*dudz(k) +                                &
+                           dwdy(k)*dvdz(k) )
           IF ( def < 0.0_wp )  def = 0.0_wp
 …
           flag  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
           tend(k,j,i) = tend(k,j,i) + km(k,j,i) * def * flag
+          tend_temp(k) = km(k,j,i) * def * flag
        ENDDO
+!        ELSE
+!
+!           DO  k = nzb+1, nzt
+! !
+! !--          Production term according to Kato and Launder (1993)
+!              def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
+!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 +              &
+!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
+!                                       dvdz(k) * dwdy(k) +                      &
+!                                       dwdx(k) * dudz(k) )       )              &
+!                        * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
+!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 -              &
+!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
+!                                       dvdz(k) * dwdy(k) +                      &
+!                                       dwdx(k) * dudz(k) )       )              &
+!                        )
+!
+!              IF ( def < 0.0_wp )  def = 0.0_wp
+!
+!              flag  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
+!
+!              tend_temp(k) = km(k,j,i) * def * flag
+!
+!           ENDDO
+!
+!        ENDIF
+    ELSE  ! not constant_flux_layer
+!        IF ( .NOT. rans_tke_e )  THEN
+!
+!--       Calculate TKE production by shear. Here, no additional
+!--       wall-bounded code is considered.
+!--       Why?
+          DO  k = nzb+1, nzt
+             dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
+             dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                &
+                                    u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
+             dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                &
+                                    u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
+             dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                &
+                                    v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
+             dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
+             dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                &
+                                    v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
+             dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                &
+                                    w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
+             dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                &
+                                    w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
+             dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
+             def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +         &
+                              dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +         &
+                              dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +         &
+.0_wp * ( dvdx(k)*dudy(k) +                                &
+                              dwdx(k)*dudz(k) +                                &
+                              dwdy(k)*dvdz(k) )
+             IF ( def < 0.0_wp )  def = 0.0_wp
+             flag        = MERGE( 1.0_wp, 0.0_wp,                              &
+                                  BTEST( wall_flags_0(k,j,i), 29 ) )
+             tend_temp(k) = km(k,j,i) * def * flag
+          ENDDO
+!        ELSE
+!
+!           DO  k = nzb+1, nzt
+!
+!              dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
+!              dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                &
+!                                     u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
+!              dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                &
+!                                     u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
+!
+!              dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                &
+!                                     v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
+!              dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
+!              dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                &
+!                                     v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
+!
+!              dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                &
+!                                     w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
+!              dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                &
+!                                     w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
+!              dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
+! !
+! !--          Production term according to Kato and Launder (1993)
+!              def = SQRT( ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
+!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 +              &
+!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
+!                                       dvdz(k) * dwdy(k) +                      &
+!                                       dwdx(k) * dudz(k) )       )              &
+!                        * ( dudy(k)**2 + dvdz(k)**2 + dwdx(k)**2 +              &
+!                            dudz(k)**2 + dvdx(k)**2 + dwdy(k)**2 -              &
+!                            2.0_wp * ( dudy(k) * dvdx(k) +                      &
+!                                       dvdz(k) * dwdy(k) +                      &
+!                                       dwdx(k) * dudz(k) )       )              &
+!                        )
+!
+!              IF ( def < 0.0_wp )  def = 0.0_wp
+!
+!              flag        = MERGE( 1.0_wp, 0.0_wp,                              &
+!                                   BTEST( wall_flags_0(k,j,i), 29 ) )
+!              tend_temp(k) = km(k,j,i) * def * flag
+!
+!           ENDDO
+!
+!        ENDIF
+    ENDIF
+    IF ( .NOT. diss_production )  THEN
+!
+!--    Production term in case of TKE production
+       DO  k = nzb+1, nzt
+          tend(k,j,i) = tend(k,j,i) + tend_temp(k)
+       ENDDO
     ELSE
+!
+!--    Calculate TKE production by shear. Here, no additional
+!--    wall-bounded code is considered.
+!--    Why?
+!--    Production term in case of dissipation-rate production (rans_tke_e)
        DO  k = nzb+1, nzt
+          dudx(k)  =           ( u(k,j,i+1) - u(k,j,i)     ) * ddx
+          dudy(k)  = 0.25_wp * ( u(k,j+1,i) + u(k,j+1,i+1) -                   &
+                                 u(k,j-1,i) - u(k,j-1,i+1) ) * ddy
+          dudz(k)  = 0.5_wp  * ( u(k+1,j,i) + u(k+1,j,i+1) -                   &
+                                 u(k-1,j,i) - u(k-1,j,i+1) ) * dd2zu(k)
+          dvdx(k)  = 0.25_wp * ( v(k,j,i+1) + v(k,j+1,i+1) -                   &
+                                 v(k,j,i-1) - v(k,j+1,i-1) ) * ddx
+          dvdy(k)  =           ( v(k,j+1,i) - v(k,j,i)     ) * ddy
+          dvdz(k)  = 0.5_wp  * ( v(k+1,j,i) + v(k+1,j+1,i) -                   &
+                                 v(k-1,j,i) - v(k-1,j+1,i) ) * dd2zu(k)
+          dwdx(k)  = 0.25_wp * ( w(k,j,i+1) + w(k-1,j,i+1) -                   &
+                                 w(k,j,i-1) - w(k-1,j,i-1) ) * ddx
+          dwdy(k)  = 0.25_wp * ( w(k,j+1,i) + w(k-1,j+1,i) -                   &
+                                 w(k,j-1,i) - w(k-1,j-1,i) ) * ddy
+          dwdz(k)  =           ( w(k,j,i)   - w(k-1,j,i)   ) * ddzw(k)
+          def = 2.0_wp * ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 ) +            &
+                           dudy(k)**2 + dvdx(k)**2 + dwdx(k)**2   +            &
+                           dwdy(k)**2 + dudz(k)**2 + dvdz(k)**2   +            &
+.0_wp * ( dvdx(k)*dudy(k) + dwdx(k)*dudz(k) + dwdy(k)*dvdz(k) )
+!
+!--       Production term according to Kato and Launder (1993)
+!           def = SQRT( ( dudx(k) + dudy(k) + dudz(k) +                        &
+!                         dvdx(k) + dvdy(k) + dvdz(k) +                        &
+!                         dwdx(k) + dwdy(k) + dwdz(k)                          &
+!                       )**4 -                                                 &
+!                       ( dudx(k)**2 + dvdy(k)**2 + dwdz(k)**2 +               &
+!                         2.0_wp * ( dudy(k) * dvdx(k) +                       &
+!                                    dudz(k) * dwdx(k) +                       &
+!                                    dvdz(k) * dwdy(k) )                       &
+!                       )**2                                                   &
+!                     )
+          IF ( def < 0.0_wp )  def = 0.0_wp
+          flag        = MERGE( 1.0_wp, 0.0_wp,                                 &
+                               BTEST( wall_flags_0(k,j,i), 29 ) )
+          tend(k,j,i) = tend(k,j,i) + km(k,j,i) * def * flag
+          ! Standard TKE-e closure
+          tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i)               &
+                                    /( e(k,j,i) + 1.0E-20_wp )                 &
+                                    * c_1
+!           ! Production according to Koblitz (2013)
+!           tend(k,j,i) = tend(k,j,i) + tend_temp(k) * diss(k,j,i)               &
+!                                     /( e(k,j,i) + 1.0E-20_wp )                 &
+!                                     * ( c_1 + ( c_2 - c_1 )                    &
+!                                             * l_wall(k,j,i) / l_max )
+!           ! Production according to Detering and Etling (1985)
+!           !> @todo us is not correct if there are vertical walls
+!           tend(k,j,i) = tend(k,j,i) + tend_temp(k) * SQRT(e(k,j,i))            &
+!                                     * c_1 * c_0**3 / c_4 * f                   &
+!              / surf_def_h(0)%us(surf_def_h(0)%start_index(j,i))
        ENDDO
     ENDIF
 …
     REAL(wp)     ::  l              !< mixing length
     REAL(wp)     ::  ll             !< adjusted l
     REAL(wp)     ::  var_reference  !<
+    REAL(wp)     ::  var_reference  !< reference temperature
 #if defined( __nopointer )
     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
+    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
 #else
     REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !<
+    REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !< temperature
 #endif
     REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  dissipation  !< TKE dissipation
 …
        DO  j = nys, nyn
           DO  k = nzb+1, nzt
+!
+!--          Predetermine flag to mask topography
+             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
+!
 …
                 CALL mixing_length_rans( i, j, k, l, ll, var, var_reference )
+                dissipation(k,j) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
+                dissipation(k,j) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
+                diss(k,j,i) = dissipation(k,j) * flag
              ELSEIF ( rans_tke_e )  THEN
 …
              ENDIF
+!
+!--          Predetermine flag to mask topography
+             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
+             tend(k,j,i) = tend(k,j,i)                                         &
+                                     + (                                       &
+             tend(k,j,i) = tend(k,j,i) + (                                     &
+                                           (                                   &
                        ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) )     &
                      - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) )     &
                                        ) * ddx2  * flag                        &
                                      + (                                       &
+                                           ) * ddx2  * flag                    &
+                                         + (                                   &
                        ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) )     &
                      - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) )     &
                                        ) * ddy2  * flag                        &
                                      + (                                       &
+                                           ) * ddy2  * flag                    &
+                                         + (                                   &
             ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1)    &
                                                           * rho_air_zw(k)      &
           - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)      &
                                                           * rho_air_zw(k-1)    &
+                                       ) * ddzw(k) * drho_air(k) * flag        &
+                                           ) * ddzw(k) * drho_air(k)           &
+                                         ) * flag * dsig_e                     &
                           - dissipation(k,j) * flag
 …
     REAL(wp)     ::  l              !< mixing length
     REAL(wp)     ::  ll             !< adjusted l
     REAL(wp)     ::  var_reference  !<
+    REAL(wp)     ::  var_reference  !< reference temperature
 #if defined( __nopointer )
     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
+    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
 #else
     REAL(wp), DIMENSION(:,:,:), POINTER ::  var     !<
+    REAL(wp), DIMENSION(:,:,:), POINTER ::  var     !< temperature
 #endif
     REAL(wp), DIMENSION(nzb+1:nzt) ::  dissipation  !< dissipation of TKE
 …
           CALL mixing_length_rans( i, j, k, l, ll, var, var_reference  )
+          dissipation(k) = c_m**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
+          dissipation(k) = c_0**3 * e(k,j,i) * SQRT( e(k,j,i) ) / ll
+          diss(k,j,i) = dissipation(k) * flag
        ELSEIF ( rans_tke_e )  THEN
 …
+!
 !--    Calculate the tendency term
        tend(k,j,i) = tend(k,j,i)                                               &
                                     + (                                        &
+       tend(k,j,i) = tend(k,j,i) + (                                           &
+                                      (                                        &
                       ( km(k,j,i)+km(k,j,i+1) ) * ( e(k,j,i+1)-e(k,j,i) )      &
                     - ( km(k,j,i)+km(k,j,i-1) ) * ( e(k,j,i)-e(k,j,i-1) )      &
                                       ) * ddx2  * flag / sig_e                 &
+                                      ) * ddx2                                 &
                                     + (                                        &
                       ( km(k,j,i)+km(k,j+1,i) ) * ( e(k,j+1,i)-e(k,j,i) )      &
                     - ( km(k,j,i)+km(k,j-1,i) ) * ( e(k,j,i)-e(k,j-1,i) )      &
                                       ) * ddy2  * flag / sig_e                 &
+                                      ) * ddy2                                 &
                                     + (                                        &
            ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1)     &
 …
          - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)       &
                                                          * rho_air_zw(k-1)     &
+                                      ) * ddzw(k) * drho_air(k) * flag / sig_e &
+                                    - dissipation(k) * flag
+                                      ) * ddzw(k) * drho_air(k)                &
+                                   ) * flag * dsig_e                           &
+                                 - dissipation(k) * flag
     ENDDO
 …
              flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
              tend(k,j,i) = tend(k,j,i)                                         &
                                + (                                             &
+             tend(k,j,i) = tend(k,j,i) +                                       &
+                         (       (                                             &
                  ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) )     &
                - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) )     &
                                  ) * ddx2  * flag                              &
+                                 ) * ddx2                                      &
                                + (                                             &
                  ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) )     &
                - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) )     &
                                  ) * ddy2  * flag                              &
+                                 ) * ddy2                                      &
                                + (                                             &
       ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1)    &
 …
     - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k)      &
                                                     * rho_air_zw(k-1)          &
+                                 ) * ddzw(k) * drho_air(k) * flag              &
+                               - c_2 * diss(k,j,i)**2                          &
+                                     / ( e(k,j,i) + 1.0E-20_wp ) * flag
+                                 ) * ddzw(k) * drho_air(k)                     &
+                         ) * flag * dsig_diss                                  &
+                         - c_2 * diss(k,j,i)**2                                &
+                               / ( e(k,j,i) + 1.0E-20_wp ) * flag
           ENDDO
 …
     IMPLICIT NONE
+    INTEGER(iwp) ::  i              !< running index x direction
+    INTEGER(iwp) ::  j              !< running index y direction
+    INTEGER(iwp) ::  k              !< running index z direction
+    REAL(wp)     ::  flag           !< flag to mask topography
+    REAL(wp), DIMENSION(nzb+1:nzt) :: tend_temp
+    INTEGER(iwp) ::  i         !< running index x direction
+    INTEGER(iwp) ::  j         !< running index y direction
+    INTEGER(iwp) ::  k         !< running index z direction
+    REAL(wp)     ::  flag      !< flag to mask topography
+!
 …
+!
 !--    Calculate the tendency term
+       tend_temp(k) =           (                                              &
+       tend(k,j,i) =  tend(k,j,i) +                                            &
+                   (            (                                              &
                 ( km(k,j,i)+km(k,j,i+1) ) * ( diss(k,j,i+1)-diss(k,j,i) )      &
               - ( km(k,j,i)+km(k,j,i-1) ) * ( diss(k,j,i)-diss(k,j,i-1) )      &
                                 ) * ddx2  * flag / sig_diss                    &
+                                ) * ddx2                                       &
                               + (                                              &
                 ( km(k,j,i)+km(k,j+1,i) ) * ( diss(k,j+1,i)-diss(k,j,i) )      &
               - ( km(k,j,i)+km(k,j-1,i) ) * ( diss(k,j,i)-diss(k,j-1,i) )      &
                                 ) * ddy2  * flag / sig_diss                    &
+                                ) * ddy2                                       &
                               + (                                              &
      ( km(k,j,i)+km(k+1,j,i) ) * ( diss(k+1,j,i)-diss(k,j,i) ) * ddzu(k+1)     &
 …
    - ( km(k,j,i)+km(k-1,j,i) ) * ( diss(k,j,i)-diss(k-1,j,i) ) * ddzu(k)       &
                                                    * rho_air_zw(k-1)           &
+                                ) * ddzw(k) * drho_air(k) * flag / sig_diss    &
+                              - c_2 * diss(k,j,i)**2                           &
+                                    / ( e(k,j,i) + 1.0E-20_wp ) * flag
+       tend(k,j,i) = tend(k,j,i) + tend_temp(k)
+                                ) * ddzw(k) * drho_air(k)                      &
+                   ) * flag * dsig_diss                                        &
+                   - c_2 * diss(k,j,i)**2 / ( e(k,j,i) + 1.0E-20_wp ) * flag
     ENDDO
 …
 !> Computation of the turbulent diffusion coefficients for momentum and heat
 !> according to Prandtl-Kolmogorov.
+!> @todo consider non-default surfaces
 !------------------------------------------------------------------------------!
  SUBROUTINE tcm_diffusivities( var, var_reference )
 …
         ONLY:  e_min, outflow_l, outflow_n, outflow_r, outflow_s
+    USE grid_variables,                                                        &
+        ONLY:  dx, dy
     USE statistics,                                                            &
         ONLY :  rmask, sums_l_l
     USE surface_mod,                                                           &
         ONLY :  bc_h, surf_def_h
+        ONLY :  bc_h, surf_def_h, surf_def_v
     IMPLICIT NONE
     INTEGER(iwp) ::  i                   !<
     INTEGER(iwp) ::  j                   !<
     INTEGER(iwp) ::  k                   !<
     INTEGER(iwp) ::  m                   !<
     INTEGER(iwp) ::  n                   !<
     INTEGER(iwp) ::  omp_get_thread_num  !<
     INTEGER(iwp) ::  sr                  !<
     INTEGER(iwp) ::  tn                  !<
     REAL(wp)     ::  flag                !<
     REAL(wp)     ::  l                   !<
     REAL(wp)     ::  ll                  !<
     REAL(wp)     ::  var_reference       !<
+    INTEGER(iwp) ::  i                   !< loop index
+    INTEGER(iwp) ::  j                   !< loop index
+    INTEGER(iwp) ::  k                   !< loop index
+    INTEGER(iwp) ::  m                   !< loop index
+    INTEGER(iwp) ::  n                   !< loop index
+    INTEGER(iwp) ::  omp_get_thread_num  !< opemmp function to get thread number
+    INTEGER(iwp) ::  sr                  !< statistic region
+    INTEGER(iwp) ::  tn                  !< thread number
+    REAL(wp)     ::  flag                !< topography flag
+    REAL(wp)     ::  l                   !< mixing length
+    REAL(wp)     ::  ll                  !< adjusted mixing length
+    REAL(wp)     ::  var_reference       !< reference temperature
 #if defined( __nopointer )
     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !<
+    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  var  !< temperature
 #else
     REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !<
+    REAL(wp), DIMENSION(:,:,:), POINTER ::  var  !< temperature
 #endif
 …
+!
 !--             Compute diffusion coefficients for momentum and heat
                 km(k,j,i) = c_m * l * SQRT( e(k,j,i) ) * flag
+                km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag
                 kh(k,j,i) = ( 1.0_wp + 2.0_wp * l / ll ) * km(k,j,i) * flag
+!
 …
+!
 !--             Compute diffusion coefficients for momentum and heat
                 km(k,j,i) = c_m * l * SQRT( e(k,j,i) ) * flag
+                km(k,j,i) = c_0 * l * SQRT( e(k,j,i) ) * flag
                 kh(k,j,i) = km(k,j,i) / prandtl_number * flag
+!
 …
+!
 !--             Compute diffusion coefficients for momentum and heat
                 km(k,j,i) = c_mu * e(k,j,i)**2 / ( diss(k,j,i) + 1.E-10 ) * flag
+                km(k,j,i) = c_0**4 * e(k,j,i)**2 / ( diss(k,j,i) + 1.0E-30_wp ) * flag
                 kh(k,j,i) = km(k,j,i) / prandtl_number * flag
+!
+!--             Summation for averaged profile of mixing length (cf. flow_statistics)
+                DO  sr = 0, statistic_regions
+                   sums_l_l(k,sr,tn) = sums_l_l(k,sr,tn) +                     &
+                      c_0**3 * e(k,j,i) * SQRT(e(k,j,i)) /                     &
+                      ( diss(k,j,i) + 1.0E-30_wp ) * rmask(j,i,sr) * flag
+                ENDDO
              ENDDO
 …
 !-- so far vertical surfaces require usage of a Prandtl-layer where the boundary
 !-- values of the diffusivities are not needed.
     IF ( .NOT. rans_tke_e )  THEN
+!
 …
           ENDDO
        ENDDO
+!
+!--    North- and southward facing surfaces
+       DO  n = 0, 1
+          DO  m = 1, surf_def_v(n)%ns
+             i = surf_def_v(n)%i(m)
+             j = surf_def_v(n)%j(m)
+             k = surf_def_v(n)%k(m)
+             km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dy
+             kh(k,j,i) = 1.35_wp * km(k,j,i)
+          ENDDO
+       ENDDO
+!
+!--    West- and eastward facing surfaces
+       DO  n = 2, 3
+          DO  m = 1, surf_def_v(n)%ns
+             i = surf_def_v(n)%i(m)
+             j = surf_def_v(n)%j(m)
+             k = surf_def_v(n)%k(m)
+             km(k,j,i) = kappa * surf_def_v(n)%us(m) * 0.5_wp * dx
+             kh(k,j,i) = 1.35_wp * km(k,j,i)
+          ENDDO
+       ENDDO
        CALL exchange_horiz( km, nbgp )
 …
     ENDIF
+!
 !-- Model top
 …
     INTEGER(iwp) ::  j      !< loop index y direction
     INTEGER(iwp) ::  k      !< loop index z direction
     INTEGER, INTENT(IN) ::  mod_count  !<
+    INTEGER, INTENT(IN) ::  mod_count  !< flag defining where pointers point to
 #if defined( __nopointer )

palm/trunk/SOURCE/vertical_nesting_mod.f90

-                      r3066
+                      r3083
 !-- Identical timestep for coarse and fine grids
           dt_3d = MIN( dtc, dtf )
+          !> @fixme setting old_dt might be obsolete at this point
+          !>   Due to changes in timestep routine, setting of old_dt might be
+          !>   not necessary any more at this point. However, could not be
+          !>   tested so far.
+          !>   2018-05-18, gronemeier
           old_dt = dt_3d
 #endif

Note: See TracChangeset for help on using the changeset viewer.

Context Navigation

Changeset 3083 for palm

Legend:

Download in other formats: