Index: /palm/trunk/SOURCE/module_interface.f90
===================================================================
--- /palm/trunk/SOURCE/module_interface.f90	(revision 4707)
+++ /palm/trunk/SOURCE/module_interface.f90	(revision 4708)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! pass fillvalue attribute to radiation output
+!
+! 4590 2020-07-06 14:34:59Z suehring
 ! Enable mpi-io for biomet
 ! 
@@ -1578,6 +1581,6 @@
 
     IF ( .NOT. found  .AND.  radiation )  THEN
-       CALL radiation_data_output_2d(                                          &
-               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
+       CALL radiation_data_output_2d(                                                       &
+               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
             )
     ENDIF
@@ -1673,5 +1676,5 @@
 
     IF ( .NOT. found  .AND.  radiation )  THEN
-       CALL radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
+       CALL radiation_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
        resorted = .TRUE.
     ENDIF
Index: /palm/trunk/SOURCE/radiation_model_mod.f90
===================================================================
--- /palm/trunk/SOURCE/radiation_model_mod.f90	(revision 4707)
+++ /palm/trunk/SOURCE/radiation_model_mod.f90	(revision 4708)
@@ -28,4 +28,9 @@
 ! -----------------
 ! $Id$
+! - Bugfix, correct mapping of RRTMG heating rates, as well as incoming/outgoing
+!   radiation onto the topography-following grid
+! - add fill values to the output
+!
+! 4694 2020-09-23 15:09:19Z pavelkrc
 ! Bugfix for tracing of maximum radiative fluxes
 !
@@ -401,5 +406,5 @@
     USE indices,                                                               &
         ONLY:  nnx, nny, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg,   &
-               nzb, nzt, topo_top_ind
+               nzb, nzt, topo_top_ind, wall_flags_total_0
 
     USE, INTRINSIC :: iso_c_binding
@@ -4390,9 +4395,14 @@
 !
 !--          Save fluxes
-             DO k = nzb, nzt+1
-                rad_lw_in(k,:,:)  = rrtm_lwdflx(0,k)
-                rad_lw_out(k,:,:) = rrtm_lwuflx(0,k)
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   k_topo_l = topo_top_ind(j,i,0)
+                   DO k = k_topo_l, nzt+1
+                      rad_lw_in(k,j,i)  = rrtm_lwdflx(0,k-k_topo_l+k_topo)
+                      rad_lw_out(k,j,i) = rrtm_lwuflx(0,k-k_topo_l+k_topo)
+                   ENDDO
+                ENDDO
              ENDDO
-             rad_lw_in_diff(:,:) = rad_lw_in(k_topo,:,:)
+             rad_lw_in_diff(:,:) = rrtm_lwdflx(0,k_topo)
 !
 !--          Save heating rates (convert from K/d to K/h).
@@ -4404,6 +4414,6 @@
                    k_topo_l = topo_top_ind(j,i,0)
                    DO k = k_topo_l+1, nzt+1
-                      rad_lw_hr(k,j,i)     = rrtm_lwhr(0,k-k_topo_l)  * d_hours_day
-                      rad_lw_cs_hr(k,j,i)  = rrtm_lwhrc(0,k-k_topo_l) * d_hours_day
+                      rad_lw_hr(k,j,i)     = rrtm_lwhr(0,k-k_topo_l+k_topo)  * d_hours_day
+                      rad_lw_cs_hr(k,j,i)  = rrtm_lwhrc(0,k-k_topo_l+k_topo) * d_hours_day
                    ENDDO
                 ENDDO
@@ -4499,7 +4509,12 @@
 !
 !--          Save radiation fluxes for the entire depth of the model domain
-             DO k = nzb, nzt+1
-                rad_sw_in(k,:,:)  = rrtm_swdflx(0,k)
-                rad_sw_out(k,:,:) = rrtm_swuflx(0,k)
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   k_topo_l = topo_top_ind(j,i,0)
+                   DO k = k_topo_l, nzt+1
+                      rad_sw_in(k,j,i)  = rrtm_swdflx(0,k-k_topo_l+k_topo)
+                      rad_sw_out(k,j,i) = rrtm_swuflx(0,k-k_topo_l+k_topo)
+                   ENDDO
+                ENDDO
              ENDDO
 !--          Save direct and diffuse SW radiation at the surface (required by RTM)
@@ -4509,7 +4524,12 @@
 !
 !--          Save heating rates (convert from K/d to K/s)
-             DO k = nzb+1, nzt+1
-                rad_sw_hr(k,:,:)     = rrtm_swhr(0,k)  * d_hours_day
-                rad_sw_cs_hr(k,:,:)  = rrtm_swhrc(0,k) * d_hours_day
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   k_topo_l = topo_top_ind(j,i,0)
+                   DO k = k_topo_l+1, nzt+1
+                      rad_sw_hr(k,j,i)     = rrtm_swhr(0,k-k_topo_l+k_topo)  * d_hours_day
+                      rad_sw_cs_hr(k,j,i)  = rrtm_swhrc(0,k-k_topo_l+k_topo) * d_hours_day
+                   ENDDO
+                ENDDO
              ENDDO
 !
@@ -4523,5 +4543,7 @@
 !
 !--    RRTMG is called for each (j,i) grid point separately, starting at the
-!--    highest topography level. Here no RTM is used since average_radiation is false
+!--    highest topography level. Here no RTM is used since average_radiation is false.
+!--    In fact, this branch is only called for homogeneous flat terrain so that the topography-top
+!--    index is actually alwasy zero.
        ELSE
 !
@@ -4694,5 +4716,5 @@
                 ENDDO
 !
-!--             Obtain topography top index (lower bound of RRTMG)
+!--             Obtain topography top index (lower bound of RRTMG). Is actually always zero.
                 k_topo = topo_top_ind(j,i,0)
 
@@ -5829,6 +5851,8 @@
 !--    Calculate tendency based on heating rate
        DO k = nzb+1, nzt+1
-          tend(k,j,i) = tend(k,j,i) + (rad_lw_hr(k,j,i) + rad_sw_hr(k,j,i))    &
-                                         * d_exner(k) * d_seconds_hour
+          tend(k,j,i) = tend(k,j,i) + (rad_lw_hr(k,j,i) + rad_sw_hr(k,j,i))                        &
+                                         * d_exner(k) * d_seconds_hour                             &
+                                         * MERGE( 1.0_wp, 0.0_wp,                                  &
+                                                  BTEST( wall_flags_total_0(k,j,i), 0 ) )
        ENDDO
 
@@ -5848,7 +5872,4 @@
  SUBROUTINE radiation_tendency ( tend )
 
-    USE indices,                                                               &
-        ONLY:  nxl, nxr, nyn, nys
-
     IMPLICIT NONE
 
@@ -5863,7 +5884,9 @@
           DO  j = nys, nyn
              DO k = nzb+1, nzt+1
-                tend(k,j,i) = tend(k,j,i) + ( rad_lw_hr(k,j,i)                 &
-                                          +  rad_sw_hr(k,j,i) ) * d_exner(k)   &
-                                          * d_seconds_hour
+                tend(k,j,i) = tend(k,j,i) + ( rad_lw_hr(k,j,i)                                      &
+                                          + rad_sw_hr(k,j,i) ) * d_exner(k)                         &
+                                          * d_seconds_hour                                          &
+                                          * MERGE( 1.0_wp, 0.0_wp,                                  &
+                                                   BTEST( wall_flags_total_0(k,j,i), 0 ) )
              ENDDO
           ENDDO
@@ -11137,5 +11160,5 @@
 !------------------------------------------------------------------------------!
  SUBROUTINE radiation_data_output_2d( av, variable, found, grid, mode,         &
-                                      local_pf, two_d, nzb_do, nzt_do )
+                                      local_pf, two_d, nzb_do, nzt_do, fill_value )
 
     USE indices
@@ -11161,5 +11184,5 @@
     LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
 
-    REAL(wp) ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
+    REAL(wp) ::  fill_value !< value for the _FillValue attribute
 
     REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
@@ -11340,5 +11363,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_lw_in(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_lw_in(k,j,i), fill_value,                       &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11364,5 +11388,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_lw_out(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_lw_out(k,j,i), fill_value,                      &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11388,5 +11413,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_lw_cs_hr(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_lw_cs_hr(k,j,i), fill_value,                    &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11412,5 +11438,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_lw_hr(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_lw_hr(k,j,i), fill_value,                       &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11436,5 +11463,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_sw_in(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_sw_in(k,j,i), fill_value,                       &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11460,5 +11488,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_sw_out(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_sw_out(k,j,i), fill_value,                      &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11484,5 +11513,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_sw_cs_hr(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_sw_cs_hr(k,j,i), fill_value,                    &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11508,5 +11538,6 @@
                 DO  j = nys, nyn
                    DO  k = nzb_do, nzt_do
-                      local_pf(i,j,k) = rad_sw_hr(k,j,i)
+                      local_pf(i,j,k) = MERGE( rad_sw_hr(k,j,i), fill_value,                       &
+                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
                    ENDDO
                 ENDDO
@@ -11542,5 +11573,5 @@
 !> Subroutine defining 3D output variables
 !------------------------------------------------------------------------------!
- SUBROUTINE radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
+ SUBROUTINE radiation_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
 
 
@@ -11561,5 +11592,5 @@
     LOGICAL      ::  found       !<
 
-    REAL(wp)     ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
+    REAL(wp)     ::  fill_value  !< value for the _FillValue attribute
 
     REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
@@ -11626,5 +11657,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_sw_in(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_sw_in(k,j,i), fill_value,                        &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11649,5 +11681,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_sw_out(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_sw_out(k,j,i), fill_value,                       &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11672,5 +11705,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_sw_cs_hr(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_sw_cs_hr(k,j,i), fill_value,                     &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11695,5 +11729,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_sw_hr(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_sw_hr(k,j,i), fill_value,                        &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11718,5 +11753,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_lw_in(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_lw_in(k,j,i), fill_value,                        &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11741,5 +11777,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_lw_out(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_lw_out(k,j,i), fill_value,                       &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11764,5 +11801,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_lw_cs_hr(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_lw_cs_hr(k,j,i), fill_value,                     &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -11787,5 +11825,6 @@
                DO  j = nys, nyn
                   DO  k = nzb_do, nzt_do
-                     local_pf(i,j,k) = rad_lw_hr(k,j,i)
+                     local_pf(i,j,k) = MERGE( rad_lw_hr(k,j,i), fill_value,                        &
+                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
                   ENDDO
                ENDDO
@@ -12227,5 +12266,6 @@
                 DO  k = 1, mask_size_l(mid,3)
                    local_pf(i,j,k) =  to_be_resorted(mask_k(mid,k), &
-                                      mask_j(mid,j),mask_i(mid,i))
+                                                     mask_j(mid,j), &
+                                                     mask_i(mid,i))
                 ENDDO
              ENDDO
Index: /palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc.000
===================================================================
--- /palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc.000	(revision 4707)
+++ /palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc.000	(revision 4708)
@@ -1,10 +1,10 @@
 
  ******************************    --------------------------------------------
- * PALM 6.0  Rev: 4700:47     *    atmosphere - run without 1D - prerun
+ * PALM 6.0  Rev: 4707M       *    atmosphere - run without 1D - prerun
  ******************************    --------------------------------------------
 
- Date:               2020-09-27    Run:       urban_environment__gfortran_defaul
- Time:                 19:53:12    Run-No.:   00
- Run on host:        gfortran_d
+ Date:               2020-09-28    Run:       urban_environment                 
+ Time:                 19:26:44    Run-No.:   00
+ Run on host:           default
  Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
  ------------------------------------------------------------------------------
@@ -404,8 +404,8 @@
   0     47 00:00:35.62   0.7570A  -2.3767  -1.0807   1.5657  0.070  0.54  -0.971E+00     46.  0.136E+01  0.341E-01  0.990E-03  0.182E-03   60   3  17   14  14   0   17  14  17    0.000    0.000       0
   0     48 00:00:36.37   0.7570A  -2.3768  -1.0892   1.5924  0.070  0.54  -0.972E+00     46.  0.136E+01  0.347E-01  0.101E-02  0.183E-03   60   3  16   14  14   0   17  14  17    0.000    0.000       0
-  0     49 00:00:37.13   0.7570A  -2.3768  -1.0970   1.6231  0.070  0.55  -0.972E+00     48.  0.136E+01  0.354E-01  0.102E-02  0.184E-03   60   2  14   14  14   0   18  14  17    0.000    0.000       0
+  0     49 00:00:37.13   0.7570A  -2.3768  -1.0970   1.6231  0.070  0.55  -0.972E+00     48.  0.136E+01  0.354E-01  0.102E-02  0.184E-03   60   2  15   14  14   0   18  14  17    0.000    0.000       0
   0     50 00:00:37.89   0.7570A  -2.3769  -1.1036   1.6518  0.070  0.55  -0.971E+00     48.  0.136E+01  0.361E-01  0.104E-02  0.185E-03   60   2  14   14  14   0   18  14  17    0.000    0.000       0
   0     51 00:00:38.64   0.7570A  -2.3770  -1.1086   1.6603  0.070  0.55  -0.969E+00     48.  0.136E+01  0.369E-01  0.105E-02  0.185E-03   60   2  13   14  14   0   18  14  17    0.000    0.000       0
-  0     52 00:00:39.40   0.7570A  -2.3771  -1.1118   1.6674  0.070  0.56  -0.968E+00     48.  0.136E+01  0.376E-01  0.107E-02  0.186E-03   60   1  12   14  14   0   19  14  17    0.000    0.000       0
+  0     52 00:00:39.40   0.7570A  -2.3771  -1.1118   1.6673  0.070  0.56  -0.968E+00     48.  0.136E+01  0.376E-01  0.107E-02  0.186E-03   60   1  12   14  14   0   19  14  17    0.000    0.000       0
   0     53 00:00:40.16   0.7570A  -2.3771  -1.1132   1.6749  0.069  0.56  -0.966E+00     48.  0.136E+01  0.383E-01  0.108E-02  0.186E-03   60   1  11   14  14   0   19  14  17    0.000    0.000       0
   0     54 00:00:40.92   0.7570A  -2.3772  -1.1128   1.6618  0.069  0.56  -0.965E+00     48.  0.136E+01  0.390E-01  0.110E-02  0.186E-03   60   1  10   14  14   0   19  14  17    0.000    0.000       0
@@ -427,5 +427,5 @@
   0     70 00:00:53.03   0.7570A  -2.3784  -1.2147   1.4712  0.069  0.61  -0.973E+00     52.  0.136E+01  0.497E-01  0.131E-02  0.192E-03   60  16  16    8  14   5    6  14   3    0.000    0.000       0
   0     71 00:00:53.78   0.7570A  -2.3785  -1.2310   1.4862  0.069  0.62  -0.970E+00     54.  0.136E+01  0.503E-01  0.132E-02  0.192E-03   60  16  16   22  12  10    7  14   3    0.000    0.000       0
-  0     72 00:00:54.54   0.7570A  -2.3785  -1.2425   1.5064  0.069  0.62  -0.967E+00     54.  0.136E+01  0.508E-01  0.133E-02  0.193E-03   60  15  14   22  12  10    7  14   3    0.000    0.000       0
+  0     72 00:00:54.54   0.7570A  -2.3786  -1.2425   1.5064  0.069  0.62  -0.967E+00     54.  0.136E+01  0.508E-01  0.133E-02  0.193E-03   60  15  14   22  12  10    7  14   3    0.000    0.000       0
   0     73 00:00:55.30   0.7570A  -2.3786  -1.2454   1.5225  0.069  0.63  -0.963E+00     54.  0.136E+01  0.514E-01  0.135E-02  0.193E-03   60  15  14   22  12  10    7  14   3    0.000    0.000       0
   0     74 00:00:56.06   0.7570A  -2.3787  -1.2396   1.5349  0.069  0.63  -0.959E+00     54.  0.136E+01  0.520E-01  0.136E-02  0.194E-03   60  15  13   22  12  10    7  14   3    0.000    0.000       0
Index: /palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc.000
===================================================================
--- /palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc.000	(revision 4707)
+++ /palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc.000	(revision 4708)
@@ -36,10 +36,10 @@
   1    120 00:01:30.84   0.7560A  -2.3821  -1.3511   1.4801  0.069  0.66  -0.926E+00     44.  0.134E+01  0.646E-01  0.183E-02  0.181E-03   60   1  12    8  13   5    7  14   3    0.000    0.000       0
   1    121 00:01:31.60   0.7560A  -2.3822  -1.3448   1.4772  0.069  0.66  -0.906E+00     44.  0.134E+01  0.646E-01  0.184E-02  0.179E-03   60   0  10    8  13   5    8  14   3    0.000    0.000       0
-  1    122 00:01:32.35   0.7560A  -2.3822  -1.3391   1.4747  0.069  0.66  -0.902E+00     44.  0.134E+01  0.646E-01  0.185E-02  0.176E-03   60   0   9    8  13   5    8  14   3    0.000    0.000       0
+  1    122 00:01:32.35   0.7560A  -2.3823  -1.3391   1.4747  0.069  0.66  -0.902E+00     44.  0.134E+01  0.646E-01  0.185E-02  0.176E-03   60   0   9    8  13   5    8  14   3    0.000    0.000       0
   1    123 00:01:33.11   0.7560A  -2.3823  -1.3340   1.4713  0.069  0.65  -0.901E+00     44.  0.134E+01  0.646E-01  0.186E-02  0.174E-03   60   0   9    8  13   5    8  14   3    0.000    0.000       0
   1    124 00:01:33.87   0.7560A  -2.3824  -1.3292   1.4673  0.069  0.66  -0.902E+00     46.  0.134E+01  0.645E-01  0.187E-02  0.172E-03   60   0   8    8  13   5    8  14   3    0.000    0.000       0
   1    125 00:01:34.62   0.7560A  -2.3825  -1.3244   1.4625  0.069  0.65  -0.905E+00     46.  0.134E+01  0.644E-01  0.188E-02  0.170E-03   60  19   7    8  13   5    8  14   3    0.000    0.000       0
   1    126 00:01:35.38   0.7550A  -2.3825  -1.3188   1.4570  0.069  0.65  -0.912E+00     46.  0.134E+01  0.644E-01  0.188E-02  0.168E-03   60  19   6    8  13   5    8  14   3    0.000    0.000       0
-  1    127 00:01:36.13   0.7550A  -2.3826  -1.3120   1.4509  0.069  0.65  -0.929E+00     46.  0.134E+01  0.643E-01  0.189E-02  0.166E-03   60  19   6    8  13   5    8  14   3    0.000    0.000       0
+  1    127 00:01:36.13   0.7550A  -2.3826  -1.3119   1.4509  0.069  0.65  -0.929E+00     46.  0.134E+01  0.643E-01  0.189E-02  0.166E-03   60  19   6    8  13   5    8  14   3    0.000    0.000       0
   1    128 00:01:36.89   0.7550A  -2.3827  -1.3034   1.4448  0.069  0.64  -0.932E+00     46.  0.134E+01  0.642E-01  0.190E-02  0.164E-03   60  19   5    8  13   5    7  14   3    0.000    0.000       0
   1    129 00:01:37.64   0.7550A  -2.3828  -1.2933   1.4391  0.069  0.64  -0.960E+00     46.  0.134E+01  0.640E-01  0.190E-02  0.163E-03   60  18   3    8  13   5    7  14   3    0.000    0.000       0
@@ -57,5 +57,5 @@
   1    141 00:01:46.70   0.7550A  -2.3836  -1.2522   1.4337  0.068  0.62  -0.990E+00     48.  0.133E+01  0.614E-01  0.195E-02  0.165E-03   60  14  11   22   0  10    8  13   3    0.000    0.000       0
   1    142 00:01:47.46   0.7550A  -2.3837  -1.2456   1.4341  0.068  0.62  -0.937E+00     48.  0.133E+01  0.611E-01  0.195E-02  0.165E-03   60  14  11   22   0   9    8  13   3    0.000    0.000       0
-  1    143 00:01:48.21   0.7550A  -2.3838  -1.2475   1.4321  0.068  0.62  -0.927E+00     48.  0.133E+01  0.608E-01  0.194E-02  0.166E-03   60  13   9   22   0   9    8  13   3    0.000    0.000       0
+  1    143 00:01:48.21   0.7550A  -2.3838  -1.2475   1.4321  0.068  0.62  -0.927E+00     48.  0.133E+01  0.608E-01  0.194E-02  0.166E-03   60  14  10   22   0   9    8  13   3    0.000    0.000       0
   1    144 00:01:48.97   0.7550A  -2.3839  -1.2648   1.4275  0.067  0.62  -0.916E+00     48.  0.133E+01  0.604E-01  0.194E-02  0.166E-03   60  13   9   22  19   9    8  13   3    0.000    0.000       0
   1    145 00:01:49.72   0.7550A  -2.3839  -1.2658   1.4337  0.067  0.63  -0.911E+00     48.  0.133E+01  0.601E-01  0.194E-02  0.166E-03   60  13   8   22  19   9    9  13   3    0.000    0.000       0
@@ -73,8 +73,8 @@
   1    157 00:01:58.78   0.7550A  -2.3848  -1.2302   1.4154  0.066  0.66  -0.878E+00     50.  0.132E+01  0.557E-01  0.189E-02  0.160E-03   60  10  19   23  16   5   20  14  16    0.000    0.000       0
   1    158 00:01:59.54   0.7550A  -2.3849  -1.2246   1.4282  0.066  0.66  -0.876E+00     50.  0.132E+01  0.553E-01  0.189E-02  0.159E-03   60   9  17   22  17   6   20  14  16    0.000    0.000       0
-  1    159 00:02:00.29   0.7550A  -2.3849  -1.2143   1.4402  0.066  0.67  -0.879E+00     50.  0.132E+01  0.549E-01  0.189E-02  0.158E-03   60   9  16   21  17   7   20  14  16    0.000    0.000       0
+  1    159 00:02:00.29   0.7550A  -2.3850  -1.2143   1.4402  0.066  0.67  -0.879E+00     50.  0.132E+01  0.549E-01  0.189E-02  0.158E-03   60   9  16   21  17   7   20  14  16    0.000    0.000       0
   1    160 00:02:01.05   0.7550A  -2.3850  -1.2267   1.4485  0.066  0.67  -0.876E+00     50.  0.132E+01  0.545E-01  0.188E-02  0.156E-03   60   9  16   21  17   7   20  14  16    0.000    0.000       0
-  1    161 00:02:01.80   0.7550A  -2.3851  -1.2264   1.4512  0.066  0.68  -0.871E+00     52.  0.131E+01  0.541E-01  0.188E-02  0.155E-03   60   9  16   21  17   7   20  14  16    0.000    0.000       0
-  1    162 00:02:02.56   0.7550A  -2.3852  -1.2129   1.4474  0.066  0.68  -0.869E+00     52.  0.131E+01  0.537E-01  0.187E-02  0.153E-03   60   8  14   21  17   7   20  14  16    0.000    0.000       0
+  1    161 00:02:01.80   0.7550A  -2.3851  -1.2264   1.4512  0.066  0.68  -0.871E+00     52.  0.131E+01  0.541E-01  0.188E-02  0.155E-03   60   9  15   21  17   7   20  14  16    0.000    0.000       0
+  1    162 00:02:02.56   0.7550A  -2.3852  -1.2130   1.4474  0.066  0.68  -0.869E+00     52.  0.131E+01  0.537E-01  0.187E-02  0.153E-03   60   8  14   21  17   7   20  14  16    0.000    0.000       0
   1    163 00:02:03.31   0.7550A  -2.3852  -1.1858   1.4383  0.066  0.69  -0.867E+00     52.  0.131E+01  0.533E-01  0.187E-02  0.152E-03   60   8  13   21  17   7   20  14  17    0.000    0.000       0
   1    164 00:02:04.07   0.7550A  -2.3853  -1.1919   1.4420  0.066  0.69  -0.865E+00     52.  0.131E+01  0.530E-01  0.186E-02  0.150E-03   60   8  13   21  16   6   20  14  17    0.000    0.000       0
@@ -82,5 +82,5 @@
   1    166 00:02:05.58   0.7550A  -2.3854  -1.1813   1.4446  0.066  0.70  -0.862E+00     52.  0.131E+01  0.522E-01  0.184E-02  0.146E-03   60   7  10   21  16   6   22  14  16    0.000    0.000       0
   1    167 00:02:06.33   0.7550A  -2.3855  -1.1855   1.4438  0.066  0.70  -0.861E+00     52.  0.131E+01  0.518E-01  0.184E-02  0.145E-03   60   7  10   20  16   7   20  14  17    0.000    0.000       0
-  1    168 00:02:07.09   0.7550A  -2.3856  -1.1984   1.4424  0.065  0.70  -0.860E+00     52.  0.131E+01  0.515E-01  0.183E-02  0.144E-03   60   6   8   20  16   7   20  14  17    0.000    0.000       0
+  1    168 00:02:07.09   0.7550A  -2.3856  -1.1985   1.4424  0.065  0.70  -0.860E+00     52.  0.131E+01  0.515E-01  0.183E-02  0.144E-03   60   6   8   20  16   7   20  14  17    0.000    0.000       0
   1    169 00:02:07.84   0.7550A  -2.3857  -1.1996   1.4399  0.065  0.71  -0.859E+00     52.  0.131E+01  0.511E-01  0.183E-02  0.143E-03   60   6   7   20  16   7   20  14  17    0.000    0.000       0
   1    170 00:02:08.60   0.7540A  -2.3857  -1.1894   1.4363  0.065  0.71  -0.859E+00     52.  0.131E+01  0.508E-01  0.183E-02  0.143E-03   60   6   7   20  16   7   23  14  16    0.000    0.000       0
@@ -90,5 +90,5 @@
   1    174 00:02:11.61   0.7540A  -2.3860  -1.2886   1.3953  0.065  0.76  -0.858E+00     60.  0.131E+01  0.495E-01  0.181E-02  0.142E-03   60  18   2   26  12  16   20  14  17    0.000    0.000       0
   1    175 00:02:12.37   0.7540A  -2.3861  -1.3028   1.3791  0.064  0.76  -0.864E+00     60.  0.130E+01  0.492E-01  0.181E-02  0.143E-03   60   4   1   26  12  16   20  14  17    0.000    0.000       0
-  1    176 00:02:13.12   0.7540A  -2.3862  -1.2817   1.3611  0.064  0.77  -0.877E+00     60.  0.130E+01  0.489E-01  0.181E-02  0.143E-03   60  17  19   26  12  16   20  14  17    0.000    0.000       0
+  1    176 00:02:13.12   0.7540A  -2.3862  -1.2818   1.3610  0.064  0.77  -0.877E+00     60.  0.130E+01  0.489E-01  0.181E-02  0.143E-03   60  17  19   26  12  16   20  14  17    0.000    0.000       0
   1    177 00:02:13.88   0.7540A  -2.3862  -1.3029   1.3395  0.064  0.77  -0.856E+00     60.  0.130E+01  0.487E-01  0.180E-02  0.143E-03   60   3  18   25  12  16   20  14  17    0.000    0.000       0
   1    178 00:02:14.63   0.7540A  -2.3863  -1.3562   1.3123  0.064  0.77  -0.851E+00     60.  0.130E+01  0.484E-01  0.180E-02  0.143E-03   60   3  18   26  11  15   20  14  17    0.000    0.000       0
@@ -96,7 +96,7 @@
   1    180 00:02:16.14   0.7540A  -2.3865  -1.3720   1.3021  0.063  0.77  -0.848E+00     60.  0.130E+01  0.479E-01  0.178E-02  0.142E-03   60   2  15   26  11  15    8  14   3    0.000    0.000       0
   1    181 00:02:16.89   0.7540A  -2.3865  -1.3415   1.3031  0.063  0.77  -0.847E+00     60.  0.130E+01  0.477E-01  0.178E-02  0.141E-03   60   2  15   26  11  15    8  14   3    0.000    0.000       0
-  1    182 00:02:17.65   0.7540A  -2.3866  -1.3268   1.3036  0.063  0.78  -0.847E+00     60.  0.130E+01  0.475E-01  0.177E-02  0.140E-03   60  15  13   26  10  14    8  14   3    0.000    0.000       0
+  1    182 00:02:17.65   0.7540A  -2.3866  -1.3268   1.3035  0.063  0.78  -0.847E+00     60.  0.130E+01  0.475E-01  0.177E-02  0.140E-03   60  15  13   26  10  14    8  14   3    0.000    0.000       0
   1    183 00:02:18.40   0.7540A  -2.3867  -1.3536   1.3035  0.063  0.78  -0.847E+00     60.  0.130E+01  0.473E-01  0.176E-02  0.138E-03   60   1  12   26  10  14    8  14   3    0.000    0.000       0
   1    184 00:02:19.15   0.7540A  -2.3868  -1.3519   1.3029  0.063  0.78  -0.848E+00     60.  0.130E+01  0.471E-01  0.175E-02  0.136E-03   60   1  12   26  10  14    8  14   3    0.000    0.000       0
   1    185 00:02:19.91   0.7540A  -2.3868  -1.3299   1.3016  0.062  0.78  -0.849E+00     60.  0.130E+01  0.470E-01  0.174E-02  0.134E-03   60   1  11   26  10  14    8  14   3    0.000    0.000       0
-  1    186 00:02:20.66   0.7540A  -2.3869  -1.2961   1.2992  0.062  0.78  -0.849E+00     60.  0.130E+01  0.468E-01  0.174E-02  0.132E-03   60   0  10   26  10  14    8  14   3    0.000    0.000       0
+  1    186 00:02:20.66   0.7540A  -2.3869  -1.2962   1.2992  0.062  0.78  -0.849E+00     60.  0.130E+01  0.468E-01  0.174E-02  0.132E-03   60   0  10   26  10  14    8  14   3    0.000    0.000       0
