Index: palm/trunk/SOURCE/bulk_cloud_model_mod.f90
===================================================================
--- palm/trunk/SOURCE/bulk_cloud_model_mod.f90	(revision 4505)
+++ palm/trunk/SOURCE/bulk_cloud_model_mod.f90	(revision 4506)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! Use correct magnus formula for liquid water temperature
+! 
+! 4502 2020-04-17 16:14:16Z schwenkel
 ! Implementation of ice microphysics
 ! 
@@ -134,5 +137,5 @@
 
     USE basic_constants_and_equations_mod,                                     &
-        ONLY:  c_p, g, lv_d_cp, lv_d_rd, l_v, magnus, magnus_ice,              &
+        ONLY:  c_p, g, lv_d_cp, lv_d_rd, l_v, magnus, magnus_ice, magnus_tl,   &
                molecular_weight_of_solute,                                     &
                molecular_weight_of_water, pi, rho_l, rho_s, r_d, r_v, vanthoff,&
@@ -1280,5 +1283,5 @@
              d_exner = exner_function_invers(hyp)
              exner = 1.0_wp / exner_function_invers(hyp)
-             hyrho  = ideal_gas_law_rho_pt(hyp, pt(:, nys, nxl) )
+             hyrho  = ideal_gas_law_rho_pt(hyp, pt_init)
 !
 !--          Compute reference density
@@ -1396,5 +1399,5 @@
              d_exner = exner_function_invers(hyp)
              exner = 1.0_wp / exner_function_invers(hyp)
-             hyrho  = ideal_gas_law_rho_pt(hyp, pt(:, nys, nxl) )
+             hyrho  = ideal_gas_law_rho_pt(hyp, pt_init)
 !
 !--          Compute reference density
@@ -6434,5 +6437,5 @@
 !
 !--    Calculate water vapor saturation pressure
-       e_s = magnus( t_l )
+       e_s = magnus_tl( t_l )
 !
 !--    Computation of saturation mixing ratio:
Index: palm/trunk/TESTS/cases/warm_air_bubble_with_bcm/MONITORING/warm_air_bubble_with_bcm_rc.000
===================================================================
--- palm/trunk/TESTS/cases/warm_air_bubble_with_bcm/MONITORING/warm_air_bubble_with_bcm_rc.000	(revision 4505)
+++ palm/trunk/TESTS/cases/warm_air_bubble_with_bcm/MONITORING/warm_air_bubble_with_bcm_rc.000	(revision 4506)
@@ -1,9 +1,9 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4325M       *    atmosphere - run without 1D - prerun
+ * PALM 6.0  Rev: 4500M       *    atmosphere - run without 1D - prerun
  ******************************    --------------------------------------------
 
- Date:               2019-12-06    Run:       warm_air_bubble_with_bcm__gfortran
- Time:                 14:27:21    Run-No.:   00
+ Date:               2020-04-21    Run:       warm_air_bubble_with_bcm__gfortran
+ Time:                 12:43:37    Run-No.:   00
  Run on host:        gfortran_d
  Number of PEs:               1    Processor grid (x,y): (   1,   1) calculated
@@ -14,5 +14,5 @@
  -----------------
 
- --> Use the Moeng_Wyngaard turbulence closure (LES mode).
+ --> Use the 1.5-order turbulence closure (LES mode).
  --> Use the boussinesq approximation for the model equations.
  --> Solve perturbation pressure via FFT using temperton-algorithm routines
@@ -214,5 +214,5 @@
  ---------------------------
 
-    Mixing length is limited to 1.80 * z
+    Mixing length is limited close to surfaces
 
 
@@ -320,76 +320,75 @@
   0     57 00:08:51.80   9.3300A   0.3737   1.7995   2.9246  0.072  0.00   0.138E-02      0.  0.230E+00  0.230E+00  0.560E-04  0.738E-14   13  16  38    1  15  17    7  20  27    0.000    0.000       0
   0     58 00:09:01.13   9.3300A   0.4133   1.8090   3.0084  0.072  0.00   0.134E-02      0.  0.243E+00  0.243E+00  0.587E-04  0.774E-14   13  16  38    1  15  17    7  20  27    0.000    0.000       0
-  0     59 00:09:10.46   9.3300A   0.4209   1.8154   3.1070  0.073  0.00   0.131E-02      0.  0.257E+00  0.257E+00  0.614E-04  0.810E-14   13  16  38    1  15  17    8  20  27    0.000    0.000       0
-  0     60 00:09:19.79   9.3300A  -0.4192  -1.8278   3.2019  0.073  0.00   0.127E-02      0.  0.271E+00  0.271E+00  0.641E-04  0.846E-14   13  16   3    2  27  21    8  20  26    0.000    0.000       0
-  0     61 00:09:29.12   9.3300A   0.4571  -1.8401   3.2932  0.074  0.00   0.124E-02      0.  0.285E+00  0.285E+00  0.668E-04  0.884E-14   13  15  38    2  27  22    8  20  26    0.000    0.000       0
-  0     62 00:09:38.46   9.3300A   0.4681  -1.8522   3.3731  0.074  0.00   0.121E-02      0.  0.300E+00  0.300E+00  0.696E-04  0.923E-14   13  15  38    2  27  22    8  20  26    0.000    0.000       0
-  0     63 00:09:47.79   9.3300A  -0.4994  -1.8605   3.4422  0.075  0.00   0.119E-02      0.  0.315E+00  0.315E+00  0.725E-04  0.964E-14   17  19  25    2  27  22    8  20  26    0.000    0.000       0
-  0     64 00:09:57.12   9.3300A   0.5292   1.8988   3.5294  0.075  0.00   0.116E-02      0.  0.330E+00  0.330E+00  0.756E-04  0.101E-13   15  19  25   15  27  32    9  20  26    0.000    0.000       0
-  0     65 00:10:06.45   9.3300A   0.5905  -1.9845   3.6050  0.075  0.00   0.114E-02      0.  0.347E+00  0.347E+00  0.789E-04  0.105E-13   15  18  25   16  14   1    9  20  26    0.000    0.000       0
-  0     66 00:10:15.78   9.3300A  -0.6320  -2.0306   3.6694  0.076  0.00   0.111E-02      0.  0.364E+00  0.364E+00  0.822E-04  0.110E-13   18  19  25   16  14   1    9  20  26    0.000    0.000       0
-  0     67 00:10:25.11   9.3300A  -0.7255   2.0707   3.7241  0.076  0.00   0.108E-02      0.  0.381E+00  0.381E+00  0.857E-04  0.115E-13   18  19  25   15  28  32    9  20  26    0.000    0.000       0
-  0     68 00:10:34.44   9.3300A  -0.7805   2.1237   3.7835  0.076  0.00   0.106E-02      0.  0.400E+00  0.400E+00  0.892E-04  0.120E-13   18  18  25   16  27   7   10  20  26    0.000    0.000       0
-  0     69 00:10:43.77   9.3300A   0.8341   2.1570   3.8420  0.076  0.00   0.104E-02      0.  0.418E+00  0.418E+00  0.928E-04  0.125E-13   16  18  25   16  27   7   10  20  26    0.000    0.000       0
-  0     70 00:10:53.10   9.3300A   0.9010   2.1577   3.8966  0.076  0.00   0.102E-02      0.  0.438E+00  0.438E+00  0.965E-04  0.130E-13   16  18  25   16  28  33   10  20  27    0.000    0.000       0
-  0     71 00:11:02.43   9.3300A  -0.9908   2.2162   3.9483  0.076  0.00   0.100E-02      0.  0.458E+00  0.458E+00  0.100E-03  0.136E-13   19  19  25   16  28   7   10  20  27    0.000    0.000       0
-  0     72 00:11:11.76   9.3300A  -1.0775   2.2420   3.9982  0.076  0.00   0.978E-03      0.  0.479E+00  0.479E+00  0.105E-03  0.141E-13   19  18  25   16  28   7   10  20  12    0.000    0.000       0
-  0     73 00:11:21.09   9.3300A  -1.1238   2.2271   4.0479  0.076  0.00   0.957E-03      0.  0.500E+00  0.500E+00  0.109E-03  0.147E-13   19  17  25   16  28   8   10  20  12    0.000    0.000       0
-  0     74 00:11:30.42   9.3300A  -1.1402   2.2306   4.1036  0.076  0.00   0.938E-03      0.  0.521E+00  0.521E+00  0.113E-03  0.153E-13   19  17  25   17  28   7   11  20  12    0.000    0.000       0
-  0     75 00:11:39.75   9.3300A  -1.2685   2.2939   4.1579  0.076  0.00   0.920E-03      0.  0.543E+00  0.543E+00  0.118E-03  0.159E-13   20  19  25   17  28   7   11  20  12    0.000    0.000       0
-  0     76 00:11:49.08   9.3300A  -1.3616   2.2689   4.2046  0.076  0.00   0.903E-03      0.  0.564E+00  0.564E+00  0.123E-03  0.165E-13   20  18  25   17  28   7   11  20  12    0.000    0.000       0
-  0     77 00:11:58.41   9.3300A   1.4149  -2.2146   4.2426  0.075  0.00   0.886E-03      0.  0.586E+00  0.586E+00  0.128E-03  0.172E-13   17  19  24   16  13   9   11  20  12    0.000    0.000       0
-  0     78 00:12:07.74   9.3300A   1.4509   2.2092   4.2785  0.075  0.00   0.870E-03      0.  0.607E+00  0.607E+00  0.133E-03  0.179E-13   17  22  24   18  28   7   12  20  12    0.000    0.000       0
-  0     79 00:12:17.07   9.3300A   1.6077   2.2534   4.3203  0.075  0.00   0.856E-03      0.  0.628E+00  0.628E+00  0.138E-03  0.185E-13   18  20  24   18  28   6   11  20  17    0.000    0.000       0
-  0     80 00:12:26.40   9.3300A   1.7852  -2.2416   4.3825  0.074  0.00   0.844E-03      0.  0.648E+00  0.648E+00  0.144E-03  0.192E-13   18  20  24   19  15   0   18  20  25    0.000    0.000       0
-  0     81 00:12:35.73   9.3300A   1.8811  -2.2962   4.5340  0.074  0.00   0.846E-03      0.  0.667E+00  0.667E+00  0.150E-03  0.199E-13   18  21  24   19  14   0   19  20  25    0.000    0.000       0
-  0     82 00:12:45.06   9.3300A   1.9184  -2.3817   4.7101  0.073  0.00   0.822E-03      0.  0.685E+00  0.685E+00  0.155E-03  0.206E-13   18  21  24   19  14   0   19  20  25    0.000    0.000       0
-  0     83 00:12:54.39   9.2700A   2.0163  -2.3672   4.8527  0.073  0.00   0.804E-03      0.  0.702E+00  0.702E+00  0.160E-03  0.212E-13   18  22  24   19  14   0   19  20  25    0.000    0.000       0
-  0     84 00:13:03.66   9.1000A   2.1154  -2.3875   4.9455  0.072  0.00   0.792E-03      0.  0.717E+00  0.717E+00  0.164E-03  0.216E-13   19  20  24   20  14   0   19  20  25    0.000    0.000       0
-  0     85 00:13:12.76   8.9900A   2.2146  -2.5041   5.0046  0.072  0.00   0.782E-03      0.  0.730E+00  0.730E+00  0.165E-03  0.218E-13   19  21  24   20  14   0   20  20  25    0.000    0.000       0
-  0     86 00:13:21.75   8.8500A   2.2848  -2.5301   5.0819  0.071  0.00   0.773E-03      0.  0.741E+00  0.741E+00  0.167E-03  0.220E-13   19  21  24   20  14   0   20  20  25    0.000    0.000       0
-  0     87 00:13:30.60   8.7700A  -2.2531   2.5764   5.2602  0.070  0.00   0.769E-03      0.  0.750E+00  0.750E+00  0.168E-03  0.222E-13   20  23  30   18  24  21   22  20  12    0.000    0.000       0
-  0     88 00:13:39.37   8.6100A   2.2886  -2.5594   5.4031  0.070  0.00   0.767E-03      0.  0.758E+00  0.757E+00  0.169E-03  0.224E-13   20  20  25   20  13   0   22  20  12    0.000    0.000       0
-  0     89 00:13:47.98   8.5300A   2.4633  -2.5616   5.4541  0.069  0.00   0.781E-03      0.  0.763E+00  0.763E+00  0.168E-03  0.223E-13   20  19  25   20  13   0   21  20  26    0.000    0.000       0
-  0     90 00:13:56.51   8.5200A   2.5249  -2.5919   5.6380  0.068  0.00   0.811E-03      0.  0.766E+00  0.766E+00  0.169E-03  0.224E-13   20  19  25   21  14   0   21  20  26    0.000    0.000       0
-  0     91 00:14:05.03   8.4200A   2.5264   2.6137   5.7696  0.068  0.00   0.864E-03      0.  0.767E+00  0.767E+00  0.171E-03  0.227E-13   20  21  25   19  24  21   21  20  26    0.000    0.000       0
-  0     92 00:14:13.45   8.3400A   2.5649   2.6404   5.8445  0.067  0.00   0.900E-03      0.  0.767E+00  0.767E+00  0.170E-03  0.227E-13   20  22  25   19  24  21   22  20  26    0.000    0.000       0
-  0     93 00:14:21.79   8.3600A   2.6069  -2.6252   6.0214  0.066  0.00   0.922E-03      0.  0.764E+00  0.764E+00  0.170E-03  0.227E-13   20  22  25   21  13   0   22  20  26    0.000    0.000       0
-  0     94 00:14:30.15   8.4400A   2.6402  -2.6226   6.0997  0.065  0.00   0.934E-03      0.  0.760E+00  0.760E+00  0.172E-03  0.230E-13   21  21  25   21  13   0   22  20  26    0.000    0.000       0
-  0     95 00:14:38.59   8.3900A   2.6507  -2.6679   6.1303  0.064  0.00   0.945E-03      0.  0.754E+00  0.754E+00  0.175E-03  0.234E-13   21  21  25   22  14   0   23  20  26    0.000    0.000       0
-  0     96 00:14:46.98   8.4000A  -2.6422  -2.6976   6.1877  0.063  0.00   0.957E-03      0.  0.747E+00  0.746E+00  0.175E-03  0.234E-13   20  19  22   22  14   0   23  20  26    0.000    0.000       0
-  0     97 00:14:55.38   8.4400A   2.6197  -2.6735   6.2025  0.063  0.00   0.972E-03      0.  0.737E+00  0.737E+00  0.176E-03  0.236E-13   21  22  25   22  15   3   23  20  27    0.000    0.000       0
-  0     98 00:15:03.82   8.3700A  -2.6472  -2.7972   6.2915  0.062  0.00   0.987E-03      0.  0.726E+00  0.726E+00  0.178E-03  0.238E-13   21  20  22   22  15   3   24  20  26    0.000    0.000       0
-  0     99 00:15:12.19   8.3800A  -2.7125  -2.8119   6.2200  0.061  0.00   0.987E-03      0.  0.714E+00  0.714E+00  0.177E-03  0.238E-13   24  20  29   22  15   3   24  20  26    0.000    0.000       0
-  0    100 00:15:20.57   8.4200A  -2.8380  -2.8368   6.2806  0.060  0.00   0.975E-03      0.  0.701E+00  0.700E+00  0.178E-03  0.240E-13   24  21  28   23  16   3   25  20  26    0.000    0.000       0
-  0    101 00:15:28.99   8.4400A  -2.9251  -2.8118   6.3188  0.059  0.00   0.964E-03      0.  0.686E+00  0.686E+00  0.178E-03  0.241E-13   24  22  28   23  15   3   25  21  26    0.000    0.000       0
-  0    102 00:15:37.43   8.3500A  -2.9653  -2.9492   6.1726  0.058  0.00   0.962E-03      0.  0.671E+00  0.670E+00  0.179E-03  0.242E-13   24  22  28   23  15   3   25  21  26    0.000    0.000       0
-  0    103 00:15:45.78   8.3200A  -2.8269  -2.9602   5.9671  0.057  0.00   0.969E-03      0.  0.654E+00  0.654E+00  0.176E-03  0.240E-13   24  22  28   23  15   3   26  20  26    0.000    0.000       0
-  0    104 00:15:54.10   8.3100A  -2.9048  -2.9405   5.9882  0.057  0.00   0.101E-02      0.  0.637E+00  0.637E+00  0.176E-03  0.239E-13   24  23  28   24  16   3   23  20  19    0.000    0.000       0
-  0    105 00:16:02.41   8.2800A  -2.8662  -2.9873   6.0499  0.056  0.00   0.126E-02      0.  0.620E+00  0.620E+00  0.176E-03  0.238E-13   24  23  28   25  16  11   23  20  19    0.000    0.000       0
-  0    106 00:16:10.69   8.2500A  -2.8878  -2.9981   6.1340  0.055  0.00   0.970E-03      0.  0.603E+00  0.602E+00  0.175E-03  0.237E-13   25  22  27   25  16  11   24  20  19    0.000    0.000       0
-  0    107 00:16:18.94   8.2600A  -2.8912   2.9599   6.1957  0.054  0.00   0.931E-03      0.  0.585E+00  0.585E+00  0.174E-03  0.236E-13   25  22  27   25  26  11   24  20  19    0.000    0.000       0
-  0    108 00:16:27.20   8.3000A  -2.8493  -2.8413   6.1933  0.054  0.00   0.910E-03      0.  0.568E+00  0.567E+00  0.175E-03  0.235E-13   25  22  27   24  15   3   24  20  19    0.000    0.000       0
-  0    109 00:16:35.50   8.3700A  -2.8947  -2.8726   6.1914  0.053  0.00   0.896E-03      0.  0.551E+00  0.550E+00  0.176E-03  0.236E-13   25  23  27   25  15  11   25  20  19    0.000    0.000       0
-  0    110 00:16:43.87   8.4500A  -2.9525  -2.8298   6.2665  0.052  0.00   0.887E-03      0.  0.533E+00  0.533E+00  0.177E-03  0.238E-13   25  23  27   25  15  11   25  20  19    0.000    0.000       0
-  0    111 00:16:52.32   8.5600A  -2.9801  -2.7141   6.2155  0.052  0.00   0.869E-03      0.  0.517E+00  0.516E+00  0.179E-03  0.240E-13   25  23  27   25  15  11   25  20  19    0.000    0.000       0
-  0    112 00:17:00.88   8.7000A  -2.9755   2.6943   6.0730  0.051  0.00   0.855E-03      0.  0.500E+00  0.500E+00  0.181E-03  0.242E-13   25  23  27   26  26  11   25  20  19    0.000    0.000       0
-  0    113 00:17:09.58   8.8400A  -2.9327   2.7209   6.1317  0.050  0.00   0.845E-03      0.  0.484E+00  0.484E+00  0.184E-03  0.247E-13   25  23  27   26  26  11   26  20  19    0.000    0.000       0
-  0    114 00:17:18.42   9.0100A  -2.8456   2.6746   6.0544  0.050  0.00   0.838E-03      0.  0.469E+00  0.469E+00  0.187E-03  0.251E-13   25  23  27   26  26  11   26  20  19    0.000    0.000       0
-  0    115 00:17:27.43   9.1600A  -2.8165  -2.7118   6.0872  0.049  0.00   0.832E-03      0.  0.455E+00  0.454E+00  0.191E-03  0.257E-13   25  24  26   25  16  25   26  20   6    0.000    0.000       0
-  0    116 00:17:36.59   9.3200A  -2.8770  -2.7160   6.1728  0.049  0.00   0.828E-03      0.  0.442E+00  0.441E+00  0.195E-03  0.262E-13   25  24  26   25  16  25   26  20   6    0.000    0.000       0
-  0    117 00:17:45.91   9.3300A  -2.9113  -2.6767   6.2011  0.048  0.00   0.823E-03      0.  0.430E+00  0.429E+00  0.199E-03  0.267E-13   25  24  26   25  16  25   26  20   6    0.000    0.000       0
-  0    118 00:17:55.24   9.3300A  -2.9222  -2.6323   6.1859  0.048  0.00   0.821E-03      0.  0.419E+00  0.418E+00  0.200E-03  0.269E-13   25  24  26   25  16  25   26  20   6    0.000    0.000       0
-  0    119 00:18:04.57   9.3300A  -2.9148  -2.5979   6.4137  0.047  0.00   0.818E-03      0.  0.410E+00  0.409E+00  0.201E-03  0.271E-13   25  24  26   25  16  25   27  20   6    0.000    0.000       0
-  0    120 00:18:13.90   9.3300A  -2.9239  -2.5688   6.5911  0.047  0.00   0.810E-03      0.  0.402E+00  0.402E+00  0.202E-03  0.272E-13   26  23  26   25  16  25   27  20   6    0.000    0.000       0
-  0    121 00:18:23.23   9.3300A  -2.9488   2.6022   6.5678  0.047  0.00   0.791E-03      0.  0.396E+00  0.395E+00  0.204E-03  0.274E-13   26  23  26   26  24  31   27  20   6    0.000    0.000       0
-  0    122 00:18:32.56   9.3300A  -2.9594   2.5707   6.3816  0.046  0.00   0.761E-03      0.  0.391E+00  0.391E+00  0.205E-03  0.277E-13   26  23  26   26  24  31   27  20   5    0.000    0.000       0
-  0    123 00:18:41.89   9.3300D  -2.9557   2.4413   6.3639  0.046  0.00   0.728E-03      0.  0.388E+00  0.387E+00  0.207E-03  0.279E-13   26  23  26   26  24  31   27  20   5    0.000    0.000       0
-  0    124 00:18:51.22   9.3300D  -3.0047  -2.3814   6.2412  0.046  0.00   0.707E-03      0.  0.386E+00  0.386E+00  0.209E-03  0.282E-13   27  19   9   26  19   7   27  20   5    0.000    0.000       0
-  0    125 00:19:00.55   9.3300D  -3.0558  -2.4018   6.2470  0.046  0.00   0.691E-03      0.  0.386E+00  0.385E+00  0.211E-03  0.285E-13   27  21   7   29  17   3   27  20   3    0.000    0.000       0
-  0    126 00:19:09.88   9.3300D  -3.1229  -2.4600   6.1559  0.045  0.00   0.689E-03      0.  0.387E+00  0.386E+00  0.213E-03  0.288E-13   27  21   7   29  17   3   27  20   3    0.000    0.000       0
-  0    127 00:19:19.21   9.3300D  -3.0820  -2.4179   6.2289  0.045  0.00   0.654E-03      0.  0.389E+00  0.388E+00  0.215E-03  0.291E-13   27  21   7   29  17   3   28  20  37    0.000    0.000       0
-  0    128 00:19:28.53   9.3300D  -2.9813   2.4661   6.2373  0.045  0.00   0.637E-03      0.  0.391E+00  0.390E+00  0.218E-03  0.293E-13   27  21   7   23  19  19   28  20  37    0.000    0.000       0
-  0    129 00:19:37.86   9.3300D  -3.0484   2.5239   6.2511  0.045  0.00   0.626E-03      0.  0.395E+00  0.394E+00  0.220E-03  0.294E-13   28  20   7   23  19  19   28  20   3    0.000    0.000       0
-  0    130 00:19:47.19   9.3300D  -3.2370   2.5417   6.2843  0.044  0.00   0.599E-03      0.  0.399E+00  0.398E+00  0.222E-03  0.296E-13   28  20   7   23  19  19   28  20   3    0.000    0.000       0
-  0    131 00:19:56.52   9.3300D  -3.3131   2.5235   6.1476  0.044  0.00   0.580E-03      0.  0.404E+00  0.403E+00  0.224E-03  0.298E-13   28  20   7   23  19  19   28  20   3    0.000    0.000       0
-  0    132 00:20:05.85   9.3300D  -3.3210   2.5312   5.9087  0.044  0.00   0.563E-03      0.  0.409E+00  0.408E+00  0.225E-03  0.299E-13   28  20   7   22  19   9   28  20   3    0.000    0.000       0
+  0     59 00:09:10.46   9.3300A   0.4208   1.8154   3.1070  0.073  0.00   0.131E-02      0.  0.257E+00  0.257E+00  0.614E-04  0.810E-14   13  16  38    1  15  17    8  20  27    0.000    0.000       0
+  0     60 00:09:19.79   9.3300A  -0.4191  -1.8278   3.2018  0.073  0.00   0.127E-02      0.  0.271E+00  0.271E+00  0.641E-04  0.846E-14   13  16   3    2  27  21    8  20  26    0.000    0.000       0
+  0     61 00:09:29.12   9.3300A   0.4569  -1.8400   3.2930  0.074  0.00   0.124E-02      0.  0.285E+00  0.285E+00  0.668E-04  0.884E-14   13  15  38    2  27  22    8  20  26    0.000    0.000       0
+  0     62 00:09:38.46   9.3300A   0.4677  -1.8521   3.3726  0.074  0.00   0.121E-02      0.  0.299E+00  0.299E+00  0.695E-04  0.922E-14   13  15  38    2  27  22    8  20  26    0.000    0.000       0
+  0     63 00:09:47.79   9.3300A  -0.4917  -1.8603   3.4413  0.075  0.00   0.119E-02      0.  0.315E+00  0.314E+00  0.724E-04  0.962E-14   17  19  25    2  27  22    8  20  26    0.000    0.000       0
+  0     64 00:09:57.12   9.3300A   0.5192   1.8995   3.5274  0.075  0.00   0.116E-02      0.  0.330E+00  0.330E+00  0.754E-04  0.100E-13   15  19  25   15  27  32    9  20  26    0.000    0.000       0
+  0     65 00:10:06.45   9.3300A   0.5791  -1.9847   3.6024  0.075  0.00   0.114E-02      0.  0.346E+00  0.346E+00  0.787E-04  0.105E-13   15  18  25   16  14   1    9  20  26    0.000    0.000       0
+  0     66 00:10:15.78   9.3300A  -0.6157  -2.0324   3.6660  0.076  0.00   0.111E-02      0.  0.363E+00  0.363E+00  0.819E-04  0.110E-13   18  19  25   16  14   1    9  20  26    0.000    0.000       0
+  0     67 00:10:25.11   9.3300A  -0.7098   2.0731   3.7198  0.076  0.00   0.108E-02      0.  0.380E+00  0.380E+00  0.853E-04  0.114E-13   18  19  25   15  28  32    9  20  26    0.000    0.000       0
+  0     68 00:10:34.44   9.3300A  -0.7664   2.1213   3.7761  0.076  0.00   0.106E-02      0.  0.398E+00  0.398E+00  0.887E-04  0.119E-13   18  19  25   16  27   7   10  20  26    0.000    0.000       0
+  0     69 00:10:43.77   9.3300A   0.8138   2.1545   3.8335  0.076  0.00   0.104E-02      0.  0.416E+00  0.416E+00  0.922E-04  0.124E-13   16  18  25   16  27   7   10  20  26    0.000    0.000       0
+  0     70 00:10:53.10   9.3300A   0.8831   2.1560   3.8868  0.076  0.00   0.102E-02      0.  0.435E+00  0.435E+00  0.959E-04  0.129E-13   16  18  25   16  28  33   10  20  27    0.000    0.000       0
+  0     71 00:11:02.43   9.3300A  -0.9638   2.2185   3.9373  0.076  0.00   0.100E-02      0.  0.455E+00  0.455E+00  0.997E-04  0.135E-13   19  19  25   16  28   7   10  20  27    0.000    0.000       0
+  0     72 00:11:11.76   9.3300A  -1.0542   2.2455   3.9858  0.076  0.00   0.978E-03      0.  0.475E+00  0.475E+00  0.104E-03  0.140E-13   19  18  25   16  28   7   10  20  12    0.000    0.000       0
+  0     73 00:11:21.09   9.3300A  -1.1024   2.2300   4.0344  0.076  0.00   0.957E-03      0.  0.495E+00  0.495E+00  0.108E-03  0.146E-13   19  18  25   16  28   8   10  20  12    0.000    0.000       0
+  0     74 00:11:30.42   9.3300A  -1.1378   2.2165   4.0847  0.076  0.00   0.938E-03      0.  0.516E+00  0.516E+00  0.112E-03  0.151E-13   19  17  25   17  28   7   11  20  12    0.000    0.000       0
+  0     75 00:11:39.75   9.3300A  -1.2291   2.2852   4.1378  0.076  0.00   0.920E-03      0.  0.537E+00  0.536E+00  0.116E-03  0.157E-13   20  19  25   17  28   7   11  20  12    0.000    0.000       0
+  0     76 00:11:49.08   9.3300A  -1.3308   2.2688   4.1835  0.076  0.00   0.903E-03      0.  0.557E+00  0.557E+00  0.121E-03  0.163E-13   20  18  25   17  28   7   11  20  12    0.000    0.000       0
+  0     77 00:11:58.41   9.3300A  -1.3838   2.2159   4.2204  0.075  0.00   0.887E-03      0.  0.578E+00  0.578E+00  0.126E-03  0.170E-13   20  18  25   16  28  10   11  20  12    0.000    0.000       0
+  0     78 00:12:07.74   9.3300A   1.4164   2.2009   4.2484  0.075  0.00   0.871E-03      0.  0.598E+00  0.598E+00  0.131E-03  0.176E-13   17  22  24   16  28  10   11  20  12    0.000    0.000       0
+  0     79 00:12:17.07   9.3300A   1.5394   2.2426   4.2951  0.075  0.00   0.857E-03      0.  0.618E+00  0.618E+00  0.136E-03  0.183E-13   18  20  24   18  28   6   11  20  17    0.000    0.000       0
+  0     80 00:12:26.40   9.3300A   1.7339  -2.2246   4.3378  0.074  0.00   0.845E-03      0.  0.637E+00  0.637E+00  0.141E-03  0.189E-13   18  20  24   19  15   0   11  20  17    0.000    0.000       0
+  0     81 00:12:35.73   9.3300A   1.8350  -2.2516   4.4188  0.074  0.00   0.851E-03      0.  0.656E+00  0.656E+00  0.147E-03  0.196E-13   18  21  24   19  14   0   19  20  25    0.000    0.000       0
+  0     82 00:12:45.06   9.3300A   1.8921  -2.3592   4.5957  0.073  0.00   0.821E-03      0.  0.673E+00  0.673E+00  0.152E-03  0.202E-13   18  21  24   19  14   0   19  20  25    0.000    0.000       0
+  0     83 00:12:54.39   9.3300A   1.9746  -2.3856   4.7437  0.073  0.00   0.805E-03      0.  0.689E+00  0.689E+00  0.157E-03  0.208E-13   18  22  24   19  14   0   19  20  25    0.000    0.000       0
+  0     84 00:13:03.72   9.2800A   2.0514  -2.3567   4.8498  0.072  0.00   0.792E-03      0.  0.703E+00  0.703E+00  0.162E-03  0.214E-13   19  20  24   20  15   0   19  20  25    0.000    0.000       0
+  0     85 00:13:13.00   9.1900A   2.1654  -2.4441   4.8983  0.072  0.00   0.782E-03      0.  0.716E+00  0.716E+00  0.165E-03  0.218E-13   19  21  24   20  14   0   20  20  25    0.000    0.000       0
+  0     86 00:13:22.19   9.0100A   2.2631  -2.5161   4.9934  0.071  0.00   0.773E-03      0.  0.727E+00  0.726E+00  0.167E-03  0.221E-13   19  21  24   20  14   0   20  20  25    0.000    0.000       0
+  0     87 00:13:31.20   8.9600A   2.2148   2.5406   5.1131  0.070  0.00   0.769E-03      0.  0.735E+00  0.735E+00  0.167E-03  0.221E-13   19  21  24   18  24  21   22  20  12    0.000    0.000       0
+  0     88 00:13:40.16   8.7900A   2.2426   2.5424   5.3017  0.070  0.00   0.767E-03      0.  0.742E+00  0.742E+00  0.169E-03  0.223E-13   19  22  25   18  24  21   22  20  12    0.000    0.000       0
+  0     89 00:13:48.95   8.6800A   2.3976  -2.5531   5.3515  0.069  0.00   0.783E-03      0.  0.747E+00  0.747E+00  0.168E-03  0.223E-13   20  19  25   20  13   0   22  20  12    0.000    0.000       0
+  0     90 00:13:57.63   8.6500A   2.4819  -2.5387   5.4991  0.068  0.00   0.817E-03      0.  0.750E+00  0.749E+00  0.168E-03  0.223E-13   20  19  25   21  14   0   21  20  26    0.000    0.000       0
+  0     91 00:14:06.28   8.6000A   2.5052   2.5688   5.6500  0.068  0.00   0.868E-03      0.  0.750E+00  0.750E+00  0.170E-03  0.225E-13   20  21  25   19  24  21   21  20  26    0.000    0.000       0
+  0     92 00:14:14.88   8.4800A   2.5443   2.6361   5.7333  0.067  0.00   0.903E-03      0.  0.749E+00  0.749E+00  0.170E-03  0.226E-13   20  22  25   19  24  21   21  20  26    0.000    0.000       0
+  0     93 00:14:23.36   8.4700A   2.5995  -2.5705   5.8872  0.066  0.00   0.923E-03      0.  0.747E+00  0.746E+00  0.169E-03  0.225E-13   20  22  25    7  28  16   22  20  26    0.000    0.000       0
+  0     94 00:14:31.83   8.5800A   2.6140  -2.5933   5.9981  0.065  0.00   0.935E-03      0.  0.742E+00  0.742E+00  0.171E-03  0.226E-13   21  21  25   21  13   0   22  20  26    0.000    0.000       0
+  0     95 00:14:40.41   8.5400A   2.6437  -2.5854   5.9981  0.064  0.00   0.947E-03      0.  0.736E+00  0.735E+00  0.174E-03  0.231E-13   21  21  25   21  15   3   22  20  26    0.000    0.000       0
+  0     96 00:14:48.95   8.5300A   2.6117  -2.6340   6.0887  0.064  0.00   0.961E-03      0.  0.728E+00  0.727E+00  0.174E-03  0.232E-13   21  21  25   22  14   0   23  20  26    0.000    0.000       0
+  0     97 00:14:57.48   8.5700A   2.5977  -2.6476   6.1143  0.063  0.00   0.978E-03      0.  0.718E+00  0.718E+00  0.175E-03  0.233E-13   21  22  25   22  14   0   23  20  27    0.000    0.000       0
+  0     98 00:15:06.05   8.5100A   2.6203  -2.7146   6.1741  0.062  0.00   0.985E-03      0.  0.707E+00  0.707E+00  0.177E-03  0.235E-13   21  22  25   22  15   3   24  20  26    0.000    0.000       0
+  0     99 00:15:14.56   8.5000A  -2.6579  -2.7759   6.1592  0.061  0.00   0.974E-03      0.  0.695E+00  0.694E+00  0.176E-03  0.235E-13   21  20  22   22  15   3   24  20  26    0.000    0.000       0
+  0    100 00:15:23.06   8.5400A  -2.7358  -2.7779   6.1081  0.060  0.00   0.963E-03      0.  0.681E+00  0.681E+00  0.176E-03  0.236E-13   24  21  28   23  16   3   25  20  26    0.000    0.000       0
+  0    101 00:15:31.60   8.5600A  -2.8233  -2.7825   6.1950  0.059  0.00   0.956E-03      0.  0.666E+00  0.666E+00  0.177E-03  0.238E-13   21  19  22   23  16   3   25  20  26    0.000    0.000       0
+  0    102 00:15:40.16   8.4700A  -2.9314   2.8503   6.1310  0.059  0.00   0.959E-03      0.  0.651E+00  0.650E+00  0.177E-03  0.239E-13   24  22  28   24  26  11   25  21  26    0.000    0.000       0
+  0    103 00:15:48.63   8.4300A  -2.8678  -2.8992   5.9199  0.058  0.00   0.974E-03      0.  0.634E+00  0.634E+00  0.175E-03  0.236E-13   24  22  28   23  15   3   25  19  26    0.000    0.000       0
+  0    104 00:15:57.06   8.4200A  -2.8272   2.8728   5.9107  0.057  0.00   0.104E-02      0.  0.618E+00  0.617E+00  0.174E-03  0.235E-13   24  23  28   24  26  12   23  20  19    0.000    0.000       0
+  0    105 00:16:05.48   8.4200A  -2.8546  -2.8737   5.9813  0.056  0.00   0.100E-02      0.  0.600E+00  0.600E+00  0.174E-03  0.235E-13   24  23  28   24  16   3   23  20  19    0.000    0.000       0
+  0    106 00:16:13.90   8.3800A  -2.7891  -2.9537   6.0319  0.056  0.00   0.941E-03      0.  0.583E+00  0.583E+00  0.174E-03  0.235E-13   25  22  27   25  16  11   23  20  19    0.000    0.000       0
+  0    107 00:16:22.28   8.3800A  -2.8342  -2.9148   6.1201  0.055  0.00   0.916E-03      0.  0.566E+00  0.565E+00  0.173E-03  0.233E-13   25  22  27   25  16  11   24  20  19    0.000    0.000       0
+  0    108 00:16:30.66   8.4200A  -2.8189   2.8537   6.1339  0.054  0.00   0.899E-03      0.  0.549E+00  0.548E+00  0.173E-03  0.233E-13   25  22  27   24  27  12   24  20  19    0.000    0.000       0
+  0    109 00:16:39.08   8.4900A  -2.8294   2.8399   6.0967  0.053  0.00   0.886E-03      0.  0.531E+00  0.531E+00  0.174E-03  0.234E-13   23  19  21   25  26  12   24  20  19    0.000    0.000       0
+  0    110 00:16:47.57   8.5700A  -2.9392   2.8297   6.1658  0.053  0.00   0.921E-03      0.  0.515E+00  0.514E+00  0.176E-03  0.236E-13   23  19  21   25  26  12   25  20  19    0.000    0.000       0
+  0    111 00:16:56.14   8.6900A  -2.9403   2.7591   6.1352  0.052  0.00   0.861E-03      0.  0.498E+00  0.498E+00  0.177E-03  0.238E-13   23  19  21   24  24  26   25  20  19    0.000    0.000       0
+  0    112 00:17:04.83   8.8200A  -2.9020   2.7213   5.9950  0.052  0.00   0.847E-03      0.  0.482E+00  0.482E+00  0.180E-03  0.241E-13   25  23  27   24  24  26   25  20  19    0.000    0.000       0
+  0    113 00:17:13.65   8.9700A  -2.8897   2.6260   5.9909  0.051  0.00   0.838E-03      0.  0.467E+00  0.467E+00  0.183E-03  0.245E-13   25  23  27   24  24  26   26  20  19    0.000    0.000       0
+  0    114 00:17:22.62   9.1500A  -2.8396   2.6918   5.9467  0.051  0.00   0.832E-03      0.  0.452E+00  0.452E+00  0.187E-03  0.250E-13   25  23  27   26  24  25   26  20  18    0.000    0.000       0
+  0    115 00:17:31.77   9.3200A  -2.8834  -2.6582   5.9365  0.050  0.00   0.829E-03      0.  0.439E+00  0.438E+00  0.191E-03  0.256E-13   24  19  20   25  16  25   26  20   6    0.000    0.000       0
+  0    116 00:17:41.09   9.3300A  -2.9117  -2.7174   6.0594  0.049  0.00   0.826E-03      0.  0.426E+00  0.425E+00  0.195E-03  0.262E-13   24  19  20   25  16  25   26  20   6    0.000    0.000       0
+  0    117 00:17:50.42   9.3300A  -2.8642  -2.7002   6.1010  0.049  0.00   0.818E-03      0.  0.415E+00  0.414E+00  0.197E-03  0.264E-13   24  19  20   25  16  25   26  20   6    0.000    0.000       0
+  0    118 00:17:59.75   9.3300A  -2.8588  -2.6500   6.0935  0.049  0.00   0.810E-03      0.  0.405E+00  0.404E+00  0.197E-03  0.265E-13   25  24  26   25  16  25   26  20   6    0.000    0.000       0
+  0    119 00:18:09.08   9.3300A  -2.8639  -2.5954   6.1245  0.048  0.00   0.800E-03      0.  0.396E+00  0.396E+00  0.198E-03  0.267E-13   25  24  26   25  16  25   27  20   6    0.000    0.000       0
+  0    120 00:18:18.41   9.3300A  -2.8476  -2.5419   6.3824  0.048  0.00   0.779E-03      0.  0.390E+00  0.389E+00  0.199E-03  0.269E-13   25  24  26   25  16  25   27  20   6    0.000    0.000       0
+  0    121 00:18:27.74   9.3300A  -2.8630   2.5542   6.4496  0.047  0.00   0.759E-03      0.  0.384E+00  0.384E+00  0.201E-03  0.271E-13   26  23  26   26  24  31   27  20   6    0.000    0.000       0
+  0    122 00:18:37.07   9.3300A  -2.8870   2.6105   6.3005  0.047  0.00   0.739E-03      0.  0.380E+00  0.380E+00  0.202E-03  0.274E-13   26  23  26   26  24  31   27  20   6    0.000    0.000       0
+  0    123 00:18:46.40   9.3300A  -2.8942   2.5677   6.2484  0.047  0.00   0.707E-03      0.  0.378E+00  0.377E+00  0.204E-03  0.276E-13   26  23  26   26  24  31   27  20   5    0.000    0.000       0
+  0    124 00:18:55.73   9.3300A  -2.8859   2.4580   6.2096  0.047  0.00   0.685E-03      0.  0.376E+00  0.376E+00  0.206E-03  0.279E-13   26  23  26   26  24  31   27  20   5    0.000    0.000       0
+  0    125 00:19:05.06   9.3300A  -2.9881  -2.3346   6.0990  0.046  0.00   0.668E-03      0.  0.376E+00  0.376E+00  0.208E-03  0.282E-13   27  19   9   26  19   7   27  20   3    0.000    0.000       0
+  0    126 00:19:14.39   9.3300D  -3.0709   2.3826   6.1200  0.046  0.00   0.654E-03      0.  0.378E+00  0.377E+00  0.210E-03  0.284E-13   27  21   7   22  19  19   27  20   3    0.000    0.000       0
+  0    127 00:19:23.72   9.3300D  -3.1153   2.4675   5.9803  0.046  0.00   0.639E-03      0.  0.380E+00  0.379E+00  0.213E-03  0.287E-13   27  21   7   22  19  19   27  20   3    0.000    0.000       0
+  0    128 00:19:33.04   9.3300D  -3.0710   2.4876   6.0516  0.046  0.00   0.626E-03      0.  0.383E+00  0.382E+00  0.215E-03  0.289E-13   27  21   7   22  19  19   28  20  37    0.000    0.000       0
+  0    129 00:19:42.37   9.3300D  -3.0103   2.4561   6.0153  0.045  0.00   0.626E-03      0.  0.386E+00  0.385E+00  0.217E-03  0.290E-13   26  23  25   22  19  19   28  20  37    0.000    0.000       0
+  0    130 00:19:51.70   9.3300D  -3.0660  -2.4443   6.0937  0.045  0.00   0.601E-03      0.  0.391E+00  0.390E+00  0.219E-03  0.292E-13   26  23  25   26  17  33   28  20   3    0.000    0.000       0
+  0    131 00:20:01.03   9.3300D  -3.1203   2.5171   6.1200  0.045  0.00   0.586E-03      0.  0.395E+00  0.394E+00  0.221E-03  0.293E-13   28  20   7   22  19   9   28  20   3    0.000    0.000       0
