Index: /palm/trunk/SOURCE/poismg_mod.f90
===================================================================
--- /palm/trunk/SOURCE/poismg_mod.f90	(revision 4831)
+++ /palm/trunk/SOURCE/poismg_mod.f90	(revision 4832)
@@ -25,4 +25,8 @@
 ! -----------------
 ! $Id$
+! bugfix in redblack algorithm: lower i,j indices need to start alternatively with even or odd
+! value on the coarsest grid level, if the subdomain has an uneven number of gridpoints along x/y
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! File re-formatted to follow the PALM coding standard
 !
@@ -742,7 +746,11 @@
     INTEGER(iwp) ::  l      !< grid level
     INTEGER(iwp) ::  n      !< loop variable GauÃ-Seidel iterations
-
-
-    LOGICAL ::  unroll  !< flag indicating whether loop unrolling is possible
+    INTEGER(iwp) ::  save_nxl_mg  !< to save nxl_mg on coarsest level 1
+    INTEGER(iwp) ::  save_nys_mg  !< to save nys_mg on coarsest level 1
+
+
+    LOGICAL ::  adjust_lower_i_index  !< adjust lower limit of i loop in case of odd number of grid points
+    LOGICAL ::  adjust_lower_j_index  !< adjust lower limit of j loop in case of odd number of grid points
+    LOGICAL ::  unroll                !< flag indicating whether loop unrolling is possible
 
     REAL(wp), DIMENSION(nzb:nzt_mg(grid_level)+1,nys_mg(grid_level)-1:nyn_mg(grid_level)+1,        &
@@ -751,8 +759,27 @@
                         nxl_mg(grid_level)-1:nxr_mg(grid_level)+1) ::  p_mg  !< perturbation pressure
 
+
     l = grid_level
 
     unroll = ( MOD( nyn_mg(l)-nys_mg(l)+1, 4 ) == 0  .AND.  MOD( nxr_mg(l)-nxl_mg(l)+1, 2 ) == 0 )
 
+!
+!-- The red/black decomposition requires that on the lower i,j indices need to start alternatively with an
+!-- even or odd value on the coarsest grid level, depending on the core-id, and if the subdomain has an
+!-- uneven number of gridpoints along x/y. Set the respective steering switches here.
+    IF ( l == 1  .AND.  MOD( myidx, 2 ) /= 0  .AND.  MOD( nxl_mg(l) - nxr_mg(l), 2 ) == 0 )  THEN
+       adjust_lower_i_index = .TRUE.
+       save_nxl_mg = nxl_mg(1)
+    ELSE
+       adjust_lower_i_index = .FALSE.
+    ENDIF
+    IF ( l == 1  .AND.  MOD( myidy, 2 ) /= 0  .AND.  MOD( nyn_mg(l) - nys_mg(l), 2 ) == 0 )  THEN
+       adjust_lower_j_index = .TRUE.
+       save_nys_mg = nys_mg(l)
+    ELSE
+       adjust_lower_j_index = .FALSE.
+    ENDIF
+
+
     DO  n = 1, ngsrb
 
@@ -762,4 +789,16 @@
 
              CALL cpu_log( log_point_s(36), 'redblack_no_unroll_f', 'start' )
+
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg + 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg - 1
+                ELSE
+                  nys_mg(l) = save_nys_mg + 1
+                ENDIF
+             ENDIF
 !
 !--          Without unrolling of loops, no cache optimization
@@ -784,4 +823,16 @@
              ENDDO
 
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg - 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg + 1
+                ELSE
+                  nys_mg(l) = save_nys_mg - 1
+                ENDIF
+             ENDIF
+
              !$OMP DO
              DO  i = nxl_mg(l)+1, nxr_mg(l), 2
@@ -803,4 +854,16 @@
              ENDDO
 
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg + 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg + 1
+                ELSE
+                  nys_mg(l) = save_nys_mg - 1
+                ENDIF
+             ENDIF
+
              !$OMP DO
              DO  i = nxl_mg(l), nxr_mg(l), 2
@@ -821,4 +884,16 @@
                 ENDDO
              ENDDO
+
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg - 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg - 1
+                ELSE
+                  nys_mg(l) = save_nys_mg + 1
+                ENDIF
+             ENDIF
 
              !$OMP DO
@@ -1013,4 +1088,15 @@
     ENDDO
 
+!
+!-- Reset lower index limits to their standard values (may happen on coarsest levels only)
+    IF ( adjust_lower_i_index )  THEN
+       nxl_mg(l) = save_nxl_mg
+    ENDIF
+
+    IF ( adjust_lower_j_index )  THEN
+       nys_mg(l) = save_nys_mg
+    ENDIF
+
+
  END SUBROUTINE redblack
 
Index: /palm/trunk/SOURCE/poismg_noopt_mod.f90
===================================================================
--- /palm/trunk/SOURCE/poismg_noopt_mod.f90	(revision 4831)
+++ /palm/trunk/SOURCE/poismg_noopt_mod.f90	(revision 4832)
@@ -25,6 +25,9 @@
 ! -----------------
 ! $Id$
+! bugfix in redblack algorithm: lower i,j indices need to start alternatively with even or odd
+! value on the coarsest grid level, if the subdomain has an uneven number of gridpoints along x/y
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! File re-formatted to follow the PALM coding standard
-!
 !
 ! 4457 2020-03-11 14:20:43Z raasch
@@ -801,14 +804,18 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) :: color  !<
-    INTEGER(iwp) :: i      !<
-    INTEGER(iwp) :: ic     !<
-    INTEGER(iwp) :: j      !<
-    INTEGER(iwp) :: jc     !<
-    INTEGER(iwp) :: jj     !<
-    INTEGER(iwp) :: k      !<
-    INTEGER(iwp) :: l      !<
-    INTEGER(iwp) :: n      !<
-
+    INTEGER(iwp) ::  color  !<
+    INTEGER(iwp) ::  i      !<
+    INTEGER(iwp) ::  ic     !<
+    INTEGER(iwp) ::  j      !<
+    INTEGER(iwp) ::  jc     !<
+    INTEGER(iwp) ::  jj     !<
+    INTEGER(iwp) ::  k      !<
+    INTEGER(iwp) ::  l      !<
+    INTEGER(iwp) ::  n      !<
+    INTEGER(iwp) ::  save_nxl_mg  !< to save nxl_mg on coarsest level 1
+    INTEGER(iwp) ::  save_nys_mg  !< to save nys_mg on coarsest level 1
+
+    LOGICAL ::  adjust_lower_i_index  !< adjust lower limit of i loop in case of odd number of grid points
+    LOGICAL ::  adjust_lower_j_index  !< adjust lower limit of j loop in case of odd number of grid points
     LOGICAL ::  unroll  !<
 
@@ -855,4 +862,22 @@
                MOD( nxr_mg(l)-nxl_mg(l)+1, 2 ) == 0 )
 
+!
+!-- The red/black decomposition requires that on the lower i,j indices need to start alternatively with an
+!-- even or odd value on the coarsest grid level, depending on the core-id, and if the subdomain has an
+!-- uneven number of gridpoints along x/y. Set the respective steering switches here.
+    IF ( l == 1  .AND.  MOD( myidx, 2 ) /= 0  .AND.  MOD( nxl_mg(l) - nxr_mg(l), 2 ) == 0 )  THEN
+       adjust_lower_i_index = .TRUE.
+       save_nxl_mg = nxl_mg(1)
+    ELSE
+       adjust_lower_i_index = .FALSE.
+    ENDIF
+    IF ( l == 1  .AND.  MOD( myidy, 2 ) /= 0  .AND.  MOD( nyn_mg(l) - nys_mg(l), 2 ) == 0 )  THEN
+       adjust_lower_j_index = .TRUE.
+       save_nys_mg = nys_mg(l)
+    ELSE
+       adjust_lower_j_index = .FALSE.
+    ENDIF
+
+
     DO  n = 1, ngsrb
 
@@ -863,4 +888,15 @@
              CALL cpu_log( log_point_s(36), 'redblack_no_unroll_noopt', 'start' )
 
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg + 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg - 1
+                ELSE
+                  nys_mg(l) = save_nys_mg + 1
+                ENDIF
+             ENDIF
 !
 !--          Without unrolling of loops, no cache optimization
@@ -897,4 +933,16 @@
                 ENDDO
              ENDDO
+
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg - 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg + 1
+                ELSE
+                  nys_mg(l) = save_nys_mg - 1
+                ENDIF
+             ENDIF
 
              DO  i = nxl_mg(l)+1, nxr_mg(l), 2
@@ -924,4 +972,16 @@
              ENDDO
 
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg + 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg + 1
+                ELSE
+                  nys_mg(l) = save_nys_mg - 1
+                ENDIF
+             ENDIF
+
              DO  i = nxl_mg(l), nxr_mg(l), 2
                 DO  j = nys_mg(l) + (color-1), nyn_mg(l), 2
@@ -949,4 +1009,16 @@
                 ENDDO
              ENDDO
+
+             IF ( adjust_lower_i_index )  THEN
+                nxl_mg(l) = save_nxl_mg - 1
+             ENDIF
+
+             IF ( adjust_lower_j_index )  THEN
+                IF ( color == 1 )  THEN
+                  nys_mg(l) = save_nys_mg - 1
+                ELSE
+                  nys_mg(l) = save_nys_mg + 1
+                ENDIF
+             ENDIF
 
              DO  i = nxl_mg(l)+1, nxr_mg(l), 2
@@ -1215,5 +1287,15 @@
 
 !
-!--    Set pressure within topography and at the topography surfaces
+!-- Reset lower index limits to their standard values (may happen on coarsest levels only)
+    IF ( adjust_lower_i_index )  THEN
+       nxl_mg(l) = save_nxl_mg
+    ENDIF
+
+    IF ( adjust_lower_j_index )  THEN
+       nys_mg(l) = save_nys_mg
+    ENDIF
+
+!
+!-- Set pressure within topography and at the topography surfaces
 !$OMP PARALLEL PRIVATE (i,j,k,wall_left,wall_north,wall_right,wall_south,wall_top,wall_total)
 !$OMP DO
Index: /palm/trunk/SOURCE/restart_data_mpi_io_mod.f90
===================================================================
--- /palm/trunk/SOURCE/restart_data_mpi_io_mod.f90	(revision 4831)
+++ /palm/trunk/SOURCE/restart_data_mpi_io_mod.f90	(revision 4832)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! some debug output flushed
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! interface for 3d logical arrays added
 !
@@ -2355,5 +2358,8 @@
                                       INT( (iog%nx+1), KIND = rd_offset_kind ) * wp
 
-    write(9,*) 'array_position real3d ',trim(name),' ',array_position
+    IF ( debug_output )  THEN
+       WRITE (9,*) 'array_position real3d ', TRIM( name ), ' ', array_position
+       FLUSH( 9 )
+    ENDIF
 
  END SUBROUTINE wrd_mpi_io_real_3d
Index: /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc.000
===================================================================
--- /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc.000	(revision 4831)
+++ /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc.000	(revision 4832)
@@ -1,9 +1,9 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4325M       *    atmosphere - nested run without 1D - prerun
+ * PALM 6.0  Rev: 4828M       *    atmosphere - nested run without 1D - prerun
  ******************************    --------------------------------------------
 
- Date:               2019-12-06    Run:       cbl_particle_nested__gfortran_defa
- Time:                 11:03:12    Run-No.:   00
+ Date:               2021-01-05    Run:       cbl_particle_nested__gfortran_defa
+ Time:                 16:37:42    Run-No.:   00
  Run on host:        gfortran_d
  Number of PEs:               2    Processor grid (x,y): (   1,   2) calculated
@@ -25,5 +25,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 multigrid method (w-cycle)
@@ -102,5 +102,5 @@
 
        z_mo =  25.00 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
-       Rif value range:     -20.00 <= rif <= 20.00
+       zeta value range:     -20.00 <= zeta <= 20.00
 
 
@@ -195,5 +195,5 @@
  ---------------------------
 
-    Mixing length is limited to 1.80 * z
+    Mixing length is limited close to surfaces
 
 
@@ -250,40 +250,39 @@
 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      0 00:00:00.00  18.9000A   1.1899D  1.2022D  0.1180  0.000  1.40   0.000E+00    800.  0.953E+00  0.103E-02  0.628E-03  0.725E-06    9   9  23    4   5  25    5  26   7    0.000    0.000       2
-  0      1 00:00:18.90   0.3520A -33.4250 -31.9911   6.8868  0.165  0.00  -0.858E+00      0.  0.194E+01  0.553E+00  0.156E-02  0.143E-04    2  20  16    2  16  20    4  15  20    0.000    0.000       2
-  0      2 00:00:19.25   0.5030A  20.8643  21.2183  -5.0981  0.120  1.34  -0.112E+01    700.  0.146E+01  0.238E+00  0.307E-03  0.101E-04    2  15  21    2  21  15    4  16  21    0.000    0.000       2
-  0      3 00:00:19.75   0.5570A  16.5277  16.6867  -5.3950  0.075  0.56  -0.184E+01     50.  0.114E+01  0.868E-01  0.242E-03  0.383E-05    2  15  21    2  21  15    4  16  21    0.000    0.000       2
-  0      4 00:00:20.31   0.7190A  14.3469  14.4897   5.4556  0.054  0.70  -0.271E+01    100.  0.106E+01  0.545E-01  0.127E-03  0.265E-05    2  15  21    2  21  15    1  20  15    0.000    0.000       2
-  0      5 00:00:21.03   0.8020A  14.0461  14.3200   5.4279  0.045  0.00  -0.325E+01      0.  0.104E+01  0.463E-01  0.747E-04  0.412E-06    2  15  21    2  21  15    1  20  15    0.000    0.000       2
-  0      6 00:00:21.83   0.8380A  14.5514  14.9741  -5.6408  0.042  0.00  -0.360E+01      0.  0.103E+01  0.449E-01  0.480E-04  0.107E-05    2  15  21    2  21  15    4  20  15    0.000    0.000       2
-  0      7 00:00:22.67   0.9250A  14.7003  15.4108  -5.2982  0.041  1.40  -0.381E+01    800.  0.103E+01  0.448E-01  0.360E-04  0.101E-05    2  15  21    3  21  15    4  20  15    0.000    0.000       2
-  0      8 00:00:23.59   1.0500A  14.2090  15.2585  -4.6465  0.040  0.56  -0.387E+01     50.  0.102E+01  0.440E-01  0.295E-04  0.212E-06    2  15  21    3  21  15    4  20  15    0.000    0.000       2
-  0      9 00:00:24.64   1.2200A  13.1540  14.1910  -4.7985  0.040  0.56  -0.386E+01     50.  0.102E+01  0.419E-01  0.286E-04  0.462E-06    2  15  21    3  21  15    3  16  20    0.000    0.000       2
-  0     10 00:00:25.86   1.4200A  11.7065  12.3453  -4.9405  0.040  0.70  -0.381E+01    100.  0.102E+01  0.390E-01  0.309E-04  0.559E-06    2  15  21    3  21  15    3  16  20    0.000    0.000       2
-  0     11 00:00:27.28   1.5700A   9.9801  10.1277  -4.8820  0.040  0.00  -0.372E+01      0.  0.102E+01  0.360E-01  0.332E-04  0.217E-06    2  15  21    2  21  15    3  16  20    0.000    0.000       2
-  0     12 00:00:28.85   1.6800A   8.1590   9.5687  -5.0578  0.041  0.00  -0.361E+01      0.  0.101E+01  0.334E-01  0.345E-04  0.155E-06    3  15  22    3  22  15    3  21  14    0.000    0.000       2
-  0     13 00:00:30.53   1.7600A   8.5870   9.4354  -5.8093  0.041  0.00  -0.350E+01      0.  0.101E+01  0.313E-01  0.356E-04  0.243E-06    3  15  22    3  22  15    3  21  14    0.000    0.000       2
-  0     14 00:00:32.29   1.8300A   8.4095   8.7247  -5.8527  0.042  0.70  -0.338E+01    100.  0.101E+01  0.296E-01  0.362E-04  0.130E-06    3  15  22    3  22  15    3  21  14    0.000    0.000       2
-  0     15 00:00:34.12   1.9000A   7.8212   7.7814  -5.3674  0.043  0.70  -0.327E+01    100.  0.101E+01  0.283E-01  0.363E-04  0.815E-07    3  15  22    3  22  15    3  21  14    0.000    0.000       2
-  0     16 00:00:36.02   1.9400A   7.0132   7.7538  -4.5972  0.044  0.70  -0.317E+01    100.  0.101E+01  0.275E-01  0.364E-04  0.113E-06    3  15  22    3  23  14    3  21  14    0.000    0.000       2
-  0     17 00:00:37.96   1.9900A   6.8381   8.0461  -3.8596  0.045  0.00  -0.307E+01      0.  0.102E+01  0.269E-01  0.363E-04  0.849E-07    3  14  23    3  23  14    2  17  20    0.000    0.000       2
-  0     18 00:00:39.95   2.0800A   7.0176   7.9419  -3.7616  0.046  0.00  -0.298E+01      0.  0.102E+01  0.265E-01  0.366E-04  0.799E-07    3  14  23    3  23  14    2  17  20    0.000    0.000       2
-  0     19 00:00:42.03   2.1400A   6.8958   7.5216  -3.5548  0.046  0.00  -0.289E+01      0.  0.102E+01  0.262E-01  0.372E-04  0.919E-07    3  14  23    3  23  14    3  22  15    0.000    0.000       2
-  0     20 00:00:44.17   2.2600A   6.5505   6.8964  -3.3458  0.047  0.00  -0.281E+01      0.  0.102E+01  0.260E-01  0.369E-04  0.884E-07    3  14  23    3  23  14    3  15  22    0.000    0.000       2
-  0     21 00:00:46.43   2.3200A   6.0519   6.6696  -3.0536  0.048  0.00  -0.273E+01      0.  0.102E+01  0.258E-01  0.372E-04  0.881E-07    3  14  23    3  24  14    3  15  22    0.000    0.000       2
-  0     22 00:00:48.75   2.4400A   5.8086   6.4602  -2.6422  0.049  0.00  -0.265E+01      0.  0.102E+01  0.256E-01  0.366E-04  0.913E-07    3  14  24    3  24  14    3  22  15    0.000    0.000       2
-  0     23 00:00:51.19   2.5500A   5.6371   6.1203  -2.5675  0.050  0.00  -0.258E+01      0.  0.102E+01  0.254E-01  0.366E-04  0.973E-07    3  14  24    3  24  14    3  23  15    0.000    0.000       2
-  0     24 00:00:53.74   2.6900A   5.3790   5.7125  -2.5507  0.051  0.00  -0.251E+01      0.  0.103E+01  0.252E-01  0.364E-04  0.101E-06    3  14  24    3  24  14    3  23  15    0.000    0.000       2
-  0     25 00:00:56.43   2.8400A   5.0888   5.5039  -2.6659  0.052  0.00  -0.244E+01      0.  0.103E+01  0.251E-01  0.367E-04  0.105E-06    3  14  24    3  25  14    2  16  22    0.000    0.000       2
-  0     26 00:00:59.27   3.0300A   4.7656   5.3942  -2.7711  0.053  1.40  -0.238E+01    800.  0.103E+01  0.249E-01  0.373E-04  0.106E-06    3  14  24    3  25  14    2  16  22    0.000    0.000       2
-  0     27 00:01:02.30   3.2100A   4.5773   5.1553  -2.7018  0.054  1.40  -0.231E+01    800.  0.103E+01  0.248E-01  0.385E-04  0.108E-06    3  14  25    3  25  14    2  16  22    0.000    0.000       2
-  0     28 00:01:05.51   3.4400A   4.5620   4.8481  -2.5241  0.055  1.40  -0.225E+01    800.  0.103E+01  0.247E-01  0.402E-04  0.109E-06    3  14  25    3  25  14    3  21  17    0.000    0.000       2
-  0     29 00:01:08.95   3.6600A   4.4420   4.5061  -2.4893  0.057  1.40  -0.219E+01    800.  0.103E+01  0.246E-01  0.425E-04  0.113E-06    3  14  25    3  25  14    3  21  17    0.000    0.000       2
-  0     30 00:01:12.61   3.9100A   4.2531   4.2750  -2.4640  0.058  1.40  -0.214E+01    800.  0.104E+01  0.245E-01  0.450E-04  0.117E-06    3  14  25    3  25  13    3  21  17    0.000    0.000       2
-  0     31 00:01:16.52   4.2300A   4.0178   4.1781  -2.4574  0.059  1.40  -0.209E+01    800.  0.104E+01  0.244E-01  0.480E-04  0.122E-06    3  14  25    3  26  14    3  21  17    0.000    0.000       2
-  0     32 00:01:20.75   4.4400A   3.9913   4.0190  -2.4518  0.060  1.40  -0.204E+01    800.  0.104E+01  0.243E-01  0.522E-04  0.131E-06    3  13  25    3  26  14    3  21  17    0.000    0.000       2
-  0     33 00:01:25.19   4.4900A   3.9102   3.8156  -2.4159  0.061  1.40  -0.199E+01    800.  0.104E+01  0.242E-01  0.555E-04  0.133E-06    3  13  25    3  26  14    3  21  17    0.000    0.000       2
-  0     34 00:01:29.68   4.6400A   3.7848   3.7235  -2.3380  0.062  1.40  -0.195E+01    800.  0.104E+01  0.240E-01  0.566E-04  0.129E-06    3  13  25    3  26  13    3  21  17    0.000    0.000       2
-  0     35 00:01:34.32   4.7800A   3.6185   3.6904  -2.2629  0.063  1.40  -0.191E+01    800.  0.104E+01  0.239E-01  0.583E-04  0.132E-06    3  13  25    3  26  13    3  17  21    0.000    0.000       2
-  0     36 00:01:39.10   5.0100A   3.5946   3.6118  -2.1993  0.064  1.40  -0.187E+01    800.  0.104E+01  0.239E-01  0.591E-04  0.135E-06    3  13  26    3  26  13    3  17  21    0.000    0.000       2
-  0     37 00:01:44.11   5.2000A   3.5345   3.4963  -2.1306  0.064  1.40  -0.185E+01    800.  0.104E+01  0.238E-01  0.608E-04  0.143E-06    3  13  26    3  26  13    3  17  21    0.000    0.000       2
+  0      0 00:00:00.00  18.9000A   1.1899D  1.2022D  0.1181  0.000  1.40   0.000E+00    800.  0.953E+00  0.103E-02  0.628E-03  0.719E-06    9   9  23    4   5  25    5  26   7    0.000    0.000       2
+  0      1 00:00:18.90   0.3890A -30.1784 -28.9672   6.4826  0.163  0.00  -0.866E+00      0.  0.190E+01  0.521E+00  0.151E-02  0.132E-04    2  20  16    2  16  20    4  15  20    0.000    0.000       2
+  0      2 00:00:19.28   0.5580A  19.3826  19.7050  -4.8216  0.119  1.40  -0.113E+01    800.  0.144E+01  0.221E+00  0.297E-03  0.912E-05    2  15  21    2  21  15    4  21  17    0.000    0.000       2
+  0      3 00:00:19.84   0.6240A  15.2675  15.4757  -4.9109  0.074  0.56  -0.186E+01     50.  0.113E+01  0.778E-01  0.234E-03  0.322E-05    2  15  21    2  21  15    4  16  21    0.000    0.000       2
+  0      4 00:00:20.47   0.7910A  13.2768  13.5177   4.9788  0.053  0.88  -0.269E+01    200.  0.105E+01  0.480E-01  0.122E-03  0.187E-05    2  15  21    2  21  15    1  20  15    0.000    0.000       2
+  0      5 00:00:21.26   0.8830A  13.0456  13.4191   4.9272  0.045  0.00  -0.325E+01      0.  0.103E+01  0.409E-01  0.715E-04  0.231E-06    2  15  21    2  21  15    1  20  15    0.000    0.000       2
+  0      6 00:00:22.14   0.9350A  13.4582  13.9466  -5.1436  0.042  0.00  -0.357E+01      0.  0.102E+01  0.397E-01  0.464E-04  0.815E-06    2  15  21    2  21  15    4  20  15    0.000    0.000       2
+  0      7 00:00:23.07   1.0300A  13.4941  14.3254  -4.8069  0.040  1.06  -0.375E+01    350.  0.102E+01  0.395E-01  0.347E-04  0.653E-06    2  15  21    3  21  15    4  20  15    0.000    0.000       2
+  0      8 00:00:24.10   1.1800A  12.9295  14.0145  -4.1814  0.040  0.56  -0.382E+01     50.  0.102E+01  0.386E-01  0.290E-04  0.557E-07    2  15  21    3  21  15    4  20  15    0.000    0.000       2
+  0      9 00:00:25.28   1.3800A  11.9921  12.8651  -4.3074  0.040  0.56  -0.381E+01     50.  0.102E+01  0.367E-01  0.292E-04  0.382E-06    3  15  21    3  21  15    3  16  20    0.000    0.000       2
+  0     10 00:00:26.66   1.6000A  10.5094  11.0438  -4.3995  0.040  0.70  -0.375E+01    100.  0.101E+01  0.342E-01  0.318E-04  0.369E-06    3  15  21    3  21  15    3  16  20    0.000    0.000       2
+  0     11 00:00:28.26   1.7200A   8.8553   9.0416  -4.4723  0.040  0.00  -0.366E+01      0.  0.101E+01  0.318E-01  0.342E-04  0.972E-07    2  15  21    2  21  15    3  21  15    0.000    0.000       2
+  0     12 00:00:29.98   1.8000A   8.2386   9.4383  -4.6579  0.041  0.00  -0.355E+01      0.  0.101E+01  0.297E-01  0.352E-04  0.144E-06    3  15  22    3  22  15    3  21  15    0.000    0.000       2
+  0     13 00:00:31.78   1.9100A   8.4369   9.1081  -5.1521  0.042  0.00  -0.343E+01      0.  0.101E+01  0.281E-01  0.357E-04  0.158E-06    3  15  22    3  22  15    3  21  14    0.000    0.000       2
+  0     14 00:00:33.70   1.9300A   8.0525   8.2807  -5.0689  0.042  0.70  -0.331E+01    100.  0.101E+01  0.267E-01  0.365E-04  0.893E-07    3  15  22    3  22  15    3  21  14    0.000    0.000       2
+  0     15 00:00:35.63   2.0600A   7.3409   7.3208  -4.5692  0.043  0.70  -0.319E+01    100.  0.101E+01  0.258E-01  0.360E-04  0.834E-07    3  15  22    3  22  15    3  21  14    0.000    0.000       2
+  0     16 00:00:37.69   2.0700A   6.4622   7.3564  -3.8332  0.044  0.70  -0.309E+01    100.  0.101E+01  0.252E-01  0.369E-04  0.897E-07    3  15  22    3  23  14    3  21  14    0.000    0.000       2
+  0     17 00:00:39.76   2.1800A   6.6219   7.5684  -3.4207  0.045  0.00  -0.299E+01      0.  0.101E+01  0.249E-01  0.366E-04  0.782E-07    3  14  23    3  23  14    2  17  20    0.000    0.000       2
+  0     18 00:00:41.94   2.2300A   6.7280   7.4219  -3.4818  0.046  0.00  -0.290E+01      0.  0.102E+01  0.247E-01  0.376E-04  0.847E-07    3  14  23    3  23  14    3  22  15    0.000    0.000       2
+  0     19 00:00:44.17   2.3400A   6.5568   6.9927  -3.3627  0.047  0.00  -0.283E+01      0.  0.102E+01  0.246E-01  0.375E-04  0.868E-07    3  14  23    3  23  14    3  15  22    0.000    0.000       2
+  0     20 00:00:46.51   2.4300A   6.1719   6.4077  -3.1450  0.048  0.00  -0.275E+01      0.  0.102E+01  0.245E-01  0.375E-04  0.875E-07    3  14  23    3  24  14    3  15  22    0.000    0.000       2
+  0     21 00:00:48.94   2.5500A   5.7065   6.4463  -2.7396  0.049  0.00  -0.267E+01      0.  0.102E+01  0.244E-01  0.371E-04  0.897E-07    3  14  24    3  24  14    3  15  22    0.000    0.000       2
+  0     22 00:00:51.49   2.6800A   5.7107   6.2693  -2.4251  0.050  0.00  -0.259E+01      0.  0.102E+01  0.243E-01  0.370E-04  0.953E-07    3  14  24    3  24  14    3  23  15    0.000    0.000       2
+  0     23 00:00:54.17   2.7800A   5.5680   5.9650  -2.4976  0.051  0.00  -0.252E+01      0.  0.102E+01  0.242E-01  0.370E-04  0.100E-06    3  14  24    3  24  14    2  16  22    0.000    0.000       2
+  0     24 00:00:56.95   2.9300A   5.3542   5.5793  -2.8245  0.052  1.40  -0.245E+01    800.  0.103E+01  0.241E-01  0.369E-04  0.102E-06    3  14  24    3  24  14    2  16  22    0.000    0.000       2
+  0     25 00:00:59.88   3.0700A   5.0772   5.4271  -2.9557  0.053  1.40  -0.238E+01    800.  0.103E+01  0.240E-01  0.375E-04  0.105E-06    3  14  24    3  25  14    2  16  22    0.000    0.000       2
+  0     26 00:01:02.95   3.2500A   4.7274   5.3021  -2.8887  0.054  1.40  -0.232E+01    800.  0.103E+01  0.239E-01  0.383E-04  0.106E-06    3  14  24    3  25  14    2  16  22    0.000    0.000       2
+  0     27 00:01:06.20   3.4400A   4.6758   5.0658  -2.6503  0.055  1.40  -0.226E+01    800.  0.103E+01  0.238E-01  0.400E-04  0.109E-06    3  14  25    3  25  14    2  16  22    0.000    0.000       2
+  0     28 00:01:09.64   3.6600A   4.6093   4.7694  -2.5041  0.056  1.40  -0.220E+01    800.  0.103E+01  0.238E-01  0.418E-04  0.111E-06    3  14  25    3  25  14    2  22  16    0.000    0.000       2
+  0     29 00:01:13.29   3.8900A   4.4540   4.4476  -2.4318  0.057  1.40  -0.215E+01    800.  0.103E+01  0.237E-01  0.444E-04  0.114E-06    3  14  25    3  25  14    3  21  17    0.000    0.000       2
+  0     30 00:01:17.18   4.1500A   4.2397   4.3180  -2.4000  0.058  1.40  -0.209E+01    800.  0.104E+01  0.237E-01  0.473E-04  0.118E-06    3  14  25    3  26  14    3  21  17    0.000    0.000       2
+  0     31 00:01:21.34   4.4400A   3.9912   4.2193  -2.3666  0.060  1.40  -0.204E+01    800.  0.104E+01  0.236E-01  0.511E-04  0.123E-06    3  14  25    3  26  14    3  21  17    0.000    0.000       2
+  0     32 00:01:25.78   4.5100A   3.7749   4.0471  -2.3100  0.061  1.40  -0.200E+01    800.  0.104E+01  0.235E-01  0.551E-04  0.129E-06    3  14  26    3  26  14    3  21  17    0.000    0.000       2
+  0     33 00:01:30.29   4.7100A   3.7265   3.8455  -2.3274  0.062  1.40  -0.195E+01    800.  0.104E+01  0.234E-01  0.564E-04  0.125E-06    3  14  26    3  26  14    3  17  21    0.000    0.000       2
+  0     34 00:01:35.00   4.8500A   3.6350   3.6402  -2.3003  0.063  1.40  -0.192E+01    800.  0.104E+01  0.233E-01  0.585E-04  0.130E-06    3  14  26    3  26  14    3  17  21    0.000    0.000       2
+  0     35 00:01:39.84   5.0800A   3.5268   3.5421  -2.2304  0.063  1.40  -0.188E+01    800.  0.104E+01  0.233E-01  0.593E-04  0.133E-06    3  14  26    3  27  14    3  17  21    0.000    0.000       2
+  0     36 00:01:44.92   5.2900A   3.4104   3.4631  -2.1461  0.064  1.40  -0.185E+01    800.  0.104E+01  0.232E-01  0.611E-04  0.142E-06    3  14  26    3  27  14    3  17  21    0.000    0.000       2
Index: /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc_N02.000
===================================================================
--- /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc_N02.000	(revision 4831)
+++ /palm/trunk/TESTS/cases/cbl_particle_nested/MONITORING/cbl_particle_nested_rc_N02.000	(revision 4832)
@@ -1,9 +1,9 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4325M       *    atmosphere - nested run without 1D - prerun
+ * PALM 6.0  Rev: 4828M       *    atmosphere - nested run without 1D - prerun
  ******************************    --------------------------------------------
 
- Date:               2019-12-06    Run:       cbl_particle_nested__gfortran_defa
- Time:                 11:03:12    Run-No.:   00
+ Date:               2021-01-05    Run:       cbl_particle_nested__gfortran_defa
+ Time:                 16:37:42    Run-No.:   00
  Run on host:        gfortran_d
  Number of PEs:               2    Processor grid (x,y): (   1,   2) calculated
@@ -25,5 +25,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 multigrid method (w-cycle)
@@ -108,5 +108,5 @@
 
        z_mo =  12.50 m   z0 = 0.1000 m   z0h = 0.10000 m   kappa = 0.40
-       Rif value range:     -20.00 <= rif <= 20.00
+       zeta value range:     -20.00 <= zeta <= 20.00
 
 
@@ -204,5 +204,5 @@
  ---------------------------
 
-    Mixing length is limited to 1.80 * z
+    Mixing length is limited close to surfaces
 
 
@@ -259,40 +259,39 @@
 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      0 00:00:00.00  18.9000A   1.1899D  1.1775D  0.1180  0.000  0.00   0.000E+00    800.  0.923E+00  0.374E-02  0.481E-03  0.119E-03   17   2  30   19  37  14    9  36  37    0.000    0.000       2
-  0      1 00:00:18.90   0.3520A -63.9848 -61.1453  17.1772  0.334  0.00  -0.329E+00      0.  0.865E+01  0.202E+01  0.333E-01  0.532E-02    4  25  16    4  16  25    8  15  25    0.000    0.000       2
-  0      2 00:00:19.25   0.5030A -44.7558 -43.0697 -13.1112  0.212  0.00  -0.529E+00     25.  0.434E+01  0.128E+01  0.125E-02  0.366E-03    4  25  16    4  16  25    8  16  26    0.000    0.000       2
-  0      3 00:00:19.75   0.5570A -40.4039 -38.6866 -12.5327  0.134  0.00  -0.101E+01     25.  0.248E+01  0.962E+00  0.174E-02  0.635E-03    4  25  15    4  15  25    8  16  26    0.000    0.000       2
-  0      4 00:00:20.31   0.7190D -30.9683 -29.9633 -11.4752  0.110  0.68  -0.138E+01    175.  0.201E+01  0.792E+00  0.760E-03  0.172E-03    4  25  15    4  15  25    8  26  16    0.000    0.000       2
-  0      5 00:00:21.03   0.8020D  25.1050  25.3409 -11.0068  0.102  0.00  -0.149E+01      0.  0.180E+01  0.662E+00  0.519E-03  0.426E-04    5  15  25    5  25  15    8  17  26    0.000    0.000       2
-  0      6 00:00:21.83   0.8380D  21.0110  21.2174 -10.2671  0.098  0.00  -0.154E+01     25.  0.167E+01  0.568E+00  0.434E-03  0.163E-04    4  15  26    4  26  15    8  17  26    0.000    0.000       2
-  0      7 00:00:22.67   0.9250D  19.9158  19.9818 -10.5432  0.096  0.00  -0.156E+01     25.  0.159E+01  0.501E+00  0.389E-03  0.123E-04    4  15  26    4  26  15    7  25  14    0.000    0.000       2
-  0      8 00:00:23.59   1.0500D  17.5472  17.4174 -10.4681  0.095  0.00  -0.158E+01     25.  0.153E+01  0.451E+00  0.372E-03  0.117E-04    4  15  26    4  26  15    7  25  14    0.000    0.000       2
-  0      9 00:00:24.64   1.2200D  14.4740  15.3214  -9.5317  0.094  0.00  -0.159E+01     25.  0.149E+01  0.411E+00  0.376E-03  0.117E-04    3  15  26    5  26  14    7  25  14    0.000    0.000       2
-  0     10 00:00:25.86   1.4200D  13.4505  14.7448  -9.9304  0.093  0.00  -0.170E+01     25.  0.146E+01  0.378E+00  0.396E-03  0.124E-04    5  14  26    5  26  14    6  17  25    0.000    0.000       2
-  0     11 00:00:27.28   1.5700A  12.8853  14.3517  -9.8167  0.093  0.00  -0.156E+01     25.  0.143E+01  0.349E+00  0.425E-03  0.132E-04    4  14  27    5  27  14    6  17  25    0.000    0.000       2
-  0     12 00:00:28.85   1.6800A  12.2611  13.3843  -9.1695  0.092  0.00  -0.154E+01     25.  0.141E+01  0.325E+00  0.439E-03  0.134E-04    4  14  27    5  28  14    6  17  25    0.000    0.000       2
-  0     13 00:00:30.53   1.7600A  11.6471  12.8129  -8.3230  0.092  0.00  -0.153E+01     25.  0.140E+01  0.304E+00  0.445E-03  0.129E-04    5  14  28    5  28  14    6  17  25    0.000    0.000       2
-  0     14 00:00:32.29   1.8300A  11.0514  12.3084  -7.7091  0.093  0.00  -0.155E+01     25.  0.139E+01  0.287E+00  0.442E-03  0.122E-04    5  14  28    5  29  14    5  18  25    0.000    0.000       2
-  0     15 00:00:34.12   1.9000A  10.7391  11.8181  -7.0408  0.093  0.00  -0.151E+01     25.  0.139E+01  0.273E+00  0.436E-03  0.115E-04    5  14  29    5  29  13    5  18  25    0.000    0.000       2
-  0     16 00:00:36.02   1.9400A  10.4235  11.6045  -6.6832  0.094  0.00  -0.147E+01     25.  0.138E+01  0.260E+00  0.430E-03  0.109E-04    5  13  29    5  30  13    5  16  26    0.000    0.000       2
-  0     17 00:00:37.96   1.9900A  10.3760  11.2840  -6.2488  0.095  0.00  -0.141E+01     25.  0.138E+01  0.250E+00  0.419E-03  0.101E-04    5  13  30    5  30  13    5  16  26    0.000    0.000       2
-  0     18 00:00:39.95   2.0800A  10.0936  10.8156  -5.7274  0.096  0.00  -0.136E+01     25.  0.138E+01  0.241E+00  0.410E-03  0.947E-05    5  13  30    5  31  13    5  16  26    0.000    0.000       2
-  0     19 00:00:42.03   2.1400A   9.7042  10.5336  -5.1659  0.097  0.00  -0.132E+01     25.  0.139E+01  0.233E+00  0.409E-03  0.899E-05    5  13  31    5  31  13    5  16  26    0.000    0.000       2
-  0     20 00:00:44.17   2.2600A   9.5063   9.9595  -5.0962  0.098  0.00  -0.129E+01     25.  0.139E+01  0.226E+00  0.403E-03  0.848E-05    5  13  31    5  31  13    4  28  15    0.000    0.000       2
-  0     21 00:00:46.43   2.3200A   8.9465   9.6992  -4.9576  0.099  0.00  -0.126E+01     25.  0.139E+01  0.219E+00  0.406E-03  0.821E-05    5  13  31    5  32  13    4  29  15    0.000    0.000       2
-  0     22 00:00:48.75   2.4400A   8.7719   9.2136  -4.9280  0.100  0.00  -0.123E+01     25.  0.140E+01  0.213E+00  0.399E-03  0.783E-05    5  13  32    5  32  13    4  29  15    0.000    0.000       2
-  0     23 00:00:51.19   2.5500A   8.3179   8.8253  -4.6153  0.101  0.00  -0.121E+01     25.  0.140E+01  0.208E+00  0.402E-03  0.773E-05    5  13  32    5  33  13    4  16  28    0.000    0.000       2
-  0     24 00:00:53.74   2.6900A   7.9621   8.3550  -4.7407  0.102  0.00  -0.119E+01     25.  0.141E+01  0.202E+00  0.400E-03  0.769E-05    5  13  33    5  33  13    4  16  28    0.000    0.000       2
-  0     25 00:00:56.43   2.8400A   7.5721   7.9157  -4.4733  0.104  0.00  -0.116E+01     25.  0.141E+01  0.197E+00  0.403E-03  0.779E-05    5  13  33    5  34  13    4  16  28    0.000    0.000       2
-  0     26 00:00:59.27   3.0300A   7.1362   7.4279  -4.1833  0.105  0.00  -0.114E+01     25.  0.142E+01  0.192E+00  0.405E-03  0.795E-05    5  13  34    5  34  13    3  16  28    0.000    0.000       2
-  0     27 00:01:02.30   3.2100A   6.7743   7.0069  -4.0566  0.106  0.00  -0.112E+01     25.  0.142E+01  0.187E+00  0.410E-03  0.825E-05    5  13  34    5  35  13    3  16  29    0.000    0.000       2
-  0     28 00:01:05.51   3.4400A   6.3246   6.5444  -3.8265  0.107  0.00  -0.110E+01     25.  0.143E+01  0.182E+00  0.414E-03  0.859E-05    5  13  35    5  36  13    3  16  29    0.000    0.000       2
-  0     29 00:01:08.95   3.6600A   5.9737   6.1447  -3.7167  0.108  0.00  -0.109E+01     25.  0.143E+01  0.177E+00  0.424E-03  0.902E-05    5  13  35    5  36  13    3  16  30    0.000    0.000       2
-  0     30 00:01:12.61   3.9100A   5.5897   5.7604  -3.4834  0.109  0.00  -0.107E+01     25.  0.144E+01  0.172E+00  0.435E-03  0.943E-05    5  12  35    5  37  13    3  30  16    0.000    0.000       2
-  0     31 00:01:16.52   4.2300A   5.2636   5.3240  -3.4369  0.111  0.00  -0.105E+01     25.  0.144E+01  0.167E+00  0.449E-03  0.987E-05    5  12  36    5  37  13    3  31  16    0.000    0.000       2
-  0     32 00:01:20.75   4.4400A   4.9538   5.0678  -3.1629  0.112  0.00  -0.103E+01     25.  0.145E+01  0.162E+00  0.470E-03  0.104E-04    5  12  36    6  35  11    3  32  16    0.000    0.000       2
-  0     33 00:01:25.19   4.4900A   4.6056   5.0137  -3.0940  0.113  0.00  -0.102E+01     25.  0.145E+01  0.157E+00  0.480E-03  0.106E-04    6  11  35    6  36  12    3  32  16    0.000    0.000       2
-  0     34 00:01:29.68   4.6400A   4.5376   4.8542  -3.0885  0.114  0.00  -0.101E+01     25.  0.145E+01  0.152E+00  0.477E-03  0.103E-04    6  12  36    6  37  12    6  26  19    0.000    0.000       2
-  0     35 00:01:34.32   4.7800A   4.5189   4.7050  -3.0195  0.115  0.00  -0.100E+01     25.  0.146E+01  0.148E+00  0.486E-03  0.102E-04    6  11  36    6  37  11    6  26  19    0.000    0.000       2
-  0     36 00:01:39.10   5.0100A   4.3825   4.4916  -2.9043  0.116  0.00  -0.993E+00     25.  0.146E+01  0.144E+00  0.498E-03  0.100E-04    6  12  37    6  37  11    6  26  19    0.000    0.000       2
-  0     37 00:01:44.11   5.2000A   4.3249   4.2997  -2.7609  0.116  0.00  -0.985E+00     25.  0.145E+01  0.140E+00  0.520E-03  0.999E-05    6  11  37    6  38  11    6  26  19    0.000    0.000       2
+  0      0 00:00:00.00  18.9000A   1.1899D  1.1775D  0.1181  0.000  0.00   0.000E+00    800.  0.923E+00  0.374E-02  0.481E-03  0.119E-03   17   2  30   19  37  14    9  36  37    0.000    0.000       2
+  0      1 00:00:18.90   0.3890A -57.8201 -55.8650  15.9007  0.329  0.00  -0.333E+00      0.  0.830E+01  0.179E+01  0.323E-01  0.516E-02    4  25  16    4  16  25    8  15  25    0.000    0.000       2
+  0      2 00:00:19.28   0.5580A -40.3509 -39.1933 -11.7860  0.208  0.00  -0.534E+00     25.  0.413E+01  0.110E+01  0.122E-02  0.354E-03    4  25  16    4  16  25    8  16  26    0.000    0.000       2
+  0      3 00:00:19.84   0.6240A -36.0364 -34.9102 -11.2993  0.131  0.00  -0.101E+01     25.  0.233E+01  0.824E+00  0.169E-02  0.617E-03    4  25  15    4  15  25    8  16  26    0.000    0.000       2
+  0      4 00:00:20.47   0.7910D -27.2677  26.7024 -10.5584  0.107  0.54  -0.142E+01    100.  0.190E+01  0.676E+00  0.741E-03  0.166E-03    4  25  15    5  25  15    8  26  16    0.000    0.000       2
+  0      5 00:00:21.26   0.8830D  22.8674  23.0801  -9.9843  0.099  0.00  -0.153E+01      0.  0.171E+01  0.568E+00  0.497E-03  0.383E-04    5  15  25    5  25  15    8  17  26    0.000    0.000       2
+  0      6 00:00:22.14   0.9350D  19.5891  19.9522   9.4314  0.095  0.00  -0.158E+01     25.  0.160E+01  0.489E+00  0.418E-03  0.147E-04    4  15  26    4  26  15    2  25  14    0.000    0.000       2
+  0      7 00:00:23.07   1.0300D  18.4335  18.6360  -9.4652  0.093  0.00  -0.160E+01     25.  0.152E+01  0.433E+00  0.379E-03  0.115E-04    4  15  26    4  26  15    7  25  14    0.000    0.000       2
+  0      8 00:00:24.10   1.1800D  16.1376  16.1405  -9.3660  0.092  0.00  -0.163E+01     25.  0.147E+01  0.390E+00  0.365E-03  0.110E-04    4  15  26    4  26  15    7  25  14    0.000    0.000       2
+  0      9 00:00:25.28   1.3800D  13.3625  13.6811  -8.7904  0.091  0.00  -0.161E+01     25.  0.144E+01  0.357E+00  0.376E-03  0.114E-04    3  15  26    5  26  14    7  17  25    0.000    0.000       2
+  0     10 00:00:26.66   1.6000D  11.9854  13.3309  -9.0287  0.091  0.00  -0.160E+01     25.  0.141E+01  0.328E+00  0.401E-03  0.122E-04    5  14  26    5  27  14    6  17  25    0.000    0.000       2
+  0     11 00:00:28.26   1.7200A  11.9411  13.0588  -8.8949  0.090  0.00  -0.157E+01     25.  0.139E+01  0.304E+00  0.431E-03  0.130E-04    5  14  27    5  27  14    6  17  25    0.000    0.000       2
+  0     12 00:00:29.98   1.8000A  11.1461  12.5302  -8.3169  0.090  0.00  -0.156E+01     25.  0.138E+01  0.285E+00  0.437E-03  0.129E-04    5  14  27    5  28  14    6  17  25    0.000    0.000       2
+  0     13 00:00:31.78   1.9100A  10.9720  11.7712  -7.5901  0.091  0.00  -0.156E+01     25.  0.137E+01  0.268E+00  0.434E-03  0.123E-04    5  14  28    5  28  14    6  17  25    0.000    0.000       2
+  0     14 00:00:33.70   1.9300A  10.2189  11.6300  -6.9607  0.091  0.00  -0.155E+01     25.  0.137E+01  0.254E+00  0.438E-03  0.118E-04    5  14  29    5  29  14    5  18  25    0.000    0.000       2
+  0     15 00:00:35.63   2.0600A  10.2722  10.9393  -6.4954  0.092  0.00  -0.149E+01     25.  0.136E+01  0.243E+00  0.422E-03  0.110E-04    5  14  29    5  29  14    5  16  26    0.000    0.000       2
+  0     16 00:00:37.69   2.0700A   9.6689  10.8704  -6.1447  0.093  0.00  -0.142E+01     25.  0.136E+01  0.233E+00  0.429E-03  0.107E-04    5  14  29    5  30  14    5  16  26    0.000    0.000       2
+  0     17 00:00:39.76   2.1800A   9.7005  10.3106  -5.6711  0.095  0.00  -0.137E+01     25.  0.137E+01  0.225E+00  0.412E-03  0.995E-05    5  14  30    5  30  14    5  16  26    0.000    0.000       2
+  0     18 00:00:41.94   2.2300A   9.2412  10.0857  -5.1245  0.096  0.00  -0.133E+01     25.  0.137E+01  0.217E+00  0.414E-03  0.960E-05    5  14  30    5  31  14    5  16  26    0.000    0.000       2
+  0     19 00:00:44.17   2.3400A   9.0825   9.6022  -4.6131  0.097  0.00  -0.129E+01     25.  0.137E+01  0.211E+00  0.404E-03  0.903E-05    5  13  31    5  31  14    5  26  16    0.000    0.000       2
+  0     20 00:00:46.51   2.4300A   8.7019   9.2422  -4.4119  0.098  0.00  -0.126E+01     25.  0.138E+01  0.205E+00  0.406E-03  0.872E-05    5  13  31    5  32  14    5  16  27    0.000    0.000       2
+  0     21 00:00:48.94   2.5500A   8.3839   8.8332  -4.5063  0.099  0.00  -0.123E+01     25.  0.138E+01  0.200E+00  0.403E-03  0.839E-05    5  13  32    5  32  13    4  29  15    0.000    0.000       2
+  0     22 00:00:51.49   2.6800A   8.0947   8.3805  -4.5658  0.100  0.00  -0.122E+01     25.  0.139E+01  0.195E+00  0.404E-03  0.822E-05    5  13  32    5  33  13    4  16  28    0.000    0.000       2
+  0     23 00:00:54.17   2.7800A   7.6473   8.0862  -4.6960  0.101  0.00  -0.119E+01     25.  0.139E+01  0.190E+00  0.405E-03  0.816E-05    5  13  33    5  33  13    4  16  28    0.000    0.000       2
+  0     24 00:00:56.95   2.9300A   7.4158   7.6777  -4.4018  0.103  0.00  -0.120E+01     25.  0.140E+01  0.186E+00  0.401E-03  0.809E-05    5  13  33    5  33  13    4  16  28    0.000    0.000       2
+  0     25 00:00:59.88   3.0700A   6.9984   7.3225  -4.2708  0.104  0.00  -0.115E+01     25.  0.141E+01  0.181E+00  0.402E-03  0.823E-05    5  13  33    5  34  13    4  16  29    0.000    0.000       2
+  0     26 00:01:02.95   3.2500A   6.7146   6.9298  -4.0033  0.105  0.00  -0.113E+01     25.  0.141E+01  0.177E+00  0.401E-03  0.840E-05    5  13  34    5  34  13    4  29  16    0.000    0.000       2
+  0     27 00:01:06.20   3.4400A   6.3365   6.5474  -3.7209  0.106  0.00  -0.111E+01     25.  0.142E+01  0.172E+00  0.405E-03  0.871E-05    5  13  34    5  35  13    4  30  16    0.000    0.000       2
+  0     28 00:01:09.64   3.6600A   6.0087   6.1558  -3.3984  0.107  0.00  -0.109E+01     25.  0.142E+01  0.168E+00  0.412E-03  0.901E-05    5  13  35    5  35  13    4  30  16    0.000    0.000       2
+  0     29 00:01:13.29   3.8900A   5.6558   5.7870  -3.1821  0.109  0.00  -0.107E+01     25.  0.143E+01  0.164E+00  0.424E-03  0.939E-05    5  13  35    5  36  13    4  31  16    0.000    0.000       2
+  0     30 00:01:17.18   4.1500A   5.3170   5.4199  -3.0883  0.110  0.00  -0.105E+01     25.  0.143E+01  0.159E+00  0.437E-03  0.977E-05    5  13  36    5  37  13    3  31  16    0.000    0.000       2
+  0     31 00:01:21.34   4.4400A   4.9940   5.0629  -3.0153  0.111  0.00  -0.104E+01     25.  0.144E+01  0.154E+00  0.452E-03  0.101E-04    5  13  36    6  35  12    6  26  19    0.000    0.000       2
+  0     32 00:01:25.78   4.5100A   4.6893   4.9884  -3.0669  0.112  0.00  -0.102E+01     25.  0.144E+01  0.150E+00  0.472E-03  0.105E-04    5  13  36    6  36  12    6  26  19    0.000    0.000       2
+  0     33 00:01:30.29   4.7100A   4.5572   4.7773  -3.0449  0.113  0.00  -0.101E+01     25.  0.144E+01  0.146E+00  0.472E-03  0.102E-04    6  12  36    6  36  12    6  26  19    0.000    0.000       2
+  0     34 00:01:35.00   4.8500A   4.3750   4.6349  -2.9550  0.114  0.00  -0.100E+01     25.  0.145E+01  0.142E+00  0.487E-03  0.102E-04    6  12  36    6  37  12    6  26  19    0.000    0.000       2
+  0     35 00:01:39.84   5.0800A   4.3144   4.4301  -2.8259  0.115  0.00  -0.997E+00     25.  0.145E+01  0.138E+00  0.499E-03  0.100E-04    6  12  37    6  37  12    6  26  19    0.000    0.000       2
+  0     36 00:01:44.92   5.2900A   4.1248   4.2526  -2.6795  0.116  0.00  -0.990E+00     25.  0.145E+01  0.135E+00  0.521E-03  0.995E-05    6  12  37    6  38  12    6  26  19    0.000    0.000       2
