Index: palm/trunk/SOURCE/Makefile
===================================================================
--- palm/trunk/SOURCE/Makefile	(revision 4101)
+++ palm/trunk/SOURCE/Makefile	(revision 4102)
@@ -25,4 +25,7 @@
 # -----------------
 # $Id$
+# Remove dependency on pmc_interface for boundary_conds
+# 
+# 4070 2019-07-03 13:51:40Z gronemeier
 # Add new data output modules
 #
@@ -802,10 +805,8 @@
 	radiation_model_mod.o
 boundary_conds.o: \
-	basic_constants_and_equations_mod.o \
 	bulk_cloud_model_mod.o \
 	chemistry_model_mod.o \
 	mod_kinds.o \
 	modules.o \
-	pmc_interface_mod.o \
 	salsa_mod.o \
 	surface_mod.o \
Index: palm/trunk/SOURCE/boundary_conds.f90
===================================================================
--- palm/trunk/SOURCE/boundary_conds.f90	(revision 4101)
+++ palm/trunk/SOURCE/boundary_conds.f90	(revision 4102)
@@ -25,4 +25,11 @@
 ! -----------------
 ! $Id$
+! - Revise setting for boundary conditions at horizontal walls, use the offset
+!   index that belongs to the data structure instead of pre-calculating 
+!   the offset index for each facing.
+! - Set boundary conditions for bulk-cloud quantities also at downward facing
+!   surfaces
+! 
+! 4087 2019-07-11 11:35:23Z gronemeier
 ! Add comment
 ! 
@@ -224,12 +231,9 @@
     USE arrays_3d,                                                             &
         ONLY:  c_u, c_u_m, c_u_m_l, c_v, c_v_m, c_v_m_l, c_w, c_w_m, c_w_m_l,  &
-               diss, diss_p, dzu, e_p, nc_p, nr_p, pt, pt_init, pt_p, q,       &
+               dzu, nc_p, nr_p, pt, pt_init, pt_p, q,                          &
                q_p, qc_p, qr_p, s, s_p, sa, sa_p, u, u_init, u_m_l, u_m_n,     &
                u_m_r, u_m_s, u_p, v, v_init, v_m_l, v_m_n, v_m_r, v_m_s, v_p,  &
                w, w_p, w_m_l, w_m_n, w_m_r, w_m_s
 
-    USE basic_constants_and_equations_mod,                                     &
-        ONLY:  kappa
-
     USE bulk_cloud_model_mod,                                                  &
         ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert
@@ -239,12 +243,12 @@
 
     USE control_parameters,                                                    &
-        ONLY:  air_chemistry, bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,  &
+        ONLY:  air_chemistry, bc_dirichlet_l,                                  &
                bc_dirichlet_s, bc_radiation_l, bc_radiation_n, bc_radiation_r, &
                bc_radiation_s, bc_pt_t_val, bc_q_t_val, bc_s_t_val,            &
-               child_domain, constant_diffusion, coupling_mode, dt_3d,         &
+               child_domain, coupling_mode, dt_3d,                             &
                humidity, ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_s_b,        &
                ibc_s_t, ibc_uv_b, ibc_uv_t, intermediate_timestep_count,       &
                nesting_offline, neutral, nudging, ocean_mode, passive_scalar,  &
-               rans_mode, rans_tke_e, tsc, salsa, use_cmax
+               tsc, salsa, use_cmax
 
     USE grid_variables,                                                        &
@@ -262,5 +266,5 @@
 
     USE pmc_interface,                                                         &
-        ONLY : nesting_mode, rans_mode_parent
+        ONLY : nesting_mode
  
     USE salsa_mod,                                                             &
@@ -268,9 +272,8 @@
 
     USE surface_mod,                                                           &
-        ONLY :  bc_h, surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v,          &
-                surf_usm_h, surf_usm_v
+        ONLY :  bc_h
 
     USE turbulence_closure_mod,                                                &
-        ONLY:  c_0
+        ONLY:  tcm_boundary_conds
 
     IMPLICIT NONE
@@ -279,5 +282,4 @@
     INTEGER(iwp) ::  j  !< grid index y direction
     INTEGER(iwp) ::  k  !< grid index z direction
-    INTEGER(iwp) ::  kb !< variable to set respective boundary value, depends on facing. 
     INTEGER(iwp) ::  l  !< running index boundary type, for up- and downward-facing walls
     INTEGER(iwp) ::  m  !< running index surface elements
@@ -297,8 +299,4 @@
 !-- of downward-facing surfaces. 
     DO  l = 0, 1
-!
-!--    Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--    for downward-facing surfaces at topography bottom (k+1). 
-       kb = MERGE( -1, 1, l == 0 )
        !$OMP PARALLEL DO PRIVATE( i, j, k )
        !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
@@ -308,5 +306,5 @@
           j = bc_h(l)%j(m)
           k = bc_h(l)%k(m)
-          w_p(k+kb,j,i) = 0.0_wp
+          w_p(k+bc_h(l)%koff,j,i) = 0.0_wp
        ENDDO
     ENDDO
@@ -341,8 +339,4 @@
        IF ( ibc_pt_b == 0 )  THEN
           DO  l = 0, 1
-!      
-!--          Set kb, for upward-facing surfaces value at topography top (k-1)
-!--          is set, for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              DO  m = 1, bc_h(l)%ns
@@ -350,5 +344,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                pt_p(k+kb,j,i) = pt(k+kb,j,i)
+                pt_p(k+bc_h(l)%koff,j,i) = pt(k+bc_h(l)%koff,j,i)
              ENDDO
           ENDDO
@@ -357,8 +351,4 @@
        ELSEIF ( ibc_pt_b == 1 )  THEN
           DO  l = 0, 1
-!      
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) 
-!--          is set, for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
@@ -368,5 +358,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                pt_p(k+kb,j,i) = pt_p(k,j,i)
+                pt_p(k+bc_h(l)%koff,j,i) = pt_p(k,j,i)
              ENDDO
           ENDDO
@@ -391,192 +381,4 @@
        ENDIF
     ENDIF
-
-!
-!-- Boundary conditions for TKE.
-!-- Generally Neumann conditions with de/dz=0 are assumed.
-    IF ( .NOT. constant_diffusion )  THEN
-
-       IF ( .NOT. rans_mode )  THEN
-          DO  l = 0, 1
-!
-!--         Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--         for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
-             !$OMP PARALLEL DO PRIVATE( i, j, k )
-             !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
-             !$ACC PRESENT(bc_h, e_p)
-             DO  m = 1, bc_h(l)%ns
-                i = bc_h(l)%i(m)            
-                j = bc_h(l)%j(m)
-                k = bc_h(l)%k(m)
-                e_p(k+kb,j,i) = e_p(k,j,i)
-             ENDDO
-          ENDDO
-       ELSE
-!
-!--       Use wall function within constant-flux layer
-!--       Note, grid points listed in bc_h are not included in any calculations in RANS mode and
-!--       are therefore not set here.
-!
-!--       Upward-facing surfaces
-!--       Default surfaces
-          DO  m = 1, surf_def_h(0)%ns
-             i = surf_def_h(0)%i(m)
-             j = surf_def_h(0)%j(m)
-             k = surf_def_h(0)%k(m)
-             e_p(k,j,i) = ( surf_def_h(0)%us(m) / c_0 )**2
-          ENDDO
-!
-!--       Natural surfaces
-          DO  m = 1, surf_lsm_h%ns
-             i = surf_lsm_h%i(m)
-             j = surf_lsm_h%j(m)
-             k = surf_lsm_h%k(m)
-             e_p(k,j,i) = ( surf_lsm_h%us(m) / c_0 )**2
-          ENDDO
-!
-!--       Urban surfaces
-          DO  m = 1, surf_usm_h%ns
-             i = surf_usm_h%i(m)
-             j = surf_usm_h%j(m)
-             k = surf_usm_h%k(m)
-             e_p(k,j,i) = ( surf_usm_h%us(m) / c_0 )**2
-          ENDDO
-!
-!--       Vertical surfaces
-          DO  l = 0, 3
-!
-!--          Default surfaces
-             DO  m = 1, surf_def_v(l)%ns
-                i = surf_def_v(l)%i(m)
-                j = surf_def_v(l)%j(m)
-                k = surf_def_v(l)%k(m)
-                e_p(k,j,i) = ( surf_def_v(l)%us(m) / c_0 )**2
-             ENDDO
-!
-!--          Natural surfaces
-             DO  m = 1, surf_lsm_v(l)%ns
-                i = surf_lsm_v(l)%i(m)
-                j = surf_lsm_v(l)%j(m)
-                k = surf_lsm_v(l)%k(m)
-                e_p(k,j,i) = ( surf_lsm_v(l)%us(m) / c_0 )**2
-             ENDDO
-!
-!--          Urban surfaces
-             DO  m = 1, surf_usm_v(l)%ns
-                i = surf_usm_v(l)%i(m)
-                j = surf_usm_v(l)%j(m)
-                k = surf_usm_v(l)%k(m)
-                e_p(k,j,i) = ( surf_usm_v(l)%us(m) / c_0 )**2
-             ENDDO
-          ENDDO
-       ENDIF
-
-       IF ( .NOT. child_domain )  THEN
-          !$ACC KERNELS PRESENT(e_p)
-          e_p(nzt+1,:,:) = e_p(nzt,:,:)
-          !$ACC END KERNELS
-!
-!--    Nesting case: if parent operates in RANS mode and child in LES mode,
-!--    no TKE is transfered. This case, set Neumann conditions at lateral and 
-!--    top child boundaries. 
-!--    If not ( both either in RANS or in LES mode ), TKE boundary condition
-!--    is treated in the nesting. 
-       ELSE
-
-          IF ( rans_mode_parent  .AND.  .NOT. rans_mode )  THEN
-
-             e_p(nzt+1,:,:) = e_p(nzt,:,:)
-             IF ( bc_dirichlet_l )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
-             IF ( bc_dirichlet_r )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
-             IF ( bc_dirichlet_s )  e_p(:,nys-1,:) = e_p(:,nys,:)
-             IF ( bc_dirichlet_n )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
-
-          ENDIF
-       ENDIF
-    ENDIF
-
-!
-!-- Boundary conditions for TKE dissipation rate. 
-    IF ( rans_tke_e )  THEN
-!
-!--    Use wall function within constant-flux layer
-!--    Upward-facing surfaces
-!--    Default surfaces
-       DO  m = 1, surf_def_h(0)%ns
-          i = surf_def_h(0)%i(m)
-          j = surf_def_h(0)%j(m)
-          k = surf_def_h(0)%k(m)
-          diss_p(k,j,i) = surf_def_h(0)%us(m)**3          &
-                        / ( kappa * surf_def_h(0)%z_mo(m) )
-       ENDDO
-!
-!--    Natural surfaces
-       DO  m = 1, surf_lsm_h%ns
-          i = surf_lsm_h%i(m)
-          j = surf_lsm_h%j(m)
-          k = surf_lsm_h%k(m)
-          diss_p(k,j,i) = surf_lsm_h%us(m)**3          &
-                        / ( kappa * surf_lsm_h%z_mo(m) )
-       ENDDO
-!
-!--    Urban surfaces
-       DO  m = 1, surf_usm_h%ns
-          i = surf_usm_h%i(m)
-          j = surf_usm_h%j(m)
-          k = surf_usm_h%k(m)
-          diss_p(k,j,i) = surf_usm_h%us(m)**3          &
-                        / ( kappa * surf_usm_h%z_mo(m) )
-       ENDDO
-!
-!--    Vertical surfaces
-       DO  l = 0, 3
-!
-!--       Default surfaces
-          DO  m = 1, surf_def_v(l)%ns
-             i = surf_def_v(l)%i(m)
-             j = surf_def_v(l)%j(m)
-             k = surf_def_v(l)%k(m)
-             diss_p(k,j,i) = surf_def_v(l)%us(m)**3          &
-                           / ( kappa * surf_def_v(l)%z_mo(m) )
-          ENDDO
-!
-!--       Natural surfaces
-          DO  m = 1, surf_lsm_v(l)%ns
-             i = surf_lsm_v(l)%i(m)
-             j = surf_lsm_v(l)%j(m)
-             k = surf_lsm_v(l)%k(m)
-             diss_p(k,j,i) = surf_lsm_v(l)%us(m)**3          &
-                           / ( kappa * surf_lsm_v(l)%z_mo(m) )
-          ENDDO
-!
-!--       Urban surfaces
-          DO  m = 1, surf_usm_v(l)%ns
-             i = surf_usm_v(l)%i(m)
-             j = surf_usm_v(l)%j(m)
-             k = surf_usm_v(l)%k(m)
-             diss_p(k,j,i) = surf_usm_v(l)%us(m)**3          &
-                           / ( kappa * surf_usm_v(l)%z_mo(m) )
-          ENDDO
-       ENDDO
-!
-!--    Limit change of diss to be between -90% and +100%. Also, set an absolute
-!--    minimum value
-       DO  i = nxl, nxr
-          DO  j = nys, nyn
-             DO  k = nzb, nzt+1
-                diss_p(k,j,i) = MAX( MIN( diss_p(k,j,i),          &
-                                          2.0_wp * diss(k,j,i) ), &
-                                     0.1_wp * diss(k,j,i),        &
-                                     0.0001_wp )
-             ENDDO
-          ENDDO
-       ENDDO
-
-       IF ( .NOT. child_domain )  THEN
-          diss_p(nzt+1,:,:) = diss_p(nzt,:,:)
-       ENDIF
-    ENDIF
-
 !
 !-- Boundary conditions for salinity
@@ -586,8 +388,4 @@
 !--    given.
        DO  l = 0, 1
-!
-!--       Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--       for downward-facing surfaces at topography bottom (k+1). 
-          kb = MERGE( -1, 1, l == 0 )
           !$OMP PARALLEL DO PRIVATE( i, j, k )
           DO  m = 1, bc_h(l)%ns
@@ -595,5 +393,5 @@
              j = bc_h(l)%j(m)
              k = bc_h(l)%k(m)
-             sa_p(k+kb,j,i) = sa_p(k,j,i)
+             sa_p(k+bc_h(l)%koff,j,i) = sa_p(k,j,i)
           ENDDO
        ENDDO
@@ -620,8 +418,4 @@
 
           DO  l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--          for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              DO  m = 1, bc_h(l)%ns
@@ -629,5 +423,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                q_p(k+kb,j,i) = q(k+kb,j,i)
+                q_p(k+bc_h(l)%koff,j,i) = q(k+bc_h(l)%koff,j,i)
              ENDDO
           ENDDO
@@ -636,8 +430,4 @@
           
           DO  l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--          for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              DO  m = 1, bc_h(l)%ns
@@ -645,5 +435,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                q_p(k+kb,j,i) = q_p(k,j,i)
+                q_p(k+bc_h(l)%koff,j,i) = q_p(k,j,i)
              ENDDO
           ENDDO
@@ -662,12 +452,14 @@
 !--       Run loop over all non-natural and natural walls. Note, in wall-datatype
 !--       the k coordinate belongs to the atmospheric grid point, therefore, set
-!--       qr_p and nr_p at k-1
+!--       qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls
+          DO  l = 0, 1
           !$OMP PARALLEL DO PRIVATE( i, j, k )
-          DO  m = 1, bc_h(0)%ns
-             i = bc_h(0)%i(m)            
-             j = bc_h(0)%j(m)
-             k = bc_h(0)%k(m)
-             qc_p(k-1,j,i) = 0.0_wp
-             nc_p(k-1,j,i) = 0.0_wp
+             DO  m = 1, bc_h(l)%ns
+                i = bc_h(l)%i(m)            
+                j = bc_h(l)%j(m)
+                k = bc_h(l)%k(m)
+                qc_p(k+bc_h(l)%koff,j,i) = 0.0_wp
+                nc_p(k+bc_h(l)%koff,j,i) = 0.0_wp
+             ENDDO
           ENDDO
 !
@@ -683,12 +475,14 @@
 !--       Run loop over all non-natural and natural walls. Note, in wall-datatype
 !--       the k coordinate belongs to the atmospheric grid point, therefore, set
-!--       qr_p and nr_p at k-1
+!--       qr_p and nr_p at upward (k-1) and downward-facing (k+1) walls
+          DO  l = 0, 1
           !$OMP PARALLEL DO PRIVATE( i, j, k )
-          DO  m = 1, bc_h(0)%ns
-             i = bc_h(0)%i(m)            
-             j = bc_h(0)%j(m)
-             k = bc_h(0)%k(m)
-             qr_p(k-1,j,i) = 0.0_wp
-             nr_p(k-1,j,i) = 0.0_wp
+             DO  m = 1, bc_h(l)%ns
+                i = bc_h(l)%i(m)            
+                j = bc_h(l)%j(m)
+                k = bc_h(l)%k(m)
+                qr_p(k+bc_h(l)%koff,j,i) = 0.0_wp
+                nr_p(k+bc_h(l)%koff,j,i) = 0.0_wp
+             ENDDO
           ENDDO
 !
@@ -711,8 +505,4 @@
           
           DO  l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--          for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              DO  m = 1, bc_h(l)%ns
@@ -720,5 +510,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                s_p(k+kb,j,i) = s(k+kb,j,i)
+                s_p(k+bc_h(l)%koff,j,i) = s(k+bc_h(l)%koff,j,i)
              ENDDO
           ENDDO
@@ -727,8 +517,4 @@
           
           DO  l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--          for downward-facing surfaces at topography bottom (k+1). 
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL DO PRIVATE( i, j, k )
              DO  m = 1, bc_h(l)%ns
@@ -736,5 +522,5 @@
                 j = bc_h(l)%j(m)
                 k = bc_h(l)%k(m)
-                s_p(k+kb,j,i) = s_p(k,j,i)
+                s_p(k+bc_h(l)%koff,j,i) = s_p(k,j,i)
              ENDDO
           ENDDO
@@ -750,5 +536,8 @@
        ENDIF
 
-    ENDIF    
+    ENDIF  
+!
+!-- Set boundary conditions for subgrid TKE and dissipation (RANS mode)
+    CALL tcm_boundary_conds
 !
 !-- Top/bottom boundary conditions for chemical species
@@ -760,15 +549,8 @@
 !-- version) these levels are handled as a prognostic level, boundary values
 !-- have to be restored here.
-!-- For the SGS-TKE, Neumann boundary conditions are used at the inflow.
     IF ( bc_dirichlet_s )  THEN
        v_p(:,nys,:) = v_p(:,nys-1,:)
-       IF ( .NOT. constant_diffusion ) e_p(:,nys-1,:) = e_p(:,nys,:)
-    ELSEIF ( bc_dirichlet_n )  THEN
-       IF ( .NOT. constant_diffusion ) e_p(:,nyn+1,:) = e_p(:,nyn,:)
     ELSEIF ( bc_dirichlet_l ) THEN
        u_p(:,:,nxl) = u_p(:,:,nxl-1)
-       IF ( .NOT. constant_diffusion ) e_p(:,:,nxl-1) = e_p(:,:,nxl)
-    ELSEIF ( bc_dirichlet_r )  THEN
-       IF ( .NOT. constant_diffusion ) e_p(:,:,nxr+1) = e_p(:,:,nxr)
     ENDIF
 
@@ -777,4 +559,6 @@
 !-- in case of nest boundaries. This must not be done in case of vertical nesting 
 !-- mode as in that case the lateral boundaries are actually cyclic. 
+!-- Lateral oundary conditions for TKE and dissipation are set 
+!-- in tcm_boundary_conds.
     IF ( nesting_mode /= 'vertical'  .OR.  nesting_offline )  THEN
        IF ( bc_dirichlet_s )  THEN
@@ -787,9 +571,9 @@
 
 !
-!-- Lateral boundary conditions for scalar quantities at the outflow
+!-- Lateral boundary conditions for scalar quantities at the outflow.
+!-- Lateral oundary conditions for TKE and dissipation are set 
+!-- in tcm_boundary_conds.
     IF ( bc_radiation_s )  THEN
        pt_p(:,nys-1,:)     = pt_p(:,nys,:)
-       IF ( .NOT. constant_diffusion )  e_p(:,nys-1,:) = e_p(:,nys,:)
-       IF ( rans_tke_e )  diss_p(:,nys-1,:) = diss_p(:,nys,:)
        IF ( humidity )  THEN
           q_p(:,nys-1,:) = q_p(:,nys,:)
@@ -806,6 +590,4 @@
     ELSEIF ( bc_radiation_n )  THEN
        pt_p(:,nyn+1,:)     = pt_p(:,nyn,:)
-       IF ( .NOT. constant_diffusion )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
-       IF ( rans_tke_e )  diss_p(:,nyn+1,:) = diss_p(:,nyn,:)
        IF ( humidity )  THEN
           q_p(:,nyn+1,:) = q_p(:,nyn,:)
@@ -822,6 +604,4 @@
     ELSEIF ( bc_radiation_l )  THEN
        pt_p(:,:,nxl-1)     = pt_p(:,:,nxl)
-       IF ( .NOT. constant_diffusion )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
-       IF ( rans_tke_e )  diss_p(:,:,nxl-1) = diss_p(:,:,nxl)
        IF ( humidity )  THEN
           q_p(:,:,nxl-1) = q_p(:,:,nxl)
@@ -838,6 +618,4 @@
     ELSEIF ( bc_radiation_r )  THEN
        pt_p(:,:,nxr+1)     = pt_p(:,:,nxr)
-       IF ( .NOT. constant_diffusion )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
-       IF ( rans_tke_e )  diss_p(:,:,nxr+1) = diss_p(:,:,nxr)
        IF ( humidity )  THEN
           q_p(:,:,nxr+1) = q_p(:,:,nxr)
Index: palm/trunk/SOURCE/chemistry_model_mod.f90
===================================================================
--- palm/trunk/SOURCE/chemistry_model_mod.f90	(revision 4101)
+++ palm/trunk/SOURCE/chemistry_model_mod.f90	(revision 4102)
@@ -27,4 +27,8 @@
 ! -----------------
 ! $Id$
+! Slightly revise setting of boundary conditions at horizontal walls, use 
+! data-structure offset index instead of pre-calculate it for each facing
+! 
+! 4080 2019-07-09 18:17:37Z suehring
 ! Restore accidantly removed limitation to positive values
 ! 
@@ -804,5 +808,4 @@
     INTEGER(iwp) ::  j                            !< grid index y direction.
     INTEGER(iwp) ::  k                            !< grid index z direction.
-    INTEGER(iwp) ::  kb                           !< variable to set respective boundary value, depends on facing. 
     INTEGER(iwp) ::  l                            !< running index boundary type, for up- and downward-facing walls.
     INTEGER(iwp) ::  m                            !< running index surface elements.
@@ -839,15 +842,8 @@
        CASE ( 'set_bc_bottomtop' )                   
 !
-!--          Bottom boundary condtions for chemical species      
+!--       Boundary condtions for chemical species at horizontal walls      
           DO  lsp = 1, nspec                                                      
-             IF ( ibc_cs_b == 0 )  THEN                   
+             IF ( ibc_cs_b == 0 )  THEN
                 DO  l = 0, 1 
-!
-!--                Set index kb: For upward-facing surfaces (l=0), kb=-1, i.e.
-!--                the chem_species(nspec)%conc_p value at the topography top (k-1)
-!--                is set; for downward-facing surfaces (l=1), kb=1, i.e. the 
-!--                value at the topography bottom (k+1) is set. 
-
-                   kb = MERGE( -1, 1, l == 0 )
                    !$OMP PARALLEL DO PRIVATE( i, j, k )
                    DO  m = 1, bc_h(l)%ns
@@ -855,5 +851,6 @@
                        j = bc_h(l)%j(m)
                        k = bc_h(l)%k(m)
-                      chem_species(lsp)%conc_p(k+kb,j,i) = chem_species(lsp)%conc(k+kb,j,i) 
+                      chem_species(lsp)%conc_p(k+bc_h(l)%koff,j,i) =           &
+                                      chem_species(lsp)%conc(k+bc_h(l)%koff,j,i) 
                    ENDDO                                        
                 ENDDO                                       
@@ -863,5 +860,4 @@
 !--             in boundary_conds there is som extra loop over m here for passive tracer
                 DO  l = 0, 1
-                   kb = MERGE( -1, 1, l == 0 )
                    !$OMP PARALLEL DO PRIVATE( i, j, k )                                           
                    DO m = 1, bc_h(l)%ns
@@ -869,5 +865,6 @@
                       j = bc_h(l)%j(m)
                       k = bc_h(l)%k(m)
-                      chem_species(lsp)%conc_p(k+kb,j,i) = chem_species(lsp)%conc_p(k,j,i)
+                      chem_species(lsp)%conc_p(k+bc_h(l)%koff,j,i) =           &
+                                         chem_species(lsp)%conc_p(k,j,i)
 
                    ENDDO
Index: palm/trunk/SOURCE/salsa_mod.f90
===================================================================
--- palm/trunk/SOURCE/salsa_mod.f90	(revision 4101)
+++ palm/trunk/SOURCE/salsa_mod.f90	(revision 4102)
@@ -26,4 +26,8 @@
 ! -----------------
 ! $Id$
+! Slightly revise setting of boundary conditions at horizontal walls, use 
+! data-structure offset index instead of pre-calculate it for each facing
+! 
+! 4079 2019-07-09 18:04:41Z suehring
 ! Application of monotonic flux limiter for the vertical scalar advection 
 ! up to the topography top (only for the cache-optimized version at the 
@@ -8017,5 +8021,4 @@
     INTEGER(iwp) ::  j    !< grid index y direction
     INTEGER(iwp) ::  k    !< grid index y direction
-    INTEGER(iwp) ::  kb   !< variable to set respective boundary value, depends on facing.
     INTEGER(iwp) ::  l    !< running index boundary type, for up- and downward-facing walls
     INTEGER(iwp) ::  m    !< running index surface elements
@@ -8030,8 +8033,4 @@
 !--       belongs to the atmospheric grid point, therefore, set s_p at k-1
           DO  l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is 
-!--          set, for downward-facing surfaces at topography bottom (k+1)
-             kb = MERGE ( -1, 1, l == 0 )
              !$OMP PARALLEL PRIVATE( ib, ic, icc, ig, i, j, k )
              !$OMP DO
@@ -8043,13 +8042,16 @@
 
                 DO  ib = 1, nbins_aerosol
-                   aerosol_number(ib)%conc_p(k+kb,j,i) = aerosol_number(ib)%conc(k+kb,j,i)
+                   aerosol_number(ib)%conc_p(k+bc_h(l)%koff,j,i) =             &
+                                    aerosol_number(ib)%conc(k+bc_h(l)%koff,j,i)
                    DO  ic = 1, ncomponents_mass
                       icc = ( ic - 1 ) * nbins_aerosol + ib
-                      aerosol_mass(icc)%conc_p(k+kb,j,i) = aerosol_mass(icc)%conc(k+kb,j,i)
+                      aerosol_mass(icc)%conc_p(k+bc_h(l)%koff,j,i) =           &
+                                    aerosol_mass(icc)%conc(k+bc_h(l)%koff,j,i)
                    ENDDO
                 ENDDO
                 IF ( .NOT. salsa_gases_from_chem )  THEN
                    DO  ig = 1, ngases_salsa
-                      salsa_gas(ig)%conc_p(k+kb,j,i) = salsa_gas(ig)%conc(k+kb,j,i)
+                      salsa_gas(ig)%conc_p(k+bc_h(l)%koff,j,i) =               &
+                                    salsa_gas(ig)%conc(k+bc_h(l)%koff,j,i)
                    ENDDO
                 ENDIF
@@ -8063,8 +8065,4 @@
 
           DO l = 0, 1
-!
-!--          Set kb, for upward-facing surfaces value at topography top (k-1) is 
-!--          set, for downward-facing surfaces at topography bottom (k+1)
-             kb = MERGE( -1, 1, l == 0 )
              !$OMP PARALLEL PRIVATE( ib, ic, icc, ig, i, j, k )
              !$OMP DO
@@ -8076,13 +8074,16 @@
 
                 DO  ib = 1, nbins_aerosol
-                   aerosol_number(ib)%conc_p(k+kb,j,i) = aerosol_number(ib)%conc_p(k,j,i)
+                   aerosol_number(ib)%conc_p(k+bc_h(l)%koff,j,i) =             &
+                                               aerosol_number(ib)%conc_p(k,j,i)
                    DO  ic = 1, ncomponents_mass
                       icc = ( ic - 1 ) * nbins_aerosol + ib
-                      aerosol_mass(icc)%conc_p(k+kb,j,i) = aerosol_mass(icc)%conc_p(k,j,i)
+                      aerosol_mass(icc)%conc_p(k+bc_h(l)%koff,j,i) =           &
+                                               aerosol_mass(icc)%conc_p(k,j,i)
                    ENDDO
                 ENDDO
                 IF ( .NOT. salsa_gases_from_chem ) THEN
                    DO  ig = 1, ngases_salsa
-                      salsa_gas(ig)%conc_p(k+kb,j,i) = salsa_gas(ig)%conc_p(k,j,i)
+                      salsa_gas(ig)%conc_p(k+bc_h(l)%koff,j,i) =               &
+                                               salsa_gas(ig)%conc_p(k,j,i)
                    ENDDO
                 ENDIF
Index: palm/trunk/SOURCE/surface_mod.f90
===================================================================
--- palm/trunk/SOURCE/surface_mod.f90	(revision 4101)
+++ palm/trunk/SOURCE/surface_mod.f90	(revision 4102)
@@ -26,4 +26,8 @@
 ! -----------------
 ! $Id$
+! - Revise initialization of the boundary data structure
+! - Add new data structure to set boundary conditions at vertical walls
+! 
+! 3943 2019-05-02 09:50:41Z maronga
 ! Removed qsws_eb as it is no longer needed.
 ! 
@@ -306,10 +310,12 @@
 !-- are applied 
     TYPE bc_type
-
-       INTEGER(iwp) :: ns                                  !< number of surface elements on the PE
-
-       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  i       !< x-index linking to the PALM 3D-grid  
-       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  j       !< y-index linking to the PALM 3D-grid    
-       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  k       !< z-index linking to the PALM 3D-grid   
+       INTEGER(iwp) :: ioff !< offset value in x-direction, used to determine index of surface element
+       INTEGER(iwp) :: joff !< offset value in y-direction, used to determine index of surface element
+       INTEGER(iwp) :: koff !< offset value in z-direction, used to determine index of surface element
+       INTEGER(iwp) :: ns   !< number of surface elements on the PE
+
+       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  i !< x-index linking to the PALM 3D-grid  
+       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  j !< y-index linking to the PALM 3D-grid    
+       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  k !< z-index linking to the PALM 3D-grid   
 
        INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: start_index !< start index within surface data type for given (j,i) 
@@ -600,4 +606,5 @@
 
     TYPE (bc_type), DIMENSION(0:1)           ::  bc_h        !< boundary condition data type, horizontal upward- and downward facing surfaces
+    TYPE (bc_type), DIMENSION(0:3)           ::  bc_v        !< boundary condition data type, vertical surfaces
 
     TYPE (surf_type), DIMENSION(0:2), TARGET ::  surf_def_h  !< horizontal default surfaces (Up, Down, and Top)
@@ -668,6 +675,6 @@
 !
 !-- Public variables
-    PUBLIC bc_h, ind_pav_green, ind_veg_wall, ind_wat_win, ns_h_on_file, ns_v_on_file, surf_def_h, &
-           surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, surf_usm_v, surf_type,                  &
+    PUBLIC bc_h, bc_v, ind_pav_green, ind_veg_wall, ind_wat_win, ns_h_on_file, ns_v_on_file,       &
+           surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, surf_usm_v, surf_type,      &
            vertical_surfaces_exist, surf_bulk_cloud_model, surf_microphysics_morrison,             &
            surf_microphysics_seifert
@@ -687,5 +694,6 @@
 ! Description:
 ! ------------
-!> Initialize data type for setting boundary conditions at horizontal surfaces. 
+!> Initialize data type for setting boundary conditions at horizontal and  
+!> vertical surfaces. 
 !------------------------------------------------------------------------------!
     SUBROUTINE init_bc
@@ -696,92 +704,175 @@
        INTEGER(iwp) ::  j         !< loop index along y-direction
        INTEGER(iwp) ::  k         !< loop index along y-direction
-
-       INTEGER(iwp), DIMENSION(0:1) ::  num_h         !< number of surfaces
-       INTEGER(iwp), DIMENSION(0:1) ::  num_h_kji     !< number of surfaces for (j,i)-grid point
-       INTEGER(iwp), DIMENSION(0:1) ::  start_index_h !< local start index
-
-!
-!--    First of all, count the number of upward- and downward-facing surfaces
-       num_h = 0
-       DO  i = nxlg, nxrg
-          DO  j = nysg, nyng
-             DO  k = nzb+1, nzt
-!
-!--             Check if current gridpoint belongs to the atmosphere
-                IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
-!
-!--                Upward-facing
-                   IF ( .NOT. BTEST( wall_flags_0(k-1,j,i), 0 ) )              &
-                      num_h(0) = num_h(0) + 1
-!
-!--                Downward-facing
-                   IF ( .NOT. BTEST( wall_flags_0(k+1,j,i), 0 ) )              &
-                      num_h(1) = num_h(1) + 1
-                ENDIF
+       INTEGER(iwp) ::  l         !< running index for differently aligned surfaces
+
+       INTEGER(iwp), DIMENSION(0:1) ::  num_h         !< number of horizontal surfaces on subdomain
+       INTEGER(iwp), DIMENSION(0:1) ::  num_h_kji     !< number of horizontal surfaces at (j,i)-grid point
+       INTEGER(iwp), DIMENSION(0:1) ::  start_index_h !< local start index of horizontal surface elements
+       
+       INTEGER(iwp), DIMENSION(0:3) ::  num_v         !< number of vertical surfaces on subdomain
+       INTEGER(iwp), DIMENSION(0:3) ::  num_v_kji     !< number of vertical surfaces at (j,i)-grid point
+       INTEGER(iwp), DIMENSION(0:3) ::  start_index_v !< local start index of vertical surface elements
+!
+!--    Set offset indices, i.e. index difference between surface element and 
+!--    surface-bounded grid point.
+!--    Horizontal surfaces - no horizontal offsets
+       bc_h(:)%ioff = 0
+       bc_h(:)%joff = 0
+!
+!--    Horizontal surfaces, upward facing (0) and downward facing (1)
+       bc_h(0)%koff   = -1
+       bc_h(1)%koff   = 1
+!
+!--    Vertical surfaces - no vertical offset
+       bc_v(0:3)%koff = 0
+!
+!--    North- and southward facing - no offset in x
+       bc_v(0:1)%ioff = 0
+!
+!--    Northward facing offset in y
+       bc_v(0)%joff = -1
+!
+!--    Southward facing offset in y
+       bc_v(1)%joff = 1
+!
+!--    East- and westward facing - no offset in y
+       bc_v(2:3)%joff = 0
+!
+!--    Eastward facing offset in x
+       bc_v(2)%ioff = -1
+!
+!--    Westward facing offset in y
+       bc_v(3)%ioff = 1
+!
+!--    Initialize data structure for horizontal surfaces, i.e. count the number
+!--    of surface elements, allocate and initialize the respective index arrays, 
+!--    and set the respective start and end indices at each (j,i)-location.
+       DO  l = 0, 1
+!
+!--       Count the number of upward- and downward-facing surfaces on subdomain
+          num_h(l) = 0
+          DO  i = nxlg, nxrg
+             DO  j = nysg, nyng
+                DO  k = nzb+1, nzt
+!         
+!--                Check if current gridpoint belongs to the atmosphere
+                   IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
+                      IF ( .NOT. BTEST( wall_flags_0(k+bc_h(l)%koff,           &
+                                                     j+bc_h(l)%joff,           &
+                                                     i+bc_h(l)%ioff), 0 ) )    &
+                         num_h(l) = num_h(l) + 1
+                   ENDIF
+                ENDDO
+             ENDDO
+          ENDDO
+!         
+!--       Save the number of horizontal surface elements
+          bc_h(l)%ns = num_h(l)
+!
+!--       ALLOCATE arrays for horizontal surfaces
+          ALLOCATE( bc_h(l)%i(1:bc_h(l)%ns) )
+          ALLOCATE( bc_h(l)%j(1:bc_h(l)%ns) )
+          ALLOCATE( bc_h(l)%k(1:bc_h(l)%ns) )
+          ALLOCATE( bc_h(l)%start_index(nysg:nyng,nxlg:nxrg) )
+          ALLOCATE( bc_h(l)%end_index(nysg:nyng,nxlg:nxrg)   )
+          bc_h(l)%start_index = 1
+          bc_h(l)%end_index   = 0
+          
+          num_h(l)         = 1
+          start_index_h(l) = 1
+          DO  i = nxlg, nxrg
+             DO  j = nysg, nyng
+          
+                num_h_kji(l) = 0
+                DO  k = nzb+1, nzt
+!         
+!--                Check if current gridpoint belongs to the atmosphere
+                   IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
+!         
+!--                   Upward-facing
+                      IF ( .NOT. BTEST( wall_flags_0(k+bc_h(l)%koff,           &
+                                                     j+bc_h(l)%joff,           &
+                                                     i+bc_h(l)%ioff), 0 )      &
+                         )  THEN
+                         bc_h(l)%i(num_h(l)) = i
+                         bc_h(l)%j(num_h(l)) = j
+                         bc_h(l)%k(num_h(l)) = k
+                         num_h_kji(l)        = num_h_kji(l) + 1
+                         num_h(l)            = num_h(l) + 1
+                      ENDIF
+                   ENDIF
+                ENDDO
+                bc_h(l)%start_index(j,i) = start_index_h(l)
+                bc_h(l)%end_index(j,i)   = bc_h(l)%start_index(j,i) +          &
+                                           num_h_kji(l) - 1
+                start_index_h(l)         = bc_h(l)%end_index(j,i) + 1
              ENDDO
           ENDDO
        ENDDO
-!
-!--    Save the number of surface elements
-       bc_h(0)%ns = num_h(0)
-       bc_h(1)%ns = num_h(1)
-!
-!--    ALLOCATE data type variables
-!--    Upward facing
-       ALLOCATE( bc_h(0)%i(1:bc_h(0)%ns) )
-       ALLOCATE( bc_h(0)%j(1:bc_h(0)%ns) )
-       ALLOCATE( bc_h(0)%k(1:bc_h(0)%ns) )
-       ALLOCATE( bc_h(0)%start_index(nysg:nyng,nxlg:nxrg) )
-       ALLOCATE( bc_h(0)%end_index(nysg:nyng,nxlg:nxrg)   )
-       bc_h(0)%start_index = 1
-       bc_h(0)%end_index   = 0
-!
-!--    Downward facing
-       ALLOCATE( bc_h(1)%i(1:bc_h(1)%ns) )
-       ALLOCATE( bc_h(1)%j(1:bc_h(1)%ns) )
-       ALLOCATE( bc_h(1)%k(1:bc_h(1)%ns) )
-       ALLOCATE( bc_h(1)%start_index(nysg:nyng,nxlg:nxrg) )
-       ALLOCATE( bc_h(1)%end_index(nysg:nyng,nxlg:nxrg)   )
-       bc_h(1)%start_index = 1
-       bc_h(1)%end_index   = 0
-!
-!--    Store the respective indices on data type
-       num_h(0:1)         = 1
-       start_index_h(0:1) = 1
-       DO  i = nxlg, nxrg
-          DO  j = nysg, nyng
-
-             num_h_kji(0:1) = 0
-             DO  k = nzb+1, nzt
-!
-!--             Check if current gridpoint belongs to the atmosphere
-                IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
-!
-!--                Upward-facing
-                   IF ( .NOT. BTEST( wall_flags_0(k-1,j,i), 0 ) )  THEN
-                      bc_h(0)%i(num_h(0)) = i
-                      bc_h(0)%j(num_h(0)) = j
-                      bc_h(0)%k(num_h(0)) = k
-                      num_h_kji(0)        = num_h_kji(0) + 1
-                      num_h(0)            = num_h(0) + 1
+
+!
+!--    Initialize data structure for vertical surfaces, i.e. count the number
+!--    of surface elements, allocate and initialize the respective index arrays, 
+!--    and set the respective start and end indices at each (j,i)-location.
+       DO  l = 0, 3
+!
+!--       Count the number of upward- and downward-facing surfaces on subdomain
+          num_v(l) = 0
+          DO  i = nxlg, nxrg
+             DO  j = nysg, nyng
+                DO  k = nzb+1, nzt
+!         
+!--                Check if current gridpoint belongs to the atmosphere
+                   IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
+                      IF ( .NOT. BTEST( wall_flags_0(k+bc_v(l)%koff,           &
+                                                     j+bc_v(l)%joff,           &
+                                                     i+bc_v(l)%ioff), 0 ) )    &
+                         num_v(l) = num_v(l) + 1
                    ENDIF
-!
-!--                Downward-facing
-                   IF ( .NOT. BTEST( wall_flags_0(k+1,j,i), 0 ) )  THEN
-                      bc_h(1)%i(num_h(1)) = i
-                      bc_h(1)%j(num_h(1)) = j
-                      bc_h(1)%k(num_h(1)) = k
-                      num_h_kji(1)        = num_h_kji(1) + 1
-                      num_h(1)            = num_h(1) + 1
+                ENDDO
+             ENDDO
+          ENDDO
+!         
+!--       Save the number of horizontal surface elements
+          bc_v(l)%ns = num_v(l)
+!
+!--       ALLOCATE arrays for horizontal surfaces
+          ALLOCATE( bc_v(l)%i(1:bc_v(l)%ns) )
+          ALLOCATE( bc_v(l)%j(1:bc_v(l)%ns) )
+          ALLOCATE( bc_v(l)%k(1:bc_v(l)%ns) )
+          ALLOCATE( bc_v(l)%start_index(nysg:nyng,nxlg:nxrg) )
+          ALLOCATE( bc_v(l)%end_index(nysg:nyng,nxlg:nxrg)   )
+          bc_v(l)%start_index = 1
+          bc_v(l)%end_index   = 0
+          
+          num_v(l)         = 1
+          start_index_v(l) = 1
+          DO  i = nxlg, nxrg
+             DO  j = nysg, nyng
+          
+                num_v_kji(l) = 0
+                DO  k = nzb+1, nzt
+!         
+!--                Check if current gridpoint belongs to the atmosphere
+                   IF ( BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
+!         
+!--                   Upward-facing
+                      IF ( .NOT. BTEST( wall_flags_0(k+bc_v(l)%koff,           &
+                                                     j+bc_v(l)%joff,           &
+                                                     i+bc_v(l)%ioff), 0 )      &
+                         )  THEN
+                         bc_v(l)%i(num_v(l)) = i
+                         bc_v(l)%j(num_v(l)) = j
+                         bc_v(l)%k(num_v(l)) = k
+                         num_v_kji(l)        = num_v_kji(l) + 1
+                         num_v(l)            = num_v(l) + 1
+                      ENDIF
                    ENDIF
-                ENDIF
+                ENDDO
+                bc_v(l)%start_index(j,i) = start_index_v(l)
+                bc_v(l)%end_index(j,i)   = bc_v(l)%start_index(j,i) +          &
+                                           num_v_kji(l) - 1
+                start_index_v(l)         = bc_v(l)%end_index(j,i) + 1
              ENDDO
-             bc_h(0)%start_index(j,i) = start_index_h(0)
-             bc_h(0)%end_index(j,i)   = bc_h(0)%start_index(j,i) + num_h_kji(0) - 1
-             start_index_h(0)         = bc_h(0)%end_index(j,i) + 1
-
-             bc_h(1)%start_index(j,i) = start_index_h(1)
-             bc_h(1)%end_index(j,i)   = bc_h(1)%start_index(j,i) + num_h_kji(1) - 1
-             start_index_h(1)         = bc_h(1)%end_index(j,i) + 1
           ENDDO
        ENDDO
Index: palm/trunk/SOURCE/turbulence_closure_mod.f90
===================================================================
--- palm/trunk/SOURCE/turbulence_closure_mod.f90	(revision 4101)
+++ palm/trunk/SOURCE/turbulence_closure_mod.f90	(revision 4102)
@@ -25,4 +25,13 @@
 ! -----------------
 ! $Id$
+! - Modularize setting of boundary conditions for TKE and dissipation
+! - Neumann boundary condition for TKE at model top is set also in child domain
+! - Revise setting of Neumann boundary conditions at non-cyclic lateral 
+!   boundaries
+! - Bugfix, set Neumann boundary condition for TKE at vertical wall instead of 
+!   an implicit Dirichlet boundary condition which implied a sink of TKE
+!   at vertical walls
+! 
+! 4048 2019-06-21 21:00:21Z knoop
 ! write out preprocessor directives; remove tailing whitespaces
 !
@@ -192,58 +201,77 @@
 !> @todo Check for random disturbances
 !> @note <Enter notes on the module>
-!------------------------------------------------------------------------------!
+!-----------------------------------------------------------------------------!
  MODULE turbulence_closure_mod
 
 
-    USE arrays_3d,                                                             &
-        ONLY:  diss, diss_1, diss_2, diss_3, diss_p, dzu, e, e_1, e_2, e_3,    &
-               e_p, kh, km, mean_inflow_profiles, prho, pt, tdiss_m,           &
+    USE arrays_3d,                                                            &
+        ONLY:  diss, diss_1, diss_2, diss_3, diss_p, dzu, e, e_1, e_2, e_3,   &
+               e_p, kh, km, mean_inflow_profiles, prho, pt, tdiss_m,          &
                te_m, tend, u, v, vpt, w
 
-    USE basic_constants_and_equations_mod,                                     &
+    USE basic_constants_and_equations_mod,                                    &
         ONLY:  g, kappa, lv_d_cp, lv_d_rd, rd_d_rv
 
-    USE control_parameters,                                                    &
-        ONLY:  constant_diffusion, dt_3d, e_init, humidity,                    &
-               initializing_actions, intermediate_timestep_count,              &
-               intermediate_timestep_count_max, km_constant,                   &
-               les_dynamic, les_mw, ocean_mode, plant_canopy, prandtl_number,  &
-               pt_reference, rans_mode, rans_tke_e, rans_tke_l,                &
-               timestep_scheme, turbulence_closure,                            &
-               turbulent_inflow, use_upstream_for_tke, vpt_reference,          &
+    USE control_parameters,                                                   &
+        ONLY:  bc_dirichlet_l,                                                &
+               bc_dirichlet_n,                                                &
+               bc_dirichlet_r,                                                &
+               bc_dirichlet_s,                                                &
+               bc_radiation_l,                                                &
+               bc_radiation_n,                                                &
+               bc_radiation_r,                                                &
+               bc_radiation_s,                                                &
+               child_domain,                                                  &
+               constant_diffusion, dt_3d, e_init, humidity,                   &
+               initializing_actions, intermediate_timestep_count,             &
+               intermediate_timestep_count_max, km_constant,                  &
+               les_dynamic, les_mw, ocean_mode, plant_canopy, prandtl_number, &
+               pt_reference, rans_mode, rans_tke_e, rans_tke_l,               &
+               timestep_scheme, turbulence_closure,                           &
+               turbulent_inflow, use_upstream_for_tke, vpt_reference,         &
                ws_scheme_sca, current_timestep_number
 
-    USE advec_ws,                                                              &
+    USE advec_ws,                                                             &
         ONLY:  advec_s_ws
 
-    USE advec_s_bc_mod,                                                        &
+    USE advec_s_bc_mod,                                                       &
         ONLY:  advec_s_bc
 
-    USE advec_s_pw_mod,                                                        &
+    USE advec_s_pw_mod,                                                       &
         ONLY:  advec_s_pw
 
-    USE advec_s_up_mod,                                                        &
+    USE advec_s_up_mod,                                                       &
         ONLY:  advec_s_up
 
-    USE cpulog,                                                                &
+    USE cpulog,                                                               &
         ONLY:  cpu_log, log_point_s
 
-    USE indices,                                                               &
-        ONLY:  nbgp, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,     &
+    USE indices,                                                              &
+        ONLY:  nbgp, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,    &
                wall_flags_0
 
     USE kinds
 
-    USE ocean_mod,                                                             &
+    USE ocean_mod,                                                            &
         ONLY:  prho_reference
 
     USE pegrid
 
-    USE plant_canopy_model_mod,                                                &
+    USE plant_canopy_model_mod,                                               &
         ONLY:  pcm_tendency
 
-    USE statistics,                                                            &
+    USE statistics,                                                           &
         ONLY:  hom, hom_sum, statistic_regions
-
+        
+    USE surface_mod,                                                          &
+        ONLY:  bc_h,                                                          &
+               bc_v,                                                          &
+               get_topography_top_index_ji,                                   &
+               surf_def_h,                                                    &
+               surf_def_v,                                                    &
+               surf_lsm_h,                                                    &
+               surf_lsm_v,                                                    &
+               surf_usm_h,                                                    &
+               surf_usm_v
 
     IMPLICIT NONE
@@ -271,5 +299,6 @@
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  l_wall !< near-wall mixing length
 
-
+!
+!-- Public variables
     PUBLIC c_0, rans_const_c, rans_const_sigma
 
@@ -277,22 +306,28 @@
 
     PRIVATE
-
-
-    PUBLIC &
-       tcm_check_parameters, &
-       tcm_check_data_output, &
-       tcm_define_netcdf_grid, &
-       tcm_init_arrays, &
-       tcm_init, &
-       tcm_actions, &
-       tcm_prognostic_equations, &
-       tcm_swap_timelevel, &
-       tcm_3d_data_averaging, &
-       tcm_data_output_2d, &
-       tcm_data_output_3d, &
+!
+!-- Public subroutines
+    PUBLIC                                                                     &
+       tcm_boundary_conds,                                                     &
+       tcm_check_parameters,                                                   &
+       tcm_check_data_output,                                                  &
+       tcm_define_netcdf_grid,                                                 &
+       tcm_init_arrays,                                                        &
+       tcm_init,                                                               &
+       tcm_actions,                                                            &
+       tcm_prognostic_equations,                                               &
+       tcm_swap_timelevel,                                                     &
+       tcm_3d_data_averaging,                                                  &
+       tcm_data_output_2d,                                                     &
+       tcm_data_output_3d,                                                     &
        tcm_diffusivities
 
 !
 !-- PALM interfaces:
+!-- Boundary conditions for subgrid TKE and dissipation
+    INTERFACE tcm_boundary_conds
+       MODULE PROCEDURE tcm_boundary_conds
+    END INTERFACE tcm_boundary_conds
+!
 !-- Input parameter checks to be done in check_parameters
     INTERFACE tcm_check_parameters
@@ -371,4 +406,249 @@
  CONTAINS
 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Check parameters routine for turbulence closure module.
+!------------------------------------------------------------------------------!
+ SUBROUTINE tcm_boundary_conds
+
+    USE pmc_interface,                                                         &
+        ONLY : rans_mode_parent
+ 
+    IMPLICIT NONE
+
+    INTEGER(iwp) ::  i  !< grid index x direction
+    INTEGER(iwp) ::  j  !< grid index y direction
+    INTEGER(iwp) ::  k  !< grid index z direction
+    INTEGER(iwp) ::  l  !< running index boundary type, for up- and downward-facing walls
+    INTEGER(iwp) ::  m  !< running index surface elements
+!
+!-- Boundary conditions for TKE.
+    IF ( .NOT. constant_diffusion )  THEN
+!
+!--    In LES mode, Neumann conditions with de/x_i=0 are assumed at solid walls.
+!--    Note, only TKE is prognostic in this case and dissipation is only 
+!--    a diagnostic quantity. 
+       IF ( .NOT. rans_mode )  THEN
+!
+!--       Horizontal walls, upward- and downward-facing
+          DO  l = 0, 1
+             !$OMP PARALLEL DO PRIVATE( i, j, k )
+             !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+             !$ACC PRESENT(bc_h, e_p)
+             DO  m = 1, bc_h(l)%ns
+                i = bc_h(l)%i(m)            
+                j = bc_h(l)%j(m)
+                k = bc_h(l)%k(m)
+                e_p(k+bc_h(l)%koff,j,i) = e_p(k,j,i)
+             ENDDO
+          ENDDO
+!
+!--       Vertical walls
+          DO  l = 0, 3
+             !$OMP PARALLEL DO PRIVATE( i, j, k )
+             !$ACC PARALLEL LOOP PRIVATE(i, j, k) &
+             !$ACC PRESENT(bc_v, e_p)
+             DO  m = 1, bc_v(l)%ns
+                i = bc_v(l)%i(m)            
+                j = bc_v(l)%j(m)
+                k = bc_v(l)%k(m)
+                e_p(k,j+bc_v(l)%joff,i+bc_v(l)%ioff) = e_p(k,j,i)
+             ENDDO
+          ENDDO
+!
+!--    In RANS mode, wall function is used as boundary condition for TKE
+       ELSE
+!
+!--       Use wall function within constant-flux layer
+!--       Note, grid points listed in bc_h are not included in any calculations in RANS mode and
+!--       are therefore not set here.
+!
+!--       Upward-facing surfaces
+!--       Default surfaces
+          DO  m = 1, surf_def_h(0)%ns
+             i = surf_def_h(0)%i(m)
+             j = surf_def_h(0)%j(m)
+             k = surf_def_h(0)%k(m)
+             e_p(k,j,i) = ( surf_def_h(0)%us(m) / c_0 )**2
+          ENDDO
+!
+!--       Natural surfaces
+          DO  m = 1, surf_lsm_h%ns
+             i = surf_lsm_h%i(m)
+             j = surf_lsm_h%j(m)
+             k = surf_lsm_h%k(m)
+             e_p(k,j,i) = ( surf_lsm_h%us(m) / c_0 )**2
+          ENDDO
+!
+!--       Urban surfaces
+          DO  m = 1, surf_usm_h%ns
+             i = surf_usm_h%i(m)
+             j = surf_usm_h%j(m)
+             k = surf_usm_h%k(m)
+             e_p(k,j,i) = ( surf_usm_h%us(m) / c_0 )**2
+          ENDDO
+!
+!--       Vertical surfaces
+          DO  l = 0, 3
+!
+!--          Default surfaces
+             DO  m = 1, surf_def_v(l)%ns
+                i = surf_def_v(l)%i(m)
+                j = surf_def_v(l)%j(m)
+                k = surf_def_v(l)%k(m)
+                e_p(k,j,i) = ( surf_def_v(l)%us(m) / c_0 )**2
+             ENDDO
+!
+!--          Natural surfaces
+             DO  m = 1, surf_lsm_v(l)%ns
+                i = surf_lsm_v(l)%i(m)
+                j = surf_lsm_v(l)%j(m)
+                k = surf_lsm_v(l)%k(m)
+                e_p(k,j,i) = ( surf_lsm_v(l)%us(m) / c_0 )**2
+             ENDDO
+!
+!--          Urban surfaces
+             DO  m = 1, surf_usm_v(l)%ns
+                i = surf_usm_v(l)%i(m)
+                j = surf_usm_v(l)%j(m)
+                k = surf_usm_v(l)%k(m)
+                e_p(k,j,i) = ( surf_usm_v(l)%us(m) / c_0 )**2
+             ENDDO
+          ENDDO
+       ENDIF
+!
+!--    Set Neumann boundary condition for TKE at model top. Do this also
+!--    in case of a nested run.
+       !$ACC KERNELS PRESENT(e_p)
+       e_p(nzt+1,:,:) = e_p(nzt,:,:)
+       !$ACC END KERNELS
+!
+!--    Nesting case: if parent operates in RANS mode and child in LES mode,
+!--    no TKE is transfered. This case, set Neumann conditions at lateral and 
+!--    top child boundaries. 
+!--    If not ( both either in RANS or in LES mode ), TKE boundary condition
+!--    is treated in the nesting. 
+       If ( child_domain )  THEN
+          IF ( rans_mode_parent  .AND.  .NOT. rans_mode )  THEN
+
+             e_p(nzt+1,:,:) = e_p(nzt,:,:)
+             IF ( bc_dirichlet_l )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
+             IF ( bc_dirichlet_r )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
+             IF ( bc_dirichlet_s )  e_p(:,nys-1,:) = e_p(:,nys,:)
+             IF ( bc_dirichlet_n )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
+
+          ENDIF
+       ENDIF
+!
+!--    At in- and outflow boundaries also set Neumann boundary conditions 
+!--    for the SGS-TKE. An exception is made for the child domain if 
+!--    both parent and child operate in RANS mode. This case no 
+!--    lateral Neumann boundary conditions will be set but Dirichlet
+!--    conditions will be set in the nesting.
+       IF ( .NOT. child_domain  .AND.  .NOT. rans_mode_parent  .AND.           &
+            .NOT. rans_mode )  THEN
+          IF ( bc_dirichlet_s  .OR.  bc_radiation_s )  THEN
+             e_p(:,nys-1,:) = e_p(:,nys,:)
+             IF ( rans_tke_e )  diss_p(:,nys-1,:) = diss_p(:,nys,:)
+          ENDIF
+          IF ( bc_dirichlet_n  .OR.  bc_radiation_n )  THEN
+             e_p(:,nyn+1,:) = e_p(:,nyn,:)
+             IF ( rans_tke_e )  diss_p(:,nyn+1,:) = diss_p(:,nyn,:)  
+          ENDIF
+          IF ( bc_dirichlet_l  .OR.  bc_radiation_l )  THEN
+             e_p(:,:,nxl-1) = e_p(:,:,nxl)
+             IF ( rans_tke_e )  diss_p(:,nyn+1,:) = diss_p(:,nyn,:)  
+          ENDIF
+          IF ( bc_dirichlet_r  .OR.  bc_radiation_r )  THEN
+             e_p(:,:,nxr+1) = e_p(:,:,nxr)
+             IF ( rans_tke_e )  diss_p(:,nyn+1,:) = diss_p(:,nyn,:) 
+          ENDIF
+       ENDIF
+    ENDIF
+
+!
+!-- Boundary conditions for TKE dissipation rate in RANS mode. 
+    IF ( rans_tke_e )  THEN
+!
+!--    Use wall function within constant-flux layer
+!--    Upward-facing surfaces
+!--    Default surfaces
+       DO  m = 1, surf_def_h(0)%ns
+          i = surf_def_h(0)%i(m)
+          j = surf_def_h(0)%j(m)
+          k = surf_def_h(0)%k(m)
+          diss_p(k,j,i) = surf_def_h(0)%us(m)**3          &
+                        / ( kappa * surf_def_h(0)%z_mo(m) )
+       ENDDO
+!
+!--    Natural surfaces
+       DO  m = 1, surf_lsm_h%ns
+          i = surf_lsm_h%i(m)
+          j = surf_lsm_h%j(m)
+          k = surf_lsm_h%k(m)
+          diss_p(k,j,i) = surf_lsm_h%us(m)**3          &
+                        / ( kappa * surf_lsm_h%z_mo(m) )
+       ENDDO
+!
+!--    Urban surfaces
+       DO  m = 1, surf_usm_h%ns
+          i = surf_usm_h%i(m)
+          j = surf_usm_h%j(m)
+          k = surf_usm_h%k(m)
+          diss_p(k,j,i) = surf_usm_h%us(m)**3          &
+                        / ( kappa * surf_usm_h%z_mo(m) )
+       ENDDO
+!
+!--    Vertical surfaces
+       DO  l = 0, 3
+!
+!--       Default surfaces
+          DO  m = 1, surf_def_v(l)%ns
+             i = surf_def_v(l)%i(m)
+             j = surf_def_v(l)%j(m)
+             k = surf_def_v(l)%k(m)
+             diss_p(k,j,i) = surf_def_v(l)%us(m)**3          &
+                           / ( kappa * surf_def_v(l)%z_mo(m) )
+          ENDDO
+!
+!--       Natural surfaces
+          DO  m = 1, surf_lsm_v(l)%ns
+             i = surf_lsm_v(l)%i(m)
+             j = surf_lsm_v(l)%j(m)
+             k = surf_lsm_v(l)%k(m)
+             diss_p(k,j,i) = surf_lsm_v(l)%us(m)**3          &
+                           / ( kappa * surf_lsm_v(l)%z_mo(m) )
+          ENDDO
+!
+!--       Urban surfaces
+          DO  m = 1, surf_usm_v(l)%ns
+             i = surf_usm_v(l)%i(m)
+             j = surf_usm_v(l)%j(m)
+             k = surf_usm_v(l)%k(m)
+             diss_p(k,j,i) = surf_usm_v(l)%us(m)**3          &
+                           / ( kappa * surf_usm_v(l)%z_mo(m) )
+          ENDDO
+       ENDDO
+!
+!--    Limit change of diss to be between -90% and +100%. Also, set an absolute
+!--    minimum value
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             DO  k = nzb, nzt+1
+                diss_p(k,j,i) = MAX( MIN( diss_p(k,j,i),          &
+                                          2.0_wp * diss(k,j,i) ), &
+                                     0.1_wp * diss(k,j,i),        &
+                                     0.0001_wp )
+             ENDDO
+          ENDDO
+       ENDDO
+
+       
+       diss_p(nzt+1,:,:) = diss_p(nzt,:,:)
+     
+    ENDIF
+ END SUBROUTINE tcm_boundary_conds 
+ 
 !------------------------------------------------------------------------------!
 ! Description:
@@ -926,7 +1206,4 @@
     USE model_1d_mod,                                                          &
         ONLY:  e1d, kh1d, km1d
-
-    USE surface_mod,                                                           &
-        ONLY:  get_topography_top_index_ji
 
     IMPLICIT NONE
@@ -1713,7 +1990,4 @@
         ONLY:  phi_m
 
-    USE surface_mod,                                                           &
-        ONLY :  surf_def_h, surf_lsm_h, surf_usm_h
-
     IMPLICIT NONE
 
@@ -2418,8 +2692,4 @@
     USE bulk_cloud_model_mod,                                                  &
         ONLY:  bulk_cloud_model
-
-    USE surface_mod,                                                           &
-        ONLY :  surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h,    &
-                surf_usm_v
 
     IMPLICIT NONE
@@ -3122,8 +3392,4 @@
         ONLY:  bulk_cloud_model
 
-    USE surface_mod,                                                           &
-        ONLY :  surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h,    &
-                surf_usm_v
-
     IMPLICIT NONE
 
@@ -3767,7 +4033,4 @@
         ONLY:  use_sgs_for_particles, wang_kernel
 
-    USE surface_mod,                                                           &
-       ONLY :  bc_h
-
     IMPLICIT NONE
 
@@ -3929,5 +4192,7 @@
 
 !
-!-- Neumann boundary condition for dissipation diss(nzb,:,:) = diss(nzb+1,:,:)
+!-- Neumann boundary condition for dissipation diss(nzb,:,:) = diss(nzb+1,:,:).
+!-- Note, bc cannot be set in tcm_boundary conditions as the dissipation 
+!-- in LES mode is only a diagnostic quantity. 
     IF ( .NOT. rans_tke_e .AND. ( use_sgs_for_particles  .OR.                  &
          wang_kernel  .OR.  collision_turbulence  ) )  THEN
@@ -3973,7 +4238,4 @@
     USE particle_attributes,                                                   &
         ONLY:  use_sgs_for_particles, wang_kernel
-
-    USE surface_mod,                                                           &
-       ONLY :  bc_h
 
     IMPLICIT NONE
@@ -4062,4 +4324,6 @@
 !--    For each surface type determine start and end index (in case of elevated
 !--    topography several up/downward facing surfaces may exist.
+!--    Note, bc cannot be set in tcm_boundary conditions as the dissipation 
+!--    in LES mode is only a diagnostic quantity. 
        surf_s = bc_h(0)%start_index(j,i)
        surf_e = bc_h(0)%end_index(j,i)
@@ -4322,8 +4586,4 @@
         ONLY:  phi_m
 
-    USE surface_mod,                                                           &
-        ONLY :  bc_h, surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v,          &
-                surf_usm_h, surf_usm_v
-
     INTEGER(iwp) ::  i          !< loop index
     INTEGER(iwp) ::  j          !< loop index
Index: palm/trunk/SOURCE/vertical_nesting_mod.f90
===================================================================
--- palm/trunk/SOURCE/vertical_nesting_mod.f90	(revision 4101)
+++ palm/trunk/SOURCE/vertical_nesting_mod.f90	(revision 4102)
@@ -26,4 +26,8 @@
 ! -----------------
 ! $Id$
+! - Slightly revise setting of boundary conditions at horizontal walls, use 
+!   data-structure offset index instead of pre-calculate it for each facing
+! 
+! 4101 2019-07-17 15:14:26Z gronemeier
 ! remove old_dt
 ! 
@@ -2711,6 +2715,5 @@
         INTEGER(iwp)                          ::  i
         INTEGER(iwp)                          ::  j
-        INTEGER(iwp)                          ::  k
-        INTEGER(iwp)                          ::  kb !< variable to set respective boundary value, depends on facing. 
+        INTEGER(iwp)                          ::  k 
         INTEGER(iwp)                          ::  l  !< running index boundary type, for up- and downward-facing walls
         INTEGER(iwp)                          ::  m  !< running index surface elements
@@ -2860,15 +2863,11 @@
             IF ( ibc_pt_b == 1 )  THEN 
                DO  l = 0, 1 
-!
-!--       Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--       for downward-facing surfaces at topography bottom (k+1). 
-                    kb = MERGE( -1, 1, l == 0 )
-                    DO  m = 1, bc_h(l)%ns
-                       i = bc_h(l)%i(m)     
-                       j = bc_h(l)%j(m)
-                       k = bc_h(l)%k(m)
-                       pt(k+kb,j,i) = pt(k,j,i)
-                    ENDDO
-                 ENDDO
+                  DO  m = 1, bc_h(l)%ns
+                     i = bc_h(l)%i(m)     
+                     j = bc_h(l)%j(m)
+                     k = bc_h(l)%k(m)
+                     pt(k+bc_h(l)%koff,j,i) = pt(k,j,i)
+                  ENDDO
+               ENDDO
             ENDIF
    
Index: palm/trunk/TESTS/cases/rans_tkee/MONITORING/rans_tkee_rc
===================================================================
--- palm/trunk/TESTS/cases/rans_tkee/MONITORING/rans_tkee_rc	(revision 4101)
+++ palm/trunk/TESTS/cases/rans_tkee/MONITORING/rans_tkee_rc	(revision 4102)
@@ -11,9 +11,9 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4100M       *    atmosphere - run with 1D - prerun
+ * PALM 6.0  Rev: 3915M       *    atmosphere - run with 1D - prerun
  ******************************    --------------------------------------------
 
  Date:               2019-07-17    Run:       rans_tkee__gfortran_default__4    
- Time:                 16:59:43    Run-No.:   00
+ Time:                 17:23:23    Run-No.:   00
  Run on host:        gfortran_d
  Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
@@ -220,40 +220,28 @@
 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   2.6700A   3.3592   0.8728   0.0000  0.127  0.00   0.000E+00      0.  0.442E+01  0.442E-02  0.000E+00  0.000E+00   20  29  29    6  29  29    0  29  29    0.000    0.000       0
-  0      1 00:00:02.67   2.6600A   3.3803   1.5756   1.1290  0.151  0.00  -0.406E-02     10.  0.444E+01  0.437E-02  0.202E-03  0.798E-04   19   7  19    2  14  13    2  13   1    0.000    0.000       2
-  0      2 00:00:05.33   0.2370D   3.3829   1.4756   1.2399  0.147  0.00  -0.468E-02     10.  0.444E+01  0.432E-02  0.137E-03  0.457E-04   19   6  21    2  14  13    2  13   1    0.000    0.000       2
-  0      3 00:00:05.56   0.5360D   3.3841   1.5585   1.3626  0.145  0.00  -0.487E-02      0.  0.444E+01  0.422E-02  0.512E-04  0.118E-04   19   6  22    2  17   0    2  13   1    0.000    0.000       2
-  0      4 00:00:06.10   0.7840D   3.3844   1.5478   1.3603  0.145  0.00  -0.497E-02      0.  0.444E+01  0.408E-02  0.412E-04  0.145E-04   19   5  23    2  17   0    2  13   1    0.000    0.000       2
-  0      5 00:00:06.88   1.0200D   3.3843   1.5248   1.3393  0.144  0.00  -0.510E-02      0.  0.444E+01  0.401E-02  0.459E-04  0.160E-04   19   5  25    2  17   0    2  13   1    0.000    0.000       2
-  0      6 00:00:07.90   1.2500D   3.3842   1.5006   1.3115  0.143  0.00  -0.528E-02      0.  0.444E+01  0.397E-02  0.521E-04  0.121E-04   19   3  27    2  17   0    2  13   1    0.000    0.000       2
-  0      7 00:00:09.15   1.2700D   3.3842   1.4774   1.2804  0.141  0.00  -0.554E-02      0.  0.444E+01  0.393E-02  0.643E-04  0.175E-04   19   0   0    2  17   0    2  13   1    0.000    0.000       2
-  0      8 00:00:10.42   1.3200D   3.3840   1.4655   1.2613  0.140  0.00  -0.558E-02      0.  0.444E+01  0.387E-02  0.621E-04  0.164E-04   19   2   0    2  17   0    2  13   1    0.000    0.000       2
-  0      9 00:00:11.74   1.4300D   3.3835   1.4555   1.2678  0.139  0.00  -0.563E-02      0.  0.444E+01  0.380E-02  0.562E-04  0.155E-04   19   1   0    2  17   0    1  12   2    0.000    0.000       2
-  0     10 00:00:13.17   1.5700D   3.3829   1.4445   1.2852  0.138  0.00  -0.566E-02      0.  0.444E+01  0.376E-02  0.568E-04  0.170E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
-  0     11 00:00:14.74   1.7600D   3.3824   1.4336   1.2960  0.137  0.00  -0.573E-02      0.  0.444E+01  0.371E-02  0.587E-04  0.172E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
-  0     12 00:00:16.50   2.0100D   3.3820   1.4223   1.2996  0.136  0.00  -0.590E-02      0.  0.444E+01  0.367E-02  0.625E-04  0.183E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
-  0     13 00:00:18.51   1.9500D   3.3816   1.4105   1.2953  0.135  0.00  -0.602E-02      0.  0.444E+01  0.362E-02  0.674E-04  0.202E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
-  0     14 00:00:20.46   1.8700D   3.3814   1.4124   1.2963  0.134  0.00  -0.644E-02      0.  0.444E+01  0.359E-02  0.635E-04  0.194E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
-  0     15 00:00:22.33   1.8000D   3.3814   1.4174   1.2972  0.133  0.00  -0.626E-02      0.  0.444E+01  0.356E-02  0.594E-04  0.185E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
-  0     16 00:00:24.13   1.7400D   3.3814   1.4235   1.2974  0.133  0.00  -0.622E-02      0.  0.444E+01  0.353E-02  0.562E-04  0.174E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     17 00:00:25.87   1.6900D   3.3815   1.4301   1.2973  0.132  0.00  -0.630E-02      0.  0.444E+01  0.351E-02  0.536E-04  0.165E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     18 00:00:27.56   1.6400D   3.3815   1.4371   1.2970  0.131  0.00  -0.639E-02      0.  0.444E+01  0.349E-02  0.513E-04  0.159E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     19 00:00:29.20   1.6000D   3.3815   1.4448   1.2972  0.131  0.00  -0.659E-02      0.  0.444E+01  0.347E-02  0.491E-04  0.152E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     20 00:00:30.80   1.5500D   3.3816   1.4526   1.2975  0.130  0.00  -0.661E-02      0.  0.444E+01  0.345E-02  0.471E-04  0.147E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     21 00:00:32.35   1.4700D   3.3816   1.4611   1.2985  0.129  0.00  -0.654E-02      0.  0.444E+01  0.343E-02  0.451E-04  0.141E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     22 00:00:33.82   1.4000D   3.3817   1.4711   1.3011  0.129  0.00  -0.649E-02      0.  0.444E+01  0.342E-02  0.426E-04  0.134E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
-  0     23 00:00:35.22   1.3400D   3.3817   1.4808   1.3040  0.129  0.00  -0.646E-02      0.  0.444E+01  0.341E-02  0.403E-04  0.126E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     24 00:00:36.56   1.2800D   3.3818   1.4899   1.3070  0.128  0.00  -0.645E-02      0.  0.444E+01  0.340E-02  0.383E-04  0.120E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     25 00:00:37.84   1.2400D   3.3818   1.4988   1.3101  0.128  0.00  -0.647E-02      0.  0.444E+01  0.339E-02  0.363E-04  0.115E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     26 00:00:39.08   1.1900D   3.3819   1.5066   1.3128  0.127  0.00  -0.650E-02      0.  0.444E+01  0.338E-02  0.349E-04  0.110E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     27 00:00:40.27   1.1500D   3.3819   1.5146   1.3157  0.127  0.00  -0.654E-02      0.  0.444E+01  0.337E-02  0.333E-04  0.106E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     28 00:00:41.42   1.1100D   3.3820   1.5219   1.3184  0.127  0.00  -0.659E-02      0.  0.444E+01  0.336E-02  0.319E-04  0.102E-04   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     29 00:00:42.53   1.0800D   3.3820   1.5289   1.3211  0.127  0.00  -0.663E-02      0.  0.444E+01  0.335E-02  0.306E-04  0.979E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     30 00:00:43.61   1.0500D   3.3821   1.5352   1.3234  0.126  0.00  -0.667E-02      0.  0.444E+01  0.334E-02  0.296E-04  0.949E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     31 00:00:44.66   1.0200D   3.3822   1.5412   1.3257  0.126  0.00  -0.667E-02      0.  0.444E+01  0.334E-02  0.286E-04  0.921E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     32 00:00:45.68   0.9930D   3.3822   1.5470   1.3278  0.126  0.00  -0.664E-02      0.  0.444E+01  0.333E-02  0.276E-04  0.893E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     33 00:00:46.68   0.9680D   3.3823   1.5525   1.3299  0.126  0.00  -0.661E-02      0.  0.444E+01  0.332E-02  0.267E-04  0.868E-05   19   8  23    2  17   0    1  12   2    0.000    0.000       2
-  0     34 00:00:47.64   0.9450D   3.3823   1.5576   1.3318  0.126  0.00  -0.659E-02      0.  0.444E+01  0.332E-02  0.259E-04  0.844E-05   19   8  24    2  17   0    1  12   2    0.000    0.000       2
-  0     35 00:00:48.59   0.9240D   3.3824   1.5625   1.3335  0.125  0.00  -0.657E-02      0.  0.444E+01  0.331E-02  0.252E-04  0.823E-05   19   8  24    2  17   0    1  12   2    0.000    0.000       2
-  0     36 00:00:49.51   0.9040D   3.3825   1.5670   1.3351  0.125  0.00  -0.655E-02      0.  0.444E+01  0.331E-02  0.245E-04  0.803E-05   19   9  24    2  17   0    1  12   2    0.000    0.000       2
-  0     37 00:00:50.42   0.8850D   3.3826   1.5714   1.3366  0.125  0.00  -0.654E-02      0.  0.444E+01  0.330E-02  0.239E-04  0.785E-05   19   9  24    2  17   0    1  12   2    0.000    0.000       2
+  0      0 00:00:00.00   0.0200A   3.3592   0.8728   0.0000  0.127  0.00   0.000E+00      0.  0.442E+01  0.442E-02  0.000E+00  0.000E+00   20  29  29    6  29  29    0  29  29    0.000    0.000       0
+  0      1 00:00:00.02   0.0400A   3.3802   1.7430  -1.2281  0.153  0.00  -0.390E-02     10.  0.444E+01  0.460E-02  0.155E-03  0.701E-04   19   7  19    2  14  13    2  26  29    0.000    0.000       2
+  0      2 00:00:00.06   0.0800A   3.3832   1.7620  -1.3655  0.151  0.20  -0.421E-02     10.  0.444E+01  0.493E-02  0.874E-04  0.359E-04   19   6  21    2  14  13    2  26  29    0.000    0.000       2
+  0      3 00:00:00.14   0.1600A   3.3844   1.7515   1.3704  0.150  0.00  -0.428E-02      0.  0.444E+01  0.496E-02  0.977E-04  0.323E-04   19   5  22    2  14  13    2  13   1    0.000    0.000       2
+  0      4 00:00:00.30   0.3200A   3.3845   1.7361   1.3743  0.150  0.00  -0.429E-02      0.  0.444E+01  0.456E-02  0.521E-04  0.687E-05   19   5  24    2  14  13    2  13   1    0.000    0.000       2
+  0      5 00:00:00.62   0.6400A   3.3838   1.7102   1.3741  0.149  0.00  -0.432E-02      0.  0.444E+01  0.437E-02  0.385E-04  0.278E-04   19   4  26    2  14  13    2  13   1    0.000    0.000       2
+  0      6 00:00:01.26   1.2800A   3.3833   1.6639   1.3660  0.148  0.00  -0.437E-02      0.  0.444E+01  0.430E-02  0.704E-04  0.332E-04   19   1   0    2  14  13    2  13   1    0.000    0.000       2
+  0      7 00:00:02.54   2.5600A   3.3827   1.5852   1.3342  0.147  0.39  -0.449E-02     30.  0.444E+01  0.420E-02  0.618E-04  0.165E-04   19   1   0    2  14  13    2  13   1    0.000    0.000       2
+  0      8 00:00:05.09   2.6600A   3.3826   1.4669   1.2489  0.145  0.00  -0.484E-02     20.  0.444E+01  0.421E-02  0.144E-03  0.499E-04   20   3  31    2  14  13    2  13   1    0.000    0.000       2
+  0      9 00:00:07.75   2.6200D   3.3825   1.4248   1.2091  0.142  0.00  -0.540E-02      0.  0.444E+01  0.399E-02  0.137E-03  0.424E-04   20   3   0    2  17   0    2  13   1    0.000    0.000       2
+  0     10 00:00:10.37   2.6600A   3.3822   1.4058   1.2176  0.140  0.00  -0.536E-02      0.  0.444E+01  0.380E-02  0.971E-04  0.324E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
+  0     11 00:00:13.03   2.6600A   3.3818   1.3920   1.2409  0.138  0.00  -0.554E-02      0.  0.444E+01  0.372E-02  0.912E-04  0.361E-04   20   2   0    2  17   0    1  12   2    0.000    0.000       2
+  0     12 00:00:15.69   2.6600A   3.3818   1.3841   1.2506  0.137  0.00  -0.573E-02      0.  0.444E+01  0.366E-02  0.867E-04  0.300E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
+  0     13 00:00:18.35   2.6600A   3.3820   1.3790   1.2520  0.135  0.00  -0.581E-02      0.  0.444E+01  0.362E-02  0.864E-04  0.282E-04   19   6  23    2  17   0    1  12   2    0.000    0.000       2
+  0     14 00:00:21.01   2.6600A   3.3822   1.3757   1.2484  0.134  0.00  -0.617E-02      0.  0.444E+01  0.358E-02  0.862E-04  0.277E-04   19   6  23    2  17   0    1  12   2    0.000    0.000       2
+  0     15 00:00:23.67   2.6600A   3.3824   1.3733   1.2423  0.133  0.00  -0.644E-02      0.  0.444E+01  0.353E-02  0.800E-04  0.249E-04   19   6  24    2  17   0    1  12   2    0.000    0.000       2
+  0     16 00:00:26.33   2.6600A   3.3825   1.3716   1.2356  0.132  0.00  -0.633E-02      0.  0.444E+01  0.348E-02  0.755E-04  0.250E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
+  0     17 00:00:29.00   2.6600A   3.3827   1.3703   1.2296  0.131  0.00  -0.648E-02      0.  0.444E+01  0.345E-02  0.734E-04  0.242E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
+  0     18 00:00:31.66   2.6600A   3.3829   1.3692   1.2245  0.130  0.00  -0.650E-02      0.  0.444E+01  0.342E-02  0.719E-04  0.232E-04   19   7  23    2  17   0    1  12   2    0.000    0.000       2
+  0     19 00:00:34.32   2.6600A   3.3830   1.3683   1.2206  0.129  0.00  -0.661E-02      0.  0.444E+01  0.339E-02  0.709E-04  0.228E-04   19   7  24    2  17   0    1  12   2    0.000    0.000       2
+  0     20 00:00:36.98   2.6600A   3.3832   1.3675   1.2176  0.128  0.00  -0.648E-02      0.  0.444E+01  0.336E-02  0.692E-04  0.224E-04   19   7  24    2  17   0    1  12   2    0.000    0.000       2
+  0     21 00:00:39.64   2.6600A   3.3834   1.3668   1.2153  0.128  0.00  -0.649E-02      0.  0.444E+01  0.334E-02  0.674E-04  0.221E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
+  0     22 00:00:42.29   2.6600A   3.3836   1.3662   1.2136  0.127  0.00  -0.658E-02      0.  0.444E+01  0.331E-02  0.661E-04  0.218E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
+  0     23 00:00:44.95   2.6600A   3.3837   1.3657   1.2122  0.127  0.00  -0.664E-02      0.  0.444E+01  0.329E-02  0.651E-04  0.215E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
+  0     24 00:00:47.61   2.6600A   3.3839   1.3652   1.2111  0.126  0.00  -0.661E-02      0.  0.444E+01  0.327E-02  0.642E-04  0.212E-04   19   8  24    2  17   0    1  12   2    0.000    0.000       2
+  0     25 00:00:50.27   2.6600A   3.3841   1.3647   1.2101  0.126  0.00  -0.654E-02      0.  0.444E+01  0.325E-02  0.634E-04  0.210E-04   19   8  25    2  17   0    1  12   2    0.000    0.000       2
Index: palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc
===================================================================
--- palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc	(revision 4101)
+++ palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc	(revision 4102)
@@ -1,9 +1,9 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4082M       *    atmosphere - run without 1D - prerun
+ * PALM 6.0  Rev: 3915M       *    atmosphere - run without 1D - prerun
  ******************************    --------------------------------------------
 
- Date:               2019-07-10    Run:       urban_environment__gfortran_defaul
- Time:                 17:38:31    Run-No.:   00
+ Date:               2019-07-17    Run:       urban_environment__gfortran_defaul
+ Time:                 17:22:42    Run-No.:   00
  Run on host:        gfortran_d
  Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
@@ -356,5 +356,5 @@
   0      9 00:00:06.82   0.7580A  -2.3738  -1.3772   0.8161  0.080  0.25  -0.913E+00     18.  0.137E+01  0.146E-01  0.760E-03  0.196E-03   60  16   9   11  14   0   20  17   0    0.000    0.000       0
   0     10 00:00:07.58   0.7580A  -2.3739  -1.3551   0.8090  0.079  0.25  -0.914E+00     18.  0.137E+01  0.149E-01  0.761E-03  0.196E-03   60  15   8   11  14   0   20  17   0    0.000    0.000       0
-  0     11 00:00:08.33   0.7580A  -2.3739  -1.3370   0.8029  0.079  0.25  -0.915E+00     18.  0.137E+01  0.153E-01  0.764E-03  0.196E-03   60  15   7   11  14   0   20  17   0    0.000    0.000       0
+  0     11 00:00:08.33   0.7580A  -2.3739  -1.3371   0.8029  0.079  0.25  -0.915E+00     18.  0.137E+01  0.153E-01  0.764E-03  0.196E-03   60  15   7   11  14   0   20  17   0    0.000    0.000       0
   0     12 00:00:09.09   0.7580A  -2.3740  -1.3232   0.7984  0.078  0.26  -0.916E+00     18.  0.136E+01  0.156E-01  0.768E-03  0.195E-03   60  15   6   16  14   0   20  17   0    0.000    0.000       0
   0     13 00:00:09.85   0.7580A  -2.3741  -1.3182   0.7956  0.078  0.26  -0.917E+00     18.  0.136E+01  0.160E-01  0.773E-03  0.195E-03   60  14   6   16  14   0   20  17   0    0.000    0.000       0
@@ -362,79 +362,79 @@
   0     15 00:00:11.37   0.7580A  -2.3743  -1.2976   0.7936  0.077  0.26  -0.922E+00     18.  0.136E+01  0.168E-01  0.783E-03  0.194E-03   60  13   4   16  14   0   20  17   0    0.000    0.000       0
   0     16 00:00:12.12   0.7580A  -2.3743  -1.2822   0.8042  0.076  0.27  -0.921E+00     18.  0.136E+01  0.172E-01  0.786E-03  0.193E-03   60  13   3   16  14   0   11  14  18    0.000    0.000       0
-  0     17 00:00:12.88   0.7580A  -2.3744  -1.2633   0.8294  0.076  0.27  -0.919E+00     18.  0.136E+01  0.176E-01  0.788E-03  0.192E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
-  0     18 00:00:13.64   0.7580A  -2.3745  -1.2574   0.8604  0.076  0.27  -0.917E+00     18.  0.136E+01  0.180E-01  0.792E-03  0.191E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
+  0     17 00:00:12.88   0.7580A  -2.3744  -1.2633   0.8295  0.076  0.27  -0.919E+00     18.  0.136E+01  0.176E-01  0.788E-03  0.192E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
+  0     18 00:00:13.64   0.7580A  -2.3745  -1.2574   0.8605  0.076  0.27  -0.917E+00     18.  0.136E+01  0.180E-01  0.792E-03  0.191E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
   0     19 00:00:14.40   0.7580A  -2.3746  -1.2531   0.8936  0.075  0.27  -0.915E+00     18.  0.136E+01  0.184E-01  0.796E-03  0.190E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
-  0     20 00:00:15.16   0.7580A  -2.3746  -1.2498   0.9251  0.075  0.28  -0.914E+00     18.  0.136E+01  0.188E-01  0.799E-03  0.188E-03   60  12   0   11  14   0   10  14  18    0.000    0.000       0
+  0     20 00:00:15.16   0.7580A  -2.3746  -1.2498   0.9251  0.075  0.28  -0.915E+00     18.  0.136E+01  0.188E-01  0.799E-03  0.188E-03   60  12   0   11  14   0   10  14  18    0.000    0.000       0
   0     21 00:00:15.91   0.7580A  -2.3747  -1.2472   0.9531  0.075  0.28  -0.914E+00     18.  0.136E+01  0.192E-01  0.802E-03  0.187E-03   60  11  19   11  14   0   10  14  18    0.000    0.000       0
   0     22 00:00:16.67   0.7580A  -2.3748  -1.2453   0.9762  0.074  0.40  -0.913E+00     40.  0.136E+01  0.196E-01  0.806E-03  0.185E-03   60  11  18   11  14   0   10  14  18    0.000    0.000       0
-  0     23 00:00:17.43   0.7580A  -2.3749  -1.2439   0.9937  0.074  0.40  -0.912E+00     40.  0.136E+01  0.200E-01  0.809E-03  0.184E-03   60  10  17   11  14   0   10  14  18    0.000    0.000       0
-  0     24 00:00:18.19   0.7580A  -2.3749  -1.2430   1.0055  0.074  0.40  -0.911E+00     40.  0.136E+01  0.205E-01  0.812E-03  0.182E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
+  0     23 00:00:17.43   0.7580A  -2.3749  -1.2439   0.9936  0.074  0.40  -0.912E+00     40.  0.136E+01  0.200E-01  0.809E-03  0.184E-03   60  10  17   11  14   0   10  14  18    0.000    0.000       0
+  0     24 00:00:18.19   0.7580A  -2.3749  -1.2430   1.0055  0.074  0.40  -0.912E+00     40.  0.136E+01  0.205E-01  0.812E-03  0.182E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
   0     25 00:00:18.95   0.7580A  -2.3750  -1.2425   1.0123  0.073  0.41  -0.911E+00     40.  0.136E+01  0.209E-01  0.815E-03  0.181E-03   60  10  16   11  14   0   10  14  18    0.000    0.000       0
   0     26 00:00:19.70   0.7580A  -2.3751  -1.2422   1.0163  0.073  0.41  -0.910E+00     40.  0.136E+01  0.213E-01  0.819E-03  0.180E-03   60  10  15   11  14   0   11  14  18    0.000    0.000       0
-  0     27 00:00:20.46   0.7580A  -2.3752  -1.2420   1.0281  0.073  0.41  -0.910E+00     40.  0.136E+01  0.217E-01  0.823E-03  0.179E-03   60   9  14   11  14   0   11  14  18    0.000    0.000       0
-  0     28 00:00:21.22   0.7580A  -2.3752  -1.2417   1.0345  0.073  0.42  -0.909E+00     40.  0.136E+01  0.222E-01  0.827E-03  0.178E-03   60   9  13   11  14   0   11  14  18    0.000    0.000       0
-  0     29 00:00:21.98   0.7580A  -2.3753  -1.2416   1.0356  0.072  0.42  -0.909E+00     40.  0.136E+01  0.226E-01  0.831E-03  0.177E-03   60   9  12   11  14   0   11  14  18    0.000    0.000       0
-  0     30 00:00:22.74   0.7580A  -2.3754  -1.2414   1.0322  0.072  0.42  -0.908E+00     40.  0.136E+01  0.230E-01  0.835E-03  0.176E-03   60   8  11   11  14   0   11  14  18    0.000    0.000       0
-  0     31 00:00:23.49   0.7580A  -2.3755  -1.2414   1.0249  0.072  0.42  -0.908E+00     40.  0.136E+01  0.235E-01  0.839E-03  0.175E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
-  0     32 00:00:24.25   0.7580A  -2.3755  -1.2414   1.0153  0.072  0.43  -0.908E+00     40.  0.136E+01  0.239E-01  0.843E-03  0.174E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
-  0     33 00:00:25.01   0.7580A  -2.3756  -1.2413   1.0047  0.071  0.43  -0.908E+00     40.  0.136E+01  0.243E-01  0.848E-03  0.173E-03   60   7   9   11  14   0   11  14  18    0.000    0.000       0
+  0     27 00:00:20.46   0.7580A  -2.3752  -1.2419   1.0281  0.073  0.41  -0.910E+00     40.  0.136E+01  0.217E-01  0.823E-03  0.179E-03   60   9  14   11  14   0   11  14  18    0.000    0.000       0
+  0     28 00:00:21.22   0.7580A  -2.3752  -1.2417   1.0344  0.073  0.42  -0.909E+00     40.  0.136E+01  0.222E-01  0.827E-03  0.178E-03   60   9  13   11  14   0   11  14  18    0.000    0.000       0
+  0     29 00:00:21.98   0.7580A  -2.3753  -1.2415   1.0356  0.072  0.42  -0.909E+00     40.  0.136E+01  0.226E-01  0.831E-03  0.177E-03   60   9  12   11  14   0   11  14  18    0.000    0.000       0
+  0     30 00:00:22.74   0.7580A  -2.3754  -1.2414   1.0321  0.072  0.42  -0.909E+00     40.  0.136E+01  0.230E-01  0.835E-03  0.176E-03   60   8  11   11  14   0   11  14  18    0.000    0.000       0
+  0     31 00:00:23.49   0.7580A  -2.3755  -1.2413   1.0248  0.072  0.42  -0.908E+00     40.  0.136E+01  0.234E-01  0.839E-03  0.175E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
+  0     32 00:00:24.25   0.7580A  -2.3755  -1.2413   1.0152  0.072  0.43  -0.908E+00     40.  0.136E+01  0.239E-01  0.843E-03  0.174E-03   60   8  10   11  14   0   11  14  18    0.000    0.000       0
+  0     33 00:00:25.01   0.7580A  -2.3756  -1.2413   1.0046  0.071  0.43  -0.908E+00     40.  0.136E+01  0.243E-01  0.848E-03  0.173E-03   60   7   9   11  14   0   11  14  18    0.000    0.000       0
   0     34 00:00:25.77   0.7580A  -2.3757  -1.2410   0.9942  0.071  0.43  -0.908E+00     40.  0.136E+01  0.247E-01  0.852E-03  0.173E-03   60   7   8   11  14   0   10  14  18    0.000    0.000       0
-  0     35 00:00:26.53   0.7580A  -2.3758  -1.2403   1.0107  0.071  0.44  -0.908E+00     40.  0.136E+01  0.252E-01  0.856E-03  0.172E-03   60   7   7   11  14   0   11  14  17    0.000    0.000       0
-  0     36 00:00:27.28   0.7580A  -2.3759  -1.2388   1.0283  0.071  0.44  -0.908E+00     40.  0.136E+01  0.256E-01  0.860E-03  0.171E-03   60   6   6   11  14   0   11  14  17    0.000    0.000       0
-  0     37 00:00:28.04   0.7580A  -2.3759  -1.2364   1.0443  0.071  0.44  -0.908E+00     40.  0.136E+01  0.261E-01  0.863E-03  0.171E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
-  0     38 00:00:28.80   0.7580A  -2.3760  -1.2328   1.0568  0.071  0.45  -0.908E+00     40.  0.136E+01  0.265E-01  0.867E-03  0.170E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
-  0     39 00:00:29.56   0.7580A  -2.3761  -1.2280   1.0647  0.070  0.45  -0.908E+00     40.  0.136E+01  0.269E-01  0.871E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
-  0     40 00:00:30.32   0.7580A  -2.3762  -1.2218   1.0674  0.070  0.45  -0.908E+00     40.  0.136E+01  0.274E-01  0.874E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
-  0     41 00:00:31.07   0.7580A  -2.3762  -1.2204   1.0649  0.070  0.46  -0.908E+00     40.  0.136E+01  0.278E-01  0.878E-03  0.168E-03   60   5   2   11  13   0   11  14  17    0.000    0.000       0
-  0     42 00:00:31.83   0.7570A  -2.3763  -1.2190   1.0762  0.070  0.46  -0.908E+00     40.  0.136E+01  0.283E-01  0.882E-03  0.168E-03   60   4   1   11  13   0   12  14  16    0.000    0.000       0
-  0     43 00:00:32.59   0.7570A  -2.3764  -1.2156   1.0984  0.070  0.46  -0.909E+00     40.  0.136E+01  0.287E-01  0.885E-03  0.168E-03   60   4   0   11  13   0   12  14  16    0.000    0.000       0
-  0     44 00:00:33.35   0.7570A  -2.3765  -1.2102   1.1184  0.070  0.47  -0.909E+00     40.  0.136E+01  0.292E-01  0.891E-03  0.168E-03   60   4  19   11  13   0   12  14  16    0.000    0.000       0
-  0     45 00:00:34.10   0.7570A  -2.3765  -1.2026   1.1342  0.070  0.47  -0.910E+00     40.  0.136E+01  0.296E-01  0.896E-03  0.168E-03   60   3  18   11  13   0   12  14  16    0.000    0.000       0
-  0     46 00:00:34.86   0.7570A  -2.3766  -1.1928   1.1438  0.070  0.47  -0.911E+00     40.  0.136E+01  0.301E-01  0.903E-03  0.168E-03   60   3  17   11  13   0   12  14  16    0.000    0.000       0
-  0     47 00:00:35.62   0.7570A  -2.3767  -1.1810   1.1535  0.070  0.48  -0.911E+00     40.  0.136E+01  0.305E-01  0.909E-03  0.168E-03   60   3  17   11  13   0   13  14  16    0.000    0.000       0
-  0     48 00:00:36.37   0.7570A  -2.3768  -1.1670   1.1658  0.070  0.48  -0.912E+00     40.  0.136E+01  0.310E-01  0.917E-03  0.168E-03   60   3  16   11  13   0   13  14  16    0.000    0.000       0
-  0     49 00:00:37.13   0.7570A  -2.3768  -1.1510   1.1795  0.069  0.48  -0.913E+00     40.  0.136E+01  0.315E-01  0.925E-03  0.169E-03   60   2  15   11  13   0   16  14  17    0.000    0.000       0
-  0     50 00:00:37.89   0.7570A  -2.3769  -1.1333   1.1993  0.069  0.49  -0.913E+00     40.  0.136E+01  0.320E-01  0.933E-03  0.169E-03   60   2  14   11  13   0   16  14  17    0.000    0.000       0
-  0     51 00:00:38.64   0.7570A  -2.3770  -1.1139   1.2169  0.069  0.49  -0.912E+00     40.  0.136E+01  0.324E-01  0.942E-03  0.169E-03   60   2  13   11  13   0   17  14  17    0.000    0.000       0
-  0     52 00:00:39.40   0.7570A  -2.3771  -1.0929   1.2376  0.069  0.49  -0.911E+00     40.  0.136E+01  0.329E-01  0.951E-03  0.170E-03   60   1  12   11  13   0   17  14  17    0.000    0.000       0
-  0     53 00:00:40.16   0.7570A  -2.3771  -1.0756   1.2530  0.069  0.50  -0.911E+00     40.  0.136E+01  0.334E-01  0.960E-03  0.170E-03   60   1  11   13  14   0   17  14  17    0.000    0.000       0
-  0     54 00:00:40.92   0.7570A  -2.3772  -1.0760   1.2625  0.069  0.50  -0.910E+00     40.  0.136E+01  0.339E-01  0.970E-03  0.171E-03   60   1  10   13  14   0   17  14  17    0.000    0.000       0
-  0     55 00:00:41.67   0.7570A  -2.3773  -1.0746   1.2666  0.069  0.50  -0.910E+00     40.  0.136E+01  0.344E-01  0.980E-03  0.172E-03   60   1  10   13  14   0   17  14  17    0.000    0.000       0
-  0     56 00:00:42.43   0.7570A  -2.3774  -1.0713   1.2691  0.069  0.51  -0.910E+00     40.  0.136E+01  0.349E-01  0.990E-03  0.172E-03   60   0   8   13  14   0   18  14  17    0.000    0.000       0
-  0     57 00:00:43.19   0.7570A  -2.3774  -1.0664   1.2730  0.069  0.51  -0.909E+00     40.  0.136E+01  0.354E-01  0.100E-02  0.173E-03   60   0   8   13  14   0   18  14  17    0.000    0.000       0
-  0     58 00:00:43.94   0.7570A  -2.3775  -1.0603   1.2697  0.069  0.55  -0.909E+00     48.  0.136E+01  0.358E-01  0.101E-02  0.173E-03   60   0   7   13  14   0   18  14  17    0.000    0.000       0
-  0     59 00:00:44.70   0.7570A  -2.3776  -1.0550   1.2800  0.069  0.56  -0.909E+00     48.  0.136E+01  0.363E-01  0.102E-02  0.174E-03   60  19   6   14  14   0   19  14  18    0.000    0.000       0
-  0     60 00:00:45.46   0.7570A  -2.3777  -1.0550   1.2932  0.069  0.56  -0.909E+00     48.  0.136E+01  0.368E-01  0.103E-02  0.175E-03   60  19   5   14  14   0   19  14  18    0.000    0.000       0
-  0     61 00:00:46.21   0.7570A  -2.3777  -1.0561   1.2996  0.069  0.56  -0.909E+00     48.  0.136E+01  0.373E-01  0.104E-02  0.175E-03   60  19   5   14  14   0   19  14  18    0.000    0.000       0
-  0     62 00:00:46.97   0.7570A  -2.3778  -1.0583   1.2988  0.069  0.53  -0.909E+00     40.  0.136E+01  0.378E-01  0.104E-02  0.176E-03   60  18   3   14  14   0   19  14  18    0.000    0.000       0
-  0     63 00:00:47.73   0.7570A  -2.3779  -1.0611   1.2908  0.069  0.53  -0.910E+00     40.  0.136E+01  0.383E-01  0.105E-02  0.176E-03   60  18   3   14  14   0   19  14  18    0.000    0.000       0
-  0     64 00:00:48.49   0.7570A  -2.3780  -1.0641   1.2764  0.069  0.53  -0.910E+00     40.  0.136E+01  0.388E-01  0.106E-02  0.177E-03   60  18   2   14  14   0   19  14  18    0.000    0.000       0
-  0     65 00:00:49.24   0.7570A  -2.3780  -1.0667   1.2565  0.069  0.53  -0.912E+00     40.  0.136E+01  0.393E-01  0.107E-02  0.177E-03   60  18   1   15  14   0   19  14  18    0.000    0.000       0
-  0     66 00:00:50.00   0.7570A  -2.3781  -1.0750   1.2323  0.069  0.53  -0.913E+00     40.  0.135E+01  0.398E-01  0.108E-02  0.178E-03   60  17   0   15  14   0   19  14  18    0.000    0.000       0
-  0     67 00:00:50.76   0.7570A  -2.3782  -1.0801   1.2051  0.069  0.53  -0.916E+00     40.  0.135E+01  0.403E-01  0.109E-02  0.178E-03   60  17  19   15  14   0   19  14  18    0.000    0.000       0
-  0     68 00:00:51.51   0.7570A  -2.3783  -1.0824   1.1761  0.069  0.53  -0.922E+00     40.  0.135E+01  0.408E-01  0.110E-02  0.178E-03   60  17  19   15  13   0   19  14  18    0.000    0.000       0
-  0     69 00:00:52.27   0.7570A  -2.3783  -1.0899   1.1464  0.069  0.53  -0.933E+00     40.  0.135E+01  0.413E-01  0.110E-02  0.179E-03   60  16  17   15  13   0   19  14  18    0.000    0.000       0
-  0     70 00:00:53.03   0.7570A  -2.3784  -1.0938   1.1212  0.069  0.53  -0.974E+00     40.  0.135E+01  0.418E-01  0.111E-02  0.179E-03   60  16  16   15  13   0    6  13   2    0.000    0.000       0
-  0     71 00:00:53.78   0.7570A  -2.3785  -1.0940   1.1570  0.069  0.53  -0.104E+01     40.  0.135E+01  0.423E-01  0.112E-02  0.180E-03   60  16  16   15  13   0    6  13   2    0.000    0.000       0
-  0     72 00:00:54.54   0.7570A  -2.3786  -1.0909   1.1893  0.069  0.53  -0.947E+00     40.  0.135E+01  0.427E-01  0.113E-02  0.180E-03   60  15  14   15  13   0    6  13   2    0.000    0.000       0
-  0     73 00:00:55.30   0.7570A  -2.3786  -1.0854   1.2176  0.069  0.53  -0.935E+00     40.  0.135E+01  0.432E-01  0.114E-02  0.180E-03   60  15  14   15  13   0    6  13   2    0.000    0.000       0
-  0     74 00:00:56.06   0.7570A  -2.3787  -1.0783   1.2412  0.069  0.54  -0.938E+00     40.  0.135E+01  0.437E-01  0.115E-02  0.181E-03   60  15  13   15  13   0    6  13   2    0.000    0.000       0
-  0     75 00:00:56.81   0.7570A  -2.3788  -1.0704   1.2599  0.069  0.54  -0.949E+00     40.  0.135E+01  0.441E-01  0.115E-02  0.181E-03   60  14  12   15  13   0    6  13   2    0.000    0.000       0
-  0     76 00:00:57.57   0.7570A  -2.3789  -1.0622   1.2927  0.069  0.54  -0.934E+00     40.  0.135E+01  0.446E-01  0.116E-02  0.182E-03   60  14  11   15  13   0    6  14   3    0.000    0.000       0
-  0     77 00:00:58.33   0.7570A  -2.3789  -1.0538   1.3263  0.070  0.54  -0.920E+00     40.  0.135E+01  0.450E-01  0.117E-02  0.182E-03   60  14  10   15  13   0    6  14   3    0.000    0.000       0
-  0     78 00:00:59.08   0.7570A  -2.3790  -1.0455   1.3553  0.070  0.54  -0.913E+00     40.  0.135E+01  0.454E-01  0.118E-02  0.182E-03   60  14  10   15  13   0    6  14   3    0.000    0.000       0
-  0     79 00:00:59.84   0.7570A  -2.3791  -1.0371   1.3795  0.070  0.54  -0.909E+00     40.  0.135E+01  0.459E-01  0.119E-02  0.183E-03   60  13   8   15  13   0    6  14   3    0.000    0.000       0
-  0     80 00:01:00.60   0.7570A  -2.3791  -1.0287   1.3988  0.070  0.54  -0.907E+00     40.  0.135E+01  0.463E-01  0.120E-02  0.183E-03   60  13   8   15  13   0    6  14   3    0.000    0.000       0
-  0     81 00:01:01.35   0.7570A  -2.3792  -1.0357   1.4134  0.070  0.54  -0.906E+00     40.  0.135E+01  0.467E-01  0.121E-02  0.183E-03   60  13   7   18  13   0    6  14   3    0.000    0.000       0
-  0     82 00:01:02.11   0.7570A  -2.3793  -1.0477   1.4235  0.070  0.54  -0.906E+00     40.  0.135E+01  0.470E-01  0.122E-02  0.183E-03   60  13   6   18  13   0    6  14   3    0.000    0.000       0
-  0     83 00:01:02.87   0.7570A  -2.3794  -1.0575   1.4295  0.070  0.54  -0.911E+00     40.  0.135E+01  0.474E-01  0.122E-02  0.183E-03   60  12   5   18  13   0    6  14   3    0.000    0.000       0
-  0     84 00:01:03.63   0.7560A  -2.3794  -1.0640   1.4321  0.070  0.55  -0.919E+00     40.  0.135E+01  0.478E-01  0.123E-02  0.183E-03   60  12   4   18  13   0    7  14   3    0.000    0.000       0
-  0     85 00:01:04.38   0.7560A  -2.3795  -1.0669   1.4539  0.070  0.55  -0.908E+00     40.  0.135E+01  0.481E-01  0.123E-02  0.183E-03   60  12   4   18  13   0    7  14   3    0.000    0.000       0
-  0     86 00:01:05.14   0.7560A  -2.3796  -1.0663   1.4705  0.070  0.55  -0.900E+00     40.  0.135E+01  0.485E-01  0.124E-02  0.183E-03   60  11   2   18  13   0    7  14   3    0.000    0.000       0
-  0     87 00:01:05.89   0.7560A  -2.3797  -1.0652   1.4821  0.070  0.56  -0.895E+00     40.  0.135E+01  0.488E-01  0.125E-02  0.183E-03   60  11   1   17  12   0    7  14   3    0.000    0.000       0
-  0     88 00:01:06.65   0.7560A  -2.3797  -1.0716   1.4893  0.070  0.56  -0.892E+00     40.  0.135E+01  0.491E-01  0.125E-02  0.183E-03   60  11   1   17  12   0    7  14   3    0.000    0.000       0
-  0     89 00:01:07.41   0.7560A  -2.3798  -1.0748   1.4927  0.070  0.56  -0.890E+00     40.  0.135E+01  0.495E-01  0.126E-02  0.183E-03   60  10  19   17  12   0    7  14   3    0.000    0.000       0
-  0     90 00:01:08.16   0.7560A  -2.3799  -1.0747   1.4929  0.070  0.57  -0.889E+00     40.  0.135E+01  0.498E-01  0.126E-02  0.183E-03   60  10  19   17  12   0    7  14   3    0.000    0.000       0
-  0     91 00:01:08.92   0.7560A  -2.3800   1.0782   1.4907  0.071  0.57  -0.887E+00     40.  0.135E+01  0.500E-01  0.127E-02  0.184E-03   60  10  18    1   2  14    7  14   3    0.000    0.000       0
-  0     92 00:01:09.67   0.7560A  -2.3800   1.0850   1.4868  0.071  0.60  -0.886E+00     46.  0.135E+01  0.503E-01  0.127E-02  0.184E-03   60   9  16    1   2  14    7  14   3    0.000    0.000       0
-  0     93 00:01:10.43   0.7560A  -2.3801   1.0912   1.4817  0.071  0.61  -0.885E+00     46.  0.135E+01  0.506E-01  0.128E-02  0.184E-03   60   9  16    1   2  14    7  14   3    0.000    0.000       0
+  0     35 00:00:26.53   0.7580A  -2.3758  -1.2402   1.0107  0.071  0.44  -0.908E+00     40.  0.136E+01  0.252E-01  0.856E-03  0.172E-03   60   7   7   11  14   0   11  14  17    0.000    0.000       0
+  0     36 00:00:27.28   0.7580A  -2.3759  -1.2387   1.0283  0.071  0.44  -0.908E+00     40.  0.136E+01  0.256E-01  0.860E-03  0.171E-03   60   6   6   11  14   0   11  14  17    0.000    0.000       0
+  0     37 00:00:28.04   0.7580A  -2.3759  -1.2363   1.0443  0.071  0.44  -0.908E+00     40.  0.136E+01  0.261E-01  0.863E-03  0.171E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
+  0     38 00:00:28.80   0.7580A  -2.3760  -1.2327   1.0568  0.071  0.45  -0.908E+00     40.  0.136E+01  0.265E-01  0.866E-03  0.170E-03   60   6   5   11  14   0   11  14  17    0.000    0.000       0
+  0     39 00:00:29.56   0.7580A  -2.3761  -1.2279   1.0647  0.070  0.45  -0.908E+00     40.  0.136E+01  0.269E-01  0.870E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
+  0     40 00:00:30.32   0.7580A  -2.3762  -1.2217   1.0673  0.070  0.45  -0.908E+00     40.  0.136E+01  0.274E-01  0.874E-03  0.169E-03   60   5   3   11  14   0   11  14  17    0.000    0.000       0
+  0     41 00:00:31.07   0.7580A  -2.3762  -1.2203   1.0648  0.070  0.46  -0.908E+00     40.  0.136E+01  0.278E-01  0.878E-03  0.168E-03   60   5   2   11  13   0   11  14  17    0.000    0.000       0
+  0     42 00:00:31.83   0.7570A  -2.3763  -1.2189   1.0761  0.070  0.46  -0.908E+00     40.  0.136E+01  0.283E-01  0.882E-03  0.168E-03   60   4   1   11  13   0   12  14  16    0.000    0.000       0
+  0     43 00:00:32.59   0.7570A  -2.3764  -1.2156   1.0983  0.070  0.46  -0.909E+00     40.  0.136E+01  0.287E-01  0.885E-03  0.168E-03   60   4   0   11  13   0   12  14  16    0.000    0.000       0
+  0     44 00:00:33.35   0.7570A  -2.3765  -1.2101   1.1183  0.070  0.47  -0.909E+00     40.  0.136E+01  0.292E-01  0.890E-03  0.168E-03   60   4  19   11  13   0   12  14  16    0.000    0.000       0
+  0     45 00:00:34.10   0.7570A  -2.3765  -1.2025   1.1340  0.070  0.47  -0.910E+00     40.  0.136E+01  0.296E-01  0.896E-03  0.168E-03   60   3  18   11  13   0   12  14  16    0.000    0.000       0
+  0     46 00:00:34.86   0.7570A  -2.3766  -1.1928   1.1437  0.070  0.47  -0.911E+00     40.  0.136E+01  0.301E-01  0.902E-03  0.168E-03   60   3  17   11  13   0   12  14  16    0.000    0.000       0
+  0     47 00:00:35.62   0.7570A  -2.3767  -1.1809   1.1533  0.070  0.48  -0.912E+00     40.  0.136E+01  0.305E-01  0.909E-03  0.168E-03   60   3  17   11  13   0   13  14  16    0.000    0.000       0
+  0     48 00:00:36.37   0.7570A  -2.3768  -1.1669   1.1656  0.069  0.48  -0.913E+00     40.  0.136E+01  0.310E-01  0.916E-03  0.168E-03   60   3  16   11  13   0   13  14  16    0.000    0.000       0
+  0     49 00:00:37.13   0.7570A  -2.3768  -1.1509   1.1794  0.069  0.48  -0.913E+00     40.  0.136E+01  0.315E-01  0.924E-03  0.169E-03   60   2  15   11  13   0   16  14  17    0.000    0.000       0
+  0     50 00:00:37.89   0.7570A  -2.3769  -1.1332   1.1992  0.069  0.49  -0.913E+00     40.  0.136E+01  0.319E-01  0.933E-03  0.169E-03   60   2  14   11  13   0   16  14  17    0.000    0.000       0
+  0     51 00:00:38.64   0.7570A  -2.3770  -1.1138   1.2167  0.069  0.49  -0.912E+00     40.  0.136E+01  0.324E-01  0.942E-03  0.169E-03   60   2  13   11  13   0   17  14  17    0.000    0.000       0
+  0     52 00:00:39.40   0.7570A  -2.3771  -1.0927   1.2374  0.069  0.49  -0.911E+00     40.  0.136E+01  0.329E-01  0.951E-03  0.170E-03   60   1  12   11  13   0   17  14  17    0.000    0.000       0
+  0     53 00:00:40.16   0.7570A  -2.3771  -1.0756   1.2528  0.069  0.50  -0.911E+00     40.  0.136E+01  0.334E-01  0.960E-03  0.170E-03   60   1  11   13  14   0   17  14  17    0.000    0.000       0
+  0     54 00:00:40.92   0.7570A  -2.3772  -1.0760   1.2624  0.069  0.50  -0.910E+00     40.  0.136E+01  0.339E-01  0.970E-03  0.171E-03   60   1  10   13  14   0   17  14  17    0.000    0.000       0
+  0     55 00:00:41.67   0.7570A  -2.3773  -1.0746   1.2665  0.069  0.50  -0.910E+00     40.  0.136E+01  0.344E-01  0.980E-03  0.171E-03   60   1  10   13  14   0   17  14  17    0.000    0.000       0
+  0     56 00:00:42.43   0.7570A  -2.3774  -1.0713   1.2687  0.069  0.51  -0.910E+00     40.  0.136E+01  0.348E-01  0.990E-03  0.172E-03   60   0   8   13  14   0   18  14  17    0.000    0.000       0
+  0     57 00:00:43.19   0.7570A  -2.3774  -1.0664   1.2727  0.069  0.51  -0.909E+00     40.  0.136E+01  0.353E-01  0.100E-02  0.173E-03   60   0   8   13  14   0   18  14  17    0.000    0.000       0
+  0     58 00:00:43.94   0.7570A  -2.3775  -1.0603   1.2694  0.069  0.55  -0.909E+00     48.  0.136E+01  0.358E-01  0.101E-02  0.173E-03   60   0   7   13  14   0   18  14  17    0.000    0.000       0
+  0     59 00:00:44.70   0.7570A  -2.3776  -1.0549   1.2782  0.069  0.56  -0.909E+00     48.  0.136E+01  0.363E-01  0.102E-02  0.174E-03   60  19   6   14  14   0   19  14  18    0.000    0.000       0
+  0     60 00:00:45.46   0.7570A  -2.3777  -1.0549   1.2915  0.069  0.56  -0.910E+00     48.  0.136E+01  0.368E-01  0.103E-02  0.175E-03   60  19   5   14  14   0   19  14  18    0.000    0.000       0
+  0     61 00:00:46.21   0.7570A  -2.3777  -1.0560   1.2979  0.069  0.56  -0.909E+00     48.  0.136E+01  0.373E-01  0.104E-02  0.175E-03   60  19   5   14  14   0   19  14  18    0.000    0.000       0
+  0     62 00:00:46.97   0.7570A  -2.3778  -1.0582   1.2972  0.069  0.53  -0.909E+00     40.  0.136E+01  0.378E-01  0.104E-02  0.176E-03   60  18   3   14  14   0   19  14  18    0.000    0.000       0
+  0     63 00:00:47.73   0.7570A  -2.3779  -1.0611   1.2893  0.069  0.53  -0.910E+00     40.  0.136E+01  0.383E-01  0.105E-02  0.176E-03   60  18   3   14  14   0   19  14  18    0.000    0.000       0
+  0     64 00:00:48.49   0.7570A  -2.3780  -1.0641   1.2750  0.069  0.53  -0.910E+00     40.  0.136E+01  0.388E-01  0.106E-02  0.177E-03   60  18   2   14  14   0   19  14  18    0.000    0.000       0
+  0     65 00:00:49.24   0.7570A  -2.3780  -1.0666   1.2552  0.069  0.53  -0.912E+00     40.  0.136E+01  0.393E-01  0.107E-02  0.177E-03   60  18   1   15  14   0   19  14  18    0.000    0.000       0
+  0     66 00:00:50.00   0.7570A  -2.3781  -1.0749   1.2310  0.069  0.53  -0.913E+00     40.  0.135E+01  0.398E-01  0.108E-02  0.178E-03   60  17   0   15  14   0   19  14  18    0.000    0.000       0
+  0     67 00:00:50.76   0.7570A  -2.3782  -1.0801   1.2039  0.069  0.53  -0.916E+00     40.  0.135E+01  0.403E-01  0.109E-02  0.178E-03   60  17  19   15  14   0   19  14  18    0.000    0.000       0
+  0     68 00:00:51.51   0.7570A  -2.3783  -1.0823   1.1751  0.069  0.53  -0.921E+00     40.  0.135E+01  0.408E-01  0.110E-02  0.178E-03   60  17  19   15  13   0   19  14  18    0.000    0.000       0
+  0     69 00:00:52.27   0.7570A  -2.3783  -1.0898   1.1455  0.069  0.53  -0.932E+00     40.  0.135E+01  0.413E-01  0.110E-02  0.179E-03   60  16  17   15  13   0   19  14  18    0.000    0.000       0
+  0     70 00:00:53.03   0.7570A  -2.3784  -1.0937   1.1201  0.069  0.53  -0.969E+00     40.  0.135E+01  0.418E-01  0.111E-02  0.179E-03   60  16  16   15  13   0    6  13   2    0.000    0.000       0
+  0     71 00:00:53.78   0.7570A  -2.3785  -1.0939   1.1556  0.069  0.53  -0.106E+01     40.  0.135E+01  0.423E-01  0.112E-02  0.180E-03   60  16  16   15  13   0    6  13   2    0.000    0.000       0
+  0     72 00:00:54.54   0.7570A  -2.3786  -1.0908   1.1876  0.069  0.53  -0.949E+00     40.  0.135E+01  0.427E-01  0.113E-02  0.180E-03   60  15  14   15  13   0    6  13   2    0.000    0.000       0
+  0     73 00:00:55.30   0.7570A  -2.3786  -1.0853   1.2156  0.069  0.54  -0.935E+00     40.  0.135E+01  0.432E-01  0.114E-02  0.180E-03   60  15  14   15  13   0    6  13   2    0.000    0.000       0
+  0     74 00:00:56.06   0.7570A  -2.3787  -1.0782   1.2390  0.069  0.54  -0.938E+00     40.  0.135E+01  0.437E-01  0.114E-02  0.181E-03   60  15  13   15  13   0    6  13   2    0.000    0.000       0
+  0     75 00:00:56.81   0.7570A  -2.3788  -1.0703   1.2575  0.069  0.54  -0.949E+00     40.  0.135E+01  0.441E-01  0.115E-02  0.181E-03   60  14  12   15  13   0    6  13   2    0.000    0.000       0
+  0     76 00:00:57.57   0.7570A  -2.3789  -1.0620   1.2906  0.069  0.54  -0.934E+00     40.  0.135E+01  0.446E-01  0.116E-02  0.182E-03   60  14  11   15  13   0    6  14   3    0.000    0.000       0
+  0     77 00:00:58.33   0.7570A  -2.3789  -1.0537   1.3239  0.070  0.54  -0.920E+00     40.  0.135E+01  0.450E-01  0.117E-02  0.182E-03   60  14  10   15  13   0    6  14   3    0.000    0.000       0
+  0     78 00:00:59.08   0.7570A  -2.3790  -1.0453   1.3527  0.070  0.54  -0.913E+00     40.  0.135E+01  0.454E-01  0.118E-02  0.182E-03   60  14  10   15  13   0    6  14   3    0.000    0.000       0
+  0     79 00:00:59.84   0.7570A  -2.3791  -1.0370   1.3767  0.070  0.54  -0.909E+00     40.  0.135E+01  0.458E-01  0.119E-02  0.183E-03   60  13   8   15  13   0    6  14   3    0.000    0.000       0
+  0     80 00:01:00.60   0.7570A  -2.3791  -1.0286   1.3959  0.070  0.54  -0.907E+00     40.  0.135E+01  0.462E-01  0.120E-02  0.183E-03   60  13   8   15  13   0    6  14   3    0.000    0.000       0
+  0     81 00:01:01.35   0.7570A  -2.3792  -1.0356   1.4104  0.070  0.54  -0.906E+00     40.  0.135E+01  0.466E-01  0.121E-02  0.183E-03   60  13   7   18  13   0    6  14   3    0.000    0.000       0
+  0     82 00:01:02.11   0.7570A  -2.3793  -1.0476   1.4206  0.070  0.54  -0.906E+00     40.  0.135E+01  0.470E-01  0.121E-02  0.183E-03   60  13   6   18  13   0    6  14   3    0.000    0.000       0
+  0     83 00:01:02.87   0.7570A  -2.3794  -1.0574   1.4267  0.070  0.55  -0.911E+00     40.  0.135E+01  0.474E-01  0.122E-02  0.183E-03   60  12   5   18  13   0    6  14   3    0.000    0.000       0
+  0     84 00:01:03.63   0.7560A  -2.3794  -1.0639   1.4300  0.070  0.55  -0.919E+00     40.  0.135E+01  0.478E-01  0.123E-02  0.183E-03   60  12   4   18  13   0    7  14   3    0.000    0.000       0
+  0     85 00:01:04.38   0.7560A  -2.3795  -1.0669   1.4516  0.070  0.55  -0.909E+00     40.  0.135E+01  0.481E-01  0.123E-02  0.183E-03   60  12   4   18  13   0    7  14   3    0.000    0.000       0
+  0     86 00:01:05.14   0.7560A  -2.3796  -1.0663   1.4680  0.070  0.55  -0.900E+00     40.  0.135E+01  0.485E-01  0.124E-02  0.183E-03   60  11   2   18  13   0    7  14   3    0.000    0.000       0
+  0     87 00:01:05.89   0.7560A  -2.3797  -1.0650   1.4796  0.070  0.56  -0.895E+00     40.  0.135E+01  0.488E-01  0.124E-02  0.183E-03   60  11   1   17  12   0    7  14   3    0.000    0.000       0
+  0     88 00:01:06.65   0.7560A  -2.3797  -1.0715   1.4868  0.070  0.56  -0.892E+00     40.  0.135E+01  0.491E-01  0.125E-02  0.183E-03   60  11   1   17  12   0    7  14   3    0.000    0.000       0
+  0     89 00:01:07.41   0.7560A  -2.3798  -1.0746   1.4901  0.070  0.56  -0.890E+00     40.  0.135E+01  0.494E-01  0.126E-02  0.183E-03   60  10  19   17  12   0    7  14   3    0.000    0.000       0
+  0     90 00:01:08.16   0.7560A  -2.3799  -1.0745   1.4904  0.070  0.57  -0.889E+00     40.  0.135E+01  0.497E-01  0.126E-02  0.183E-03   60  10  19   17  12   0    7  14   3    0.000    0.000       0
+  0     91 00:01:08.92   0.7560A  -2.3800   1.0764   1.4883  0.070  0.57  -0.887E+00     40.  0.135E+01  0.500E-01  0.127E-02  0.184E-03   60  10  18    1   2  14    7  14   3    0.000    0.000       0
+  0     92 00:01:09.67   0.7560A  -2.3800   1.0833   1.4844  0.071  0.60  -0.886E+00     46.  0.135E+01  0.503E-01  0.127E-02  0.184E-03   60   9  16    1   2  14    7  14   3    0.000    0.000       0
+  0     93 00:01:10.43   0.7560A  -2.3801   1.0895   1.4794  0.071  0.61  -0.885E+00     46.  0.135E+01  0.506E-01  0.128E-02  0.184E-03   60   9  16    1   2  14    7  14   3    0.000    0.000       0
Index: palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc
===================================================================
--- palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc	(revision 4101)
+++ palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc	(revision 4102)
@@ -7,96 +7,96 @@
 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
 ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
-  1     93 00:01:10.43   0.7560A  -2.3801   1.0912   1.4817  0.071  0.11  -0.885E+00     46.  0.135E+01  0.506E-01  0.000E+00  0.000E+00   60   9  16    1   2  14    7  14   3    0.000    0.000       0
-  1     94 00:01:11.19   0.7560A  -2.3802   1.0968   1.4761  0.071  0.61  -0.884E+00     46.  0.135E+01  0.509E-01  0.128E-02  0.185E-03   60   9  15    1   2  14    7  14   3    0.000    0.000       0
-  1     95 00:01:11.94   0.7560A  -2.3803   1.1019   1.4704  0.071  0.62  -0.883E+00     46.  0.135E+01  0.511E-01  0.129E-02  0.185E-03   60   8  13    1   2  14    7  14   3    0.000    0.000       0
-  1     96 00:01:12.70   0.7560A  -2.3803   1.1066   1.4650  0.071  0.62  -0.883E+00     46.  0.135E+01  0.514E-01  0.129E-02  0.185E-03   60   8  13    1   2  14    7  14   3    0.000    0.000       0
-  1     97 00:01:13.45   0.7560A  -2.3804   1.1108   1.4601  0.071  0.63  -0.882E+00     46.  0.135E+01  0.516E-01  0.130E-02  0.186E-03   60   8  12    1   2  14    7  14   3    0.000    0.000       0
-  1     98 00:01:14.21   0.7560A  -2.3805   1.1146   1.4557  0.071  0.63  -0.882E+00     46.  0.135E+01  0.519E-01  0.130E-02  0.186E-03   60   8  12    1   2  14    7  14   3    0.000    0.000       0
-  1     99 00:01:14.97   0.7560A  -2.3806   1.1180   1.4520  0.071  0.63  -0.881E+00     46.  0.135E+01  0.521E-01  0.131E-02  0.187E-03   60   7  10    1   2  14    7  14   3    0.000    0.000       0
-  1    100 00:01:15.72   0.7560A  -2.3806   1.1209   1.4490  0.071  0.64  -0.881E+00     46.  0.135E+01  0.524E-01  0.131E-02  0.187E-03   60   7  10    1   2  14    7  14   3    0.000    0.000       0
-  1    101 00:01:16.48   0.7560A  -2.3807   1.1235   1.4468  0.071  0.65  -0.881E+00     48.  0.135E+01  0.526E-01  0.132E-02  0.188E-03   60   7   9    1   2  14    7  14   3    0.000    0.000       0
-  1    102 00:01:17.23   0.7560A  -2.3808  -1.1326   1.4453  0.071  0.66  -0.881E+00     48.  0.135E+01  0.528E-01  0.132E-02  0.188E-03   60   6   7   13  13   0    7  14   3    0.000    0.000       0
-  1    103 00:01:17.99   0.7560A  -2.3809  -1.1492   1.4445  0.071  0.66  -0.882E+00     48.  0.135E+01  0.530E-01  0.133E-02  0.189E-03   60   6   7   14  14   0    7  14   3    0.000    0.000       0
-  1    104 00:01:18.75   0.7560A  -2.3809  -1.1669   1.4444  0.071  0.66  -0.883E+00     48.  0.135E+01  0.532E-01  0.134E-02  0.190E-03   60   6   6   14  14   0    7  14   3    0.000    0.000       0
-  1    105 00:01:19.50   0.7560A  -2.3810  -1.1848   1.4447  0.071  0.68  -0.885E+00     50.  0.135E+01  0.534E-01  0.134E-02  0.190E-03   60   6   6   14  14   0    7  14   3    0.000    0.000       0
-  1    106 00:01:20.26   0.7560A  -2.3811  -1.2038   1.4454  0.071  0.68  -0.890E+00     50.  0.135E+01  0.536E-01  0.135E-02  0.191E-03   60   5   4   14  13   0    7  14   3    0.000    0.000       0
-  1    107 00:01:21.01   0.7560A  -2.3811  -1.2249   1.4465  0.071  0.68  -0.902E+00     50.  0.135E+01  0.538E-01  0.136E-02  0.192E-03   60   5   4   14  13   0    7  14   3    0.000    0.000       0
-  1    108 00:01:21.77   0.7560A  -2.3812  -1.2464   1.4477  0.071  0.69  -0.890E+00     50.  0.135E+01  0.540E-01  0.136E-02  0.193E-03   60   5   3   14  13   0    7  14   3    0.000    0.000       0
-  1    109 00:01:22.53   0.7560A  -2.3813  -1.2675   1.4489  0.071  0.69  -0.885E+00     50.  0.135E+01  0.542E-01  0.137E-02  0.193E-03   60   4   1   14  13   0    7  14   3    0.000    0.000       0
-  1    110 00:01:23.28   0.7560A  -2.3814  -1.2873   1.4501  0.071  0.69  -0.884E+00     50.  0.134E+01  0.543E-01  0.138E-02  0.194E-03   60   4   1   14  13   0    7  14   3    0.000    0.000       0
-  1    111 00:01:24.04   0.7560A  -2.3814  -1.3050   1.4512  0.071  0.69  -0.883E+00     50.  0.134E+01  0.545E-01  0.139E-02  0.194E-03   60   4   0   14  13   0    7  14   3    0.000    0.000       0
-  1    112 00:01:24.79   0.7560A  -2.3815  -1.3196   1.4522  0.071  0.70  -0.883E+00     50.  0.134E+01  0.547E-01  0.139E-02  0.194E-03   60   3  19   14  13   0    7  14   3    0.000    0.000       0
-  1    113 00:01:25.55   0.7560A  -2.3816  -1.3309   1.4530  0.071  0.70  -0.883E+00     50.  0.134E+01  0.548E-01  0.140E-02  0.194E-03   60   3  18   14  13   0    7  14   3    0.000    0.000       0
-  1    114 00:01:26.31   0.7560A  -2.3817  -1.3388   1.4537  0.071  0.71  -0.884E+00     52.  0.134E+01  0.550E-01  0.141E-02  0.194E-03   60   3  17   14  13   0    7  14   3    0.000    0.000       0
-  1    115 00:01:27.06   0.7560A  -2.3817  -1.3431   1.4543  0.071  0.71  -0.884E+00     52.  0.134E+01  0.551E-01  0.142E-02  0.194E-03   60   2  16   14  13   0    7  14   3    0.000    0.000       0
-  1    116 00:01:27.82   0.7560A  -2.3818  -1.3439   1.4564  0.071  0.71  -0.885E+00     52.  0.134E+01  0.552E-01  0.142E-02  0.194E-03   60   2  15   14  13   0    8  14   3    0.000    0.000       0
-  1    117 00:01:28.57   0.7560A  -2.3819  -1.3412   1.4598  0.071  0.71  -0.886E+00     52.  0.134E+01  0.553E-01  0.143E-02  0.193E-03   60   2  14   14  13   0    8  14   3    0.000    0.000       0
-  1    118 00:01:29.33   0.7560A  -2.3819  -1.3350   1.4627  0.071  0.71  -0.887E+00     52.  0.134E+01  0.554E-01  0.143E-02  0.192E-03   60   1  13   14  13   0    8  14   3    0.000    0.000       0
-  1    119 00:01:30.09   0.7560A  -2.3820  -1.3253   1.4650  0.071  0.64  -0.887E+00     40.  0.134E+01  0.555E-01  0.144E-02  0.191E-03   60   1  12   14  13   0    8  14   3    0.000    0.000       0
-  1    120 00:01:30.84   0.7560A  -2.3821  -1.3125   1.4668  0.071  0.64  -0.887E+00     40.  0.134E+01  0.556E-01  0.145E-02  0.190E-03   60   1  12   14  13   0    8  14   3    0.000    0.000       0
-  1    121 00:01:31.60   0.7560A  -2.3822  -1.2967   1.4680  0.071  0.64  -0.888E+00     40.  0.134E+01  0.557E-01  0.146E-02  0.189E-03   60   0  10   14  13   0    8  14   3    0.000    0.000       0
-  1    122 00:01:32.35   0.7560A  -2.3822  -1.2830   1.4685  0.071  0.64  -0.889E+00     40.  0.134E+01  0.558E-01  0.146E-02  0.187E-03   60   0   9   13  13   0    8  14   3    0.000    0.000       0
-  1    123 00:01:33.11   0.7560A  -2.3823  -1.2692   1.4685  0.071  0.63  -0.890E+00     40.  0.134E+01  0.558E-01  0.147E-02  0.186E-03   60   0   9   13  13   0    8  14   3    0.000    0.000       0
-  1    124 00:01:33.87   0.7560A  -2.3824  -1.2787   1.4680  0.071  0.63  -0.892E+00     40.  0.134E+01  0.559E-01  0.148E-02  0.184E-03   60   0   8    8  14   5    8  14   3    0.000    0.000       0
-  1    125 00:01:34.62   0.7560A  -2.3825  -1.2876   1.4669  0.071  0.63  -0.895E+00     40.  0.134E+01  0.559E-01  0.149E-02  0.183E-03   60  19   7    8  14   5    8  14   3    0.000    0.000       0
-  1    126 00:01:35.38   0.7550A  -2.3825  -1.2951   1.4654  0.071  0.64  -0.897E+00     42.  0.134E+01  0.560E-01  0.149E-02  0.181E-03   60  19   6    8  14   5    8  14   3    0.000    0.000       0
-  1    127 00:01:36.13   0.7550A  -2.3826  -1.3016   1.4634  0.071  0.63  -0.899E+00     42.  0.134E+01  0.560E-01  0.150E-02  0.179E-03   60  19   6    8  14   5    8  14   3    0.000    0.000       0
-  1    128 00:01:36.89   0.7550A  -2.3827  -1.3075   1.4613  0.071  0.62  -0.902E+00     40.  0.134E+01  0.560E-01  0.151E-02  0.178E-03   60  18   4    8  14   5    8  14   3    0.000    0.000       0
-  1    129 00:01:37.64   0.7550A  -2.3828  -1.3132   1.4592  0.070  0.62  -0.911E+00     40.  0.134E+01  0.560E-01  0.152E-02  0.176E-03   60  18   3    8  14   5    8  14   3    0.000    0.000       0
-  1    130 00:01:38.40   0.7550A  -2.3828  -1.3183   1.4574  0.070  0.61  -0.916E+00     40.  0.134E+01  0.560E-01  0.153E-02  0.174E-03   60  18   3    8  14   5    8  14   3    0.000    0.000       0
-  1    131 00:01:39.15   0.7550A  -2.3829  -1.3222   1.4560  0.070  0.61  -0.915E+00     40.  0.134E+01  0.559E-01  0.154E-02  0.173E-03   60  18   2    8  14   5    8  14   3    0.000    0.000       0
-  1    132 00:01:39.91   0.7550A  -2.3830  -1.3240   1.4559  0.070  0.61  -0.915E+00     40.  0.134E+01  0.559E-01  0.156E-02  0.171E-03   60  17   1    8  14   5    7  14   3    0.000    0.000       0
-  1    133 00:01:40.66   0.7550A  -2.3831  -1.3225   1.4575  0.070  0.60  -0.912E+00     40.  0.133E+01  0.558E-01  0.157E-02  0.170E-03   60  17   0    8  14   5    7  14   3    0.000    0.000       0
-  1    134 00:01:41.42   0.7550A  -2.3831  -1.3164   1.4593  0.070  0.60  -0.909E+00     40.  0.133E+01  0.557E-01  0.158E-02  0.169E-03   60  17   0    8  14   5    7  14   3    0.000    0.000       0
-  1    135 00:01:42.17   0.7550A  -2.3832  -1.3049   1.4608  0.070  0.61  -0.908E+00     42.  0.133E+01  0.556E-01  0.159E-02  0.168E-03   60  16  18    8  14   5    7  14   3    0.000    0.000       0
-  1    136 00:01:42.93   0.7550A  -2.3833  -1.2873   1.4618  0.070  0.61  -0.908E+00     42.  0.133E+01  0.555E-01  0.160E-02  0.167E-03   60  16  17    8  14   5    7  14   3    0.000    0.000       0
-  1    137 00:01:43.68   0.7550A  -2.3833  -1.2642   1.4638  0.070  0.62  -0.909E+00     44.  0.133E+01  0.554E-01  0.161E-02  0.166E-03   60  16  17    8  14   5    8  13   3    0.000    0.000       0
-  1    138 00:01:44.44   0.7550A  -2.3834  -1.2366   1.4694  0.070  0.62  -0.910E+00     44.  0.133E+01  0.552E-01  0.162E-02  0.165E-03   60  15  15    8  14   5    8  13   3    0.000    0.000       0
-  1    139 00:01:45.19   0.7550A  -2.3835  -1.2220   1.4730  0.070  0.61  -0.916E+00     44.  0.133E+01  0.550E-01  0.163E-02  0.164E-03   60  15  14    6  13   5    8  13   3    0.000    0.000       0
-  1    140 00:01:45.95   0.7550A  -2.3836  -1.2184   1.4743  0.070  0.61  -0.917E+00     44.  0.133E+01  0.549E-01  0.164E-02  0.163E-03   60  15  14    6  13   5    8  13   3    0.000    0.000       0
-  1    141 00:01:46.70   0.7550A  -2.3836  -1.2098   1.4726  0.069  0.61  -0.912E+00     44.  0.133E+01  0.546E-01  0.165E-02  0.162E-03   60  14  12    6  13   5    8  13   3    0.000    0.000       0
-  1    142 00:01:47.46   0.7550A  -2.3837  -1.2086   1.4679  0.069  0.61  -0.912E+00     44.  0.133E+01  0.544E-01  0.165E-02  0.161E-03   60  14  11   17   8  19    8  13   3    0.000    0.000       0
-  1    143 00:01:48.21   0.7550A  -2.3838  -1.2028   1.4643  0.069  0.60  -0.913E+00     44.  0.133E+01  0.542E-01  0.166E-02  0.160E-03   60  14  11   17   8  19    8  14   3    0.000    0.000       0
-  1    144 00:01:48.97   0.7550A  -2.3839  -1.1926   1.4664  0.069  0.60  -0.915E+00     44.  0.133E+01  0.539E-01  0.167E-02  0.159E-03   60  13   9   17   8  19    8  14   3    0.000    0.000       0
-  1    145 00:01:49.72   0.7550A  -2.3839  -1.1785   1.4689  0.069  0.60  -0.917E+00     44.  0.133E+01  0.537E-01  0.167E-02  0.158E-03   60  13   8   17   8  19    8  14   3    0.000    0.000       0
-  1    146 00:01:50.48   0.7550A  -2.3840  -1.1750   1.4719  0.069  0.61  -0.917E+00     46.  0.133E+01  0.534E-01  0.168E-02  0.157E-03   60  13   8   16   8  19    8  14   3    0.000    0.000       0
-  1    147 00:01:51.23   0.7550A  -2.3841  -1.1698   1.4753  0.069  0.61  -0.917E+00     46.  0.133E+01  0.531E-01  0.168E-02  0.157E-03   60  13   7   16   8  19    8  14   3    0.000    0.000       0
-  1    148 00:01:51.99   0.7550A  -2.3841  -1.1592   1.4790  0.068  0.60  -0.918E+00     46.  0.133E+01  0.528E-01  0.168E-02  0.156E-03   60  12   6   16   8  19    8  14   3    0.000    0.000       0
-  1    149 00:01:52.74   0.7550A  -2.3842  -1.1440   1.4830  0.068  0.60  -0.920E+00     46.  0.133E+01  0.525E-01  0.169E-02  0.156E-03   60  12   5   15   8  19    8  14   3    0.000    0.000       0
-  1    150 00:01:53.50   0.7550A  -2.3843  -1.1376   1.4870  0.068  0.60  -0.922E+00     46.  0.132E+01  0.522E-01  0.169E-02  0.155E-03   60  12   5   17   7  18    8  14   3    0.000    0.000       0
-  1    151 00:01:54.25   0.7550A  -2.3844  -1.1399   1.4908  0.068  0.60  -0.925E+00     46.  0.132E+01  0.519E-01  0.169E-02  0.155E-03   60  11   3   17   7  18    8  14   3    0.000    0.000       0
-  1    152 00:01:55.01   0.7550A  -2.3844  -1.1357   1.4942  0.068  0.61  -0.928E+00     48.  0.132E+01  0.515E-01  0.170E-02  0.155E-03   60  11   2   17   7  18    8  14   3    0.000    0.000       0
-  1    153 00:01:55.76   0.7550A  -2.3845  -1.1241   1.4971  0.067  0.61  -0.929E+00     48.  0.132E+01  0.512E-01  0.170E-02  0.155E-03   60  11   2   17   7  18    8  14   3    0.000    0.000       0
-  1    154 00:01:56.52   0.7550A  -2.3846  -1.1043   1.4991  0.067  0.61  -0.930E+00     48.  0.132E+01  0.508E-01  0.170E-02  0.155E-03   60  11   1   17   7  18    8  14   3    0.000    0.000       0
-  1    155 00:01:57.27   0.7550A  -2.3847  -1.1019   1.4997  0.067  0.61  -0.931E+00     48.  0.132E+01  0.505E-01  0.170E-02  0.155E-03   60  10   0   18   7  17    8  14   3    0.000    0.000       0
-  1    156 00:01:58.03   0.7550A  -2.3847  -1.0973   1.4989  0.067  0.61  -0.933E+00     48.  0.132E+01  0.502E-01  0.169E-02  0.155E-03   60  10  19   18   7  17    8  14   3    0.000    0.000       0
-  1    157 00:01:58.78   0.7550A  -2.3848  -1.1023   1.4962  0.067  0.61  -0.935E+00     48.  0.132E+01  0.498E-01  0.169E-02  0.155E-03   60  10  19   16   5  18    8  14   3    0.000    0.000       0
-  1    158 00:01:59.54   0.7550A  -2.3849  -1.1089   1.4915  0.067  0.61  -0.939E+00     48.  0.132E+01  0.495E-01  0.169E-02  0.155E-03   60   9  17   16   5  18    8  14   3    0.000    0.000       0
-  1    159 00:02:00.29   0.7550A  -2.3849  -1.1068   1.4844  0.066  0.61  -0.942E+00     48.  0.132E+01  0.491E-01  0.169E-02  0.155E-03   60   9  16   16   5  18    8  14   3    0.000    0.000       0
-  1    160 00:02:01.05   0.7550A  -2.3850  -1.0961   1.4745  0.066  0.62  -0.944E+00     48.  0.132E+01  0.488E-01  0.169E-02  0.155E-03   60   9  16   16   5  18    8  14   3    0.000    0.000       0
-  1    161 00:02:01.80   0.7550A  -2.3851  -1.0773   1.4617  0.066  0.62  -0.947E+00     48.  0.132E+01  0.484E-01  0.169E-02  0.154E-03   60   8  14   16   5  18    8  14   3    0.000    0.000       0
-  1    162 00:02:02.56   0.7550A  -2.3852  -1.0860   1.4457  0.066  0.62  -0.948E+00     48.  0.132E+01  0.481E-01  0.169E-02  0.154E-03   60   8  13   16   5  17    8  14   3    0.000    0.000       0
-  1    163 00:02:03.31   0.7550A  -2.3852  -1.1032   1.4263  0.066  0.62  -0.951E+00     48.  0.132E+01  0.477E-01  0.168E-02  0.153E-03   60   8  13   16   5  17    8  14   3    0.000    0.000       0
-  1    164 00:02:04.07   0.7550A  -2.3853  -1.1080   1.4101  0.066  0.63  -0.957E+00     48.  0.132E+01  0.474E-01  0.168E-02  0.152E-03   60   7  11   16   5  17    9  14   3    0.000    0.000       0
-  1    165 00:02:04.82   0.7550A  -2.3854  -1.0998   1.4185  0.066  0.63  -0.962E+00     48.  0.132E+01  0.470E-01  0.167E-02  0.151E-03   60   7  10   16   5  17    8  14   2    0.000    0.000       0
-  1    166 00:02:05.58   0.7550A  -2.3854  -1.0790   1.4267  0.066  0.63  -0.962E+00     48.  0.131E+01  0.467E-01  0.167E-02  0.150E-03   60   7  10   16   5  17    8  14   2    0.000    0.000       0
-  1    167 00:02:06.33   0.7550A  -2.3855  -1.0532   1.4328  0.065  0.64  -0.959E+00     48.  0.131E+01  0.463E-01  0.166E-02  0.149E-03   60   6   8   16   5  16    8  14   2    0.000    0.000       0
-  1    168 00:02:07.09   0.7550A  -2.3856  -1.0680   1.4363  0.065  0.64  -0.958E+00     48.  0.131E+01  0.460E-01  0.166E-02  0.147E-03   60   6   7   16   5  16    8  14   2    0.000    0.000       0
-  1    169 00:02:07.84   0.7550A  -2.3857  -1.0709   1.4369  0.065  0.64  -0.957E+00     48.  0.131E+01  0.456E-01  0.165E-02  0.144E-03   60   6   7   21  11  18    8  14   2    0.000    0.000       0
-  1    170 00:02:08.60   0.7540A  -2.3857  -1.0774   1.4344  0.065  0.64  -0.958E+00     48.  0.131E+01  0.453E-01  0.164E-02  0.142E-03   60   6   7   22  11  17    8  14   2    0.000    0.000       0
-  1    171 00:02:09.35   0.7540A  -2.3858  -1.0936   1.4288  0.065  0.65  -0.960E+00     48.  0.131E+01  0.450E-01  0.163E-02  0.139E-03   60   5   5   22  11  17    8  14   2    0.000    0.000       0
-  1    172 00:02:10.11   0.7540A  -2.3859  -1.1017   1.4205  0.065  0.65  -0.963E+00     48.  0.131E+01  0.446E-01  0.162E-02  0.136E-03   60   5   4   22  11  17    8  14   2    0.000    0.000       0
-  1    173 00:02:10.86   0.7540A  -2.3860  -1.1152   1.4096  0.065  0.65  -0.968E+00     48.  0.131E+01  0.443E-01  0.161E-02  0.133E-03   60   5   4   21  10  17    8  14   2    0.000    0.000       0
-  1    174 00:02:11.61   0.7540A  -2.3860  -1.1221   1.3966  0.064  0.65  -0.975E+00     48.  0.131E+01  0.439E-01  0.160E-02  0.130E-03   60   4   2   21  10  17    8  14   2    0.000    0.000       0
-  1    175 00:02:12.37   0.7540A  -2.3861  -1.1173   1.3876  0.064  0.67  -0.977E+00     50.  0.131E+01  0.436E-01  0.159E-02  0.127E-03   60   4   1   21  10  17    9  14   2    0.000    0.000       0
-  1    176 00:02:13.12   0.7540A  -2.3862  -1.1054   1.3955  0.064  0.67  -0.977E+00     50.  0.131E+01  0.433E-01  0.158E-02  0.124E-03   60   4   1   22  10  16    9  14   2    0.000    0.000       0
-  1    177 00:02:13.88   0.7540A  -2.3862  -1.1083   1.3996  0.064  0.67  -0.983E+00     50.  0.131E+01  0.429E-01  0.157E-02  0.122E-03   60   4   1   22  10  16    9  14   2    0.000    0.000       0
-  1    178 00:02:14.63   0.7540A  -2.3863  -1.0999   1.4001  0.064  0.67  -0.100E+01     50.  0.131E+01  0.426E-01  0.156E-02  0.119E-03   60   3  19   22  10  16    9  14   2    0.000    0.000       0
-  1    179 00:02:15.38   0.7540A  -2.3864  -1.0791   1.3975  0.063  0.67  -0.103E+01     50.  0.131E+01  0.423E-01  0.156E-02  0.117E-03   60   3  18   21   9  16    9  14   2    0.000    0.000       0
-  1    180 00:02:16.14   0.7540A  -2.3865  -1.0799   1.3925  0.063  0.68  -0.993E+00     50.  0.131E+01  0.420E-01  0.155E-02  0.115E-03   60   3  18   21   9  16    9  14   2    0.000    0.000       0
-  1    181 00:02:16.89   0.7540A  -2.3865  -1.0667   1.3857  0.063  0.68  -0.101E+01     50.  0.131E+01  0.417E-01  0.154E-02  0.113E-03   60   2  16   21   9  16    9  14   2    0.000    0.000       0
-  1    182 00:02:17.65   0.7540A  -2.3866  -1.0471   1.3775  0.063  0.68  -0.983E+00     50.  0.131E+01  0.414E-01  0.154E-02  0.112E-03   60   2  15   22   9  15    9  14   2    0.000    0.000       0
-  1    183 00:02:18.40   0.7540A  -2.3867  -1.0440   1.3686  0.063  0.68  -0.964E+00     50.  0.131E+01  0.411E-01  0.153E-02  0.111E-03   60   2  15   22   9  15    9  14   2    0.000    0.000       0
-  1    184 00:02:19.15   0.7540A  -2.3867  -1.0475   1.3589  0.063  0.69  -0.956E+00     52.  0.130E+01  0.408E-01  0.153E-02  0.110E-03   60   1  13   18   1  12    9  14   2    0.000    0.000       0
-  1    185 00:02:19.91   0.7540A  -2.3868  -1.0550   1.3488  0.063  0.69  -0.951E+00     52.  0.130E+01  0.406E-01  0.152E-02  0.109E-03   60   1  12   18   1  12    9  14   2    0.000    0.000       0
-  1    186 00:02:20.66   0.7540A  -2.3869  -1.0554   1.3380  0.063  0.69  -0.946E+00     52.  0.130E+01  0.403E-01  0.152E-02  0.108E-03   60   1  12   18   1  12    9  14   2    0.000    0.000       0
+  1     93 00:01:10.43   0.7560A  -2.3801   1.0895   1.4794  0.071  0.11  -0.885E+00     46.  0.135E+01  0.506E-01  0.000E+00  0.000E+00   60   9  16    1   2  14    7  14   3    0.000    0.000       0
+  1     94 00:01:11.19   0.7560A  -2.3802   1.0952   1.4738  0.071  0.61  -0.884E+00     46.  0.135E+01  0.509E-01  0.128E-02  0.185E-03   60   9  15    1   2  14    7  14   3    0.000    0.000       0
+  1     95 00:01:11.94   0.7560A  -2.3803   1.1004   1.4683  0.071  0.62  -0.883E+00     46.  0.135E+01  0.511E-01  0.129E-02  0.185E-03   60   8  13    1   2  14    7  14   3    0.000    0.000       0
+  1     96 00:01:12.70   0.7560A  -2.3803   1.1051   1.4629  0.071  0.62  -0.883E+00     46.  0.135E+01  0.514E-01  0.129E-02  0.185E-03   60   8  13    1   2  14    7  14   3    0.000    0.000       0
+  1     97 00:01:13.45   0.7560A  -2.3804   1.1093   1.4580  0.071  0.63  -0.882E+00     46.  0.135E+01  0.516E-01  0.130E-02  0.186E-03   60   8  12    1   2  14    7  14   3    0.000    0.000       0
+  1     98 00:01:14.21   0.7560A  -2.3805   1.1132   1.4537  0.071  0.63  -0.882E+00     46.  0.135E+01  0.519E-01  0.130E-02  0.186E-03   60   8  12    1   2  14    7  14   3    0.000    0.000       0
+  1     99 00:01:14.97   0.7560A  -2.3806   1.1166   1.4501  0.071  0.63  -0.882E+00     46.  0.135E+01  0.521E-01  0.131E-02  0.187E-03   60   7  10    1   2  14    7  14   3    0.000    0.000       0
+  1    100 00:01:15.72   0.7560A  -2.3806   1.1196   1.4471  0.071  0.64  -0.881E+00     46.  0.135E+01  0.523E-01  0.131E-02  0.187E-03   60   7  10    1   2  14    7  14   3    0.000    0.000       0
+  1    101 00:01:16.48   0.7560A  -2.3807   1.1223   1.4449  0.071  0.65  -0.881E+00     48.  0.135E+01  0.526E-01  0.132E-02  0.188E-03   60   7   9    1   2  14    7  14   3    0.000    0.000       0
+  1    102 00:01:17.23   0.7560A  -2.3808  -1.1319   1.4435  0.071  0.66  -0.882E+00     48.  0.135E+01  0.528E-01  0.132E-02  0.188E-03   60   6   7   13  13   0    7  14   3    0.000    0.000       0
+  1    103 00:01:17.99   0.7560A  -2.3809  -1.1487   1.4427  0.071  0.66  -0.882E+00     48.  0.135E+01  0.530E-01  0.133E-02  0.189E-03   60   6   7   14  14   0    7  14   3    0.000    0.000       0
+  1    104 00:01:18.75   0.7560A  -2.3809  -1.1665   1.4425  0.071  0.66  -0.883E+00     48.  0.135E+01  0.532E-01  0.134E-02  0.190E-03   60   6   6   14  14   0    7  14   3    0.000    0.000       0
+  1    105 00:01:19.50   0.7560A  -2.3810  -1.1843   1.4429  0.071  0.67  -0.885E+00     48.  0.135E+01  0.534E-01  0.134E-02  0.190E-03   60   6   6   14  14   0    7  14   3    0.000    0.000       0
+  1    106 00:01:20.26   0.7560A  -2.3811  -1.2034   1.4436  0.071  0.68  -0.890E+00     50.  0.135E+01  0.536E-01  0.135E-02  0.191E-03   60   5   4   14  13   0    7  14   3    0.000    0.000       0
+  1    107 00:01:21.01   0.7560A  -2.3811  -1.2245   1.4446  0.071  0.68  -0.900E+00     50.  0.135E+01  0.538E-01  0.136E-02  0.192E-03   60   5   4   14  13   0    7  14   3    0.000    0.000       0
+  1    108 00:01:21.77   0.7560A  -2.3812  -1.2460   1.4458  0.071  0.69  -0.890E+00     50.  0.135E+01  0.540E-01  0.136E-02  0.193E-03   60   5   3   14  13   0    7  14   3    0.000    0.000       0
+  1    109 00:01:22.53   0.7560A  -2.3813  -1.2671   1.4471  0.071  0.69  -0.885E+00     50.  0.135E+01  0.541E-01  0.137E-02  0.193E-03   60   4   1   14  13   0    7  14   3    0.000    0.000       0
+  1    110 00:01:23.28   0.7560A  -2.3814  -1.2870   1.4483  0.071  0.69  -0.884E+00     50.  0.134E+01  0.543E-01  0.138E-02  0.194E-03   60   4   1   14  13   0    7  14   3    0.000    0.000       0
+  1    111 00:01:24.04   0.7560A  -2.3814  -1.3046   1.4494  0.071  0.69  -0.883E+00     50.  0.134E+01  0.545E-01  0.138E-02  0.194E-03   60   4   0   14  13   0    7  14   3    0.000    0.000       0
+  1    112 00:01:24.79   0.7560A  -2.3815  -1.3193   1.4504  0.071  0.70  -0.883E+00     50.  0.134E+01  0.546E-01  0.139E-02  0.194E-03   60   3  19   14  13   0    7  14   3    0.000    0.000       0
+  1    113 00:01:25.55   0.7560A  -2.3816  -1.3307   1.4512  0.071  0.70  -0.883E+00     50.  0.134E+01  0.548E-01  0.140E-02  0.194E-03   60   3  18   14  13   0    7  14   3    0.000    0.000       0
+  1    114 00:01:26.31   0.7560A  -2.3817  -1.3385   1.4519  0.071  0.71  -0.884E+00     52.  0.134E+01  0.549E-01  0.141E-02  0.194E-03   60   3  17   14  13   0    7  14   3    0.000    0.000       0
+  1    115 00:01:27.06   0.7560A  -2.3817  -1.3428   1.4525  0.071  0.71  -0.884E+00     52.  0.134E+01  0.550E-01  0.141E-02  0.194E-03   60   2  16   14  13   0    7  14   3    0.000    0.000       0
+  1    116 00:01:27.82   0.7560A  -2.3818  -1.3436   1.4542  0.071  0.71  -0.885E+00     52.  0.134E+01  0.552E-01  0.142E-02  0.194E-03   60   2  15   14  13   0    8  14   3    0.000    0.000       0
+  1    117 00:01:28.57   0.7560A  -2.3819  -1.3409   1.4575  0.071  0.71  -0.886E+00     52.  0.134E+01  0.553E-01  0.143E-02  0.193E-03   60   2  14   14  13   0    8  14   3    0.000    0.000       0
+  1    118 00:01:29.33   0.7560A  -2.3819  -1.3347   1.4604  0.071  0.71  -0.887E+00     52.  0.134E+01  0.554E-01  0.143E-02  0.192E-03   60   1  13   14  13   0    8  14   3    0.000    0.000       0
+  1    119 00:01:30.09   0.7560A  -2.3820  -1.3251   1.4628  0.071  0.64  -0.887E+00     40.  0.134E+01  0.555E-01  0.144E-02  0.191E-03   60   1  12   14  13   0    8  14   3    0.000    0.000       0
+  1    120 00:01:30.84   0.7560A  -2.3821  -1.3123   1.4646  0.071  0.64  -0.887E+00     40.  0.134E+01  0.556E-01  0.145E-02  0.190E-03   60   1  12   14  13   0    8  14   3    0.000    0.000       0
+  1    121 00:01:31.60   0.7560A  -2.3822  -1.2966   1.4658  0.071  0.64  -0.888E+00     40.  0.134E+01  0.557E-01  0.145E-02  0.189E-03   60   0  10   14  13   0    8  14   3    0.000    0.000       0
+  1    122 00:01:32.35   0.7560A  -2.3822  -1.2829   1.4664  0.071  0.64  -0.889E+00     40.  0.134E+01  0.557E-01  0.146E-02  0.187E-03   60   0   9   13  13   0    8  14   3    0.000    0.000       0
+  1    123 00:01:33.11   0.7560A  -2.3823  -1.2691   1.4664  0.071  0.63  -0.890E+00     40.  0.134E+01  0.558E-01  0.147E-02  0.186E-03   60   0   9   13  13   0    8  14   3    0.000    0.000       0
+  1    124 00:01:33.87   0.7560A  -2.3824  -1.2773   1.4659  0.071  0.63  -0.892E+00     40.  0.134E+01  0.558E-01  0.148E-02  0.184E-03   60   0   8    8  14   5    8  14   3    0.000    0.000       0
+  1    125 00:01:34.62   0.7560A  -2.3825  -1.2862   1.4649  0.071  0.63  -0.895E+00     40.  0.134E+01  0.559E-01  0.148E-02  0.183E-03   60  19   7    8  14   5    8  14   3    0.000    0.000       0
+  1    126 00:01:35.38   0.7550A  -2.3825  -1.2937   1.4634  0.071  0.64  -0.897E+00     42.  0.134E+01  0.559E-01  0.149E-02  0.181E-03   60  19   6    8  14   5    8  14   3    0.000    0.000       0
+  1    127 00:01:36.13   0.7550A  -2.3826  -1.3002   1.4615  0.071  0.63  -0.899E+00     42.  0.134E+01  0.559E-01  0.150E-02  0.179E-03   60  19   6    8  14   5    8  14   3    0.000    0.000       0
+  1    128 00:01:36.89   0.7550A  -2.3827  -1.3062   1.4594  0.071  0.62  -0.902E+00     40.  0.134E+01  0.559E-01  0.151E-02  0.177E-03   60  18   4    8  14   5    8  14   3    0.000    0.000       0
+  1    129 00:01:37.64   0.7550A  -2.3828  -1.3119   1.4574  0.070  0.62  -0.911E+00     40.  0.134E+01  0.559E-01  0.152E-02  0.176E-03   60  18   3    8  14   5    8  14   3    0.000    0.000       0
+  1    130 00:01:38.40   0.7550A  -2.3828  -1.3170   1.4556  0.070  0.61  -0.916E+00     40.  0.134E+01  0.559E-01  0.153E-02  0.174E-03   60  18   3    8  14   5    8  14   3    0.000    0.000       0
+  1    131 00:01:39.15   0.7550A  -2.3829  -1.3209   1.4542  0.070  0.61  -0.915E+00     40.  0.134E+01  0.559E-01  0.154E-02  0.173E-03   60  18   2    8  14   5    8  14   3    0.000    0.000       0
+  1    132 00:01:39.91   0.7550A  -2.3830  -1.3227   1.4543  0.070  0.61  -0.915E+00     40.  0.134E+01  0.558E-01  0.155E-02  0.171E-03   60  17   1    8  14   5    7  14   3    0.000    0.000       0
+  1    133 00:01:40.66   0.7550A  -2.3831  -1.3213   1.4559  0.070  0.60  -0.912E+00     40.  0.133E+01  0.558E-01  0.156E-02  0.170E-03   60  17   0    8  14   5    7  14   3    0.000    0.000       0
+  1    134 00:01:41.42   0.7550A  -2.3831  -1.3153   1.4577  0.070  0.60  -0.910E+00     40.  0.133E+01  0.557E-01  0.158E-02  0.169E-03   60  17   0    8  14   5    7  14   3    0.000    0.000       0
+  1    135 00:01:42.17   0.7550A  -2.3832  -1.3038   1.4592  0.070  0.61  -0.908E+00     42.  0.133E+01  0.556E-01  0.159E-02  0.168E-03   60  16  18    8  14   5    7  14   3    0.000    0.000       0
+  1    136 00:01:42.93   0.7550A  -2.3833  -1.2864   1.4601  0.070  0.61  -0.908E+00     42.  0.133E+01  0.554E-01  0.160E-02  0.167E-03   60  16  17    8  14   5    7  14   3    0.000    0.000       0
+  1    137 00:01:43.68   0.7550A  -2.3833  -1.2634   1.4626  0.070  0.62  -0.909E+00     44.  0.133E+01  0.553E-01  0.161E-02  0.166E-03   60  16  17    8  14   5    8  13   3    0.000    0.000       0
+  1    138 00:01:44.44   0.7550A  -2.3834  -1.2359   1.4682  0.070  0.62  -0.910E+00     44.  0.133E+01  0.552E-01  0.162E-02  0.165E-03   60  15  15    8  14   5    8  13   3    0.000    0.000       0
+  1    139 00:01:45.19   0.7550A  -2.3835  -1.2203   1.4717  0.070  0.61  -0.915E+00     44.  0.133E+01  0.550E-01  0.163E-02  0.164E-03   60  15  14    6  13   5    8  13   3    0.000    0.000       0
+  1    140 00:01:45.95   0.7550A  -2.3836  -1.2167   1.4729  0.070  0.61  -0.917E+00     44.  0.133E+01  0.548E-01  0.164E-02  0.163E-03   60  15  14    6  13   5    8  13   3    0.000    0.000       0
+  1    141 00:01:46.70   0.7550A  -2.3836  -1.2094   1.4712  0.069  0.61  -0.912E+00     44.  0.133E+01  0.546E-01  0.164E-02  0.162E-03   60  14  12   17   8  19    8  13   3    0.000    0.000       0
+  1    142 00:01:47.46   0.7550A  -2.3837  -1.2085   1.4664  0.069  0.61  -0.912E+00     44.  0.133E+01  0.544E-01  0.165E-02  0.161E-03   60  14  11   17   8  19    8  13   3    0.000    0.000       0
+  1    143 00:01:48.21   0.7550A  -2.3838  -1.2026   1.4624  0.069  0.60  -0.913E+00     44.  0.133E+01  0.541E-01  0.166E-02  0.160E-03   60  14  11   17   8  19    8  14   3    0.000    0.000       0
+  1    144 00:01:48.97   0.7550A  -2.3839  -1.1923   1.4644  0.069  0.60  -0.915E+00     44.  0.133E+01  0.539E-01  0.166E-02  0.159E-03   60  13   9   17   8  19    8  14   3    0.000    0.000       0
+  1    145 00:01:49.72   0.7550A  -2.3839  -1.1783   1.4669  0.069  0.60  -0.917E+00     44.  0.133E+01  0.536E-01  0.167E-02  0.158E-03   60  13   8   17   8  19    8  14   3    0.000    0.000       0
+  1    146 00:01:50.48   0.7550A  -2.3840  -1.1749   1.4698  0.069  0.61  -0.917E+00     46.  0.133E+01  0.533E-01  0.167E-02  0.157E-03   60  13   8   16   8  19    8  14   3    0.000    0.000       0
+  1    147 00:01:51.23   0.7550A  -2.3841  -1.1696   1.4732  0.069  0.61  -0.917E+00     46.  0.133E+01  0.531E-01  0.168E-02  0.157E-03   60  13   7   16   8  19    8  14   3    0.000    0.000       0
+  1    148 00:01:51.99   0.7550A  -2.3841  -1.1589   1.4768  0.068  0.60  -0.918E+00     46.  0.133E+01  0.528E-01  0.168E-02  0.156E-03   60  12   6   16   8  19    8  14   3    0.000    0.000       0
+  1    149 00:01:52.74   0.7550A  -2.3842  -1.1438   1.4807  0.068  0.60  -0.920E+00     46.  0.133E+01  0.524E-01  0.168E-02  0.156E-03   60  12   5   15   8  19    8  14   3    0.000    0.000       0
+  1    150 00:01:53.50   0.7550A  -2.3843  -1.1373   1.4846  0.068  0.60  -0.922E+00     46.  0.132E+01  0.521E-01  0.169E-02  0.155E-03   60  12   5   17   7  18    8  14   3    0.000    0.000       0
+  1    151 00:01:54.25   0.7550A  -2.3844  -1.1397   1.4884  0.068  0.60  -0.925E+00     46.  0.132E+01  0.518E-01  0.169E-02  0.155E-03   60  11   3   17   7  18    8  14   3    0.000    0.000       0
+  1    152 00:01:55.01   0.7550A  -2.3844  -1.1355   1.4917  0.068  0.61  -0.928E+00     48.  0.132E+01  0.515E-01  0.169E-02  0.155E-03   60  11   2   17   7  18    8  14   3    0.000    0.000       0
+  1    153 00:01:55.76   0.7550A  -2.3845  -1.1239   1.4945  0.067  0.61  -0.929E+00     48.  0.132E+01  0.511E-01  0.170E-02  0.155E-03   60  11   2   17   7  18    8  14   3    0.000    0.000       0
+  1    154 00:01:56.52   0.7550A  -2.3846  -1.1042   1.4964  0.067  0.61  -0.930E+00     48.  0.132E+01  0.508E-01  0.170E-02  0.155E-03   60  11   1   17   7  18    8  14   3    0.000    0.000       0
+  1    155 00:01:57.27   0.7550A  -2.3847  -1.1017   1.4970  0.067  0.61  -0.931E+00     48.  0.132E+01  0.505E-01  0.169E-02  0.155E-03   60  10   0   18   7  17    8  14   3    0.000    0.000       0
+  1    156 00:01:58.03   0.7550A  -2.3847  -1.0972   1.4960  0.067  0.61  -0.933E+00     48.  0.132E+01  0.501E-01  0.169E-02  0.155E-03   60  10  19   18   7  17    8  14   3    0.000    0.000       0
+  1    157 00:01:58.78   0.7550A  -2.3848  -1.1025   1.4933  0.067  0.61  -0.935E+00     48.  0.132E+01  0.498E-01  0.169E-02  0.155E-03   60  10  19   16   5  18    8  14   3    0.000    0.000       0
+  1    158 00:01:59.54   0.7550A  -2.3849  -1.1090   1.4886  0.067  0.61  -0.939E+00     48.  0.132E+01  0.494E-01  0.169E-02  0.155E-03   60   9  17   16   5  18    8  14   3    0.000    0.000       0
+  1    159 00:02:00.29   0.7550A  -2.3849  -1.1069   1.4814  0.066  0.61  -0.942E+00     48.  0.132E+01  0.491E-01  0.169E-02  0.155E-03   60   9  16   16   5  18    8  14   3    0.000    0.000       0
+  1    160 00:02:01.05   0.7550A  -2.3850  -1.0961   1.4715  0.066  0.62  -0.944E+00     48.  0.132E+01  0.487E-01  0.169E-02  0.155E-03   60   9  16   16   5  18    8  14   3    0.000    0.000       0
+  1    161 00:02:01.80   0.7550A  -2.3851  -1.0772   1.4587  0.066  0.62  -0.947E+00     48.  0.132E+01  0.484E-01  0.169E-02  0.154E-03   60   8  14   16   5  18    8  14   3    0.000    0.000       0
+  1    162 00:02:02.56   0.7550A  -2.3852  -1.0857   1.4427  0.066  0.62  -0.948E+00     48.  0.132E+01  0.480E-01  0.168E-02  0.154E-03   60   8  13   16   5  17    8  14   3    0.000    0.000       0
+  1    163 00:02:03.31   0.7550A  -2.3852  -1.1029   1.4233  0.066  0.62  -0.951E+00     48.  0.132E+01  0.477E-01  0.168E-02  0.153E-03   60   8  13   16   5  17    8  14   3    0.000    0.000       0
+  1    164 00:02:04.07   0.7550A  -2.3853  -1.1078   1.4090  0.066  0.63  -0.957E+00     48.  0.132E+01  0.473E-01  0.168E-02  0.152E-03   60   7  11   16   5  17    8  14   2    0.000    0.000       0
+  1    165 00:02:04.82   0.7550A  -2.3854  -1.0995   1.4187  0.066  0.63  -0.963E+00     48.  0.132E+01  0.470E-01  0.167E-02  0.151E-03   60   7  10   16   5  17    8  14   2    0.000    0.000       0
+  1    166 00:02:05.58   0.7550A  -2.3854  -1.0787   1.4268  0.066  0.63  -0.962E+00     48.  0.131E+01  0.466E-01  0.166E-02  0.150E-03   60   7  10   16   5  17    8  14   2    0.000    0.000       0
+  1    167 00:02:06.33   0.7550A  -2.3855  -1.0526   1.4328  0.065  0.64  -0.959E+00     48.  0.131E+01  0.463E-01  0.166E-02  0.149E-03   60   6   8   16   5  16    8  14   2    0.000    0.000       0
+  1    168 00:02:07.09   0.7550A  -2.3856  -1.0675   1.4363  0.065  0.64  -0.958E+00     48.  0.131E+01  0.459E-01  0.165E-02  0.147E-03   60   6   7   16   5  16    8  14   2    0.000    0.000       0
+  1    169 00:02:07.84   0.7550A  -2.3857  -1.0679   1.4368  0.065  0.64  -0.958E+00     48.  0.131E+01  0.456E-01  0.165E-02  0.144E-03   60   6   7   21  11  18    8  14   2    0.000    0.000       0
+  1    170 00:02:08.60   0.7540A  -2.3857  -1.0746   1.4343  0.065  0.64  -0.958E+00     48.  0.131E+01  0.453E-01  0.164E-02  0.142E-03   60   6   7   22  11  17    8  14   2    0.000    0.000       0
+  1    171 00:02:09.35   0.7540A  -2.3858  -1.0912   1.4288  0.065  0.65  -0.960E+00     48.  0.131E+01  0.449E-01  0.163E-02  0.139E-03   60   5   5   22  11  17    8  14   2    0.000    0.000       0
+  1    172 00:02:10.11   0.7540A  -2.3859  -1.0997   1.4204  0.065  0.65  -0.963E+00     48.  0.131E+01  0.446E-01  0.162E-02  0.136E-03   60   5   4   22  11  17    8  14   2    0.000    0.000       0
+  1    173 00:02:10.86   0.7540A  -2.3860  -1.1130   1.4095  0.064  0.65  -0.968E+00     48.  0.131E+01  0.442E-01  0.160E-02  0.133E-03   60   5   4   21  10  17    8  14   2    0.000    0.000       0
+  1    174 00:02:11.61   0.7540A  -2.3860  -1.1206   1.3965  0.064  0.65  -0.975E+00     48.  0.131E+01  0.439E-01  0.159E-02  0.130E-03   60   4   2   21  10  17    8  14   2    0.000    0.000       0
+  1    175 00:02:12.37   0.7540A  -2.3861  -1.1167   1.3873  0.064  0.67  -0.977E+00     50.  0.131E+01  0.436E-01  0.159E-02  0.127E-03   60   4   1   21  10  17    9  14   2    0.000    0.000       0
+  1    176 00:02:13.12   0.7540A  -2.3862  -1.1043   1.3951  0.064  0.67  -0.977E+00     50.  0.131E+01  0.432E-01  0.158E-02  0.124E-03   60   4   1   22  10  16    9  14   2    0.000    0.000       0
+  1    177 00:02:13.88   0.7540A  -2.3862  -1.1081   1.3991  0.064  0.67  -0.982E+00     50.  0.131E+01  0.429E-01  0.157E-02  0.122E-03   60   4   1   22  10  16    9  14   2    0.000    0.000       0
+  1    178 00:02:14.63   0.7540A  -2.3863  -1.1008   1.3996  0.064  0.67  -0.998E+00     50.  0.131E+01  0.426E-01  0.156E-02  0.119E-03   60   3  19   22  10  16    9  14   2    0.000    0.000       0
+  1    179 00:02:15.38   0.7540A  -2.3864  -1.0795   1.3969  0.063  0.68  -0.106E+01     50.  0.131E+01  0.423E-01  0.155E-02  0.117E-03   60   3  18   22  10  16    9  14   2    0.000    0.000       0
+  1    180 00:02:16.14   0.7540A  -2.3865  -1.0806   1.3917  0.063  0.68  -0.996E+00     50.  0.131E+01  0.419E-01  0.155E-02  0.115E-03   60   3  18   21   9  16    9  14   2    0.000    0.000       0
+  1    181 00:02:16.89   0.7540A  -2.3865  -1.0686   1.3847  0.063  0.68  -0.997E+00     50.  0.131E+01  0.416E-01  0.154E-02  0.113E-03   60   2  16   21   9  16    9  14   2    0.000    0.000       0
+  1    182 00:02:17.65   0.7540A  -2.3866  -1.0480   1.3764  0.063  0.68  -0.982E+00     50.  0.131E+01  0.414E-01  0.153E-02  0.112E-03   60   2  15   22   9  15    9  14   2    0.000    0.000       0
+  1    183 00:02:18.40   0.7540A  -2.3867  -1.0460   1.3672  0.063  0.68  -0.965E+00     50.  0.131E+01  0.411E-01  0.153E-02  0.111E-03   60   2  15   22   9  15    9  14   2    0.000    0.000       0
+  1    184 00:02:19.15   0.7540A  -2.3867  -1.0472   1.3574  0.063  0.69  -0.956E+00     52.  0.130E+01  0.408E-01  0.153E-02  0.110E-03   60   1  13   18   1  12    9  14   2    0.000    0.000       0
+  1    185 00:02:19.91   0.7540A  -2.3868  -1.0554   1.3470  0.063  0.69  -0.951E+00     52.  0.130E+01  0.405E-01  0.152E-02  0.109E-03   60   1  12   18   1  12    9  14   2    0.000    0.000       0
+  1    186 00:02:20.66   0.7540A  -2.3869  -1.0564   1.3360  0.063  0.69  -0.946E+00     52.  0.130E+01  0.403E-01  0.152E-02  0.108E-03   60   1  12   18   1  12    9  14   2    0.000    0.000       0
