Changeset 4708


Ignore:
Timestamp:
Sep 28, 2020 5:42:58 PM (4 years ago)
Author:
suehring
Message:

Bugfix, correct mapping of RRTMG heating rates, as well as incoming/outgoing radiation onto the topography-following grid; add fill values to the radiation output

Location:
palm/trunk
Files:
10 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/module_interface.f90

    r4590 r4708  
    2525! -----------------
    2626! $Id$
     27! pass fillvalue attribute to radiation output
     28!
     29! 4590 2020-07-06 14:34:59Z suehring
    2730! Enable mpi-io for biomet
    2831!
     
    15781581
    15791582    IF ( .NOT. found  .AND.  radiation )  THEN
    1580        CALL radiation_data_output_2d(                                          &
    1581                av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do&
     1583       CALL radiation_data_output_2d(                                                       &
     1584               av, variable, found, grid, mode, local_pf, two_d, nzb_do, nzt_do, fill_value &
    15821585            )
    15831586    ENDIF
     
    16731676
    16741677    IF ( .NOT. found  .AND.  radiation )  THEN
    1675        CALL radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
     1678       CALL radiation_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
    16761679       resorted = .TRUE.
    16771680    ENDIF
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r4694 r4708  
    2828! -----------------
    2929! $Id$
     30! - Bugfix, correct mapping of RRTMG heating rates, as well as incoming/outgoing
     31!   radiation onto the topography-following grid
     32! - add fill values to the output
     33!
     34! 4694 2020-09-23 15:09:19Z pavelkrc
    3035! Bugfix for tracing of maximum radiative fluxes
    3136!
     
    401406    USE indices,                                                               &
    402407        ONLY:  nnx, nny, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg,   &
    403                nzb, nzt, topo_top_ind
     408               nzb, nzt, topo_top_ind, wall_flags_total_0
    404409
    405410    USE, INTRINSIC :: iso_c_binding
     
    43904395!
    43914396!--          Save fluxes
    4392              DO k = nzb, nzt+1
    4393                 rad_lw_in(k,:,:)  = rrtm_lwdflx(0,k)
    4394                 rad_lw_out(k,:,:) = rrtm_lwuflx(0,k)
     4397             DO  i = nxl, nxr
     4398                DO  j = nys, nyn
     4399                   k_topo_l = topo_top_ind(j,i,0)
     4400                   DO k = k_topo_l, nzt+1
     4401                      rad_lw_in(k,j,i)  = rrtm_lwdflx(0,k-k_topo_l+k_topo)
     4402                      rad_lw_out(k,j,i) = rrtm_lwuflx(0,k-k_topo_l+k_topo)
     4403                   ENDDO
     4404                ENDDO
    43954405             ENDDO
    4396              rad_lw_in_diff(:,:) = rad_lw_in(k_topo,:,:)
     4406             rad_lw_in_diff(:,:) = rrtm_lwdflx(0,k_topo)
    43974407!
    43984408!--          Save heating rates (convert from K/d to K/h).
     
    44044414                   k_topo_l = topo_top_ind(j,i,0)
    44054415                   DO k = k_topo_l+1, nzt+1
    4406                       rad_lw_hr(k,j,i)     = rrtm_lwhr(0,k-k_topo_l)  * d_hours_day
    4407                       rad_lw_cs_hr(k,j,i)  = rrtm_lwhrc(0,k-k_topo_l) * d_hours_day
     4416                      rad_lw_hr(k,j,i)     = rrtm_lwhr(0,k-k_topo_l+k_topo)  * d_hours_day
     4417                      rad_lw_cs_hr(k,j,i)  = rrtm_lwhrc(0,k-k_topo_l+k_topo) * d_hours_day
    44084418                   ENDDO
    44094419                ENDDO
     
    44994509!
    45004510!--          Save radiation fluxes for the entire depth of the model domain
    4501              DO k = nzb, nzt+1
    4502                 rad_sw_in(k,:,:)  = rrtm_swdflx(0,k)
    4503                 rad_sw_out(k,:,:) = rrtm_swuflx(0,k)
     4511             DO  i = nxl, nxr
     4512                DO  j = nys, nyn
     4513                   k_topo_l = topo_top_ind(j,i,0)
     4514                   DO k = k_topo_l, nzt+1
     4515                      rad_sw_in(k,j,i)  = rrtm_swdflx(0,k-k_topo_l+k_topo)
     4516                      rad_sw_out(k,j,i) = rrtm_swuflx(0,k-k_topo_l+k_topo)
     4517                   ENDDO
     4518                ENDDO
    45044519             ENDDO
    45054520!--          Save direct and diffuse SW radiation at the surface (required by RTM)
     
    45094524!
    45104525!--          Save heating rates (convert from K/d to K/s)
    4511              DO k = nzb+1, nzt+1
    4512                 rad_sw_hr(k,:,:)     = rrtm_swhr(0,k)  * d_hours_day
    4513                 rad_sw_cs_hr(k,:,:)  = rrtm_swhrc(0,k) * d_hours_day
     4526             DO  i = nxl, nxr
     4527                DO  j = nys, nyn
     4528                   k_topo_l = topo_top_ind(j,i,0)
     4529                   DO k = k_topo_l+1, nzt+1
     4530                      rad_sw_hr(k,j,i)     = rrtm_swhr(0,k-k_topo_l+k_topo)  * d_hours_day
     4531                      rad_sw_cs_hr(k,j,i)  = rrtm_swhrc(0,k-k_topo_l+k_topo) * d_hours_day
     4532                   ENDDO
     4533                ENDDO
    45144534             ENDDO
    45154535!
     
    45234543!
    45244544!--    RRTMG is called for each (j,i) grid point separately, starting at the
    4525 !--    highest topography level. Here no RTM is used since average_radiation is false
     4545!--    highest topography level. Here no RTM is used since average_radiation is false.
     4546!--    In fact, this branch is only called for homogeneous flat terrain so that the topography-top
     4547!--    index is actually alwasy zero.
    45264548       ELSE
    45274549!
     
    46944716                ENDDO
    46954717!
    4696 !--             Obtain topography top index (lower bound of RRTMG)
     4718!--             Obtain topography top index (lower bound of RRTMG). Is actually always zero.
    46974719                k_topo = topo_top_ind(j,i,0)
    46984720
     
    58295851!--    Calculate tendency based on heating rate
    58305852       DO k = nzb+1, nzt+1
    5831           tend(k,j,i) = tend(k,j,i) + (rad_lw_hr(k,j,i) + rad_sw_hr(k,j,i))    &
    5832                                          * d_exner(k) * d_seconds_hour
     5853          tend(k,j,i) = tend(k,j,i) + (rad_lw_hr(k,j,i) + rad_sw_hr(k,j,i))                        &
     5854                                         * d_exner(k) * d_seconds_hour                             &
     5855                                         * MERGE( 1.0_wp, 0.0_wp,                                  &
     5856                                                  BTEST( wall_flags_total_0(k,j,i), 0 ) )
    58335857       ENDDO
    58345858
     
    58485872 SUBROUTINE radiation_tendency ( tend )
    58495873
    5850     USE indices,                                                               &
    5851         ONLY:  nxl, nxr, nyn, nys
    5852 
    58535874    IMPLICIT NONE
    58545875
     
    58635884          DO  j = nys, nyn
    58645885             DO k = nzb+1, nzt+1
    5865                 tend(k,j,i) = tend(k,j,i) + ( rad_lw_hr(k,j,i)                 &
    5866                                           +  rad_sw_hr(k,j,i) ) * d_exner(k)   &
    5867                                           * d_seconds_hour
     5886                tend(k,j,i) = tend(k,j,i) + ( rad_lw_hr(k,j,i)                                      &
     5887                                          + rad_sw_hr(k,j,i) ) * d_exner(k)                         &
     5888                                          * d_seconds_hour                                          &
     5889                                          * MERGE( 1.0_wp, 0.0_wp,                                  &
     5890                                                   BTEST( wall_flags_total_0(k,j,i), 0 ) )
    58685891             ENDDO
    58695892          ENDDO
     
    1113711160!------------------------------------------------------------------------------!
    1113811161 SUBROUTINE radiation_data_output_2d( av, variable, found, grid, mode,         &
    11139                                       local_pf, two_d, nzb_do, nzt_do )
     11162                                      local_pf, two_d, nzb_do, nzt_do, fill_value )
    1114011163
    1114111164    USE indices
     
    1116111184    LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
    1116211185
    11163     REAL(wp) ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
     11186    REAL(wp) ::  fill_value !< value for the _FillValue attribute
    1116411187
    1116511188    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
     
    1134011363                DO  j = nys, nyn
    1134111364                   DO  k = nzb_do, nzt_do
    11342                       local_pf(i,j,k) = rad_lw_in(k,j,i)
     11365                      local_pf(i,j,k) = MERGE( rad_lw_in(k,j,i), fill_value,                       &
     11366                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1134311367                   ENDDO
    1134411368                ENDDO
     
    1136411388                DO  j = nys, nyn
    1136511389                   DO  k = nzb_do, nzt_do
    11366                       local_pf(i,j,k) = rad_lw_out(k,j,i)
     11390                      local_pf(i,j,k) = MERGE( rad_lw_out(k,j,i), fill_value,                      &
     11391                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1136711392                   ENDDO
    1136811393                ENDDO
     
    1138811413                DO  j = nys, nyn
    1138911414                   DO  k = nzb_do, nzt_do
    11390                       local_pf(i,j,k) = rad_lw_cs_hr(k,j,i)
     11415                      local_pf(i,j,k) = MERGE( rad_lw_cs_hr(k,j,i), fill_value,                    &
     11416                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1139111417                   ENDDO
    1139211418                ENDDO
     
    1141211438                DO  j = nys, nyn
    1141311439                   DO  k = nzb_do, nzt_do
    11414                       local_pf(i,j,k) = rad_lw_hr(k,j,i)
     11440                      local_pf(i,j,k) = MERGE( rad_lw_hr(k,j,i), fill_value,                       &
     11441                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1141511442                   ENDDO
    1141611443                ENDDO
     
    1143611463                DO  j = nys, nyn
    1143711464                   DO  k = nzb_do, nzt_do
    11438                       local_pf(i,j,k) = rad_sw_in(k,j,i)
     11465                      local_pf(i,j,k) = MERGE( rad_sw_in(k,j,i), fill_value,                       &
     11466                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1143911467                   ENDDO
    1144011468                ENDDO
     
    1146011488                DO  j = nys, nyn
    1146111489                   DO  k = nzb_do, nzt_do
    11462                       local_pf(i,j,k) = rad_sw_out(k,j,i)
     11490                      local_pf(i,j,k) = MERGE( rad_sw_out(k,j,i), fill_value,                      &
     11491                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1146311492                   ENDDO
    1146411493                ENDDO
     
    1148411513                DO  j = nys, nyn
    1148511514                   DO  k = nzb_do, nzt_do
    11486                       local_pf(i,j,k) = rad_sw_cs_hr(k,j,i)
     11515                      local_pf(i,j,k) = MERGE( rad_sw_cs_hr(k,j,i), fill_value,                    &
     11516                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1148711517                   ENDDO
    1148811518                ENDDO
     
    1150811538                DO  j = nys, nyn
    1150911539                   DO  k = nzb_do, nzt_do
    11510                       local_pf(i,j,k) = rad_sw_hr(k,j,i)
     11540                      local_pf(i,j,k) = MERGE( rad_sw_hr(k,j,i), fill_value,                       &
     11541                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1151111542                   ENDDO
    1151211543                ENDDO
     
    1154211573!> Subroutine defining 3D output variables
    1154311574!------------------------------------------------------------------------------!
    11544  SUBROUTINE radiation_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
     11575 SUBROUTINE radiation_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
    1154511576
    1154611577
     
    1156111592    LOGICAL      ::  found       !<
    1156211593
    11563     REAL(wp)     ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
     11594    REAL(wp)     ::  fill_value !< value for the _FillValue attribute
    1156411595
    1156511596    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<
     
    1162611657               DO  j = nys, nyn
    1162711658                  DO  k = nzb_do, nzt_do
    11628                      local_pf(i,j,k) = rad_sw_in(k,j,i)
     11659                     local_pf(i,j,k) = MERGE( rad_sw_in(k,j,i), fill_value,                        &
     11660                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1162911661                  ENDDO
    1163011662               ENDDO
     
    1164911681               DO  j = nys, nyn
    1165011682                  DO  k = nzb_do, nzt_do
    11651                      local_pf(i,j,k) = rad_sw_out(k,j,i)
     11683                     local_pf(i,j,k) = MERGE( rad_sw_out(k,j,i), fill_value,                       &
     11684                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1165211685                  ENDDO
    1165311686               ENDDO
     
    1167211705               DO  j = nys, nyn
    1167311706                  DO  k = nzb_do, nzt_do
    11674                      local_pf(i,j,k) = rad_sw_cs_hr(k,j,i)
     11707                     local_pf(i,j,k) = MERGE( rad_sw_cs_hr(k,j,i), fill_value,                     &
     11708                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1167511709                  ENDDO
    1167611710               ENDDO
     
    1169511729               DO  j = nys, nyn
    1169611730                  DO  k = nzb_do, nzt_do
    11697                      local_pf(i,j,k) = rad_sw_hr(k,j,i)
     11731                     local_pf(i,j,k) = MERGE( rad_sw_hr(k,j,i), fill_value,                        &
     11732                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1169811733                  ENDDO
    1169911734               ENDDO
     
    1171811753               DO  j = nys, nyn
    1171911754                  DO  k = nzb_do, nzt_do
    11720                      local_pf(i,j,k) = rad_lw_in(k,j,i)
     11755                     local_pf(i,j,k) = MERGE( rad_lw_in(k,j,i), fill_value,                        &
     11756                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1172111757                  ENDDO
    1172211758               ENDDO
     
    1174111777               DO  j = nys, nyn
    1174211778                  DO  k = nzb_do, nzt_do
    11743                      local_pf(i,j,k) = rad_lw_out(k,j,i)
     11779                     local_pf(i,j,k) = MERGE( rad_lw_out(k,j,i), fill_value,                       &
     11780                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1174411781                  ENDDO
    1174511782               ENDDO
     
    1176411801               DO  j = nys, nyn
    1176511802                  DO  k = nzb_do, nzt_do
    11766                      local_pf(i,j,k) = rad_lw_cs_hr(k,j,i)
     11803                     local_pf(i,j,k) = MERGE( rad_lw_cs_hr(k,j,i), fill_value,                     &
     11804                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1176711805                  ENDDO
    1176811806               ENDDO
     
    1178711825               DO  j = nys, nyn
    1178811826                  DO  k = nzb_do, nzt_do
    11789                      local_pf(i,j,k) = rad_lw_hr(k,j,i)
     11827                     local_pf(i,j,k) = MERGE( rad_lw_hr(k,j,i), fill_value,                        &
     11828                                              BTEST( wall_flags_total_0(k,j,i), 0 ) )
    1179011829                  ENDDO
    1179111830               ENDDO
     
    1222712266                DO  k = 1, mask_size_l(mid,3)
    1222812267                   local_pf(i,j,k) =  to_be_resorted(mask_k(mid,k), &
    12229                                       mask_j(mid,j),mask_i(mid,i))
     12268                                                     mask_j(mid,j), &
     12269                                                     mask_i(mid,i))
    1223012270                ENDDO
    1223112271             ENDDO
  • palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc.000

    r4702 r4708  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 6.0  Rev: 4700:47     *    atmosphere - run without 1D - prerun
     3 * PALM 6.0  Rev: 4707M       *    atmosphere - run without 1D - prerun
    44 ******************************    --------------------------------------------
    55
    6  Date:               2020-09-27    Run:       urban_environment__gfortran_defaul
    7  Time:                 19:53:12    Run-No.:   00
    8  Run on host:        gfortran_d
     6 Date:               2020-09-28    Run:       urban_environment                 
     7 Time:                 19:26:44    Run-No.:   00
     8 Run on host:           default
    99 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
    1010 ------------------------------------------------------------------------------
     
    404404  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
    405405  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
    406   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
     406  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
    407407  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
    408408  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
    409   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
     409  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
    410410  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
    411411  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
     
    427427  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
    428428  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
    429   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
     429  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
    430430  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
    431431  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
  • palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc.000

    r4702 r4708  
    3636  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
    3737  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
    38   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
     38  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
    3939  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
    4040  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
    4141  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
    4242  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
    43   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
     43  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
    4444  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
    4545  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
     
    5757  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
    5858  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
    59   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
     59  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
    6060  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
    6161  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
     
    7373  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
    7474  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
    75   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
     75  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
    7676  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
    77   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
    78   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
     77  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
     78  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
    7979  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
    8080  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
     
    8282  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
    8383  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
    84   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
     84  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
    8585  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
    8686  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
     
    9090  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
    9191  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
    92   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
     92  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
    9393  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
    9494  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
     
    9696  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
    9797  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
    98   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
     98  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
    9999  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
    100100  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
    101101  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
    102   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
     102  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
Note: See TracChangeset for help on using the changeset viewer.