Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Changeset 3449 for palm


Ignore:
Timestamp:
Oct 29, 2018 7:36:56 PM (6 years ago)
Author:
suehring
Message:

Branch resler -r 3439 re-integrated into current trunk: RTM 3.0, transpiration of plant canopy, output fixes in USM

Location:
palm/trunk
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

    r3361 r3449  
    2525! -----------------
    2626! $Id$
     27! +degc_to_k
     28!
     29! 3361 2018-10-16 20:39:37Z knoop
    2730! New module (introduced with modularization of bulk cloud physics model)
    2831!
    … …  
    4346    IMPLICIT NONE
    4447
    45     REAL(wp), PARAMETER ::  pi = 3.141592654_wp  !< PI
    46 
    4748    REAL(wp), PARAMETER ::  c_p = 1005.0_wp                           !< heat capacity of dry air (J kg-1 K-1)
     49    REAL(wp), PARAMETER ::  degc_to_k = 273.15_wp                     !< temperature (in K) of 0 deg C (K)
    4850    REAL(wp), PARAMETER ::  g = 9.81_wp                               !< gravitational acceleration (m s-2)
    4951    REAL(wp), PARAMETER ::  kappa = 0.4_wp                            !< von Karman constant
    … …  
    5557    REAL(wp), PARAMETER ::  molecular_weight_of_nh4no3 = 0.08004_wp   !< mol. m. ammonium sulfate (kg mol-1)
    5658    REAL(wp), PARAMETER ::  molecular_weight_of_water = 0.01801528_wp !< mol. m. H2O (kg mol-1)
     59    REAL(wp), PARAMETER ::  pi = 3.141592654_wp                       !< PI
    5760    REAL(wp), PARAMETER ::  rho_l = 1.0E3_wp                          !< density of water (kg m-3)
    5861    REAL(wp), PARAMETER ::  rho_nacl = 2165.0_wp                      !< density of NaCl (kg m-3)
    … …  
    6164    REAL(wp), PARAMETER ::  r_d = 287.0_wp                            !< sp. gas const. dry air (J kg-1 K-1)
    6265    REAL(wp), PARAMETER ::  r_v = 461.51_wp                           !< sp. gas const. water vapor (J kg-1 K-1)
    63     REAL(wp), PARAMETER ::  sigma_sb = 5.67E-08_wp                    !< Stefan-Boltzmann constant
     66    REAL(wp), PARAMETER ::  sigma_sb = 5.67037E-08_wp                 !< Stefan-Boltzmann constant
    6467    REAL(wp), PARAMETER ::  solar_constant = 1368.0_wp                !< solar constant at top of atmosphere
    6568    REAL(wp), PARAMETER ::  vanthoff_nacl = 2.0_wp                    !< van't Hoff factor for NaCl
    … …  
    143146!
    144147!--    Saturation vapor pressure for a specific temperature:
    145        magnus_0d =  611.2_wp * EXP( 17.62_wp * ( t - 273.15_wp ) /             &
     148       magnus_0d =  611.2_wp * EXP( 17.62_wp * ( t - degc_to_k ) /             &
    146149                                                   ( t - 29.65_wp  ) )
    147150
    … …  
    163166!
    164167!--    Saturation vapor pressure for a specific temperature:
    165        magnus_1d =  611.2_wp * EXP( 17.62_wp * ( t - 273.15_wp ) /             &
     168       magnus_1d =  611.2_wp * EXP( 17.62_wp * ( t - degc_to_k ) /             &
    166169                                               ( t - 29.65_wp  ) )
    167170

  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r3435 r3449  
    2727! -----------------
    2828! $Id$
     29! additional output - merged from branch resler
     30!
     31! 3438 2018-10-28 19:31:42Z pavelkrc
    2932! Add terrain-following masked output
    3033!
    … …  
    262265                                 !< neumann = zero gradient
    263266
    264     REAL(kind=wp), PUBLIC :: cs_time_step
     267    REAL(kind=wp), PUBLIC :: cs_time_step = 0._wp
    265268    CHARACTER(10), PUBLIC :: photolysis_scheme
    266269                                         ! 'constant',
    … …  
    410413       IF ( mode == 'allocate' )  THEN
    411414          DO lsp = 1, nspec
    412              IF (TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
     415             IF ( TRIM(variable(1:3)) == 'kc_' .AND. &
     416                  TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
    413417                chem_species(lsp)%conc_av = 0.0_wp
    414                
    415418             ENDIF
    416419          ENDDO
    … …  
    418421       ELSEIF ( mode == 'sum' )  THEN
    419422
    420 
    421423          DO lsp = 1, nspec
    422              IF (TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
     424             IF ( TRIM(variable(1:3)) == 'kc_' .AND. &
     425                  TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
    423426                DO  i = nxlg, nxrg
    424427                   DO  j = nysg, nyng
    … …  
    462465       ELSEIF ( mode == 'average' )  THEN
    463466          DO lsp = 1, nspec
    464              IF (TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
     467             IF ( TRIM(variable(1:3)) == 'kc_' .AND. &
     468                  TRIM(variable(4:)) == TRIM(chem_species(lsp)%name)) THEN
    465469                DO  i = nxlg, nxrg
    466470                   DO  j = nysg, nyng
    … …  
    811815       unit = 'illegal'
    812816
    813        spec_name = TRIM(var(4:))             !< var 1:3 is 'kc_'.
     817       spec_name = TRIM(var(4:))             !< var 1:3 is 'kc_' or 'em_'.
     818
     819       IF ( TRIM(var(1:3)) == 'em_' )  THEN
     820          DO lsp=1,nspec
     821             IF (TRIM(spec_name) == TRIM(chem_species(lsp)%name))   THEN
     822                unit = 'mol m-2 s-1'
     823             ENDIF
     824             !-- It is possible to plant PM10 and PM25 into the gasphase chemistry code
     825             !-- as passive species (e.g. 'passive' in GASPHASE_PREPROC/mechanisms):
     826             !-- set unit to micrograms per m**3 for PM10 and PM25 (PM2.5)
     827             IF (spec_name(1:2) == 'PM')   THEN
     828                unit = 'kg m-2 s-1'
     829             ENDIF
     830          ENDDO
     831
     832       ELSE
    814833
    815834          DO lsp=1,nspec
    … …  
    821840!            set unit to micrograms per m**3 for PM10 and PM25 (PM2.5)
    822841             IF (spec_name(1:2) == 'PM')   THEN 
    823             unit = 'kg m-3'
     842               unit = 'kg m-3'
    824843             ENDIF
    825844          ENDDO
    826845
    827           DO lsp=1,nphot                                                     
     846          DO lsp=1,nphot
    828847             IF (TRIM(spec_name) == TRIM(phot_frequen(lsp)%name))   THEN
    829848                unit = 'sec-1'
    830849             ENDIF
    831850          ENDDO
    832 
    833 
    834        RETURN
     851       ENDIF
     852
     853
     854      RETURN
    835855    END SUBROUTINE chem_check_data_output
    836856!
    … …  
    933953          CALL message( 'chem_check_parameters', 'CM0426', 1, 2, 0, 6, 0 )
    934954       ENDIF
     955
     956!---------------------
     957!--    chem_check_parameters is called before the array chem_species is allocated!
     958!--    temporary switch of this part of the check
     959       RETURN
     960!---------------------
    935961
    936962!--    check for initial chem species input
    … …  
    10761102       USE kinds
    10771103
     1104       USE surface_mod
     1105
     1106       !USE chem_modules, ONLY: nvar
    10781107
    10791108       IMPLICIT NONE
    … …  
    10911120
    10921121       !-- local variables
    1093 
    1094        INTEGER              ::  i, j, k, lsp
    10951122       CHARACTER(len=16)    ::  spec_name
     1123       INTEGER(iwp)         ::  i, j, k
     1124       INTEGER(iwp)         ::  m, l    !< running indices for surfaces
     1125       INTEGER(iwp)         ::  lsp  !< running index for chem spcs
    10961126
    10971127
    10981128       found = .FALSE.
     1129       IF ( .NOT. (variable(1:3) == 'kc_' .OR. variable(1:3) == 'em_' ) ) THEN
     1130          RETURN
     1131       ENDIF
    10991132
    11001133       spec_name = TRIM(variable(4:))
    11011134
    1102        DO lsp=1,nspec
    1103           IF (TRIM(spec_name) == TRIM(chem_species(lsp)%name))   THEN
     1135       IF ( variable(1:3) == 'em_' ) THEN
     1136
     1137         local_pf = 0.0_wp
     1138
     1139         DO lsp = 1, nvar  !!! cssws - nvar species, chem_species - nspec species !!!
     1140            IF ( TRIM(spec_name) == TRIM(chem_species(lsp)%name) )   THEN
     1141               ! no average for now
     1142               DO m = 1, surf_usm_h%ns
     1143                  local_pf(surf_usm_h%i(m),surf_usm_h%j(m),surf_usm_h%k(m)) = &
     1144                     local_pf(surf_usm_h%i(m),surf_usm_h%j(m),surf_usm_h%k(m)) + surf_usm_h%cssws(lsp,m)
     1145               ENDDO
     1146               DO m = 1, surf_lsm_h%ns
     1147                  local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),surf_lsm_h%k(m)) = &
     1148                     local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),surf_lsm_h%k(m)) + surf_lsm_h%cssws(lsp,m)
     1149               ENDDO
     1150               DO l = 0, 3
     1151                  DO m = 1, surf_usm_v(l)%ns
     1152                     local_pf(surf_usm_v(l)%i(m),surf_usm_v(l)%j(m),surf_usm_v(l)%k(m)) = &
     1153                        local_pf(surf_usm_v(l)%i(m),surf_usm_v(l)%j(m),surf_usm_v(l)%k(m)) + surf_usm_v(l)%cssws(lsp,m)
     1154                  ENDDO
     1155                  DO m = 1, surf_lsm_v(l)%ns
     1156                     local_pf(surf_lsm_v(l)%i(m),surf_lsm_v(l)%j(m),surf_lsm_v(l)%k(m)) = &
     1157                        local_pf(surf_lsm_v(l)%i(m),surf_lsm_v(l)%j(m),surf_lsm_v(l)%k(m)) + surf_lsm_v(l)%cssws(lsp,m)
     1158                  ENDDO
     1159               ENDDO
     1160               found = .TRUE.
     1161            ENDIF
     1162         ENDDO
     1163       ELSE
     1164         DO lsp=1,nspec
     1165            IF (TRIM(spec_name) == TRIM(chem_species(lsp)%name))   THEN
    11041166!
    11051167!--   todo: remove or replace by "CALL message" mechanism (kanani)
    11061168!              IF(myid == 0 .AND. chem_debug0 )  WRITE(6,*) 'Output of species ' // TRIM(variable)  //  &
    11071169!                                                           TRIM(chem_species(lsp)%name)       
    1108              
    1109              IF (av == 0) THEN
    1110                 DO  i = nxl, nxr
    1111                    DO  j = nys, nyn
    1112                       DO  k = nzb_do, nzt_do
    1113                           local_pf(i,j,k) = MERGE(                             &
    1114                                               chem_species(lsp)%conc(k,j,i),   &
    1115                                               REAL( fill_value, KIND = wp ),   &
    1116                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    1117                       ENDDO
    1118                    ENDDO
    1119                 ENDDO
    1120            
    1121              ELSE
    1122                 DO  i = nxl, nxr
    1123                    DO  j = nys, nyn
    1124                       DO  k = nzb_do, nzt_do
    1125                           local_pf(i,j,k) = MERGE(                             &
    1126                                               chem_species(lsp)%conc_av(k,j,i),&
    1127                                               REAL( fill_value, KIND = wp ),   &
    1128                                               BTEST( wall_flags_0(k,j,i), 0 ) )
    1129                       ENDDO
    1130                    ENDDO
    1131                 ENDDO
    1132              ENDIF
    1133 
    1134              found = .TRUE.
    1135           ENDIF
    1136        ENDDO
     1170               IF (av == 0) THEN
     1171                  DO  i = nxl, nxr
     1172                     DO  j = nys, nyn
     1173                        DO  k = nzb_do, nzt_do
     1174                            local_pf(i,j,k) = MERGE(                             &
     1175                                                chem_species(lsp)%conc(k,j,i),   &
     1176                                                REAL( fill_value, KIND = wp ),   &
     1177                                                BTEST( wall_flags_0(k,j,i), 0 ) )
     1178                        ENDDO
     1179                     ENDDO
     1180                  ENDDO
     1181
     1182               ELSE
     1183                  DO  i = nxl, nxr
     1184                     DO  j = nys, nyn
     1185                        DO  k = nzb_do, nzt_do
     1186                            local_pf(i,j,k) = MERGE(                             &
     1187                                                chem_species(lsp)%conc_av(k,j,i),&
     1188                                                REAL( fill_value, KIND = wp ),   &
     1189                                                BTEST( wall_flags_0(k,j,i), 0 ) )
     1190                        ENDDO
     1191                     ENDDO
     1192                  ENDDO
     1193               ENDIF
     1194               found = .TRUE.
     1195            ENDIF
     1196         ENDDO
     1197       ENDIF
    11371198
    11381199       RETURN
     1200
    11391201    END SUBROUTINE chem_data_output_3d
    11401202!
    … …  
    11771239
    11781240       spec_name = TRIM( variable(4:) )
    1179    !av == 0
     1241       !av == 0
    11801242
    11811243       DO lsp=1,nspec
    … …  
    12921354       found  = .TRUE.
    12931355
    1294        IF (var(1:3) == 'kc_')   THEN                   !< always the same grid for chemistry variables
     1356       IF ( var(1:3) == 'kc_' .OR. var(1:3) == 'em_' )   THEN                   !< always the same grid for chemistry variables
    12951357          grid_x = 'x'
    12961358          grid_y = 'y'
    … …  
    18731935 !--   Enable chemistry model
    18741936       air_chemistry = .TRUE.                   
    1875       GOTO 20
    1876 
    1877  10   BACKSPACE( 11 )
    1878       READ( 11 , '(A)') line
    1879       CALL parin_fail_message( 'chemistry_parameters', line )
    1880 
    1881  20   CONTINUE
     1937       GOTO 20
     1938
     1939 10    BACKSPACE( 11 )
     1940       READ( 11 , '(A)') line
     1941       CALL parin_fail_message( 'chemistry_parameters', line )
     1942
     1943 20    CONTINUE
    18821944
    18831945!

  • palm/trunk/SOURCE/plant_canopy_model_mod.f90

    r3337 r3449  
    1616!
    1717! Copyright 1997-2018 Leibniz Universitaet Hannover
     18! Copyright 2018 Institute of Computer Science of the
     19!                Czech Academy of Sciences, Prague
    1820!------------------------------------------------------------------------------!
    1921!
    … …  
    2527! -----------------
    2628! $Id$
     29! Add calculation of transpiration for resolved plant canopy (trees, shrubs)
     30! (V. Fuka, MFF UK Prague, J.Resler, ICS AS, Prague)
     31!
     32! Fix reading plant canopy over buildings
     33!
     34! 3337 2018-10-12 15:17:09Z kanani
    2735! Fix reading plant canopy over buildings
    2836!
    … …  
    214222 
    215223    USE arrays_3d,                                                             &
    216         ONLY:  dzu, dzw, e, q, s, tend, u, v, w, zu, zw
     224        ONLY:  dzu, dzw, e, exner, hyp, pt, q, s, tend, u, v, w, zu, zw
     225
     226    USE basic_constants_and_equations_mod,                                     &
     227        ONLY:  c_p, degc_to_k, l_v, lv_d_cp, r_d, rd_d_rv
     228
     229    USE control_parameters,                                                   &
     230        ONLY: humidity
    217231
    218232    USE indices,                                                               &
    … …  
    222236    USE kinds
    223237
     238    USE pegrid
     239
    224240    USE surface_mod,                                                           &
    225241        ONLY:  get_topography_top_index_ji
    … …  
    229245
    230246
    231     CHARACTER (LEN=30)   ::  canopy_mode = 'block' !< canopy coverage
    232 
    233     INTEGER(iwp) ::  pch_index = 0                               !< plant canopy height/top index
    234     INTEGER(iwp) ::  lad_vertical_gradient_level_ind(10) = -9999 !< lad-profile levels (index)
    235 
    236     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  pch_index_ji   !< local plant canopy top
    237 
    238     LOGICAL ::  calc_beta_lad_profile = .FALSE. !< switch for calc. of lad from beta func.
     247    CHARACTER (LEN=30)   ::  canopy_mode = 'block'                 !< canopy coverage
     248    LOGICAL              ::  plant_canopy_transpiration = .FALSE.  !< flag to switch calculation of transpiration and corresponding latent heat
     249                                                                   !< for resolved plant canopy inside radiation model
     250                                                                   !< (calls subroutine pcm_calc_transpiration_rate from module plant_canopy_mod)
     251
     252    INTEGER(iwp) ::  pch_index = 0                                 !< plant canopy height/top index
     253    INTEGER(iwp) ::  lad_vertical_gradient_level_ind(10) = -9999   !< lad-profile levels (index)
     254
     255    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  pch_index_ji     !< local plant canopy top
     256
     257    LOGICAL ::  calc_beta_lad_profile = .FALSE.   !< switch for calc. of lad from beta func.
    239258
    240259    REAL(wp) ::  alpha_lad = 9999999.9_wp         !< coefficient for lad calculation
    … …  
    261280    REAL(wp), DIMENSION(:), ALLOCATABLE ::  pre_lad        !< preliminary lad
    262281   
    263     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  cum_lai_hf       !< cumulative lai for heatflux calc.
    264     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  lad_s            !< lad on scalar-grid
    265     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pc_heating_rate  !< plant canopy heating rate
     282    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  cum_lai_hf             !< cumulative lai for heatflux calc.
     283    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  lad_s                  !< lad on scalar-grid
     284    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pc_heating_rate        !< plant canopy heating rate
    266285    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pc_transpiration_rate  !< plant canopy transpiration rate
     286    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pc_latent_rate         !< plant canopy latent heating rate
    267287
    268288    SAVE
    … …  
    273293!
    274294!-- Public functions
    275     PUBLIC pcm_check_data_output, pcm_check_parameters, pcm_data_output_3d,    &
    276            pcm_define_netcdf_grid, pcm_header, pcm_init, pcm_parin, pcm_tendency
     295    PUBLIC pcm_calc_transpiration_rate, pcm_check_data_output,                &
     296           pcm_check_parameters, pcm_data_output_3d, pcm_define_netcdf_grid,  &
     297           pcm_header, pcm_init, pcm_parin, pcm_tendency
    277298
    278299!
    279300!-- Public variables and constants
    280     PUBLIC pc_heating_rate, pc_transpiration_rate, canopy_mode, cthf, dt_plant_canopy, lad, lad_s,   &
    281            pch_index
    282            
     301    PUBLIC pc_heating_rate, pc_transpiration_rate, pc_latent_rate,             &
     302            canopy_mode, cthf, dt_plant_canopy, lad, lad_s, pch_index,         &
     303            plant_canopy_transpiration
     304
     305    INTERFACE pcm_calc_transpiration_rate
     306       MODULE PROCEDURE pcm_calc_transpiration_rate
     307    END INTERFACE pcm_calc_transpiration_rate
    283308
    284309    INTERFACE pcm_check_data_output
    … …  
    323348
    324349
     350
     351!------------------------------------------------------------------------------!
     352! Description:
     353! ------------
     354!> Calculation of the plant canopy transpiration rate based on the Jarvis-Stewart
     355!> with parametrizations described in Daudet et al. (1999; Agricult. and Forest
     356!> Meteorol. 97) and Ngao, Adam and Saudreau (2017;  Agricult. and Forest Meteorol
     357!> 237-238). Model functions f1-f4 were adapted from Stewart (1998; Agric.
     358!> and Forest. Meteorol. 43) instead, because they are valid for broader intervals
     359!> of values. Funcion f4 used in form present in van Wijk et al. (1998;
     360!> Tree Physiology 20).
     361!>
     362!> This subroutine is called from subroutine radiation_interaction
     363!> after the calculation of radiation in plant canopy boxes.
     364!> (arrays pcbinsw and pcbinlw).
     365!>
     366!------------------------------------------------------------------------------!
     367 SUBROUTINE pcm_calc_transpiration_rate(i, j, k, kk, pcbsw, pcblw, pcbtr, pcblh)
     368
     369     USE control_parameters,                                                   &
     370        ONLY: dz
     371
     372     USE grid_variables,                                                       &
     373        ONLY:  dx, dy
     374
     375     IMPLICIT NONE
     376!--  input parameters
     377     INTEGER(iwp), INTENT(IN)          :: i, j, k, kk        !< indices of the pc gridbox
     378     REAL(wp), INTENT(IN)              :: pcbsw              !< sw radiation in gridbox (W)
     379     REAL(wp), INTENT(IN)              :: pcblw              !< lw radiation in gridbox (W)
     380     REAL(wp), INTENT(OUT)             :: pcbtr              !< transpiration rate dq/dt (kg/kg/s)
     381     REAL(wp), INTENT(OUT)             :: pcblh              !< latent heat from transpiration dT/dt (K/s)
     382
     383!--  variables and parameters for calculation of transpiration rate
     384     REAL(wp)                          :: sat_press, sat_press_d, temp, v_lad
     385     REAL(wp)                          :: d_fact, g_b, g_s, wind_speed, evapor_rate
     386     REAL(wp)                          :: f1, f2, f3, f4, vpd, rswc, e_eq, e_imp, rad
     387     REAL(wp), PARAMETER               :: gama_psychr = 66   !< psychrometric constant (Pa/K)
     388     REAL(wp), PARAMETER               :: g_s_max = 0.01     !< maximum stomatal conductivity (m/s)
     389     REAL(wp), PARAMETER               :: m_soil = 0.4_wp    !< soil water content (needs to adjust or take from LSM)
     390     REAL(wp), PARAMETER               :: m_wilt = 0.01_wp   !< wilting point soil water content (needs to adjust or take from LSM)
     391     REAL(wp), PARAMETER               :: m_sat = 0.51_wp    !< saturation soil water content (needs to adjust or take from LSM)
     392     REAL(wp), PARAMETER               :: t2_min = 0.0_wp    !< minimal temperature for calculation of f2
     393     REAL(wp), PARAMETER               :: t2_max = 40.0_wp   !< maximal temperature for calculation of f2
     394
     395
     396!--  Temperature (deg C)
     397     temp = pt(k,j,i) * exner(k) - degc_to_k
     398!--  Coefficient for conversion of radiation to grid to radiation to unit leaves surface
     399     v_lad = 1.0_wp / ( MAX(lad_s(kk,j,i), 1.0e-10) * dx * dy * dz(1) )
     400!--  Magnus formula for the saturation pressure (see Ngao, Adam and Saudreau (2017) eq. 1)
     401!--  There are updated formulas available, kept consistent with the rest of the parametrization
     402     sat_press = 610.8_wp * exp(17.27_wp * temp/(temp + 237.3_wp))
     403!--  Saturation pressure derivative (derivative of the above)
     404     sat_press_d = sat_press * 17.27_wp * 237.3_wp / (temp + 237.3_wp)**2
     405!--  Wind speed
     406     wind_speed = SQRT( ((u(k,j,i) + u(k,j,i+1))/2.0_wp )**2 +           &
     407                        ((v(k,j,i) + v(k,j,i+1))/2.0_wp )**2 +           &
     408                        ((w(k,j,i) + w(k,j,i+1))/2.0_wp )**2 )
     409!--  Aerodynamic conductivity (Daudet et al. (1999) eq. 14
     410     g_b = 0.01_wp * wind_speed + 0.0071_wp
     411!--  Radiation flux per leaf surface unit
     412     rad = pcbsw * v_lad
     413!--  First function for calculation of stomatal conductivity (radiation dependency)
     414!--  Stewart (1988; Agric. and Forest. Meteorol. 43) eq. 17
     415     f1 = rad * (1000._wp+42.1_wp) / 1000._wp / (rad+42.1_wp)
     416!--  Second function for calculation of stomatal conductivity (temperature dependency)
     417!--  Stewart (1988; Agric. and Forest. Meteorol. 43) eq. 21
     418     f2 = MAX(t2_min, (temp-t2_min) * (t2_max-temp)**((t2_max-16.9_wp)/(16.9_wp-t2_min)) / &
     419              ((16.9_wp-t2_min) * (t2_max-16.9_wp)**((t2_max-16.9_wp)/(16.9_wp-t2_min))) )
     420!--  Water pressure deficit
     421!--  Ngao, Adam and Saudreau (2017) eq. 6 but with water vapour partial pressure
     422     vpd = max( sat_press - q(k,j,i) * hyp(k) / rd_d_rv, 0._wp )
     423!--  Third function for calculation of stomatal conductivity (water pressure deficit dependency)
     424!--  Ngao, Adam and Saudreau (2017) Table 1, limited from below according to Stewart (1988)
     425!--  The coefficients of the linear dependence should better correspond to broad-leaved trees
     426!--  than the coefficients from Stewart (1988) which correspond to conifer trees.
     427     vpd = MIN(MAX(vpd,770.0_wp),3820.0_wp)
     428     f3 = -2e-4_wp * vpd + 1.154_wp
     429!--  Fourth function for calculation of stomatal conductivity (soil moisture dependency)
     430!--  Residual soil water content
     431!--  van Wijk et al. (1998; Tree Physiology 20) eq. 7
     432!--  TODO - over LSM surface might be calculated from LSM parameters
     433     rswc = ( m_sat - m_soil ) / ( m_sat - m_wilt )
     434!--  van Wijk et al. (1998; Tree Physiology 20) eq. 5-6 (it is a reformulation of eq. 22-23 of Stewart(1988))
     435     f4 = MAX(0._wp, MIN(1.0_wp - 0.041_wp * EXP(3.2_wp * rswc), 1.0_wp - 0.041_wp))
     436!--  Stomatal conductivity
     437!--  Stewart (1988; Agric. and Forest. Meteorol. 43) eq. 12
     438!--  (notation according to Ngao, Adam and Saudreau (2017) and others)
     439     g_s = g_s_max * f1 * f2 * f3 * f4 + 1.0e-10_wp
     440!--  Decoupling factor
     441!--  Daudet et al. (1999) eq. 6
     442     d_fact = (sat_press_d / gama_psychr + 2._wp ) /                        &
     443              (sat_press_d / gama_psychr + 2._wp + 2 * g_b / g_s )
     444!--  Equilibrium evaporation rate
     445!--  Daudet et al. (1999) eq. 4
     446     e_eq = (pcbsw + pcblw) * v_lad * sat_press_d /                         &
     447                 gama_psychr /( sat_press_d / gama_psychr + 2.0_wp ) / l_v
     448!--  Imposed evaporation rate
     449!--  Daudet et al. (1999) eq. 5
     450     e_imp = r_d * pt(k,j,i) * exner(k) / hyp(k) * c_p * g_s * vpd / gama_psychr / l_v
     451!--  Evaporation rate
     452!--  Daudet et al. (1999) eq. 3
     453!--  (evaporation rate is limited to non-negative values)
     454     evapor_rate = MAX(d_fact * e_eq + ( 1.0_wp - d_fact ) * e_imp, 0.0_wp)
     455!--  Conversion of evaporation rate to q tendency in gridbox
     456!--  dq/dt = E * LAD * V_g / (rho_air * V_g)
     457     pcbtr = evapor_rate * r_d * pt(k,j,i) * exner(k) * lad_s(kk,j,i) / hyp(k)  !-- = dq/dt
     458!--  latent heat from evaporation
     459     pcblh = pcbtr * lv_d_cp  !-- = - dT/dt
     460
     461 END SUBROUTINE pcm_calc_transpiration_rate
     462
     463
    325464!------------------------------------------------------------------------------!
    326465! Description:
    … …  
    352491       CASE ( 'pcm_transpirationrate' )
    353492          unit = 'kg kg-1 s-1'
     493
     494       CASE ( 'pcm_latentrate' )
     495          unit = 'K s-1'
     496
     497       CASE ( 'pcm_bowenratio' )
     498          unit = 'K s-1'
    354499
    355500       CASE ( 'pcm_lad' )
    … …  
    503648            ENDDO
    504649         ENDIF
    505    
     650
     651       CASE ( 'pcm_latentrate' )
     652         IF ( av == 0 )  THEN
     653            DO  i = nxl, nxr
     654               DO  j = nys, nyn
     655                  IF ( pch_index_ji(j,i) /= 0 )  THEN
     656                     k_topo = get_topography_top_index_ji( j, i, 's' )
     657                     DO  k = k_topo, k_topo + pch_index_ji(j,i)
     658                        local_pf(i,j,k) = pc_latent_rate(k-k_topo,j,i)
     659                     ENDDO
     660                  ENDIF
     661               ENDDO
     662            ENDDO
     663         ENDIF
     664
     665       CASE ( 'pcm_bowenratio' )
     666         IF ( av == 0 )  THEN
     667            DO  i = nxl, nxr
     668               DO  j = nys, nyn
     669                  IF ( pch_index_ji(j,i) /= 0 )  THEN
     670                     k_topo = get_topography_top_index_ji( j, i, 's' )
     671                     DO  k = k_topo, k_topo + pch_index_ji(j,i)
     672                        IF ( pc_latent_rate(k-k_topo,j,i) /= 0._wp ) THEN
     673                           local_pf(i,j,k) = pc_heating_rate(k-k_topo,j,i) / &
     674                                             pc_latent_rate(k-k_topo,j,i)
     675                        ENDIF
     676                     ENDDO
     677                  ENDIF
     678               ENDDO
     679            ENDDO
     680         ENDIF
     681
    506682      CASE ( 'pcm_lad' )
    507683         IF ( av == 0 )  THEN
    … …  
    550726     SELECT CASE ( TRIM( var ) )
    551727
    552         CASE ( 'pcm_heatrate', 'pcm_lad', 'pcm_transpirationrate')
     728        CASE ( 'pcm_heatrate', 'pcm_lad', 'pcm_transpirationrate', 'pcm_latentrate', 'pcm_bowenratio')
    553729           grid_x = 'x'
    554730           grid_y = 'y'
    … …  
    691867
    692868       USE control_parameters,                                                 &
    693            ONLY: humidity, message_string, ocean_mode, urban_surface
     869           ONLY: message_string, ocean_mode, urban_surface
    694870
    695871       USE netcdf_data_input_mod,                                              &
    … …  
    705881       INTEGER(iwp) ::  k   !< running index
    706882       INTEGER(iwp) ::  m   !< running index
     883       INTEGER(iwp) ::  pch_index_l
    707884
    708885       REAL(wp) ::  int_bpdf        !< vertical integral for lad-profile construction
    … …  
    9031080       CALL exchange_horiz_2d_int( pch_index_ji, nys, nyn, nxl, nxr, nbgp )
    9041081
     1082!--    Calculate global pch_index value (index of top of plant canopy from ground)
     1083       pch_index = MAXVAL(pch_index_ji)
     1084!--    Exchange pch_index from all processors
     1085#if defined( __parallel )
     1086       CALL MPI_Allreduce(MPI_IN_PLACE, pch_index, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr)
     1087       IF ( ierr /= 0 )  THEN
     1088          WRITE (9,*) 'Error MPI_Allreduce pch_index:', ierr
     1089          FLUSH(9)
     1090       ENDIF
     1091#endif
     1092
     1093!--    Allocation of arrays pc_heating_rate, pc_transpiration_rate and pc_latent_rate
     1094       ALLOCATE( pc_heating_rate(0:pch_index,nysg:nyng,nxlg:nxrg) )
     1095       IF ( humidity )  THEN
     1096          ALLOCATE( pc_transpiration_rate(0:pch_index,nysg:nyng,nxlg:nxrg) )
     1097          pc_transpiration_rate = 0.0_wp
     1098          ALLOCATE( pc_latent_rate(0:pch_index,nysg:nyng,nxlg:nxrg) )
     1099          pc_latent_rate = 0.0_wp
     1100       ENDIF
     1101
    9051102!
    9061103!--    Initialization of the canopy heat source distribution due to heating
    … …  
    9151112!--    73–96). This approach has been applied e.g. by Shaw & Schumann (1992;
    9161113!--    Bound.-Layer Meteorol. 61, 47–64).
    917 !--    When using the urban surface model (USM), canopy heating (pc_heating_rate)
    918 !--    by radiation is calculated in the USM.
    919        IF ( cthf /= 0.0_wp  .AND. .NOT.  urban_surface )  THEN
    920 
    921           ALLOCATE( cum_lai_hf(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                 &
    922                     pc_heating_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     1114!--    When using the radiation_interactions, canopy heating (pc_heating_rate)
     1115!--    and plant canopy transpiration (pc_transpiration_rate, pc_latent_rate)
     1116!--    are calculated in the RTM after the calculation of radiation.
     1117!--    We cannot use variable radiation_interactions here to determine the situation
     1118!--    as it is assigned in init_3d_model after the call of pcm_init.
     1119       IF ( cthf /= 0.0_wp )  THEN
     1120
     1121          ALLOCATE( cum_lai_hf(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    9231122!
    9241123!--       Piecewise calculation of the cumulative leaf area index by vertical
    … …  
    9961195
    9971196       ENDIF
    998 !
    999 !--    Allocate transpiration rate
    1000        IF ( humidity )                                                         &
    1001           ALLOCATE( pc_transpiration_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    1002 
    10031197
    10041198
    … …  
    10281222                                  lai_beta,                                    &
    10291223                                  leaf_scalar_exch_coeff,                      &
    1030                                   leaf_surface_conc, pch_index
     1224                                  leaf_surface_conc, pch_index,                &
     1225                                  plant_canopy_transpiration
    10311226
    10321227       NAMELIST /canopy_par/      alpha_lad, beta_lad, canopy_drag_coeff,      &
    … …  
    10371232                                  lai_beta,                                    &
    10381233                                  leaf_scalar_exch_coeff,                      &
    1039                                   leaf_surface_conc, pch_index
     1234                                  leaf_surface_conc, pch_index,                &
     1235                                  plant_canopy_transpiration
    10401236
    10411237       line = ' '
    … …  
    14341630!--       potential temperature
    14351631          CASE ( 4 )
    1436              DO  i = nxl, nxr
    1437                 DO  j = nys, nyn
    1438 !
    1439 !--                Determine topography-top index on scalar-grid
    1440                    k_wall = get_topography_top_index_ji( j, i, 's' )
    1441 
    1442                    DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
    1443 
    1444                       kk = k - k_wall   !- lad arrays are defined flat
    1445                       tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i)
     1632             IF ( humidity ) THEN
     1633                DO  i = nxl, nxr
     1634                   DO  j = nys, nyn
     1635!--                   Determine topography-top index on scalar-grid
     1636                      k_wall = get_topography_top_index_ji( j, i, 's' )
     1637                      DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     1638                         kk = k - k_wall   !- lad arrays are defined flat
     1639                         tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i) - pc_latent_rate(kk,j,i)
     1640                      ENDDO
    14461641                   ENDDO
    14471642                ENDDO
    1448              ENDDO
     1643             ELSE
     1644                DO  i = nxl, nxr
     1645                   DO  j = nys, nyn
     1646!--                   Determine topography-top index on scalar-grid
     1647                      k_wall = get_topography_top_index_ji( j, i, 's' )
     1648                      DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     1649                         kk = k - k_wall   !- lad arrays are defined flat
     1650                         tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i)
     1651                      ENDDO
     1652                   ENDDO
     1653                ENDDO
     1654             ENDIF
    14491655
    14501656!
    … …  
    14601666
    14611667                      kk = k - k_wall   !- lad arrays are defined flat
    1462                       pc_transpiration_rate(kk,j,i) =  - lsec                  &
    1463                                        * lad_s(kk,j,i) *                       &
    1464                                        SQRT( ( 0.5_wp * ( u(k,j,i) +           &
    1465                                                           u(k,j,i+1) )         &
    1466                                              )**2 +                            &
    1467                                              ( 0.5_wp * ( v(k,j,i) +           &
    1468                                                           v(k,j+1,i) )         &
    1469                                              )**2 +                            &
    1470                                              ( 0.5_wp * ( w(k-1,j,i) +         &
    1471                                                           w(k,j,i) )           &
    1472                                              )**2                              &
    1473                                            ) *                                 &
    1474                                        ( q(k,j,i) - lsc )
     1668
     1669                      IF ( .NOT. plant_canopy_transpiration ) THEN
     1670                         ! pc_transpiration_rate is calculated in radiation model
     1671                         ! in case of plant_canopy_transpiration = .T.
     1672                         ! to include also the dependecy to the radiation
     1673                         ! in the plant canopy box
     1674                         pc_transpiration_rate(kk,j,i) =  - lsec                  &
     1675                                          * lad_s(kk,j,i) *                       &
     1676                                          SQRT( ( 0.5_wp * ( u(k,j,i) +           &
     1677                                                             u(k,j,i+1) )         &
     1678                                                )**2 +                            &
     1679                                                ( 0.5_wp * ( v(k,j,i) +           &
     1680                                                             v(k,j+1,i) )         &
     1681                                                )**2 +                            &
     1682                                                ( 0.5_wp * ( w(k-1,j,i) +         &
     1683                                                             w(k,j,i) )           &
     1684                                                )**2                              &
     1685                                              ) *                                 &
     1686                                          ( q(k,j,i) - lsc )
     1687                      ENDIF
    14751688
    14761689                      tend(k,j,i) = tend(k,j,i) + pc_transpiration_rate(kk,j,i)
    … …  
    17781991             k_wall = get_topography_top_index_ji( j, i, 's' )
    17791992
     1993             IF ( humidity ) THEN
     1994                DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     1995                   kk = k - k_wall  !- lad arrays are defined flat
     1996                   tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i) - pc_latent_rate(kk,j,i)
     1997                ENDDO
     1998             ELSE
     1999                DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     2000                   kk = k - k_wall  !- lad arrays are defined flat
     2001                   tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i)
     2002                ENDDO
     2003             ENDIF
     2004
     2005!
     2006!--       humidity
     2007          CASE ( 5 )
     2008!
     2009!--          Determine topography-top index on scalar grid
     2010             k_wall = get_topography_top_index_ji( j, i, 's' )
     2011
    17802012             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
    17812013                kk = k - k_wall  !- lad arrays are defined flat
    1782                 tend(k,j,i) = tend(k,j,i) + pc_heating_rate(kk,j,i)
    1783              ENDDO
    1784 
    1785 
    1786 !
    1787 !--       humidity
    1788           CASE ( 5 )
    1789 !
    1790 !--          Determine topography-top index on scalar grid
    1791              k_wall = get_topography_top_index_ji( j, i, 's' )
    1792 
    1793              DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
    1794                 kk = k - k_wall  !- lad arrays are defined flat
    1795 
    1796                 pc_transpiration_rate(kk,j,i) = - lsec                         &
    1797                                  * lad_s(kk,j,i) *                             &
    1798                                  SQRT( ( 0.5_wp * ( u(k,j,i) +                 &
    1799                                                     u(k,j,i+1) )               &
    1800                                        )**2  +                                 &
    1801                                        ( 0.5_wp * ( v(k,j,i) +                 &
    1802                                                     v(k,j+1,i) )               &
    1803                                        )**2 +                                  &
    1804                                        ( 0.5_wp * ( w(k-1,j,i) +               &
    1805                                                     w(k,j,i) )                 &
    1806                                        )**2                                    &
    1807                                      ) *                                       &
    1808                                  ( q(k,j,i) - lsc )
     2014                IF ( .NOT. plant_canopy_transpiration ) THEN
     2015                   ! pc_transpiration_rate is calculated in radiation model
     2016                   ! in case of plant_canopy_transpiration = .T.
     2017                   ! to include also the dependecy to the radiation
     2018                   ! in the plant canopy box
     2019                   pc_transpiration_rate(kk,j,i) = - lsec                         &
     2020                                    * lad_s(kk,j,i) *                             &
     2021                                    SQRT( ( 0.5_wp * ( u(k,j,i) +                 &
     2022                                                       u(k,j,i+1) )               &
     2023                                          )**2  +                                 &
     2024                                          ( 0.5_wp * ( v(k,j,i) +                 &
     2025                                                       v(k,j+1,i) )               &
     2026                                          )**2 +                                  &
     2027                                          ( 0.5_wp * ( w(k-1,j,i) +               &
     2028                                                       w(k,j,i) )                 &
     2029                                          )**2                                    &
     2030                                        ) *                                       &
     2031                                    ( q(k,j,i) - lsc )
     2032                ENDIF
    18092033
    18102034                tend(k,j,i) = tend(k,j,i) + pc_transpiration_rate(kk,j,i)

  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3435 r3449  
    2828! -----------------
    2929! $Id$
     30! New RTM version 3.0: (Pavel Krc, Jaroslav Resler, ICS, Prague)
     31!   - Interaction of plant canopy with LW radiation
     32!   - Transpiration from resolved plant canopy dependent on radiation
     33!     called from RTM
     34!
     35!
     36! 3435 2018-10-26 18:25:44Z gronemeier
    3037! - workaround: return unit=illegal in check_data_output for certain variables
    3138!   when check called from init_masks
    … …  
    499506
    500507    USE basic_constants_and_equations_mod,                                     &
    501         ONLY:  c_p, g, lv_d_cp, l_v, pi, r_d, rho_l, solar_constant,           &
     508        ONLY:  c_p, g, lv_d_cp, l_v, pi, r_d, rho_l, solar_constant,      &
    502509               barometric_formula
    503510
    … …  
    544551
    545552    USE plant_canopy_model_mod,                                                &
    546         ONLY:  lad_s, pc_heating_rate, pc_transpiration_rate
     553        ONLY:  lad_s, pc_heating_rate, pc_transpiration_rate, pc_latent_rate,  &
     554               plant_canopy_transpiration, pcm_calc_transpiration_rate
    547555
    548556    USE pegrid
    … …  
    853861    INTEGER(iwp), DIMENSION(0:nsurf_type), PARAMETER ::  jdir = (/0, 0,1,-1,0, 0,0,1,-1,0, 0/)   !< surface normal direction y indices
    854862    INTEGER(iwp), DIMENSION(0:nsurf_type), PARAMETER ::  kdir = (/1,-1,0, 0,0, 0,1,0, 0,0, 0/)   !< surface normal direction z indices
    855                                                                                           !< parameter but set in the code
     863    REAL(wp),     DIMENSION(0:nsurf_type)          :: facearea                            !< area of single face in respective
     864                                                                                          !< direction (will be calc'd)
    856865
    857866
    … …  
    887896!-- configuration parameters (they can be setup in PALM config)
    888897    LOGICAL                                        ::  raytrace_mpi_rma = .TRUE.          !< use MPI RMA to access LAD and gridsurf from remote processes during raytracing
    889     LOGICAL                                        ::  read_svf_on_init = .FALSE.         !< flag parameter indicating wheather SVFs will be read from a file at initialization
    890     LOGICAL                                        ::  write_svf_on_init = .FALSE.        !< flag parameter indicating wheather SVFs will be written out to a file
    891     LOGICAL                                        ::  rad_angular_discretization = .TRUE.!< whether to use fixed resolution discretization of view factors for
     898    LOGICAL                                        ::  rad_angular_discretization = .TRUE.!< whether to use fixed resolution discretization of view factors for
    892899                                                                                          !< reflected radiation (as opposed to all mutually visible pairs)
     900    LOGICAL                                        ::  plant_lw_interact = .TRUE.         !< whether plant canopy interacts with LW radiation (in addition to SW)
    893901    INTEGER(iwp)                                   ::  mrt_nlevels = 0                    !< number of vertical boxes above surface for which to calculate MRT
    894902    LOGICAL                                        ::  mrt_skip_roof = .TRUE.             !< do not calculate MRT above roof surfaces
    895903    LOGICAL                                        ::  mrt_include_sw = .TRUE.            !< should MRT calculation include SW radiation as well?
    896     INTEGER(iwp)                                   ::  nrefsteps = 0                      !< number of reflection steps to perform
     904    INTEGER(iwp)                                   ::  nrefsteps = 3                      !< number of reflection steps to perform
    897905    REAL(wp), PARAMETER                            ::  ext_coef = 0.6_wp                  !< extinction coefficient (a.k.a. alpha)
    898906    INTEGER(iwp), PARAMETER                        ::  rad_version_len = 10               !< length of identification string of rad version
    899     CHARACTER(rad_version_len), PARAMETER          ::  rad_version = 'RAD v. 1.1'         !< identification of version of binary svf and restart files
     907    CHARACTER(rad_version_len), PARAMETER          ::  rad_version = 'RAD v. 3.0'         !< identification of version of binary svf and restart files
    900908    INTEGER(iwp)                                   ::  raytrace_discrete_elevs = 40       !< number of discretization steps for elevation (nadir to zenith)
    901909    INTEGER(iwp)                                   ::  raytrace_discrete_azims = 80       !< number of discretization steps for azimuth (out of 360 degrees)
    … …  
    929937        INTEGER(iwp)                               :: ity               !<
    930938        INTEGER(iwp)                               :: itz               !<
    931         INTEGER(iwp)                               :: isurfs            !<
    932         REAL(wp)                                   :: rsvf              !<
     939        INTEGER(iwp)                               :: isurfs            !< Idx of source face / -1 for sky
     940        REAL(wp)                                   :: rcvf              !< Canopy view factor for faces /
     941                                                                        !< canopy sink factor for sky (-1)
    933942    END TYPE
    934943
    … …  
    976985    REAL(wp), DIMENSION(:), ALLOCATABLE            ::  surfouts         !< array of reflected sw radiation for all surfaces in i-th reflection
    977986    REAL(wp), DIMENSION(:), ALLOCATABLE            ::  surfoutl         !< array of reflected + emitted lw radiation for all surfaces in i-th reflection
     987    REAL(wp), DIMENSION(:), ALLOCATABLE            ::  surfinlg         !< global array of incoming lw radiation from plant canopy
    978988    REAL(wp), DIMENSION(:), ALLOCATABLE            ::  surfoutsw        !< array of total sw radiation outgoing from nonvirtual surfaces surfaces after all reflection
    979989    REAL(wp), DIMENSION(:), ALLOCATABLE            ::  surfoutlw        !< array of total lw radiation outgoing from nonvirtual surfaces surfaces after all reflection
    … …  
    29282938       CHARACTER (LEN=80) ::  line  !< dummy string that contains the current line of the parameter file
    29292939       
    2930        NAMELIST /radiation_par/   albedo, albedo_type, albedo_lw_dir,          &
    2931                                   albedo_lw_dif, albedo_sw_dir, albedo_sw_dif, &
    2932                                   constant_albedo, dt_radiation, emissivity,   &
    2933                                   lw_radiation, mrt_nlevels, mrt_skip_roof,    &
    2934                                   mrt_include_sw,  net_radiation,              &
    2935                                   radiation_scheme, skip_time_do_radiation,    &
    2936                                   sw_radiation, unscheduled_radiation_calls,   &
    2937                                   max_raytracing_dist, min_irrf_value,         &
    2938                                   nrefsteps, raytrace_mpi_rma,                 &
    2939                                   surface_reflections, svfnorm_report_thresh,  &
    2940                                   radiation_interactions_on,                   &
    2941                                   rad_angular_discretization,                  &
    2942                                   raytrace_discrete_azims,                     &
    2943                                   raytrace_discrete_elevs
     2940       NAMELIST /radiation_par/   albedo, albedo_lw_dif, albedo_lw_dir,         &
     2941                                  albedo_sw_dif, albedo_sw_dir, albedo_type,    &
     2942                                  constant_albedo, dt_radiation, emissivity,    &
     2943                                  lw_radiation, max_raytracing_dist,            &
     2944                                  min_irrf_value, mrt_include_sw, mrt_nlevels,  &
     2945                                  mrt_skip_roof, net_radiation, nrefsteps,      &
     2946                                  plant_lw_interact, rad_angular_discretization,&
     2947                                  radiation_interactions_on, radiation_scheme,  &
     2948                                  raytrace_discrete_azims,                      &
     2949                                  raytrace_discrete_elevs, raytrace_mpi_rma,    &
     2950                                  skip_time_do_radiation, surface_reflections,  &
     2951                                  svfnorm_report_thresh, sw_radiation,          &
     2952                                  unscheduled_radiation_calls
     2953
    29442954   
    2945        NAMELIST /radiation_parameters/   albedo, albedo_type, albedo_lw_dir,   &
    2946                                   albedo_lw_dif, albedo_sw_dir, albedo_sw_dif, &
    2947                                   constant_albedo, dt_radiation, emissivity,   &
    2948                                   lw_radiation, mrt_nlevels, mrt_skip_roof,    &
    2949                                   mrt_include_sw,  net_radiation,              &
    2950                                   radiation_scheme, skip_time_do_radiation,    &
    2951                                   sw_radiation, unscheduled_radiation_calls,   &
    2952                                   max_raytracing_dist, min_irrf_value,         &
    2953                                   nrefsteps, raytrace_mpi_rma,                 &
    2954                                   surface_reflections, svfnorm_report_thresh,  &
    2955                                   radiation_interactions_on,                   &
    2956                                   rad_angular_discretization,                  &
    2957                                   raytrace_discrete_azims,                     &
    2958                                   raytrace_discrete_elevs
     2955       NAMELIST /radiation_parameters/ albedo, albedo_lw_dif, albedo_lw_dir,    &
     2956                                  albedo_sw_dif, albedo_sw_dir, albedo_type,    &
     2957                                  constant_albedo, dt_radiation, emissivity,    &
     2958                                  lw_radiation, max_raytracing_dist,            &
     2959                                  min_irrf_value, mrt_include_sw, mrt_nlevels,  &
     2960                                  mrt_skip_roof, net_radiation, nrefsteps,      &
     2961                                  plant_lw_interact, rad_angular_discretization,&
     2962                                  radiation_interactions_on, radiation_scheme,  &
     2963                                  raytrace_discrete_azims,                      &
     2964                                  raytrace_discrete_elevs, raytrace_mpi_rma,    &
     2965                                  skip_time_do_radiation, surface_reflections,  &
     2966                                  svfnorm_report_thresh, sw_radiation,          &
     2967                                  unscheduled_radiation_calls
    29592968   
    29602969       line = ' '
    … …  
    46484657     REAL(wp), DIMENSION(0:nsurf_type) :: costheta           !< direct irradiance factor of solar angle
    46494658     REAL(wp), DIMENSION(nzub:nzut)    :: pchf_prep          !< precalculated factor for canopy temperature tendency
    4650      REAL(wp), DIMENSION(nzub:nzut)    :: pctf_prep          !< precalculated factor for canopy transpiration tendency
    46514659     REAL(wp), PARAMETER               :: alpha = 0._wp      !< grid rotation (TODO: add to namelist or remove)
    46524660     REAL(wp)                          :: pc_box_area, pc_abs_frac, pc_abs_eff
    4653      REAL(wp), DIMENSION(0:nsurf_type) :: facearea
     4661     REAL(wp)                          :: asrc               !< area of source face
     4662     REAL(wp)                          :: pcrad              !< irradiance from plant canopy
    46544663     REAL(wp)                          :: pabsswl  = 0.0_wp  !< total absorbed SW radiation energy in local processor (W)
    46554664     REAL(wp)                          :: pabssw   = 0.0_wp  !< total absorbed SW radiation energy in all processors (W)
    … …  
    46724681         pchf_prep(:) = r_d * exner(nzub:nzut)                                 &
    46734682                     / (c_p * hyp(nzub:nzut) * dx*dy*dz(1)) !< equals to 1 / (rho * c_p * Vbox * T)
    4674          pctf_prep(:) = r_d * exner(nzub:nzut)                                 &
    4675                      / (l_v * hyp(nzub:nzut) * dx*dy*dz(1))
    46764683     ENDIF
    46774684#endif
    … …  
    47294736        mrtinlw(:) = 0._wp
    47304737     ENDIF
     4738     surfinlg(:)  = 0._wp !global
    47314739
    47324740
    … …  
    48244832!
    48254833!--        For surface-to-surface factors we calculate thermal radiation in 1st pass
    4826            surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * surfoutl(isurfsrc)
     4834           IF ( plant_lw_interact )  THEN
     4835              surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * svf(2,isvf) * surfoutl(isurfsrc)
     4836           ELSE
     4837              surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * surfoutl(isurfsrc)
     4838           ENDIF
    48274839        ENDDO
    48284840     ENDIF
    … …  
    48334845        i = surfl(ix, isurf)
    48344846        surfinswdif(isurf) = rad_sw_in_diff(j,i) * skyvft(isurf)
    4835         surfinlwdif(isurf) = rad_lw_in_diff(j,i) * skyvf(isurf)
     4847        IF ( plant_lw_interact )  THEN
     4848           surfinlwdif(isurf) = rad_lw_in_diff(j,i) * skyvft(isurf)
     4849        ELSE
     4850           surfinlwdif(isurf) = rad_lw_in_diff(j,i) * skyvf(isurf)
     4851        ENDIF
    48364852     ENDDO
    48374853!
    … …  
    48804896         pcbinswdir(:) = 0._wp
    48814897         pcbinswdif(:) = 0._wp
    4882          pcbinlw(:) = 0._wp  !< will stay always 0 since we don't absorb lw anymore
    4883 !
    4884 !--         pcsf first pass
     4898         pcbinlw(:) = 0._wp
     4899!
     4900!--      pcsf first pass
    48854901         DO icsf = 1, ncsfl
    48864902             ipcgb = csfsurf(1, icsf)
    … …  
    48914907
    48924908             IF ( isurfsrc == -1 )  THEN
    4893 !--                 Diffuse rad from sky.
    4894                  pcbinswdif(ipcgb) = csf(1,icsf) * rad_sw_in_diff(j,i)
    4895 
    4896                  !--Direct rad
    4897                  IF ( zenith(0) > 0 )  THEN
    4898                     !--Estimate directed box absorption
    4899                     pc_abs_frac = 1._wp - exp(pc_abs_eff * lad_s(k,j,i))
    4900 
    4901                     !--isd has already been established, see 1)
    4902                     pcbinswdir(ipcgb) = rad_sw_in_dir(j, i) * pc_box_area &
    4903                                         * pc_abs_frac * dsitransc(ipcgb, isd)
    4904                  ENDIF
     4909!
     4910!--             Diffuse rad from sky.
     4911                pcbinswdif(ipcgb) = csf(1,icsf) * rad_sw_in_diff(j,i)
     4912!
     4913!--             Absorbed diffuse LW from sky minus emitted to sky
     4914                IF ( plant_lw_interact )  THEN
     4915                   pcbinlw(ipcgb) = csf(1,icsf)                                  &
     4916                                       * (rad_lw_in_diff(j, i)                   &
     4917                                          - sigma_sb * (pt(k,j,i)*exner(k))**4)
     4918                ENDIF
     4919!
     4920!--             Direct rad
     4921                IF ( zenith(0) > 0 )  THEN
     4922!--                Estimate directed box absorption
     4923                   pc_abs_frac = 1._wp - exp(pc_abs_eff * lad_s(k,j,i))
     4924!
     4925!--                isd has already been established, see 1)
     4926                   pcbinswdir(ipcgb) = rad_sw_in_dir(j, i) * pc_box_area &
     4927                                       * pc_abs_frac * dsitransc(ipcgb, isd)
     4928                ENDIF
     4929             ELSE
     4930                IF ( plant_lw_interact )  THEN
     4931!
     4932!--                Thermal emission from plan canopy towards respective face
     4933                   pcrad = sigma_sb * (pt(k,j,i) * exner(k))**4 * csf(1,icsf)
     4934                   surfinlg(isurfsrc) = surfinlg(isurfsrc) + pcrad
     4935!
     4936!--                Remove the flux above + absorb LW from first pass from surfaces
     4937                   asrc = facearea(surf(id, isurfsrc))
     4938                   pcbinlw(ipcgb) = pcbinlw(ipcgb)                      &
     4939                                    + (csf(1,icsf) * surfoutl(isurfsrc) & ! Absorb from first pass surf emit
     4940                                       - pcrad)                         & ! Remove emitted heatflux
     4941                                    * asrc
     4942                ENDIF
    49054943             ENDIF
    49064944         ENDDO
    … …  
    49084946         pcbinsw(:) = pcbinswdir(:) + pcbinswdif(:)
    49094947     ENDIF
     4948
     4949     IF ( plant_lw_interact )  THEN
     4950!
     4951!--     Exchange incoming lw radiation from plant canopy
     4952#if defined( __parallel )
     4953        CALL MPI_Allreduce(MPI_IN_PLACE, surfinlg, nsurf, MPI_REAL, MPI_SUM, comm2d, ierr)
     4954        IF ( ierr /= 0 )  THEN
     4955           WRITE (9,*) 'Error MPI_Allreduce:', ierr
     4956           FLUSH(9)
     4957        ENDIF
     4958        surfinl(:) = surfinl(:) + surfinlg(surfstart(myid)+1:surfstart(myid+1))
     4959#else
     4960        surfinl(:) = surfinl(:) + surfinlg(:)
     4961#endif
     4962     ENDIF
     4963
    49104964     surfins = surfinswdir + surfinswdif
    49114965     surfinl = surfinl + surfinlwdif
    … …  
    49665020             isurfsrc = svfsurf(2, isvf)
    49675021             surfins(isurf) = surfins(isurf) + svf(1,isvf) * svf(2,isvf) * surfouts(isurfsrc)
    4968              surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * surfoutl(isurfsrc)
     5022             IF ( plant_lw_interact )  THEN
     5023                surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * svf(2,isvf) * surfoutl(isurfsrc)
     5024             ELSE
     5025                surfinl(isurf) = surfinl(isurf) + svf(1,isvf) * surfoutl(isurfsrc)
     5026             ENDIF
    49695027         ENDDO
    4970 
    4971 !--         radiation absorbed by plant canopy
    4972          DO icsf = 1, ncsfl
     5028!
     5029!--      NOTE: PC absorbtion and MRT from reflected can both be done at once
     5030!--      after all reflections if we do one more MPI_ALLGATHERV on surfout.
     5031!--      Advantage: less local computation. Disadvantage: one more collective
     5032!--      MPI call.
     5033!
     5034!--      Radiation absorbed by plant canopy
     5035         DO  icsf = 1, ncsfl
    49735036             ipcgb = csfsurf(1, icsf)
    49745037             isurfsrc = csfsurf(2, icsf)
    49755038             IF ( isurfsrc == -1 )  CYCLE ! sky->face only in 1st pass, not here
    4976 
    4977              pcbinsw(ipcgb) = pcbinsw(ipcgb) + csf(1,icsf) * surfouts(isurfsrc)
     5039!
     5040!--          Calculate source surface area. If the `surf' array is removed
     5041!--          before timestepping starts (future version), then asrc must be
     5042!--          stored within `csf'
     5043             asrc = facearea(surf(id, isurfsrc))
     5044             pcbinsw(ipcgb) = pcbinsw(ipcgb) + csf(1,icsf) * surfouts(isurfsrc) * asrc
     5045             IF ( plant_lw_interact )  THEN
     5046                pcbinlw(ipcgb) = pcbinlw(ipcgb) + csf(1,icsf) * surfoutl(isurfsrc) * asrc
     5047             ENDIF
    49785048         ENDDO
    49795049!
    … …  
    49965066     IF ( npcbl > 0 )  THEN
    49975067         pc_heating_rate(:,:,:) = 0.0_wp
    4998          pc_transpiration_rate(:,:,:) = 0.0_wp
    49995068         DO ipcgb = 1, npcbl
    5000                  
    50015069             j = pcbl(iy, ipcgb)
    50025070             i = pcbl(ix, ipcgb)
    50035071             k = pcbl(iz, ipcgb)
    50045072!
    5005 !--             Following expression equals former kk = k - nzb_s_inner(j,i)
     5073!--          Following expression equals former kk = k - nzb_s_inner(j,i)
    50065074             kk = k - get_topography_top_index_ji( j, i, 's' )  !- lad arrays are defined flat
    50075075             pc_heating_rate(kk, j, i) = (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
    50085076                 * pchf_prep(k) * pt(k, j, i) !-- = dT/dt
    5009 
    5010 !             pc_transpiration_rate(kk,j,i) = 0.75_wp* (pcbinsw(ipcgb) + pcbinlw(ipcgb)) &
    5011 !                 * pctf_prep(k) * pt(k, j, i) !-- = dq/dt
    5012 
    50135077         ENDDO
     5078
     5079         IF ( humidity .AND. plant_canopy_transpiration ) THEN
     5080!--          Calculation of plant canopy transpiration rate and correspondidng latent heat rate
     5081             pc_transpiration_rate(:,:,:) = 0.0_wp
     5082             pc_latent_rate(:,:,:) = 0.0_wp
     5083             DO ipcgb = 1, npcbl
     5084                 i = pcbl(ix, ipcgb)
     5085                 j = pcbl(iy, ipcgb)
     5086                 k = pcbl(iz, ipcgb)
     5087                 kk = k - get_topography_top_index_ji( j, i, 's' )  !- lad arrays are defined flat
     5088                 CALL pcm_calc_transpiration_rate( i, j, k, kk, pcbinsw(ipcgb), pcbinlw(ipcgb), &
     5089                                                   pc_transpiration_rate(kk,j,i), pc_latent_rate(kk,j,i) )
     5090              ENDDO
     5091         ENDIF
    50145092     ENDIF
    50155093!
    … …  
    51615239!--  domain when using average_radiation in the respective radiation model
    51625240
    5163 !--  Precalculate face areas for all face directions using normal vector
    5164      DO d = 0, nsurf_type
    5165         facearea(d) = 1._wp
    5166         IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx
    5167         IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy
    5168         IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1)
    5169      ENDDO
    51705241!--  calculate horizontal area
    51715242! !!! ATTENTION!!! uniform grid is assumed here
    … …  
    54295500       IMPLICIT NONE
    54305501
    5431        INTEGER(iwp) :: i, j, k, l, m
     5502       INTEGER(iwp) :: i, j, k, l, m, d
    54325503       INTEGER(iwp) :: k_topo     !< vertical index indicating topography top for given (j,i)
    54335504       INTEGER(iwp) :: nzptl, nzubl, nzutl, isurf, ipcgb, imrt
    … …  
    54395510#endif
    54405511
     5512!
     5513!--     precalculate face areas for different face directions using normal vector
     5514        DO d = 0, nsurf_type
     5515            facearea(d) = 1._wp
     5516            IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx
     5517            IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy
     5518            IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1)
     5519        ENDDO
    54415520!
    54425521!--    Find nzub, nzut, nzu via wall_flag_0 array (nzb_s_inner will be
    … …  
    59045983       ALLOCATE( surfouts(nsurf) )
    59055984       ALLOCATE( surfoutl(nsurf) )
     5985       ALLOCATE( surfinlg(nsurf) )
    59065986       ALLOCATE( skyvf(nsurfl) )
    59075987       ALLOCATE( skyvft(nsurfl) )
    … …  
    59406020        REAL(wp)                                      :: az1, az2      !< relative azimuth of section borders
    59416021        REAL(wp)                                      :: azmid         !< ray (center) azimuth
     6022        REAL(wp)                                      :: yxlen         !< |yxdir|
    59426023        REAL(wp), DIMENSION(2)                        :: yxdir         !< y,x *unit* vector of ray direction (in grid units)
    59436024        REAL(wp), DIMENSION(:), ALLOCATABLE           :: zdirs         !< directions in z (tangent of elevation)
     6025        REAL(wp), DIMENSION(:), ALLOCATABLE           :: zcent         !< zenith angle centers
    59446026        REAL(wp), DIMENSION(:), ALLOCATABLE           :: zbdry         !< zenith angle boundaries
    59456027        REAL(wp), DIMENSION(:), ALLOCATABLE           :: vffrac        !< view factor fractions for individual rays
    … …  
    59496031        INTEGER(iwp), DIMENSION(:), ALLOCATABLE       :: itarget       !< face indices of detected obstacles
    59506032        INTEGER(iwp)                                  :: itarg0, itarg1
    5951 #if defined( __parallel )
    5952 #endif
    5953 
    5954 
    5955 
    5956         REAL(wp),     DIMENSION(0:nsurf_type)         :: facearea
     6033
    59576034        INTEGER(iwp)                                  :: udim
    59586035        INTEGER(iwp), DIMENSION(:), ALLOCATABLE,TARGET:: nzterrl_l
    … …  
    59666043        LOGICAL                                       :: visible
    59676044        REAL(wp), DIMENSION(3)                        :: sa, ta          !< real coordinates z,y,x of source and target
     6045        REAL(wp)                                      :: difvf           !< differential view factor
    59686046        REAL(wp)                                      :: transparency, rirrf, sqdist, svfsum
    59696047        INTEGER(iwp)                                  :: isurflt, isurfs, isurflt_prev
    … …  
    59836061        CALL location_message( '    calculation of SVF and CSF', .TRUE. )
    59846062        CALL radiation_write_debug_log('Start calculation of SVF and CSF')
    5985 
    5986 !--     precalculate face areas for different face directions using normal vector
    5987         DO d = 0, nsurf_type
    5988             facearea(d) = 1._wp
    5989             IF ( idir(d) == 0 ) facearea(d) = facearea(d) * dx
    5990             IF ( jdir(d) == 0 ) facearea(d) = facearea(d) * dy
    5991             IF ( kdir(d) == 0 ) facearea(d) = facearea(d) * dz(1)
    5992         ENDDO
    59936063
    59946064!--     initialize variables and temporary arrays for calculation of svf and csf
    … …  
    62196289            END SELECT
    62206290
    6221             ALLOCATE ( zdirs(1:nzn), zbdry(0:nzn), vffrac(1:nzn*naz), &
     6291            ALLOCATE ( zdirs(1:nzn), zcent(1:nzn), zbdry(0:nzn), vffrac(1:nzn*naz), &
    62226292                       ztransp(1:nzn*naz), itarget(1:nzn*naz) )   !FIXME allocate itarget only
    62236293                                                                  !in case of rad_angular_discretization
    … …  
    62256295            itarg0 = 1
    62266296            itarg1 = nzn
    6227             zdirs(:) = (/( TAN(pi/2 - (zn0+(REAL(izn,wp)-.5_wp)*zns)), izn=1, nzn )/)
     6297            zcent(:) = (/( zn0+(REAL(izn,wp)-.5_wp)*zns, izn=1, nzn )/)
    62286298            zbdry(:) = (/( zn0+REAL(izn,wp)*zns, izn=0, nzn )/)
    62296299            IF ( td == iup_u  .OR.  td == iup_l )  THEN
    … …  
    62516321!--               sum of whole vffrac equals 1, verified
    62526322               ENDIF
    6253                yxdir = (/ COS(azmid), SIN(azmid) /)
     6323               yxdir(:) = (/ COS(azmid) / dy, SIN(azmid) / dx /)
     6324               yxlen = SQRT(SUM(yxdir(:)**2))
     6325               zdirs(:) = COS(zcent(:)) / (dz(1) * yxlen * SIN(zcent(:)))
     6326               yxdir(:) = yxdir(:) / yxlen
     6327
    62546328               CALL raytrace_2d(ta, yxdir, nzn, zdirs,                        &
    62556329                                    surfstart(myid) + isurflt, facearea(td),  &
    … …  
    62576331                                    .FALSE., lowest_free_ray,                 &
    62586332                                    ztransp(itarg0:itarg1),                   &
    6259                                     itarget(itarg0:itarg1))   !FIXME unit vect in grid units + zdirs
    6260                                                               !FIXME itarget available only in
    6261                                                               !case of rad_angular_discretization
     6333                                    itarget(itarg0:itarg1))
     6334
    62626335               skyvf(isurflt) = skyvf(isurflt) + &
    62636336                                SUM(vffrac(itarg0:itarg0+lowest_free_ray-1))
    … …  
    63376410            ENDIF ! rad_angular_discretization
    63386411
    6339             DEALLOCATE ( zdirs, zbdry, vffrac, ztransp, itarget ) !FIXME itarget shall be allocated only
     6412            DEALLOCATE ( zdirs, zcent, zbdry, vffrac, ztransp, itarget ) !FIXME itarget shall be allocated only
    63406413                                                                  !in case of rad_angular_discretization
    63416414!
    … …  
    63676440                  ENDIF
    63686441                 
     6442                  difvf = dot_product((/ kdir(sd), jdir(sd), idir(sd) /), uv) & ! cosine of source normal and direction
     6443                      * dot_product((/ kdir(td), jdir(td), idir(td) /), -uv) &  ! cosine of target normal and reverse direction
     6444                      / (pi * sqdist) ! square of distance between centers
     6445!
    63696446!--               irradiance factor (our unshaded shape view factor) = view factor per differential target area * source area
    6370                   rirrf = dot_product((/ kdir(sd), jdir(sd), idir(sd) /), uv) & ! cosine of source normal and direction
    6371                       * dot_product((/ kdir(td), jdir(td), idir(td) /), -uv) &  ! cosine of target normal and reverse direction
    6372                       / (pi * sqdist) & ! square of distance between centers
    6373                       * facearea(sd)
     6447                  rirrf = difvf * facearea(sd)
    63746448
    63756449!--               reject raytracing for potentially too small view factor values
    … …  
    63806454
    63816455!--               raytrace + process plant canopy sinks within
    6382                   CALL raytrace(sa, ta, isurfs, rirrf, facearea(td), .TRUE., &
     6456                  CALL raytrace(sa, ta, isurfs, difvf, facearea(td), .TRUE., &
    63836457                                visible, transparency)
    63846458
    … …  
    64286502        nzn = raytrace_discrete_elevs / 2
    64296503        zns = pi / 2._wp / REAL(nzn, wp)
    6430         ALLOCATE ( zdirs(1:nzn), vffrac(1:nzn), ztransp(1:nzn), &
     6504        ALLOCATE ( zdirs(1:nzn), zcent(1:nzn), vffrac(1:nzn), ztransp(1:nzn), &
    64316505               itarget(1:nzn) )
    6432         zdirs(:) = (/( TAN(pi/2 - (zn0+(REAL(izn,wp)-.5_wp)*zns)), izn=1, nzn )/)
     6506        zcent(:) = (/( zn0+(REAL(izn,wp)-.5_wp)*zns, izn=1, nzn )/)
    64336507        vffrac(:) = 0._wp
    64346508
    … …  
    64406514           DO  iaz = 1, naz
    64416515              azmid = az0 + (REAL(iaz, wp) - .5_wp) * azs
    6442               yxdir = (/ COS(azmid), SIN(azmid) /)
     6516              yxdir(:) = (/ COS(azmid) / dy, SIN(azmid) / dx /)
     6517              yxlen = SQRT(SUM(yxdir(:)**2))
     6518              zdirs(:) = COS(zcent(:)) / (dz(1) * yxlen * SIN(zcent(:)))
     6519              yxdir(:) = yxdir(:) / yxlen
    64436520              CALL raytrace_2d(ta, yxdir, nzn, zdirs,                                &
    64446521                                   -999, -999._wp, vffrac, .FALSE., .FALSE., .TRUE., &
    6445                                    lowest_free_ray, ztransp, itarget) !FIXME unit vect in grid units + zdirs
     6522                                   lowest_free_ray, ztransp, itarget)
    64466523
    64476524!--           Save direct solar transparency
    … …  
    64566533           ENDDO
    64576534        ENDDO
    6458         DEALLOCATE ( zdirs, vffrac, ztransp, itarget )
     6535        DEALLOCATE ( zdirs, zcent, vffrac, ztransp, itarget )
    64596536!--
    64606537!--     Raytrace to MRT boxes
    … …  
    64696546           nzn = raytrace_discrete_elevs
    64706547           zns = pi / REAL(nzn, wp)
    6471            ALLOCATE ( zdirs(1:nzn), zbdry(0:nzn), vffrac(1:nzn*naz), vffrac0(1:nzn), &
     6548           ALLOCATE ( zdirs(1:nzn), zcent(1:nzn), zbdry(0:nzn), vffrac(1:nzn*naz), vffrac0(1:nzn), &
    64726549                      ztransp(1:nzn*naz), itarget(1:nzn*naz) )   !FIXME allocate itarget only
    64736550                                                                 !in case of rad_angular_discretization
    64746551
    6475            zdirs(:) = (/( TAN(pi/2 - (zn0+(REAL(izn,wp)-.5_wp)*zns)), izn=1, nzn )/)
     6552           zcent(:) = (/( zn0+(REAL(izn,wp)-.5_wp)*zns, izn=1, nzn )/)
    64766553           zbdry(:) = (/( zn0+REAL(izn,wp)*zns, izn=0, nzn )/)
    64776554           vffrac0(:) = (COS(zbdry(0:nzn-1)) - COS(zbdry(1:nzn))) / 2._wp / REAL(naz, wp)
    … …  
    64936570              DO  iaz = 1, naz
    64946571                 azmid = az0 + (REAL(iaz, wp) - .5_wp) * azs
    6495                  CALL raytrace_2d(ta, (/ COS(azmid), SIN(azmid) /), nzn, zdirs,  &
     6572                 yxdir(:) = (/ COS(azmid) / dy, SIN(azmid) / dx /)
     6573                 yxlen = SQRT(SUM(yxdir(:)**2))
     6574                 zdirs(:) = COS(zcent(:)) / (dz(1) * yxlen * SIN(zcent(:)))
     6575                 yxdir(:) = yxdir(:) / yxlen
     6576
     6577                 CALL raytrace_2d(ta, yxdir, nzn, zdirs,                         &
    64966578                                  -999, -999._wp, vffrac(itarg0:itarg1), .TRUE., &
    64976579                                  .FALSE., .TRUE., lowest_free_ray,              &
    64986580                                  ztransp(itarg0:itarg1),                        &
    6499                                   itarget(itarg0:itarg1))   !FIXME unit vect in grid units + zdirs
    6500                                                             !FIXME itarget available only in
    6501                                                             !case of rad_angular_discretization
     6581                                  itarget(itarg0:itarg1))
    65026582
    65036583!--              Sky view factors for MRT
    … …  
    65746654
    65756655           ENDDO ! imrt
    6576            DEALLOCATE ( zdirs, zbdry, vffrac, vffrac0, ztransp, itarget )
     6656           DEALLOCATE ( zdirs, zcent, zbdry, vffrac, vffrac0, ztransp, itarget )
    65776657!
    65786658!--        Move MRT factors to final arrays
    … …  
    67786858                    j = 0
    67796859                ENDIF
    6780 !--             fill out real values of rsvf, rtransp
    6781                 csflt(1,kcsf) = acsf(kcsf)%rsvf
     6860                csflt(1,kcsf) = acsf(kcsf)%rcvf
    67826861!--             fill out integer values of itz,ity,itx,isurfs
    67836862                kcsflt(1,kcsf) = acsf(kcsf)%itz
    … …  
    69467025!>    doesn't need to be the same in all three dimensions).
    69477026!------------------------------------------------------------------------------!
    6948     SUBROUTINE raytrace(src, targ, isrc, rirrf, atarg, create_csf, visible, transparency)
     7027    SUBROUTINE raytrace(src, targ, isrc, difvf, atarg, create_csf, visible, transparency)
    69497028        IMPLICIT NONE
    69507029
    69517030        REAL(wp), DIMENSION(3), INTENT(in)     :: src, targ    !< real coordinates z,y,x
    69527031        INTEGER(iwp), INTENT(in)               :: isrc         !< index of source face for csf
    6953         REAL(wp), INTENT(in)                   :: rirrf        !< irradiance factor for csf
     7032        REAL(wp), INTENT(in)                   :: difvf        !< differential view factor for csf
    69547033        REAL(wp), INTENT(in)                   :: atarg        !< target surface area for csf
    69557034        LOGICAL, INTENT(in)                    :: create_csf   !< whether to generate new CSFs during raytracing
    … …  
    71137192                    acsf(ncsfl)%itz = boxes(1,i)
    71147193                    acsf(ncsfl)%isurfs = isrc
    7115                     acsf(ncsfl)%rsvf = cursink*transparency*rirrf*atarg
     7194                    acsf(ncsfl)%rcvf = cursink*transparency*difvf*atarg
    71167195                ENDIF  !< create_csf
    71177196
    … …  
    75147593                           acsf(ncsfl)%itz = iz
    75157594                           acsf(ncsfl)%isurfs = iorig
    7516                            acsf(ncsfl)%rsvf = (1._wp - curtrans)*transparency(k)*aorig*vffrac(k)
     7595                           acsf(ncsfl)%rcvf = (1._wp - curtrans)*transparency(k)*vffrac(k)
    75177596                        ENDIF
    75187597
    … …  
    75717650                        acsf(ncsfl)%ity = rt2_track(1,i)
    75727651                        acsf(ncsfl)%itz = iz
    7573                         acsf(ncsfl)%isurfs = itarget(k) !if above horizon, then itarget(k)==-1, which
    7574                                                         !is also a special ID indicating sky
    7575                         acsf(ncsfl)%rsvf = (1._wp - curtrans)*transparency(k)*aorig*vffrac(k)
     7652                        IF ( itarget(k) /= -1 )  ERROR STOP !FIXME remove after test
     7653                        acsf(ncsfl)%isurfs = -1
     7654                        acsf(ncsfl)%rcvf = (1._wp - curtrans)*transparency(k)*aorig*vffrac(k)
    75767655                     ENDIF  ! create_csf
    75777656
    … …  
    76007679         INTEGER(iwp), TARGET, INTENT(out)   ::  isurfl
    76017680         INTEGER(iwp), INTENT(out)           ::  iproc
     7681#if defined( __parallel )
     7682#else
     7683         INTEGER(iwp)                        ::  target_displ  !< index of the grid in the local gridsurf array
     7684#endif
    76027685         INTEGER(iwp)                        ::  px, py        !< number of processors in x and y direction
    76037686                                                               !< before the processor in the question
    … …  
    81388221                         .AND.  acsfnew(iwrite)%isurfs == acsf(iread)%isurfs )  THEN
    81398222
    8140                     acsfnew(iwrite)%rsvf = acsfnew(iwrite)%rsvf + acsf(iread)%rsvf
     8223                    acsfnew(iwrite)%rcvf = acsfnew(iwrite)%rcvf + acsf(iread)%rcvf
    81418224!--                 advance reading index, keep writing index
    81428225                ELSE

  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3435 r3449  
    2828! -----------------
    2929! $Id$
     30! Bugfix: Fix average arrays allocations in usm_average_3d_data (J.Resler)
     31! Bugfix: Fix reading wall temperatures (J.Resler)
     32! Bugfix: Fix treating of outputs for wall temperature and sky view factors (J.Resler)
     33!
     34!
     35! 3435 2018-10-26 18:25:44Z gronemeier
    3036! Bugfix: allocate gamma_w_green_sat until nzt_wall+1
    3137!
    … …  
    411417   
    412418    USE plant_canopy_model_mod,                                                &
    413         ONLY:  pc_heating_rate, pc_transpiration_rate
     419        ONLY:  pc_heating_rate, pc_transpiration_rate, pc_latent_rate
    414420   
    415421    USE radiation_model_mod,                                                   &
    … …  
    18451851                CASE ( 'usm_rad_net' )
    18461852!--                 array of complete radiation balance
    1847                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%rad_net_av) )  THEN
    1848                         ALLOCATE( surf_usm_h%rad_net_av(1:surf_usm_h%ns) )
    1849                         surf_usm_h%rad_net_av = 0.0_wp
     1853                    IF ( l == -1 ) THEN
     1854                        IF ( .NOT.  ALLOCATED(surf_usm_h%rad_net_av) )  THEN
     1855                            ALLOCATE( surf_usm_h%rad_net_av(1:surf_usm_h%ns) )
     1856                            surf_usm_h%rad_net_av = 0.0_wp
     1857                        ENDIF
    18501858                    ELSE
    18511859                        IF ( .NOT.  ALLOCATED(surf_usm_v(l)%rad_net_av) )  THEN
    … …  
    19681976                CASE ( 'usm_rad_hf' )
    19691977!--                 array of heat flux from radiation for surfaces after i-th reflection
    1970                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%surfhf_av) )  THEN
    1971                         ALLOCATE( surf_usm_h%surfhf_av(1:surf_usm_h%ns) )
    1972                         surf_usm_h%surfhf_av = 0.0_wp
     1978                    IF ( l == -1 ) THEN
     1979                        IF ( .NOT.  ALLOCATED(surf_usm_h%surfhf_av) )  THEN
     1980                           ALLOCATE( surf_usm_h%surfhf_av(1:surf_usm_h%ns) )
     1981                           surf_usm_h%surfhf_av = 0.0_wp
     1982                        ENDIF
    19731983                    ELSE
    19741984                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%surfhf_av) )  THEN
    … …  
    19811991!--                 array of sensible heat flux from surfaces
    19821992!--                 land surfaces
    1983                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%wshf_eb_av) )  THEN
    1984                         ALLOCATE( surf_usm_h%wshf_eb_av(1:surf_usm_h%ns) )
    1985                         surf_usm_h%wshf_eb_av = 0.0_wp
     1993                    IF ( l == -1 ) THEN
     1994                       IF ( .NOT.  ALLOCATED(surf_usm_h%wshf_eb_av) )  THEN
     1995                          ALLOCATE( surf_usm_h%wshf_eb_av(1:surf_usm_h%ns) )
     1996                          surf_usm_h%wshf_eb_av = 0.0_wp
     1997                       ENDIF
    19861998                    ELSE
    19871999                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%wshf_eb_av) )  THEN
    … …  
    20362048                CASE ( 'usm_wghf' )
    20372049!--                 array of heat flux from ground (wall, roof, land)
    2038                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%wghf_eb_av) )  THEN
    2039                         ALLOCATE( surf_usm_h%wghf_eb_av(1:surf_usm_h%ns) )
    2040                         surf_usm_h%wghf_eb_av = 0.0_wp
     2050                    IF ( l == -1 ) THEN
     2051                       IF ( .NOT.  ALLOCATED(surf_usm_h%wghf_eb_av) )  THEN
     2052                           ALLOCATE( surf_usm_h%wghf_eb_av(1:surf_usm_h%ns) )
     2053                           surf_usm_h%wghf_eb_av = 0.0_wp
     2054                       ENDIF
    20412055                    ELSE
    20422056                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%wghf_eb_av) )  THEN
    … …  
    20482062                CASE ( 'usm_wghf_window' )
    20492063!--                 array of heat flux from window ground (wall, roof, land)
    2050                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%wghf_eb_window_av) )  THEN
    2051                         ALLOCATE( surf_usm_h%wghf_eb_window_av(1:surf_usm_h%ns) )
    2052                         surf_usm_h%wghf_eb_window_av = 0.0_wp
     2064                    IF ( l == -1 ) THEN
     2065                       IF ( .NOT.  ALLOCATED(surf_usm_h%wghf_eb_window_av) )  THEN
     2066                           ALLOCATE( surf_usm_h%wghf_eb_window_av(1:surf_usm_h%ns) )
     2067                           surf_usm_h%wghf_eb_window_av = 0.0_wp
     2068                       ENDIF
    20532069                    ELSE
    20542070                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%wghf_eb_window_av) )  THEN
    … …  
    20602076                CASE ( 'usm_wghf_green' )
    20612077!--                 array of heat flux from green ground (wall, roof, land)
    2062                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%wghf_eb_green_av) )  THEN
    2063                         ALLOCATE( surf_usm_h%wghf_eb_green_av(1:surf_usm_h%ns) )
    2064                         surf_usm_h%wghf_eb_green_av = 0.0_wp
     2078                    IF ( l == -1 ) THEN
     2079                       IF ( .NOT.  ALLOCATED(surf_usm_h%wghf_eb_green_av) )  THEN
     2080                           ALLOCATE( surf_usm_h%wghf_eb_green_av(1:surf_usm_h%ns) )
     2081                           surf_usm_h%wghf_eb_green_av = 0.0_wp
     2082                       ENDIF
    20652083                    ELSE
    20662084                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%wghf_eb_green_av) )  THEN
    … …  
    20722090                CASE ( 'usm_iwghf' )
    20732091!--                 array of heat flux from indoor ground (wall, roof, land)
    2074                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%iwghf_eb_av) )  THEN
    2075                         ALLOCATE( surf_usm_h%iwghf_eb_av(1:surf_usm_h%ns) )
    2076                         surf_usm_h%iwghf_eb_av = 0.0_wp
     2092                    IF ( l == -1 ) THEN
     2093                       IF ( .NOT.  ALLOCATED(surf_usm_h%iwghf_eb_av) )  THEN
     2094                           ALLOCATE( surf_usm_h%iwghf_eb_av(1:surf_usm_h%ns) )
     2095                           surf_usm_h%iwghf_eb_av = 0.0_wp
     2096                       ENDIF
    20772097                    ELSE
    20782098                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%iwghf_eb_av) )  THEN
    … …  
    20842104                CASE ( 'usm_iwghf_window' )
    20852105!--                 array of heat flux from indoor window ground (wall, roof, land)
    2086                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%iwghf_eb_window_av) )  THEN
    2087                         ALLOCATE( surf_usm_h%iwghf_eb_window_av(1:surf_usm_h%ns) )
    2088                         surf_usm_h%iwghf_eb_window_av = 0.0_wp
     2106                    IF ( l == -1 ) THEN
     2107                       IF ( .NOT.  ALLOCATED(surf_usm_h%iwghf_eb_window_av) )  THEN
     2108                           ALLOCATE( surf_usm_h%iwghf_eb_window_av(1:surf_usm_h%ns) )
     2109                           surf_usm_h%iwghf_eb_window_av = 0.0_wp
     2110                       ENDIF
    20892111                    ELSE
    20902112                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%iwghf_eb_window_av) )  THEN
    … …  
    20932115                       ENDIF
    20942116                    ENDIF
    2095                    
     2117
    20962118                CASE ( 'usm_t_surf_wall' )
    20972119!--                 surface temperature for surfaces
    2098                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_surf_wall_av) )  THEN
    2099                         ALLOCATE( surf_usm_h%t_surf_wall_av(1:surf_usm_h%ns) )
    2100                         surf_usm_h%t_surf_wall_av = 0.0_wp
     2120                    IF ( l == -1 ) THEN
     2121                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_surf_wall_av) )  THEN
     2122                           ALLOCATE( surf_usm_h%t_surf_wall_av(1:surf_usm_h%ns) )
     2123                           surf_usm_h%t_surf_wall_av = 0.0_wp
     2124                       ENDIF
    21012125                    ELSE
    21022126                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_surf_wall_av) )  THEN
    … …  
    21082132                CASE ( 'usm_t_surf_window' )
    21092133!--                 surface temperature for window surfaces
    2110                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_surf_window_av) )  THEN
    2111                         ALLOCATE( surf_usm_h%t_surf_window_av(1:surf_usm_h%ns) )
    2112                         surf_usm_h%t_surf_window_av = 0.0_wp
     2134                    IF ( l == -1 ) THEN
     2135                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_surf_window_av) )  THEN
     2136                           ALLOCATE( surf_usm_h%t_surf_window_av(1:surf_usm_h%ns) )
     2137                           surf_usm_h%t_surf_window_av = 0.0_wp
     2138                       ENDIF
    21132139                    ELSE
    21142140                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_surf_window_av) )  THEN
    … …  
    21202146                CASE ( 'usm_t_surf_green' )
    21212147!--                 surface temperature for green surfaces
    2122                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_surf_green_av) )  THEN
    2123                         ALLOCATE( surf_usm_h%t_surf_green_av(1:surf_usm_h%ns) )
    2124                         surf_usm_h%t_surf_green_av = 0.0_wp
     2148                    IF ( l == -1 ) THEN
     2149                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_surf_green_av) )  THEN
     2150                           ALLOCATE( surf_usm_h%t_surf_green_av(1:surf_usm_h%ns) )
     2151                           surf_usm_h%t_surf_green_av = 0.0_wp
     2152                       ENDIF
    21252153                    ELSE
    21262154                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_surf_green_av) )  THEN
    … …  
    21322160                CASE ( 'usm_t_surf_10cm' )
    21332161!--                 near surface temperature for whole surfaces
    2134                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_surf_10cm_av) )  THEN
    2135                         ALLOCATE( surf_usm_h%t_surf_10cm_av(1:surf_usm_h%ns) )
    2136                         surf_usm_h%t_surf_10cm_av = 0.0_wp
     2162                    IF ( l == -1 ) THEN
     2163                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_surf_10cm_av) )  THEN
     2164                           ALLOCATE( surf_usm_h%t_surf_10cm_av(1:surf_usm_h%ns) )
     2165                           surf_usm_h%t_surf_10cm_av = 0.0_wp
     2166                       ENDIF
    21372167                    ELSE
    21382168                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_surf_10cm_av) )  THEN
    … …  
    21442174                CASE ( 'usm_t_wall' )
    21452175!--                 wall temperature for iwl layer of walls and land
    2146                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_wall_av) )  THEN
    2147                         ALLOCATE( surf_usm_h%t_wall_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
    2148                         surf_usm_h%t_wall_av = 0.0_wp
     2176                    IF ( l == -1 ) THEN
     2177                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_wall_av) )  THEN
     2178                           ALLOCATE( surf_usm_h%t_wall_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
     2179                           surf_usm_h%t_wall_av = 0.0_wp
     2180                       ENDIF
    21492181                    ELSE
    21502182                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_wall_av) )  THEN
    … …  
    21562188                CASE ( 'usm_t_window' )
    21572189!--                 window temperature for iwl layer of walls and land
    2158                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_window_av) )  THEN
    2159                         ALLOCATE( surf_usm_h%t_window_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
    2160                         surf_usm_h%t_window_av = 0.0_wp
     2190                    IF ( l == -1 ) THEN
     2191                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_window_av) )  THEN
     2192                           ALLOCATE( surf_usm_h%t_window_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
     2193                           surf_usm_h%t_window_av = 0.0_wp
     2194                       ENDIF
    21612195                    ELSE
    21622196                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_window_av) )  THEN
    … …  
    21682202                CASE ( 'usm_t_green' )
    21692203!--                 green temperature for iwl layer of walls and land
    2170                     IF ( l == -1 .AND. .NOT.  ALLOCATED(surf_usm_h%t_green_av) )  THEN
    2171                         ALLOCATE( surf_usm_h%t_green_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
    2172                         surf_usm_h%t_green_av = 0.0_wp
     2204                    IF ( l == -1 ) THEN
     2205                       IF ( .NOT.  ALLOCATED(surf_usm_h%t_green_av) )  THEN
     2206                           ALLOCATE( surf_usm_h%t_green_av(nzb_wall:nzt_wall,1:surf_usm_h%ns) )
     2207                           surf_usm_h%t_green_av = 0.0_wp
     2208                       ENDIF
    21732209                    ELSE
    21742210                       IF ( .NOT.  ALLOCATED(surf_usm_v(l)%t_green_av) )  THEN
    … …  
    30913127        CHARACTER(LEN=*),INTENT(OUT)   ::  unit       !<
    30923128
    3093         INTEGER(iwp)                                  :: i,j          !< index
     3129        INTEGER(iwp)                                  :: i,j,l        !< index
     3130        CHARACTER(LEN=2)                              :: ls
    30943131        CHARACTER(LEN=varnamelength)                  :: var          !< TRIM(variable)
    3095         INTEGER(iwp), PARAMETER                       :: nl1 = 32     !< number of directional usm variables
     3132        INTEGER(iwp), PARAMETER                       :: nl1 = 31     !< number of directional usm variables
    30963133        CHARACTER(LEN=varnamelength), DIMENSION(nl1)  :: varlist1 = & !< list of directional usm variables
    30973134                  (/'usm_rad_net                   ', &
    … …  
    31193156                    'usm_surfalb                   ', &
    31203157                    'usm_surfemis                  ', &
    3121                     'usm_t_surf                    ', &
     3158                    'usm_t_surf_wall               ', &
    31223159                    'usm_t_surf_window             ', &
    31233160                    'usm_t_surf_green              ', &
    31243161                    'usm_t_green                   ', &
    31253162                    'usm_t_surf_10cm               ', &
    3126                     'usm_t_wall                    ', &
    3127                     'usm_t_window                  ', &
    3128                     'usm_t_green                   '/)
    3129 
    3130         INTEGER(iwp), PARAMETER                       :: nl2 = 9      !< number of other variables
     3163                    'usm_skyvf                     ', &
     3164                    'usm_skyvft                    '/)
     3165
     3166        INTEGER(iwp), PARAMETER                       :: nl2 = 7      !< number of other variables
    31313167        CHARACTER(LEN=varnamelength), DIMENSION(nl2)  :: varlist2 = & !< list of other usm variables
    3132                   (/'usm_skyvf                     ', &
    3133                     'usm_skyvft                    ', &
    3134                     'usm_svf                       ', &
     3168                  (/'usm_svf                       ', &
    31353169                    'usm_dif                       ', &
    31363170                    'usm_rad_pc_inlw               ', &
    … …  
    31393173                    'usm_rad_pc_inswdif            ', &
    31403174                    'usm_rad_pc_inswref            '/)
     3175
     3176        INTEGER(iwp), PARAMETER                       :: nl3 = 3      !< number of directional layer usm variables
     3177        CHARACTER(LEN=varnamelength), DIMENSION(nl3)  :: varlist3 = & !< list of directional layer usm variables
     3178                  (/'usm_t_wall                    ', &
     3179                    'usm_t_window                  ', &
     3180                    'usm_t_green                   '/)
    31413181
    31423182        INTEGER(iwp), PARAMETER                       :: nd = 5     !< number of directions
    … …  
    31583198                 EXIT
    31593199              ENDIF
     3200              IF ( lfound ) EXIT
     3201           ENDDO
     3202        ENDDO
     3203        IF ( lfound ) GOTO 10
     3204!       directional layer variables
     3205        DO i = 1, nl3
     3206           DO j = 1, nd
     3207              DO l = nzb_wall, nzt_wall
     3208                 WRITE(ls,'(A1,I1)') '_',l
     3209                 IF ( TRIM(var) == TRIM(varlist3(i))//TRIM(ls)//TRIM(dirname(j)) ) THEN
     3210                    lfound = .TRUE.
     3211                    EXIT
     3212                 ENDIF
     3213              ENDDO
    31603214              IF ( lfound ) EXIT
    31613215           ENDDO
    … …  
    45034557             var(1:9) == 'usm_qsws_'  .OR.  var(1:13) == 'usm_qsws_veg_'  .OR.              &
    45044558             var(1:13) == 'usm_qsws_liq_' .OR.                                              &
    4505              var(1:10) == 'usm_t_surf_wall'  .OR.  var(1:10) == 'usm_t_wall'  .OR.          &
     4559             var(1:15) == 'usm_t_surf_wall'  .OR.  var(1:10) == 'usm_t_wall'  .OR.          &
    45064560             var(1:17) == 'usm_t_surf_window'  .OR.  var(1:12) == 'usm_t_window'  .OR.      &
    45074561             var(1:16) == 'usm_t_surf_green'  .OR. var(1:11) == 'usm_t_green' .OR.          &
    … …  
    59105964        ENDIF
    59115965
    5912         IF ( plant_canopy )  THEN
    5913            
    5914             IF ( .NOT.  ALLOCATED( pc_heating_rate) )  THEN
    5915 !--             then pc_heating_rate is allocated in init_plant_canopy
    5916 !--             in case of cthf /= 0 => we need to allocate it for our use here
    5917                 ALLOCATE( pc_heating_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    5918 
    5919                 pc_heating_rate = 0.0_wp
    5920 
    5921             ENDIF
    5922 
    5923             IF ( .NOT.  ALLOCATED( pc_transpiration_rate) )  THEN
    5924 !--             then pc_heating_rate is allocated in init_plant_canopy
    5925 !--             in case of cthf /= 0 => we need to allocate it for our use here
    5926                 ALLOCATE( pc_transpiration_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    5927 
    5928                 pc_transpiration_rate = 0.0_wp
    5929 
    5930 
    5931             ENDIF
    5932         ENDIF
    59335966!
    59345967!--    Check for consistent initialization.
    … …  
    60926125               ENDDO
    60936126            ENDDO
    6094         ELSE
    6095 !--         If specified, replace constant wall temperatures with fully 3D values from file
    6096             IF ( read_wall_temp_3d )  CALL usm_read_wall_temperature()
    6097 !
    60986127        ENDIF
    6099        
    6100 !--   
     6128
     6129!--     If specified, replace constant wall temperatures with fully 3D values from file
     6130        IF ( read_wall_temp_3d )  CALL usm_read_wall_temperature()
     6131
     6132!--
    61016133!--     Possibly DO user-defined actions (e.g. define heterogeneous wall surface)
    61026134        CALL user_init_urban_surface
    … …  
    88708902!> x, y, d, z coordinates
    88718903!------------------------------------------------------------------------------!
    8872     PURE FUNCTION advance_surface(isurfl_start, isurfl_stop, x, y, z, d) &
    8873             result(isurfl)
    8874 
    8875         INTEGER(iwp), INTENT(in)                :: isurfl_start, isurfl_stop
     8904    PURE FUNCTION find_surface( x, y, z, d ) result(isurfl)
     8905
    88768906        INTEGER(iwp), INTENT(in)                :: x, y, z, d
    8877         INTEGER(iwp)                            :: isx, isy, isz, isd
    88788907        INTEGER(iwp)                            :: isurfl
    8879 
    8880         DO isurfl = isurfl_start, isurfl_stop
    8881             isx = surfl(ix, isurfl)
    8882             isy = surfl(iy, isurfl)
    8883             isz = surfl(iz, isurfl)
    8884             isd = surfl(id, isurfl)
    8885             IF ( isx==x .and. isy==y .and. isz==z .and. isd==d )  RETURN
    8886         ENDDO
     8908        INTEGER(iwp)                            :: isx, isy, isz
     8909
     8910        IF ( d == 0 ) THEN
     8911           DO  isurfl = 1, surf_usm_h%ns
     8912              isx = surf_usm_h%i(isurfl)
     8913              isy = surf_usm_h%j(isurfl)
     8914              isz = surf_usm_h%k(isurfl)
     8915              IF ( isx==x .and. isy==y .and. isz==z )  RETURN
     8916           ENDDO
     8917        ELSE
     8918           DO  isurfl = 1, surf_usm_v(d-1)%ns
     8919              isx = surf_usm_v(d-1)%i(isurfl)
     8920              isy = surf_usm_v(d-1)%j(isurfl)
     8921              isz = surf_usm_v(d-1)%k(isurfl)
     8922              IF ( isx==x .and. isy==y .and. isz==z )  RETURN
     8923           ENDDO
     8924        ENDIF
    88878925
    88888926!--     coordinate not found
    … …  
    89278965                        nys <= j .and. j <= nyn)  THEN  !< local processor
    89288966!--                     identify surface id
    8929                         isurfl = advance_surface(isurfl+1, nsurfl, i, j, k, d)
     8967                        isurfl = find_surface( i, j, k, d )
    89308968                        IF ( isurfl == -1 )  THEN
    89318969                            WRITE(message_string, '(a,4i5,a,i5,a)') 'Coordinates (xyzd) ', i, j, k, d, &
    … …  
    89398977                           t_surf_wall_h(isurfl) = rtsurf
    89408978                           t_wall_h(:,isurfl) = rtwall(:)
     8979                           t_window_h(:,isurfl) = rtwall(:)
     8980                           t_green_h(:,isurfl) = rtwall(:)
    89418981                        ELSE
    89428982                           t_surf_wall_v(d-1)%t(isurfl) = rtsurf
    89438983                           t_wall_v(d-1)%t(:,isurfl) = rtwall(:)
     8984                           t_window_v(d-1)%t(:,isurfl) = rtwall(:)
     8985                           t_green_v(d-1)%t(:,isurfl) = rtwall(:)
    89448986                        ENDIF
    89458987                    ENDIF
    … …  
    96639705                                          surf_usm_v(l)%rad_lw_out(m)
    96649706
    9665 if ((abs(t_surf_wall_v(l)%t(m)-276.).gt.100.).or.(abs(t_surf_window_v(l)%t(m)-276.).gt.300.) &
    9666 .or.(abs(t_surf_green_v(l)%t(m)-276.).gt.100)) then
    9667 print*, "tsurfvvvv",m,t_surf_wall_v(l)%t(m),t_surf_window_v(l)%t(m),t_surf_green_v(l)%t(m)
    9668 print*, "params", surf_usm_v(l)%emissivity(ind_veg_wall,m),lambda_surface, t_wall_v(l)%t(nzb_wall:nzt_wall,m), &
    9669                   surf_usm_v(l)%emissivity(ind_wat_win,m),lambda_surface_window,t_window_v(l)%t(nzb_wall:nzt_wall,m),t_green_v(l)%t(nzb_wall:nzt_wall,m)
    9670 print*, "dicken",surf_usm_v(l)%zw(nzb_wall:nzt_wall,m),surf_usm_v(l)%zw_window(nzb_wall:nzt_wall,m),surf_usm_v(l)%zw_green(nzb_wall:nzt_wall,m)
    9671 print*, "c",surf_usm_v(l)%c_surface_window(m),surf_usm_v(l)%c_surface(m)
    9672 !if ((abs(t_surf_v(l)%t(m)-276.).gt.10.).or.(abs(t_surf_window_v(l)%t(m)-276.).gt.10.)) then
    9673 stop
    9674 endif
    96759707!--           numerator of the prognostic equation
    96769708              coef_1 = surf_usm_v(l)%rad_net_l(m) +                            & ! coef +1 corresponds to -lwout included in calculation of radnet_l

  • palm/trunk/TESTS/cases/urban_environment/MONITORING/urban_environment_rc

    r3417 r3449  
    11
    22 ******************************    --------------------------------------------
    3  * PALM 5.0  Rev: 3401M       *    atmosphere - run initialized with COSMO data
     3 * PALM 5.0  Rev: 3448M       *    atmosphere - run initialized with COSMO data
    44 ******************************    --------------------------------------------
    55
    6  Date:                 24-10-18    Run:       urban_environment                 
    7  Time:                 12:09:47    Run-No.:   00
    8  Run on host:        imuktrunkr
     6 Date:                 29-10-18    Run:       urban_environment                 
     7 Time:                 20:15:36    Run-No.:   00
     8 Run on host:        lcmuk_para
    99 Number of PEs:               4    Processor grid (x,y): (   2,   2) calculated
    1010 ------------------------------------------------------------------------------
    … …  
    203203       Profile: u, u*2, wu, w*u*, w"u", v, v*2, wv, w*v*, w"v", w, w*2,     
    204204
    205               : pt, pt*2, wpt, w*pt*, w"pt", q, q*2, e, e*, km, l, t_soil,   
    206 
    207               : m_soil, rad_lw_in, rad_lw_out,                               
     205              : theta, theta*2, wtheta, w*theta*, w"theta", q, q*2,         
     206
     207              : e, e*, km, l, t_soil, m_soil, rad_lw_in, rad_lw_out,         
     208
     209              :                                                             
    208210       Output every                60.00 s
    209211       Time averaged over          60.00 s
    … …  
    214216
    215217
    216        XY-cross-section  Arrays:  shf*, ghf*, u*, z0*, t*, ol*, m_liq*, r_s*, r_a*, rad_net*, rad_lw_in
     218       XY-cross-section  Arrays:  shf*, ghf*, us*, z0*, t*, ol*, m_liq*, r_s*, r_a*, rad_net*, rad_lw_i
    217219       Output every                60.00 s  and at the start
    218220       Cross sections at k = /0/
    … …  
    242244
    243245
    244        Arrays:  u, v, w, pt, q, e, p, m_soil, t_soil, pcm_lad, pcm_heatrate, kc_PM10,
     246       Arrays:  u, v, w, theta, q, e, p, m_soil, t_soil, pcm_lad, pcm_heatrate, kc_PM
    245247       Output every                60.00 s  and at the start
    246248       Upper output limit at      121.00 m  (GP   61)
    … …  
    254256
    255257
    256        Arrays:  u, v, w, pt, q, e, p, m_soil, t_soil, kc_PM10, kc_PM25,              
     258       Arrays:  u, v, w, theta, q, e, p, m_soil, t_soil, kc_PM10, kc_PM25,           
    257259       Output every                60.00 s  and at the start
    258260       Time averaged over          60.00 s
    … …  
    342344RUN  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
    343345---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    344   0      0 00:00:00.00   0.7580A  -2.3731D -1.1942D  0.4177  0.070  0.00  -0.360E+00     16.  0.137E+01  0.743E-02  0.234E-01  0.446E-02   60  17  18   19   0  14   20  18   0    0.000    0.000       0
     346  0      0 00:00:00.00   0.7580A  -2.3731D -1.1942D  0.4177  0.071  0.00  -0.361E+00     16.  0.137E+01  0.743E-02  0.234E-01  0.446E-02   60  17  18   19   0  14   20  18   0    0.000    0.000       0
    345347  0      1 00:00:00.75   0.7580A  -2.3732  -1.5382   0.7410  0.086  0.33  -0.108E+01      2.  0.137E+01  0.116E-01  0.269E-02  0.892E-03   60  18  15   16  15   0   20  18   0    0.000    0.000       0
    346348  0      2 00:00:01.51   0.7580A  -2.3733  -1.6029   0.8137  0.085  0.00  -0.108E+01      2.  0.137E+01  0.125E-01  0.953E-03  0.351E-03   60  17  14   16  15   0   20  17   0    0.000    0.000       0
    … …  
    350352  0      6 00:00:04.54   0.7580A  -2.3736  -1.4555   0.8403  0.082  0.24  -0.112E+01     18.  0.137E+01  0.137E-01  0.769E-03  0.202E-03   60  17  11   17  14   0   20  17   0    0.000    0.000       0
    351353  0      7 00:00:05.30   0.7580A  -2.3736  -1.4258   0.8339  0.082  0.25  -0.112E+01     18.  0.137E+01  0.140E-01  0.767E-03  0.200E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
    352   0      8 00:00:06.06   0.7580A  -2.3737  -1.3955   0.8272  0.081  0.26  -0.112E+01     18.  0.137E+01  0.143E-01  0.768E-03  0.199E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
    353   0      9 00:00:06.82   0.7580A  -2.3738  -1.3671   0.8203  0.081  0.26  -0.111E+01     18.  0.137E+01  0.147E-01  0.769E-03  0.198E-03   60  16   9   11  14   0   20  17   0    0.000    0.000       0
    354   0     10 00:00:07.58   0.7580A  -2.3739  -1.3438   0.8138  0.080  0.27  -0.111E+01     18.  0.137E+01  0.150E-01  0.771E-03  0.198E-03   60  15   8   11  14   0   20  17   0    0.000    0.000       0
    355   0     11 00:00:08.33   0.7580A  -2.3739  -1.3267   0.8083  0.080  0.28  -0.111E+01     18.  0.137E+01  0.154E-01  0.775E-03  0.198E-03   60  15   7   16  14   0   20  17   0    0.000    0.000       0
     354  0      8 00:00:06.06   0.7580A  -2.3737  -1.3954   0.8272  0.081  0.26  -0.112E+01     18.  0.137E+01  0.143E-01  0.768E-03  0.199E-03   60  16  10   11  14   0   20  17   0    0.000    0.000       0
     355  0      9 00:00:06.82   0.7580A  -2.3738  -1.3671   0.8203  0.081  0.26  -0.112E+01     18.  0.137E+01  0.147E-01  0.769E-03  0.198E-03   60  16   9   11  14   0   20  17   0    0.000    0.000       0
     356  0     10 00:00:07.58   0.7580A  -2.3739  -1.3438   0.8139  0.080  0.27  -0.111E+01     18.  0.137E+01  0.150E-01  0.771E-03  0.198E-03   60  15   8   11  14   0   20  17   0    0.000    0.000       0
     357  0     11 00:00:08.33   0.7580A  -2.3740  -1.3267   0.8083  0.080  0.28  -0.111E+01     18.  0.137E+01  0.154E-01  0.775E-03  0.198E-03   60  15   7   16  14   0   20  17   0    0.000    0.000       0
    356358  0     12 00:00:09.09   0.7580A  -2.3740  -1.3246   0.8042  0.079  0.28  -0.111E+01     18.  0.136E+01  0.158E-01  0.780E-03  0.198E-03   60  15   6   16  14   0   20  17   0    0.000    0.000       0
    357   0     13 00:00:09.85   0.7580A  -2.3741  -1.3195   0.8018  0.079  0.29  -0.110E+01     18.  0.136E+01  0.162E-01  0.786E-03  0.197E-03   60  14   6   16  14   0   20  17   0    0.000    0.000       0
     359  0     13 00:00:09.85   0.7580A  -2.3741  -1.3195   0.8018  0.079  0.29  -0.111E+01     18.  0.136E+01  0.162E-01  0.786E-03  0.197E-03   60  14   6   16  14   0   20  17   0    0.000    0.000       0
    358360  0     14 00:00:10.61   0.7580A  -2.3742  -1.3110   0.8363  0.078  0.29  -0.110E+01     18.  0.136E+01  0.166E-01  0.791E-03  0.197E-03   60  14   5   16  14   0   11  14  18    0.000    0.000       0
    359   0     15 00:00:11.37   0.7580A  -2.3743  -1.2989   0.8854  0.078  0.37  -0.110E+01     28.  0.136E+01  0.170E-01  0.796E-03  0.197E-03   60  13   4   16  14   0   11  14  18    0.000    0.000       0
     361  0     15 00:00:11.37   0.7580A  -2.3743  -1.2989   0.8855  0.078  0.37  -0.110E+01     28.  0.136E+01  0.170E-01  0.796E-03  0.197E-03   60  13   4   16  14   0   11  14  18    0.000    0.000       0
    360362  0     16 00:00:12.12   0.7580A  -2.3743  -1.2833   0.9238  0.078  0.38  -0.110E+01     28.  0.136E+01  0.175E-01  0.799E-03  0.196E-03   60  13   3   16  14   0   11  14  18    0.000    0.000       0
    361   0     17 00:00:12.88   0.7580A  -2.3744  -1.2642   0.9549  0.077  0.44  -0.109E+01     40.  0.136E+01  0.180E-01  0.802E-03  0.195E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
     363  0     17 00:00:12.88   0.7580A  -2.3744  -1.2642   0.9550  0.077  0.45  -0.110E+01     40.  0.136E+01  0.180E-01  0.802E-03  0.195E-03   60  12   2   16  14   0   11  14  18    0.000    0.000       0
    362364  0     18 00:00:13.64   0.7580A  -2.3745  -1.2461   0.9857  0.077  0.46  -0.109E+01     40.  0.136E+01  0.185E-01  0.807E-03  0.194E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
    363   0     19 00:00:14.40   0.7580A  -2.3746  -1.2422   1.0287  0.077  0.47  -0.109E+01     40.  0.136E+01  0.190E-01  0.811E-03  0.192E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
     365  0     19 00:00:14.40   0.7580A  -2.3746  -1.2422   1.0288  0.077  0.47  -0.109E+01     40.  0.136E+01  0.190E-01  0.812E-03  0.193E-03   60  12   1   11  14   0   10  14  18    0.000    0.000       0
    364366  0     20 00:00:15.16   0.7580A  -2.3746  -1.2391   1.0686  0.077  0.48  -0.109E+01     40.  0.136E+01  0.195E-01  0.816E-03  0.191E-03   60  12   0   11  14   0   10  14  18    0.000    0.000       0
    365   0     21 00:00:15.91   0.7580A  -2.3747  -1.2366   1.1033  0.076  0.49  -0.108E+01     40.  0.136E+01  0.200E-01  0.820E-03  0.190E-03   60  11  19   11  14   0   10  14  18    0.000    0.000       0
    366   0     22 00:00:16.67   0.7580A  -2.3748  -1.2343   1.1366  0.076  0.51  -0.108E+01     40.  0.136E+01  0.206E-01  0.824E-03  0.188E-03   60  11  18   11  14   0   12  14  17    0.000    0.000       0
    367   0     23 00:00:17.43   0.7580A  -2.3749  -1.2323   1.1786  0.076  0.52  -0.108E+01     40.  0.136E+01  0.212E-01  0.829E-03  0.187E-03   60  10  17   11  14   0   12  14  17    0.000    0.000       0
    368   0     24 00:00:18.19   0.7580A  -2.3749  -1.2303   1.2093  0.076  0.53  -0.108E+01     40.  0.136E+01  0.217E-01  0.834E-03  0.186E-03   60  10  16   11  14   0   11  14  18    0.000    0.000       0
    369   0     25 00:00:18.95   0.7580A  -2.3750  -1.2286   1.2416  0.076  0.54  -0.107E+01     40.  0.136E+01  0.224E-01  0.841E-03  0.185E-03   60  10  16   11  14   0   11  14  18    0.000    0.000       0
    370   0     26 00:00:19.70   0.7580A  -2.3751  -1.2272   1.2663  0.075  0.55  -0.107E+01     40.  0.136E+01  0.230E-01  0.846E-03  0.184E-03   60  10  15   11  14   0   11  14  18    0.000    0.000       0
    371   0     27 00:00:20.46   0.7580A  -2.3752  -1.2260   1.2924  0.075  0.56  -0.107E+01     40.  0.136E+01  0.236E-01  0.854E-03  0.183E-03   60   9  14   11  14   0   13  14  17    0.000    0.000       0
    372   0     28 00:00:21.22   0.7580A  -2.3752  -1.2247   1.3070  0.075  0.56  -0.107E+01     40.  0.136E+01  0.243E-01  0.861E-03  0.183E-03   60   9  13   11  14   0   13  14  17    0.000    0.000       0
    373   0     29 00:00:21.98   0.7580A  -2.3753  -1.2230   1.3219  0.075  0.56  -0.107E+01     40.  0.136E+01  0.250E-01  0.871E-03  0.183E-03   60   9  13   11  14   0   14  14  17    0.000    0.000       0
    374   0     30 00:00:22.74   0.7580A  -2.3754  -1.2207   1.3989  0.075  0.57  -0.107E+01     40.  0.136E+01  0.257E-01  0.881E-03  0.183E-03   60   8  11   11  14   0    6  14  15    0.000    0.000       0
    375   0     31 00:00:23.49   0.7580A  -2.3755  -1.2173   1.4707  0.074  0.57  -0.107E+01     40.  0.136E+01  0.264E-01  0.892E-03  0.183E-03   60   8  11   11  14   0    6  14  15    0.000    0.000       0
    376   0     32 00:00:24.25   0.7580A  -2.3756  -1.2125   1.5578  0.074  0.57  -0.107E+01     40.  0.136E+01  0.271E-01  0.903E-03  0.184E-03   60   8  10   11  14   0    7  14  15    0.000    0.000       0
    377   0     33 00:00:25.01   0.7580A  -2.3756  -1.2059   1.6304  0.074  0.58  -0.107E+01     40.  0.136E+01  0.279E-01  0.915E-03  0.185E-03   60   7   9   11  14   0    7  14  15    0.000    0.000       0
    378   0     34 00:00:25.77   0.7580A  -2.3757  -1.1973   1.7210  0.074  0.58  -0.107E+01     40.  0.136E+01  0.287E-01  0.925E-03  0.186E-03   60   7   8   11  14   0    8  14  15    0.000    0.000       0
    379   0     35 00:00:26.53   0.7580A  -2.3758  -1.1861   1.7884  0.074  0.58  -0.107E+01     40.  0.136E+01  0.295E-01  0.937E-03  0.187E-03   60   7   7   11  14   0    9  14  15    0.000    0.000       0
    380   0     36 00:00:27.28   0.7580A  -2.3759  -1.1721   1.8841  0.074  0.61  -0.107E+01     44.  0.136E+01  0.303E-01  0.950E-03  0.189E-03   60   6   6   11  14   0    9  14  15    0.000    0.000       0
    381   0     37 00:00:28.04   0.7580A  -2.3759  -1.1675   1.9632  0.073  0.61  -0.107E+01     44.  0.136E+01  0.312E-01  0.965E-03  0.190E-03   60   6   5   12  14   0   10  14  15    0.000    0.000       0
    382   0     38 00:00:28.80   0.7580A  -2.3760  -1.1630   2.0273  0.073  0.62  -0.106E+01     46.  0.136E+01  0.321E-01  0.981E-03  0.192E-03   60   6   5   12  14   0   10  14  15    0.000    0.000       0
    383   0     39 00:00:29.56   0.7580A  -2.3761  -1.1572   2.1068  0.073  0.62  -0.106E+01     46.  0.136E+01  0.330E-01  0.997E-03  0.194E-03   60   6   4   12  14   0   11  14  15    0.000    0.000       0
    384   0     40 00:00:30.32   0.7580A  -2.3762  -1.1498   2.1299  0.073  0.62  -0.106E+01     46.  0.136E+01  0.339E-01  0.102E-02  0.195E-03   60   5   3   12  14   0   11  14  15    0.000    0.000       0
    385   0     41 00:00:31.07   0.7580A  -2.3762  -1.1406   2.1864  0.073  0.62  -0.106E+01     46.  0.136E+01  0.349E-01  0.104E-02  0.196E-03   60   5   2   12  14   0   12  14  15    0.000    0.000       0
    386   0     42 00:00:31.83   0.7570A  -2.3763  -1.1298   2.1802  0.073  0.62  -0.106E+01     46.  0.136E+01  0.359E-01  0.106E-02  0.197E-03   60   4   1   12  14   0   12  14  15    0.000    0.000       0
    387   0     43 00:00:32.59   0.7570A  -2.3764  -1.1178   2.1937  0.073  0.63  -0.106E+01     48.  0.136E+01  0.369E-01  0.108E-02  0.199E-03   60   4   0   12  14   0   13  14  15    0.000    0.000       0
    388   0     44 00:00:33.35   0.7570A  -2.3765  -1.1054   2.1759  0.073  0.63  -0.106E+01     48.  0.136E+01  0.379E-01  0.110E-02  0.200E-03   60   4  19   14  14   0   13  14  15    0.000    0.000       0
    389   0     45 00:00:34.10   0.7570A  -2.3765  -1.1159   2.1763  0.073  0.62  -0.106E+01     48.  0.137E+01  0.390E-01  0.113E-02  0.202E-03   60   4  19   14  14   0   16  14  16    0.000    0.000       0
    390   0     46 00:00:34.86   0.7570A  -2.3766  -1.1262   2.1870  0.073  0.62  -0.106E+01     48.  0.137E+01  0.400E-01  0.115E-02  0.203E-03   60   3  17   14  14   0   16  14  16    0.000    0.000       0
    391   0     47 00:00:35.62   0.7570A  -2.3767  -1.1351   2.1866  0.072  0.62  -0.106E+01     48.  0.137E+01  0.411E-01  0.117E-02  0.203E-03   60   3  17   14  14   0   17  14  16    0.000    0.000       0
    392   0     48 00:00:36.37   0.7570A  -2.3768  -1.1417   2.2127  0.072  0.62  -0.106E+01     50.  0.137E+01  0.422E-01  0.120E-02  0.203E-03   60   3  16   14  14   0   20  14  17    0.000    0.000       0
    393   0     49 00:00:37.13   0.7570A  -2.3768  -1.1457   2.2094  0.072  0.62  -0.106E+01     50.  0.137E+01  0.433E-01  0.122E-02  0.203E-03   60   2  15   14  14   0   20  14  17    0.000    0.000       0
    394   0     50 00:00:37.89   0.7570A  -2.3769  -1.1472   2.1873  0.072  0.62  -0.106E+01     50.  0.137E+01  0.444E-01  0.125E-02  0.202E-03   60   2  14   14  14   0   18  14  16    0.000    0.000       0
    395   0     51 00:00:38.64   0.7570A  -2.3770  -1.1466   2.1653  0.072  0.62  -0.106E+01     50.  0.137E+01  0.455E-01  0.127E-02  0.202E-03   60   2  13   14  14   0   18  14  16    0.000    0.000       0
    396   0     52 00:00:39.40   0.7570A  -2.3771  -1.1441   2.1348  0.072  0.63  -0.106E+01     52.  0.137E+01  0.466E-01  0.129E-02  0.202E-03   60   1  12   14  14   0   18  14  16    0.000    0.000       0
    397   0     53 00:00:40.16   0.7570A  -2.3771  -1.1614   2.1122  0.072  0.63  -0.106E+01     52.  0.137E+01  0.476E-01  0.132E-02  0.202E-03   60   1  11   22  14  13   19  14  16    0.000    0.000       0
    398   0     54 00:00:40.92   0.7570A  -2.3772  -1.2365   2.1081  0.072  0.63  -0.106E+01     52.  0.137E+01  0.487E-01  0.134E-02  0.203E-03   60   1  11   22  14  13   20  14  16    0.000    0.000       0
    399   0     55 00:00:41.67   0.7570A  -2.3773  -1.2991   2.0894  0.072  0.63  -0.106E+01     52.  0.137E+01  0.497E-01  0.136E-02  0.203E-03   60   1  10   22  14  13   20  14  16    0.000    0.000       0
    400   0     56 00:00:42.43   0.7570A  -2.3774  -1.3363   2.0876  0.072  0.64  -0.107E+01     52.  0.137E+01  0.507E-01  0.137E-02  0.203E-03   60   0   9   22  14  13   22  13  16    0.000    0.000       0
    401   0     57 00:00:43.19   0.7570A  -2.3774  -1.3594   2.0990  0.072  0.65  -0.107E+01     54.  0.137E+01  0.517E-01  0.139E-02  0.204E-03   60   0   8    8  14   5   23  13  16    0.000    0.000       0
    402   0     58 00:00:43.94   0.7570A  -2.3775  -1.3961   2.0887  0.072  0.65  -0.108E+01     54.  0.137E+01  0.527E-01  0.141E-02  0.204E-03   60   0   7    8  14   5   23  13  16    0.000    0.000       0
    403   0     59 00:00:44.70   0.7570A  -2.3776  -1.4255   2.0185  0.072  0.66  -0.108E+01     54.  0.137E+01  0.537E-01  0.142E-02  0.205E-03   60  19   6    8  14   5   23  13  16    0.000    0.000       0
    404   0     60 00:00:45.46   0.7570A  -2.3777  -1.4514   1.8960  0.072  0.66  -0.108E+01     54.  0.137E+01  0.546E-01  0.144E-02  0.206E-03   60  19   6    8  14   5   21  14  16    0.000    0.000       0
    405   0     61 00:00:46.21   0.7570A  -2.3777  -1.4753   1.8564  0.072  0.68  -0.108E+01     56.  0.137E+01  0.556E-01  0.146E-02  0.207E-03   60  19   5    8  14   5   21  14  16    0.000    0.000       0
    406   0     62 00:00:46.97   0.7570A  -2.3778  -1.4968   1.8575  0.072  0.68  -0.108E+01     56.  0.137E+01  0.565E-01  0.148E-02  0.208E-03   60  19   4    8  14   5   19  14  17    0.000    0.000       0
    407   0     63 00:00:47.73   0.7570A  -2.3779  -1.5143   1.8652  0.072  0.69  -0.108E+01     56.  0.137E+01  0.574E-01  0.151E-02  0.208E-03   60  18   3    8  14   5   19  14  17    0.000    0.000       0
    408   0     64 00:00:48.49   0.7570A  -2.3780  -1.5266   1.8703  0.072  0.69  -0.107E+01     56.  0.137E+01  0.583E-01  0.153E-02  0.209E-03   60  18   2    8  14   5   19  14  17    0.000    0.000       0
    409   0     65 00:00:49.24   0.7570A  -2.3780  -1.5329   1.8775  0.072  0.70  -0.107E+01     56.  0.137E+01  0.592E-01  0.155E-02  0.209E-03   60  18   2    8  14   5   19  14  17    0.000    0.000       0
    410   0     66 00:00:50.00   0.7570A  -2.3781  -1.5330   1.8872  0.072  0.71  -0.106E+01     58.  0.137E+01  0.601E-01  0.158E-02  0.210E-03   60  17   0    8  14   5   19  14  17    0.000    0.000       0
    411   0     67 00:00:50.76   0.7570A  -2.3782  -1.5265   1.8967  0.072  0.72  -0.106E+01     58.  0.137E+01  0.609E-01  0.160E-02  0.210E-03   60  17   0    8  14   5   19  14  17    0.000    0.000       0
    412   0     68 00:00:51.51   0.7570A  -2.3783  -1.5126   1.9036  0.072  0.72  -0.106E+01     58.  0.137E+01  0.618E-01  0.162E-02  0.211E-03   60  17  19    8  14   5   19  14  17    0.000    0.000       0
    413   0     69 00:00:52.27   0.7570A  -2.3783  -1.5277   1.9068  0.072  0.73  -0.105E+01     58.  0.137E+01  0.626E-01  0.164E-02  0.211E-03   60  16  17   23  12  10   19  14  17    0.000    0.000       0
    414   0     70 00:00:53.03   0.7570A  -2.3784  -1.5317   1.9079  0.072  0.74  -0.105E+01     58.  0.137E+01  0.634E-01  0.167E-02  0.211E-03   60  16  17   23  12  10   19  14  17    0.000    0.000       0
    415   0     71 00:00:53.78   0.7570A  -2.3785  -1.5203   1.9104  0.072  0.75  -0.105E+01     60.  0.137E+01  0.643E-01  0.169E-02  0.212E-03   60  16  16   23  12  10   19  14  17    0.000    0.000       0
    416   0     72 00:00:54.54   0.7570A  -2.3786  -1.5277   1.9276  0.072  0.76  -0.105E+01     60.  0.137E+01  0.651E-01  0.171E-02  0.213E-03   60  15  14   23  11   9   17  14  17    0.000    0.000       0
    417   0     73 00:00:55.30   0.7570A  -2.3786  -1.5235   1.9656  0.072  0.76  -0.104E+01     60.  0.137E+01  0.659E-01  0.173E-02  0.213E-03   60  15  14   23  11   9   17  14  17    0.000    0.000       0
    418   0     74 00:00:56.06   0.7570A  -2.3787  -1.5014   1.9903  0.072  0.77  -0.104E+01     60.  0.137E+01  0.666E-01  0.175E-02  0.214E-03   60  15  13   23  12   9   18  14  17    0.000    0.000       0
    419   0     75 00:00:56.81   0.7570A  -2.3788  -1.4942   2.0221  0.072  0.77  -0.104E+01     60.  0.137E+01  0.674E-01  0.177E-02  0.215E-03   60  14  11   23  12   9   18  14  17    0.000    0.000       0
    420   0     76 00:00:57.57   0.7570A  -2.3789  -1.4971   2.0358  0.072  0.78  -0.103E+01     60.  0.137E+01  0.682E-01  0.180E-02  0.217E-03   60  14  11   23  11   8   18  14  17    0.000    0.000       0
    421   0     77 00:00:58.33   0.7570A  -2.3789  -1.4822   2.0535  0.072  0.78  -0.103E+01     60.  0.137E+01  0.689E-01  0.182E-02  0.218E-03   60  14  10   23  11   8   19  14  17    0.000    0.000       0
    422   0     78 00:00:59.08   0.7570A  -2.3790  -1.4489   2.0559  0.072  0.79  -0.103E+01     60.  0.137E+01  0.697E-01  0.185E-02  0.219E-03   60  13   9   23  10   7   19  14  17    0.000    0.000       0
    423   0     79 00:00:59.84   0.7570A  -2.3791  -1.4398   2.0309  0.072  0.79  -0.102E+01     60.  0.137E+01  0.704E-01  0.187E-02  0.219E-03   60  13   8   23  11   7   19  14  17    0.000    0.000       0
    424   0     80 00:01:00.60   0.7570A  -2.3792  -1.4269   1.9909  0.072  0.80  -0.102E+01     62.  0.137E+01  0.711E-01  0.189E-02  0.220E-03   60  13   7   23  10   6   20  14  17    0.000    0.000       0
    425   0     81 00:01:01.35   0.7570A  -2.3792  -1.4439   1.9515  0.072  0.81  -0.102E+01     62.  0.137E+01  0.718E-01  0.192E-02  0.221E-03   60  13   7   23  10   6   20  14  17    0.000    0.000       0
    426   0     82 00:01:02.11   0.7570A  -2.3793  -1.4366   1.9048  0.072  0.81  -0.102E+01     62.  0.137E+01  0.725E-01  0.194E-02  0.221E-03   60  13   6   23  10   6   22  13  16    0.000    0.000       0
    427   0     83 00:01:02.87   0.7560A  -2.3794  -1.4117   1.9630  0.072  0.82  -0.102E+01     62.  0.137E+01  0.732E-01  0.197E-02  0.222E-03   60  12   5   11  12   4   22  13  16    0.000    0.000       0
    428   0     84 00:01:03.62   0.7560A  -2.3795  -1.4064   1.9978  0.072  0.82  -0.102E+01     62.  0.137E+01  0.739E-01  0.199E-02  0.224E-03   60  12   4   11  12   4   23  13  16    0.000    0.000       0
    429   0     85 00:01:04.38   0.7560A  -2.3795  -1.4246   1.9913  0.072  0.82  -0.102E+01     62.  0.137E+01  0.745E-01  0.201E-02  0.225E-03   60  12   4    8  13   5   23  13  16    0.000    0.000       0
    430   0     86 00:01:05.14   0.7560A  -2.3796  -1.4442   1.9883  0.073  0.83  -0.102E+01     62.  0.137E+01  0.751E-01  0.202E-02  0.227E-03   60  11   2    8  13   5   16  14  14    0.000    0.000       0
    431   0     87 00:01:05.89   0.7560A  -2.3797  -1.4622   1.9979  0.073  0.83  -0.102E+01     62.  0.137E+01  0.757E-01  0.204E-02  0.229E-03   60  11   2    8  13   5   16  14  14    0.000    0.000       0
    432   0     88 00:01:06.65   0.7560A  -2.3798  -1.4784   2.0064  0.073  0.83  -0.102E+01     62.  0.137E+01  0.763E-01  0.206E-02  0.231E-03   60  11   1    8  13   5   17  14  14    0.000    0.000       0
    433   0     89 00:01:07.40   0.7560A  -2.3798  -1.4928   2.0170  0.073  0.83  -0.101E+01     62.  0.137E+01  0.768E-01  0.208E-02  0.233E-03   60  11   1    8  13   5   17  14  14    0.000    0.000       0
    434   0     90 00:01:08.16   0.7560A  -2.3799  -1.5053   2.0339  0.073  0.84  -0.101E+01     62.  0.137E+01  0.774E-01  0.210E-02  0.235E-03   60  10  19    8  13   5   18  14  14    0.000    0.000       0
    435   0     91 00:01:08.92   0.7560A  -2.3800  -1.5158   2.0510  0.073  0.84  -0.101E+01     62.  0.137E+01  0.779E-01  0.212E-02  0.237E-03   60  10  18    8  13   5   18  14  14    0.000    0.000       0
    436   0     92 00:01:09.67   0.7560A  -2.3800  -1.5243   1.9579  0.073  0.84  -0.101E+01     62.  0.137E+01  0.784E-01  0.214E-02  0.240E-03   60  10  18    8  13   5   18  14  14    0.000    0.000       0
    437   0     93 00:01:10.43   0.7560A  -2.3801  -1.5308   1.9738  0.073  0.85  -0.101E+01     62.  0.137E+01  0.789E-01  0.216E-02  0.243E-03   60   9  16    8  13   5   19  14  14    0.000    0.000       0
     367  0     21 00:00:15.91   0.7580A  -2.3747  -1.2365   1.1034  0.076  0.49  -0.109E+01     40.  0.136E+01  0.200E-01  0.820E-03  0.190E-03   60  11  19   11  14   0   10  14  18    0.000    0.000       0
     368  0     22 00:00:16.67   0.7580A  -2.3748  -1.2343   1.1367  0.076  0.51  -0.108E+01     40.  0.136E+01  0.206E-01  0.825E-03  0.188E-03   60  11  18   11  14   0   12  14  17    0.000    0.000       0
     369  0     23 00:00:17.43   0.7580A  -2.3749  -1.2322   1.1787  0.076  0.52  -0.108E+01     40.  0.136E+01  0.212E-01  0.830E-03  0.187E-03   60  10  17   11  14   0   12  14  17    0.000    0.000       0
     370  0     24 00:00:18.19   0.7580A  -2.3749  -1.2303   1.2094  0.076  0.53  -0.108E+01     40.  0.136E+01  0.217E-01  0.835E-03  0.186E-03   60  10  16   11  14   0   11  14  18    0.000    0.000       0
     371  0     25 00:00:18.95   0.7580A  -2.3750  -1.2285   1.2417  0.076  0.54  -0.108E+01     40.  0.136E+01  0.224E-01  0.841E-03  0.185E-03   60  10  16   11  14   0   11  14  18    0.000    0.000       0
     372  0     26 00:00:19.70   0.7580A  -2.3751  -1.2271   1.2664  0.075  0.55  -0.107E+01     40.  0.136E+01  0.230E-01  0.847E-03  0.184E-03   60  10  15   11  14   0   11  14  18    0.000    0.000       0
     373  0     27 00:00:20.46   0.7580A  -2.3752  -1.2259   1.2926  0.075  0.56  -0.107E+01     40.  0.136E+01  0.236E-01  0.854E-03  0.183E-03   60   9  14   11  14   0   13  14  17    0.000    0.000       0
     374  0     28 00:00:21.22   0.7580A  -2.3753  -1.2246   1.3071  0.075  0.56  -0.107E+01     40.  0.136E+01  0.243E-01  0.862E-03  0.183E-03   60   9  13   11  14   0   13  14  17    0.000    0.000       0
     375  0     29 00:00:21.98   0.7580A  -2.3753  -1.2229   1.3222  0.075  0.56  -0.107E+01     40.  0.136E+01  0.250E-01  0.871E-03  0.183E-03   60   9  13   11  14   0   14  14  17    0.000    0.000       0
     376  0     30 00:00:22.74   0.7580A  -2.3754  -1.2206   1.4032  0.075  0.57  -0.107E+01     40.  0.136E+01  0.257E-01  0.882E-03  0.183E-03   60   8  11   11  14   0    6  14  15    0.000    0.000       0
     377  0     31 00:00:23.49   0.7580A  -2.3755  -1.2172   1.4746  0.074  0.57  -0.107E+01     40.  0.136E+01  0.264E-01  0.892E-03  0.183E-03   60   8  11   11  14   0    6  14  15    0.000    0.000       0
     378  0     32 00:00:24.25   0.7580A  -2.3756  -1.2124   1.5624  0.074  0.57  -0.107E+01     40.  0.136E+01  0.271E-01  0.904E-03  0.184E-03   60   8  10   11  14   0    7  14  15    0.000    0.000       0
     379  0     33 00:00:25.01   0.7580A  -2.3756  -1.2058   1.6345  0.074  0.58  -0.107E+01     40.  0.136E+01  0.279E-01  0.915E-03  0.185E-03   60   7   9   11  14   0    7  14  15    0.000    0.000       0
     380  0     34 00:00:25.77   0.7580A  -2.3757  -1.1972   1.7258  0.074  0.58  -0.107E+01     40.  0.136E+01  0.287E-01  0.925E-03  0.186E-03   60   7   8   11  14   0    8  14  15    0.000    0.000       0
     381  0     35 00:00:26.53   0.7580A  -2.3758  -1.1860   1.7935  0.074  0.58  -0.107E+01     40.  0.136E+01  0.295E-01  0.937E-03  0.187E-03   60   7   7   11  14   0    9  14  15    0.000    0.000       0
     382  0     36 00:00:27.28   0.7580A  -2.3759  -1.1720   1.8888  0.074  0.61  -0.107E+01     44.  0.136E+01  0.303E-01  0.951E-03  0.189E-03   60   6   6   11  14   0    9  14  15    0.000    0.000       0
     383  0     37 00:00:28.04   0.7580A  -2.3759  -1.1675   1.9682  0.073  0.61  -0.107E+01     44.  0.136E+01  0.312E-01  0.965E-03  0.190E-03   60   6   5   12  14   0   10  14  15    0.000    0.000       0
     384  0     38 00:00:28.80   0.7580A  -2.3760  -1.1630   2.0314  0.073  0.62  -0.107E+01     46.  0.136E+01  0.321E-01  0.982E-03  0.192E-03   60   6   5   12  14   0   10  14  15    0.000    0.000       0
     385  0     39 00:00:29.56   0.7580A  -2.3761  -1.1572   2.1114  0.073  0.63  -0.107E+01     46.  0.136E+01  0.330E-01  0.998E-03  0.194E-03   60   6   4   12  14   0   11  14  15    0.000    0.000       0
     386  0     40 00:00:30.32   0.7580A  -2.3762  -1.1498   2.1334  0.073  0.62  -0.107E+01     46.  0.136E+01  0.339E-01  0.102E-02  0.195E-03   60   5   3   12  14   0   11  14  15    0.000    0.000       0
     387  0     41 00:00:31.07   0.7580A  -2.3762  -1.1406   2.1907  0.073  0.62  -0.107E+01     46.  0.136E+01  0.349E-01  0.104E-02  0.196E-03   60   5   2   12  14   0   12  14  15    0.000    0.000       0
     388  0     42 00:00:31.83   0.7570A  -2.3763  -1.1298   2.1834  0.073  0.62  -0.106E+01     46.  0.136E+01  0.359E-01  0.106E-02  0.198E-03   60   4   1   12  14   0   12  14  15    0.000    0.000       0
     389  0     43 00:00:32.59   0.7570A  -2.3764  -1.1178   2.1983  0.073  0.63  -0.106E+01     48.  0.136E+01  0.369E-01  0.108E-02  0.199E-03   60   4   0   12  14   0   13  14  15    0.000    0.000       0
     390  0     44 00:00:33.35   0.7570A  -2.3765  -1.1055   2.1793  0.073  0.63  -0.106E+01     48.  0.136E+01  0.379E-01  0.110E-02  0.200E-03   60   4  19   14  14   0   13  14  15    0.000    0.000       0
     391  0     45 00:00:34.10   0.7570A  -2.3765  -1.1160   2.1784  0.073  0.62  -0.106E+01     48.  0.137E+01  0.390E-01  0.113E-02  0.202E-03   60   4  19   14  14   0   16  14  16    0.000    0.000       0
     392  0     46 00:00:34.86   0.7570A  -2.3766  -1.1264   2.1885  0.073  0.62  -0.106E+01     48.  0.137E+01  0.401E-01  0.115E-02  0.203E-03   60   3  17   14  14   0   16  14  16    0.000    0.000       0
     393  0     47 00:00:35.62   0.7570A  -2.3767  -1.1353   2.1889  0.072  0.62  -0.106E+01     48.  0.137E+01  0.411E-01  0.117E-02  0.203E-03   60   3  17   14  14   0   17  14  16    0.000    0.000       0
     394  0     48 00:00:36.37   0.7570A  -2.3768  -1.1419   2.2142  0.072  0.62  -0.106E+01     50.  0.137E+01  0.422E-01  0.120E-02  0.203E-03   60   3  16   14  14   0   20  14  17    0.000    0.000       0
     395  0     49 00:00:37.13   0.7570A  -2.3768  -1.1459   2.2101  0.072  0.62  -0.106E+01     50.  0.137E+01  0.433E-01  0.122E-02  0.203E-03   60   2  15   14  14   0   20  14  17    0.000    0.000       0
     396  0     50 00:00:37.89   0.7570A  -2.3769  -1.1473   2.1891  0.072  0.62  -0.106E+01     50.  0.137E+01  0.444E-01  0.125E-02  0.202E-03   60   2  14   14  14   0   18  14  16    0.000    0.000       0
     397  0     51 00:00:38.64   0.7570A  -2.3770  -1.1467   2.1666  0.072  0.62  -0.106E+01     50.  0.137E+01  0.455E-01  0.127E-02  0.202E-03   60   2  13   14  14   0   18  14  16    0.000    0.000       0
     398  0     52 00:00:39.40   0.7570A  -2.3771  -1.1443   2.1360  0.072  0.63  -0.106E+01     52.  0.137E+01  0.466E-01  0.130E-02  0.202E-03   60   1  12   14  14   0   18  14  16    0.000    0.000       0
     399  0     53 00:00:40.16   0.7570A  -2.3771  -1.1623   2.1140  0.072  0.63  -0.106E+01     52.  0.137E+01  0.476E-01  0.132E-02  0.202E-03   60   1  11   22  14  13   20  14  16    0.000    0.000       0
     400  0     54 00:00:40.92   0.7570A  -2.3772  -1.2374   2.1096  0.072  0.63  -0.107E+01     52.  0.137E+01  0.487E-01  0.134E-02  0.203E-03   60   1  11   22  14  13   20  14  16    0.000    0.000       0
     401  0     55 00:00:41.67   0.7570A  -2.3773  -1.2997   2.0905  0.072  0.64  -0.107E+01     52.  0.137E+01  0.497E-01  0.136E-02  0.203E-03   60   1  10   22  14  13   20  14  16    0.000    0.000       0
     402  0     56 00:00:42.43   0.7570A  -2.3774  -1.3365   2.0920  0.072  0.64  -0.107E+01     52.  0.137E+01  0.507E-01  0.138E-02  0.203E-03   60   0   9   22  14  13   22  13  16    0.000    0.000       0
     403  0     57 00:00:43.19   0.7570A  -2.3774  -1.3596   2.1035  0.072  0.65  -0.107E+01     54.  0.137E+01  0.518E-01  0.139E-02  0.204E-03   60   0   8    8  14   5   23  13  16    0.000    0.000       0
     404  0     58 00:00:43.94   0.7570A  -2.3775  -1.3962   2.0921  0.072  0.65  -0.108E+01     54.  0.137E+01  0.527E-01  0.141E-02  0.204E-03   60   0   7    8  14   5   23  13  16    0.000    0.000       0
     405  0     59 00:00:44.70   0.7570A  -2.3776  -1.4256   2.0198  0.072  0.66  -0.108E+01     54.  0.137E+01  0.537E-01  0.142E-02  0.205E-03   60  19   6    8  14   5   23  13  16    0.000    0.000       0
     406  0     60 00:00:45.46   0.7570A  -2.3777  -1.4515   1.8968  0.072  0.66  -0.109E+01     54.  0.137E+01  0.547E-01  0.144E-02  0.206E-03   60  19   6    8  14   5   21  14  16    0.000    0.000       0
     407  0     61 00:00:46.21   0.7570A  -2.3777  -1.4755   1.8570  0.072  0.68  -0.108E+01     56.  0.137E+01  0.556E-01  0.146E-02  0.207E-03   60  19   5    8  14   5   21  14  16    0.000    0.000       0
     408  0     62 00:00:46.97   0.7570A  -2.3778  -1.4970   1.8597  0.072  0.68  -0.108E+01     56.  0.137E+01  0.565E-01  0.149E-02  0.208E-03   60  19   4    8  14   5   19  14  17    0.000    0.000       0
     409  0     63 00:00:47.73   0.7570A  -2.3779  -1.5145   1.8675  0.072  0.69  -0.108E+01     56.  0.137E+01  0.575E-01  0.151E-02  0.208E-03   60  18   3    8  14   5   19  14  17    0.000    0.000       0
     410  0     64 00:00:48.49   0.7570A  -2.3780  -1.5269   1.8729  0.072  0.69  -0.107E+01     56.  0.137E+01  0.584E-01  0.153E-02  0.209E-03   60  18   2    8  14   5   19  14  17    0.000    0.000       0
     411  0     65 00:00:49.24   0.7570A  -2.3780  -1.5331   1.8803  0.072  0.70  -0.107E+01     56.  0.137E+01  0.592E-01  0.155E-02  0.210E-03   60  18   2    8  14   5   19  14  17    0.000    0.000       0
     412  0     66 00:00:50.00   0.7570A  -2.3781  -1.5332   1.8901  0.072  0.71  -0.107E+01     58.  0.137E+01  0.601E-01  0.158E-02  0.210E-03   60  17   0    8  14   5   19  14  17    0.000    0.000       0
     413  0     67 00:00:50.76   0.7570A  -2.3782  -1.5267   1.8994  0.072  0.72  -0.106E+01     58.  0.137E+01  0.610E-01  0.160E-02  0.210E-03   60  17   0    8  14   5   19  14  17    0.000    0.000       0
     414  0     68 00:00:51.51   0.7570A  -2.3783  -1.5128   1.9059  0.072  0.73  -0.106E+01     58.  0.137E+01  0.618E-01  0.162E-02  0.211E-03   60  17  19    8  14   5   19  14  17    0.000    0.000       0
     415  0     69 00:00:52.27   0.7570A  -2.3783  -1.5285   1.9086  0.072  0.73  -0.106E+01     58.  0.137E+01  0.627E-01  0.165E-02  0.211E-03   60  16  17   23  12  10   19  14  17    0.000    0.000       0
     416  0     70 00:00:53.03   0.7570A  -2.3784  -1.5324   1.9091  0.072  0.74  -0.105E+01     58.  0.137E+01  0.635E-01  0.167E-02  0.212E-03   60  16  17   23  12  10   19  14  17    0.000    0.000       0
     417  0     71 00:00:53.78   0.7570A  -2.3785  -1.5209   1.9111  0.072  0.75  -0.105E+01     60.  0.137E+01  0.643E-01  0.169E-02  0.212E-03   60  16  16   23  12  10   19  14  17    0.000    0.000       0
     418  0     72 00:00:54.54   0.7570A  -2.3786  -1.5281   1.9285  0.072  0.76  -0.105E+01     60.  0.137E+01  0.651E-01  0.171E-02  0.213E-03   60  15  14   23  11   9   17  14  17    0.000    0.000       0
     419  0     73 00:00:55.30   0.7570A  -2.3786  -1.5237   1.9667  0.072  0.76  -0.105E+01     60.  0.137E+01  0.659E-01  0.173E-02  0.214E-03   60  15  14   23  11   9   17  14  17    0.000    0.000       0
     420  0     74 00:00:56.06   0.7570A  -2.3787  -1.5023   1.9915  0.072  0.77  -0.104E+01     60.  0.137E+01  0.667E-01  0.175E-02  0.214E-03   60  15  13   23  12   9   18  14  17    0.000    0.000       0
     421  0     75 00:00:56.81   0.7570A  -2.3788  -1.4950   2.0234  0.072  0.77  -0.104E+01     60.  0.137E+01  0.675E-01  0.177E-02  0.215E-03   60  14  11   23  12   9   18  14  17    0.000    0.000       0
     422  0     76 00:00:57.57   0.7570A  -2.3789  -1.4977   2.0370  0.072  0.78  -0.104E+01     60.  0.137E+01  0.682E-01  0.180E-02  0.217E-03   60  14  11   23  11   8   18  14  17    0.000    0.000       0
     423  0     77 00:00:58.33   0.7570A  -2.3789  -1.4828   2.0548  0.072  0.78  -0.104E+01     60.  0.137E+01  0.690E-01  0.182E-02  0.218E-03   60  14  10   23  11   8   19  14  17    0.000    0.000       0
     424  0     78 00:00:59.08   0.7570A  -2.3790  -1.4494   2.0570  0.072  0.79  -0.103E+01     60.  0.137E+01  0.697E-01  0.185E-02  0.219E-03   60  13   9   23  10   7   19  14  17    0.000    0.000       0
     425  0     79 00:00:59.84   0.7570A  -2.3791  -1.4404   2.0315  0.072  0.79  -0.103E+01     60.  0.137E+01  0.705E-01  0.187E-02  0.219E-03   60  13   8   23  11   7   19  14  17    0.000    0.000       0
     426  0     80 00:01:00.60   0.7570A  -2.3792  -1.4277   1.9922  0.072  0.81  -0.103E+01     62.  0.137E+01  0.712E-01  0.190E-02  0.220E-03   60  13   7   23  10   6   20  14  17    0.000    0.000       0
     427  0     81 00:01:01.35   0.7570A  -2.3792  -1.4445   1.9523  0.072  0.81  -0.102E+01     62.  0.137E+01  0.719E-01  0.192E-02  0.221E-03   60  13   7   23  10   6   20  14  17    0.000    0.000       0
     428  0     82 00:01:02.11   0.7570A  -2.3793  -1.4372   1.9055  0.072  0.81  -0.102E+01     62.  0.137E+01  0.726E-01  0.195E-02  0.222E-03   60  13   6   23  10   6   22  13  16    0.000    0.000       0
     429  0     83 00:01:02.87   0.7560A  -2.3794  -1.4125   1.9619  0.072  0.82  -0.103E+01     62.  0.137E+01  0.733E-01  0.197E-02  0.223E-03   60  12   5   11  12   4   22  13  16    0.000    0.000       0
     430  0     84 00:01:03.62   0.7560A  -2.3795  -1.4072   1.9975  0.072  0.82  -0.102E+01     62.  0.137E+01  0.739E-01  0.199E-02  0.224E-03   60  12   4   11  12   4   23  13  16    0.000    0.000       0
     431  0     85 00:01:04.38   0.7560A  -2.3795  -1.4225   1.9909  0.073  0.82  -0.102E+01     62.  0.137E+01  0.746E-01  0.201E-02  0.225E-03   60  12   4    8  13   5   23  13  16    0.000    0.000       0
     432  0     86 00:01:05.14   0.7560A  -2.3796  -1.4420   1.9955  0.073  0.83  -0.102E+01     62.  0.137E+01  0.752E-01  0.203E-02  0.227E-03   60  11   2    8  13   5   16  14  14    0.000    0.000       0
     433  0     87 00:01:05.89   0.7560A  -2.3797  -1.4597   1.9970  0.073  0.83  -0.102E+01     62.  0.137E+01  0.758E-01  0.204E-02  0.229E-03   60  11   2    8  13   5   16  14  14    0.000    0.000       0
     434  0     88 00:01:06.65   0.7560A  -2.3798  -1.4756   2.0149  0.073  0.83  -0.102E+01     62.  0.137E+01  0.763E-01  0.206E-02  0.231E-03   60  11   1    8  13   5   17  14  14    0.000    0.000       0
     435  0     89 00:01:07.40   0.7560A  -2.3798  -1.4897   2.0161  0.073  0.84  -0.102E+01     62.  0.137E+01  0.769E-01  0.208E-02  0.233E-03   60  11   1    8  13   5   17  14  14    0.000    0.000       0
     436  0     90 00:01:08.16   0.7560A  -2.3799  -1.5019   2.0432  0.073  0.84  -0.102E+01     62.  0.137E+01  0.774E-01  0.210E-02  0.235E-03   60  10  19    8  13   5   18  14  14    0.000    0.000       0
     437  0     91 00:01:08.92   0.7560A  -2.3800  -1.5120   2.0499  0.073  0.84  -0.102E+01     62.  0.137E+01  0.780E-01  0.212E-02  0.238E-03   60  10  18    8  13   5   18  14  14    0.000    0.000       0
     438  0     92 00:01:09.67   0.7560A  -2.3800  -1.5201   1.9605  0.073  0.84  -0.102E+01     62.  0.137E+01  0.785E-01  0.214E-02  0.240E-03   60  10  18    8  13   5   19  14  14    0.000    0.000       0
     439  0     93 00:01:10.43   0.7560A  -2.3801  -1.5261   1.9734  0.073  0.85  -0.101E+01     62.  0.137E+01  0.789E-01  0.216E-02  0.243E-03   60   9  16    8  13   5   19  14  14    0.000    0.000       0

  • palm/trunk/TESTS/cases/urban_environment_restart/MONITORING/urban_environment_restart_rc

    r3417 r3449  
    1 
    2 
    31
    42Run-control output:
    … …  
    75RUN  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
    86---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
    9   1     93 00:01:10.43   0.7560A  -2.3801  -1.5308   1.9738  0.073  0.09  -0.101E+01     62.  0.137E+01  0.789E-01  0.000E+00  0.000E+00   60   9  16    8  13   5   19  14  14    0.000    0.000       0
    10   1     94 00:01:11.18   0.7560A  -2.3802  -1.5352   1.8803  0.073  0.85  -0.101E+01     62.  0.137E+01  0.794E-01  0.218E-02  0.245E-03   60   9  16    8  13   5   19  14  14    0.000    0.000       0
    11   1     95 00:01:11.94   0.7560A  -2.3803  -1.5376   1.7937  0.073  0.86  -0.101E+01     64.  0.136E+01  0.798E-01  0.220E-02  0.246E-03   60   9  15    8  13   5   20  14  14    0.000    0.000       0
    12   1     96 00:01:12.70   0.7560A  -2.3803  -1.5381   1.6947  0.073  0.86  -0.100E+01     64.  0.136E+01  0.803E-01  0.222E-02  0.248E-03   60   8  13    8  13   5   14  14  13    0.000    0.000       0
    13   1     97 00:01:13.45   0.7560A  -2.3804  -1.5368   1.6872  0.073  0.86  -0.100E+01     64.  0.136E+01  0.807E-01  0.224E-02  0.249E-03   60   8  13    8  13   5   15  14  13    0.000    0.000       0
    14   1     98 00:01:14.21   0.7560A  -2.3805  -1.5340   1.6688  0.073  0.87  -0.100E+01     64.  0.136E+01  0.811E-01  0.226E-02  0.250E-03   60   8  12    8  13   5   19  14  13    0.000    0.000       0
    15   1     99 00:01:14.96   0.7560A  -2.3806  -1.5301   1.6452  0.073  0.87  -0.100E+01     64.  0.136E+01  0.815E-01  0.228E-02  0.251E-03   60   8  12    8  13   5   19  14  13    0.000    0.000       0
    16   1    100 00:01:15.72   0.7560A  -2.3806  -1.5253   1.6320  0.073  0.87  -0.100E+01     64.  0.136E+01  0.818E-01  0.229E-02  0.252E-03   60   7  10    8  13   5    7  14   3    0.000    0.000       0
    17   1    101 00:01:16.48   0.7560A  -2.3807  -1.5197   1.6338  0.073  0.87  -0.100E+01     64.  0.136E+01  0.822E-01  0.231E-02  0.252E-03   60   7   9    8  13   5    7  14   3    0.000    0.000       0
    18   1    102 00:01:17.23   0.7560A  -2.3808  -1.5135   1.6351  0.073  0.87  -0.100E+01     64.  0.136E+01  0.825E-01  0.232E-02  0.252E-03   60   7   9    8  13   5    7  14   3    0.000    0.000       0
    19   1    103 00:01:17.99   0.7560A  -2.3808  -1.5117   1.6357  0.073  0.87  -0.100E+01     64.  0.136E+01  0.828E-01  0.235E-02  0.251E-03   60   6   7    7  13   5    7  14   3    0.000    0.000       0
    20   1    104 00:01:18.74   0.7560A  -2.3809  -1.5200   1.6355  0.073  0.78  -0.100E+01     48.  0.136E+01  0.831E-01  0.237E-02  0.250E-03   60   6   7    7  13   5    7  14   3    0.000    0.000       0
    21   1    105 00:01:19.50   0.7560A  -2.3810  -1.5238   1.6415  0.073  0.78  -0.100E+01     48.  0.136E+01  0.834E-01  0.239E-02  0.249E-03   60   6   6    7  13   5   14  14  18    0.000    0.000       0
    22   1    106 00:01:20.26   0.7560A  -2.3811  -1.5235   1.6530  0.073  0.78  -0.100E+01     48.  0.136E+01  0.837E-01  0.241E-02  0.247E-03   60   5   4    7  13   5   14  14  18    0.000    0.000       0
    23   1    107 00:01:21.01   0.7560A  -2.3811  -1.5194   1.6614  0.073  0.77  -0.100E+01     48.  0.136E+01  0.839E-01  0.243E-02  0.245E-03   60   5   4    7  13   5   14  14  18    0.000    0.000       0
    24   1    108 00:01:21.77   0.7560A  -2.3812  -1.5118   1.6672  0.073  0.76  -0.100E+01     46.  0.136E+01  0.842E-01  0.245E-02  0.242E-03   60   5   3    7  13   5   14  14  18    0.000    0.000       0
    25   1    109 00:01:22.52   0.7560A  -2.3813  -1.5015   1.6706  0.073  0.75  -0.100E+01     46.  0.136E+01  0.844E-01  0.247E-02  0.239E-03   60   4   1    7  13   5   14  14  18    0.000    0.000       0
    26   1    110 00:01:23.28   0.7560A  -2.3814  -1.4984   1.6719  0.073  0.75  -0.100E+01     46.  0.136E+01  0.846E-01  0.249E-02  0.236E-03   60   4   0   13  13   0   15  14  18    0.000    0.000       0
    27   1    111 00:01:24.04   0.7560A  -2.3814  -1.4920   1.6906  0.073  0.77  -0.100E+01     50.  0.136E+01  0.848E-01  0.251E-02  0.233E-03   60   3  19   13  13   0   16  14  17    0.000    0.000       0
    28   1    112 00:01:24.79   0.7560A  -2.3815  -1.4835   1.7136  0.073  0.77  -0.100E+01     50.  0.135E+01  0.849E-01  0.252E-02  0.230E-03   60   3  18   13  13   0   16  14  17    0.000    0.000       0
    29   1    113 00:01:25.55   0.7560A  -2.3816  -1.4746   1.7339  0.073  0.77  -0.100E+01     52.  0.135E+01  0.851E-01  0.254E-02  0.226E-03   60   3  17   13  13   0   16  14  17    0.000    0.000       0
    30   1    114 00:01:26.30   0.7560A  -2.3817  -1.4660   1.7536  0.073  0.77  -0.101E+01     52.  0.135E+01  0.853E-01  0.255E-02  0.223E-03   60   2  16   13  13   0   17  14  17    0.000    0.000       0
    31   1    115 00:01:27.06   0.7560A  -2.3817  -1.4677   1.7755  0.073  0.76  -0.101E+01     52.  0.135E+01  0.854E-01  0.256E-02  0.221E-03   60   2  15   12  12   0   17  14  17    0.000    0.000       0
    32   1    116 00:01:27.82   0.7560A  -2.3818  -1.4702   1.7948  0.073  0.76  -0.102E+01     52.  0.135E+01  0.855E-01  0.256E-02  0.218E-03   60   2  14   12  12   0   17  14  17    0.000    0.000       0
    33   1    117 00:01:28.57   0.7560A  -2.3819  -1.4678   1.8172  0.073  0.75  -0.102E+01     52.  0.135E+01  0.856E-01  0.257E-02  0.215E-03   60   1  13   12  12   0   18  14  17    0.000    0.000       0
    34   1    118 00:01:29.33   0.7560A  -2.3820  -1.4621   1.8367  0.072  0.75  -0.102E+01     52.  0.135E+01  0.858E-01  0.257E-02  0.213E-03   60   1  12   12  12   0   18  14  17    0.000    0.000       0
    35   1    119 00:01:30.08   0.7560A  -2.3821  -1.4546   1.8535  0.072  0.74  -0.102E+01     52.  0.135E+01  0.859E-01  0.258E-02  0.211E-03   60   1  12   12  12   0   18  14  17    0.000    0.000       0
    36   1    120 00:01:30.84   0.7560A  -2.3821  -1.4460   1.8675  0.072  0.74  -0.102E+01     52.  0.135E+01  0.859E-01  0.258E-02  0.208E-03   60   1  11   12  12   0   18  14  17    0.000    0.000       0
    37   1    121 00:01:31.60   0.7560A  -2.3822  -1.4366   1.8787  0.072  0.73  -0.102E+01     52.  0.135E+01  0.860E-01  0.259E-02  0.207E-03   60   0  10   12  12   0   18  14  17    0.000    0.000       0
    38   1    122 00:01:32.35   0.7560A  -2.3823  -1.4261   1.8871  0.072  0.73  -0.104E+01     52.  0.135E+01  0.861E-01  0.259E-02  0.205E-03   60   0   9   12  12   0   18  14  17    0.000    0.000       0
    39   1    123 00:01:33.11   0.7560A  -2.3824  -1.4141   1.8980  0.072  0.72  -0.104E+01     52.  0.135E+01  0.861E-01  0.260E-02  0.204E-03   60   0   9   12  12   0   19  14  17    0.000    0.000       0
    40   1    124 00:01:33.86   0.7560A  -2.3824  -1.4003   1.9099  0.072  0.72  -0.103E+01     52.  0.135E+01  0.862E-01  0.260E-02  0.202E-03   60   0   8   12  12   0   19  14  17    0.000    0.000       0
    41   1    125 00:01:34.62   0.7560A  -2.3825  -1.3991   1.9196  0.072  0.72  -0.103E+01     52.  0.135E+01  0.862E-01  0.261E-02  0.201E-03   60  19   7   15  13   0   19  14  17    0.000    0.000       0
    42   1    126 00:01:35.38   0.7550A  -2.3826  -1.3979   1.9273  0.073  0.71  -0.103E+01     52.  0.134E+01  0.862E-01  0.261E-02  0.200E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
    43   1    127 00:01:36.13   0.7550A  -2.3827  -1.3954   1.9331  0.073  0.72  -0.102E+01     54.  0.134E+01  0.862E-01  0.261E-02  0.198E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
    44   1    128 00:01:36.89   0.7550A  -2.3827  -1.3915   1.9369  0.073  0.72  -0.102E+01     54.  0.134E+01  0.862E-01  0.262E-02  0.197E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
    45   1    129 00:01:37.64   0.7550A  -2.3828  -1.3859   1.9383  0.073  0.71  -0.101E+01     54.  0.134E+01  0.861E-01  0.262E-02  0.196E-03   60  19   5   15  13   0   19  14  17    0.000    0.000       0
    46   1    130 00:01:38.40   0.7550A  -2.3829  -1.3783   1.9457  0.073  0.71  -0.101E+01     54.  0.134E+01  0.860E-01  0.263E-02  0.194E-03   60  18   4   15  13   0   20  14  16    0.000    0.000       0
    47   1    131 00:01:39.15   0.7550A  -2.3829  -1.3719   1.9642  0.073  0.71  -0.101E+01     54.  0.134E+01  0.859E-01  0.263E-02  0.193E-03   60  18   3   16  12   0   20  14  16    0.000    0.000       0
    48   1    132 00:01:39.91   0.7550A  -2.3830  -1.3732   1.9796  0.073  0.71  -0.100E+01     54.  0.134E+01  0.858E-01  0.263E-02  0.192E-03   60  18   3   16  12   0   20  14  16    0.000    0.000       0
    49   1    133 00:01:40.66   0.7550A  -2.3831  -1.3711   2.0032  0.073  0.71  -0.100E+01     54.  0.134E+01  0.857E-01  0.264E-02  0.192E-03   60  18   3   16  12   0   22  14  16    0.000    0.000       0
    50   1    134 00:01:41.42   0.7550A  -2.3831  -1.3633   2.0324  0.074  0.71  -0.999E+00     54.  0.134E+01  0.855E-01  0.263E-02  0.191E-03   60  18   2   16  12   0   22  14  16    0.000    0.000       0
    51   1    135 00:01:42.17   0.7550A  -2.3832  -1.3539   2.0559  0.074  0.71  -0.994E+00     54.  0.134E+01  0.854E-01  0.263E-02  0.190E-03   60  18   2   17  12   0   22  14  16    0.000    0.000       0
    52   1    136 00:01:42.93   0.7550A  -2.3833  -1.3460   2.0721  0.074  0.71  -0.990E+00     54.  0.134E+01  0.852E-01  0.263E-02  0.190E-03   60  17   0   17  12   0   22  14  16    0.000    0.000       0
    53   1    137 00:01:43.68   0.7550A  -2.3834  -1.3416   2.0802  0.074  0.71  -0.988E+00     54.  0.134E+01  0.849E-01  0.262E-02  0.189E-03   60  17   0   18  10  19   22  14  16    0.000    0.000       0
    54   1    138 00:01:44.44   0.7550A  -2.3834  -1.3600   2.0794  0.074  0.71  -0.989E+00     54.  0.134E+01  0.847E-01  0.262E-02  0.189E-03   60   5   4   18  10  19   22  14  16    0.000    0.000       0
    55   1    139 00:01:45.19   0.7550A  -2.3835  -1.3688   2.0704  0.074  0.72  -0.996E+00     56.  0.133E+01  0.845E-01  0.261E-02  0.188E-03   60   5   3   18  10  19   22  14  16    0.000    0.000       0
    56   1    140 00:01:45.95   0.7550A  -2.3836  -1.3664   2.0535  0.074  0.73  -0.100E+01     56.  0.133E+01  0.842E-01  0.261E-02  0.188E-03   60   5   3   18  10  19   22  14  16    0.000    0.000       0
    57   1    141 00:01:46.70   0.7550A  -2.3836  -1.3594   2.0300  0.075  0.73  -0.983E+00     56.  0.133E+01  0.839E-01  0.260E-02  0.188E-03   60   5   3   17   9  19   22  14  16    0.000    0.000       0
    58   1    142 00:01:47.46   0.7550A  -2.3837  -1.3455   2.0024  0.075  0.73  -0.974E+00     56.  0.133E+01  0.836E-01  0.260E-02  0.189E-03   60   4   1   17   9  19   22  14  16    0.000    0.000       0
    59   1    143 00:01:48.21   0.7550A  -2.3838  -1.3126   1.9731  0.075  0.73  -0.970E+00     56.  0.133E+01  0.833E-01  0.260E-02  0.190E-03   60   4   1   17   9  19   22  14  16    0.000    0.000       0
    60   1    144 00:01:48.97   0.7550A  -2.3839  -1.2728   1.9441  0.075  0.73  -0.968E+00     56.  0.133E+01  0.829E-01  0.259E-02  0.192E-03   60   4   1   17   8  18   21  14  16    0.000    0.000       0
    61   1    145 00:01:49.72   0.7550A  -2.3839  -1.2768   1.9274  0.075  0.73  -0.968E+00     56.  0.133E+01  0.825E-01  0.259E-02  0.194E-03   60   4   1   17   8  18   20  14  16    0.000    0.000       0
    62   1    146 00:01:50.48   0.7550A  -2.3840  -1.2650   1.9103  0.075  0.73  -0.970E+00     56.  0.133E+01  0.822E-01  0.258E-02  0.195E-03   60   4   0   17   8  18   20  14  16    0.000    0.000       0
    63   1    147 00:01:51.23   0.7550A  -2.3841  -1.2371   1.8919  0.075  0.73  -0.972E+00     56.  0.133E+01  0.818E-01  0.258E-02  0.197E-03   60   4   0   17   8  18   20  14  16    0.000    0.000       0
    64   1    148 00:01:51.99   0.7550A  -2.3842  -1.2361   1.8720  0.074  0.74  -0.976E+00     56.  0.133E+01  0.814E-01  0.256E-02  0.198E-03   60   4   0   18  14   0   20  14  16    0.000    0.000       0
    65   1    149 00:01:52.74   0.7550A  -2.3842  -1.2427   1.8504  0.074  0.74  -0.981E+00     56.  0.133E+01  0.810E-01  0.255E-02  0.199E-03   60   4   0   18  14   0   20  14  16    0.000    0.000       0
    66   1    150 00:01:53.50   0.7550A  -2.3843  -1.2456   1.8266  0.074  0.74  -0.985E+00     56.  0.132E+01  0.806E-01  0.254E-02  0.199E-03   60   3  18   18  14   0   20  14  16    0.000    0.000       0
    67   1    151 00:01:54.25   0.7550A  -2.3844  -1.2438   1.8004  0.074  0.74  -0.989E+00     58.  0.132E+01  0.801E-01  0.253E-02  0.199E-03   60   3  18   18  14   0   20  14  16    0.000    0.000       0
    68   1    152 00:01:55.01   0.7550A  -2.3845  -1.2475   1.7718  0.074  0.74  -0.997E+00     58.  0.132E+01  0.797E-01  0.252E-02  0.199E-03   60   3  18   17  14   0   20  14  16    0.000    0.000       0
    69   1    153 00:01:55.76   0.7550A  -2.3846  -1.2498   1.7411  0.073  0.74  -0.103E+01     58.  0.132E+01  0.794E-01  0.251E-02  0.198E-03   60   3  17   17  14   0   20  14  16    0.000    0.000       0
    70   1    154 00:01:56.52   0.7550A  -2.3846  -1.2453   1.7087  0.073  0.74  -0.118E+01     58.  0.132E+01  0.790E-01  0.250E-02  0.197E-03   60   3  17   15  12   4   20  14  16    0.000    0.000       0
    71   1    155 00:01:57.27   0.7550A  -2.3847  -1.2420   1.6747  0.073  0.74  -0.102E+01     58.  0.132E+01  0.786E-01  0.249E-02  0.196E-03   60   9  17   15  12   4   20  14  16    0.000    0.000       0
    72   1    156 00:01:58.03   0.7550A  -2.3848  -1.2433   1.6394  0.072  0.74  -0.102E+01     58.  0.132E+01  0.782E-01  0.248E-02  0.196E-03   60   9  17   14  10   4   20  14  16    0.000    0.000       0
    73   1    157 00:01:58.78   0.7550A  -2.3849  -1.2641   1.6024  0.072  0.74  -0.117E+01     58.  0.132E+01  0.779E-01  0.248E-02  0.195E-03   60   9  16   14  10   4   20  14  16    0.000    0.000       0
    74   1    158 00:01:59.54   0.7550A  -2.3849  -1.2830   1.5630  0.072  0.73  -0.103E+01     58.  0.132E+01  0.775E-01  0.247E-02  0.195E-03   60   9  16   14  10   4   20  14  16    0.000    0.000       0
    75   1    159 00:02:00.29   0.7550A  -2.3850  -1.2968   1.5297  0.071  0.73  -0.101E+01     58.  0.132E+01  0.772E-01  0.246E-02  0.194E-03   60   9  15   14  10   4   21  14  16    0.000    0.000       0
    76   1    160 00:02:01.05   0.7550A  -2.3851  -1.3025   1.5371  0.071  0.73  -0.101E+01     58.  0.131E+01  0.769E-01  0.246E-02  0.194E-03   60   8  14   14  10   4    6  14   2    0.000    0.000       0
    77   1    161 00:02:01.80   0.7550A  -2.3852  -1.2990   1.5702  0.071  0.74  -0.101E+01     60.  0.131E+01  0.766E-01  0.245E-02  0.193E-03   60   8  13   14  10   4    7  14   2    0.000    0.000       0
    78   1    162 00:02:02.56   0.7550A  -2.3853  -1.2983   1.6187  0.070  0.74  -0.101E+01     60.  0.131E+01  0.763E-01  0.245E-02  0.193E-03   60   8  13   14   9   4    7  14   2    0.000    0.000       0
    79   1    163 00:02:03.31   0.7550A  -2.3853  -1.3131   1.6587  0.070  0.75  -0.101E+01     60.  0.131E+01  0.760E-01  0.243E-02  0.192E-03   60   8  12   14   9   4    7  14   2    0.000    0.000       0
    80   1    164 00:02:04.07   0.7550A  -2.3854  -1.3127   1.6889  0.070  0.68  -0.101E+01     46.  0.131E+01  0.757E-01  0.242E-02  0.191E-03   60   7  11   14   9   4    7  14   2    0.000    0.000       0
    81   1    165 00:02:04.82   0.7550A  -2.3855  -1.3052   1.7100  0.069  0.68  -0.102E+01     46.  0.131E+01  0.755E-01  0.241E-02  0.191E-03   60   7  11   15  10   4    7  14   2    0.000    0.000       0
    82   1    166 00:02:05.58   0.7550A  -2.3856  -1.2859   1.7306  0.069  0.69  -0.102E+01     46.  0.131E+01  0.753E-01  0.240E-02  0.190E-03   60   7  10   15  10   4    8  14   2    0.000    0.000       0
    83   1    167 00:02:06.33   0.7550A  -2.3856  -1.2703   1.7662  0.069  0.69  -0.102E+01     46.  0.131E+01  0.750E-01  0.239E-02  0.191E-03   60   7  10   15   9   4    8  14   2    0.000    0.000       0
    84   1    168 00:02:07.09   0.7540A  -2.3857  -1.3004   1.7907  0.069  0.70  -0.102E+01     46.  0.131E+01  0.748E-01  0.237E-02  0.191E-03   60   7   9    8  14   5    8  14   2    0.000    0.000       0
    85   1    169 00:02:07.84   0.7540A  -2.3858  -1.3226   1.8063  0.068  0.70  -0.102E+01     46.  0.131E+01  0.746E-01  0.236E-02  0.192E-03   60   7   9    8  14   5    8  14   2    0.000    0.000       0
    86   1    170 00:02:08.59   0.7540A  -2.3858  -1.3359   1.8154  0.068  0.71  -0.101E+01     46.  0.131E+01  0.744E-01  0.234E-02  0.193E-03   60   6   8    8  14   5    8  14   2    0.000    0.000       0
    87   1    171 00:02:09.35   0.7540A  -2.3859  -1.3415   1.8195  0.068  0.71  -0.101E+01     46.  0.131E+01  0.743E-01  0.232E-02  0.195E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
    88   1    172 00:02:10.10   0.7540A  -2.3860  -1.3413   1.8192  0.068  0.72  -0.101E+01     46.  0.131E+01  0.741E-01  0.230E-02  0.196E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
    89   1    173 00:02:10.86   0.7540A  -2.3861  -1.3374   1.8143  0.068  0.72  -0.101E+01     46.  0.130E+01  0.740E-01  0.229E-02  0.197E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
    90   1    174 00:02:11.61   0.7540A  -2.3861  -1.3318   1.8049  0.068  0.72  -0.101E+01     46.  0.130E+01  0.738E-01  0.227E-02  0.198E-03   60   6   6    8  14   5    8  14   2    0.000    0.000       0
    91   1    175 00:02:12.36   0.7540A  -2.3862  -1.3268   1.7918  0.068  0.72  -0.101E+01     46.  0.130E+01  0.737E-01  0.225E-02  0.199E-03   60   6   6    8  14   5   12  14  17    0.000    0.000       0
    92   1    176 00:02:13.12   0.7540A  -2.3863  -1.3247   1.8233  0.068  0.72  -0.102E+01     46.  0.130E+01  0.736E-01  0.224E-02  0.200E-03   60   5   5    8  14   5   12  14  17    0.000    0.000       0
    93   1    177 00:02:13.87   0.7540A  -2.3863  -1.3489   1.8481  0.068  0.72  -0.104E+01     46.  0.130E+01  0.735E-01  0.223E-02  0.200E-03   60   5   4   13  11   4   13  14  17    0.000    0.000       0
    94   1    178 00:02:14.63   0.7540A  -2.3864  -1.3686   1.8745  0.068  0.74  -0.102E+01     48.  0.130E+01  0.734E-01  0.221E-02  0.199E-03   60   5   4   13  11   4   12  14  16    0.000    0.000       0
    95   1    179 00:02:15.38   0.7540A  -2.3865  -1.3571   1.9040  0.068  0.74  -0.101E+01     48.  0.130E+01  0.734E-01  0.220E-02  0.199E-03   60   5   4   13  11   4   13  14  16    0.000    0.000       0
    96   1    180 00:02:16.13   0.7540A  -2.3865  -1.3713   1.9493  0.068  0.74  -0.100E+01     48.  0.130E+01  0.733E-01  0.219E-02  0.198E-03   60   5   4    8  14   5   13  14  16    0.000    0.000       0
    97   1    181 00:02:16.89   0.7540A  -2.3866  -1.3933   1.9866  0.068  0.74  -0.998E+00     48.  0.130E+01  0.733E-01  0.217E-02  0.196E-03   60   5   3    8  14   5   14  14  16    0.000    0.000       0
    98   1    182 00:02:17.64   0.7540A  -2.3867  -1.4153   2.0088  0.068  0.74  -0.993E+00     48.  0.130E+01  0.733E-01  0.216E-02  0.195E-03   60   5   3    8  14   5   15  14  16    0.000    0.000       0
    99   1    183 00:02:18.40   0.7540A  -2.3868  -1.4411   2.0533  0.069  0.75  -0.989E+00     48.  0.130E+01  0.733E-01  0.215E-02  0.195E-03   60  12   3    7  14   5   15  14  16    0.000    0.000       0
    100   1    184 00:02:19.15   0.7540A  -2.3868  -1.4644   2.0732  0.069  0.75  -0.985E+00     48.  0.130E+01  0.733E-01  0.215E-02  0.194E-03   60  11   2    7  14   5   16  14  16    0.000    0.000       0
    101   1    185 00:02:19.90   0.7540A  -2.3869  -1.4798   2.0875  0.069  0.75  -0.983E+00     48.  0.130E+01  0.733E-01  0.215E-02  0.195E-03   60  11   2    7  14   5   16  14  16    0.000    0.000       0
    102   1    186 00:02:20.66   0.7540A  -2.3870  -1.4873   2.1132  0.069  0.77  -0.982E+00     50.  0.130E+01  0.734E-01  0.215E-02  0.194E-03   60  11   1    7  14   5   17  14  16    0.000    0.000       0
     7  1     93 00:01:10.43   0.7560A  -2.3801  -1.5261   1.9734  0.073  0.09  -0.101E+01     62.  0.137E+01  0.789E-01  0.000E+00  0.000E+00   60   9  16    8  13   5   19  14  14    0.000    0.000       0
     8  1     94 00:01:11.18   0.7560A  -2.3802  -1.5300   1.8775  0.073  0.85  -0.101E+01     62.  0.137E+01  0.794E-01  0.218E-02  0.245E-03   60   9  16    8  13   5   19  14  14    0.000    0.000       0
     9  1     95 00:01:11.94   0.7560A  -2.3803  -1.5320   1.7916  0.073  0.86  -0.101E+01     64.  0.136E+01  0.799E-01  0.220E-02  0.247E-03   60   9  15    8  13   5   20  14  14    0.000    0.000       0
     10  1     96 00:01:12.70   0.7560A  -2.3803  -1.5320   1.7005  0.073  0.86  -0.101E+01     64.  0.136E+01  0.803E-01  0.222E-02  0.248E-03   60   8  13    8  13   5   14  14  13    0.000    0.000       0
     11  1     97 00:01:13.45   0.7560A  -2.3804  -1.5302   1.6938  0.073  0.87  -0.101E+01     64.  0.136E+01  0.807E-01  0.224E-02  0.249E-03   60   8  13    8  13   5   15  14  13    0.000    0.000       0
     12  1     98 00:01:14.21   0.7560A  -2.3805  -1.5272   1.6685  0.073  0.87  -0.101E+01     64.  0.136E+01  0.811E-01  0.226E-02  0.251E-03   60   8  12    8  13   5   19  14  13    0.000    0.000       0
     13  1     99 00:01:14.96   0.7560A  -2.3806  -1.5231   1.6484  0.073  0.87  -0.101E+01     64.  0.136E+01  0.815E-01  0.228E-02  0.252E-03   60   8  12    8  13   5   19  14  13    0.000    0.000       0
     14  1    100 00:01:15.72   0.7560A  -2.3806  -1.5181   1.6316  0.073  0.87  -0.101E+01     64.  0.136E+01  0.819E-01  0.229E-02  0.252E-03   60   7  10    8  13   5    7  14   3    0.000    0.000       0
     15  1    101 00:01:16.48   0.7560A  -2.3807  -1.5125   1.6332  0.073  0.87  -0.101E+01     64.  0.136E+01  0.822E-01  0.231E-02  0.252E-03   60   7   9    8  13   5    7  14   3    0.000    0.000       0
     16  1    102 00:01:17.23   0.7560A  -2.3808  -1.5063   1.6342  0.073  0.87  -0.101E+01     64.  0.136E+01  0.826E-01  0.233E-02  0.252E-03   60   7   9    8  13   5    7  14   3    0.000    0.000       0
     17  1    103 00:01:17.99   0.7560A  -2.3808  -1.5115   1.6345  0.073  0.87  -0.101E+01     64.  0.136E+01  0.829E-01  0.235E-02  0.252E-03   60   6   7    7  13   5    7  14   3    0.000    0.000       0
     18  1    104 00:01:18.74   0.7560A  -2.3809  -1.5189   1.6340  0.073  0.78  -0.101E+01     48.  0.136E+01  0.832E-01  0.237E-02  0.251E-03   60   6   7    7  13   5    7  14   3    0.000    0.000       0
     19  1    105 00:01:19.50   0.7560A  -2.3810  -1.5219   1.6381  0.073  0.78  -0.100E+01     48.  0.136E+01  0.835E-01  0.239E-02  0.249E-03   60   6   6    7  13   5   14  14  18    0.000    0.000       0
     20  1    106 00:01:20.26   0.7560A  -2.3811  -1.5207   1.6497  0.073  0.78  -0.100E+01     48.  0.136E+01  0.837E-01  0.241E-02  0.247E-03   60   5   4    7  13   5   14  14  18    0.000    0.000       0
     21  1    107 00:01:21.01   0.7560A  -2.3811  -1.5156   1.6583  0.073  0.77  -0.100E+01     48.  0.136E+01  0.840E-01  0.243E-02  0.245E-03   60   5   4    7  13   5   14  14  18    0.000    0.000       0
     22  1    108 00:01:21.77   0.7560A  -2.3812  -1.5071   1.6642  0.073  0.76  -0.100E+01     46.  0.136E+01  0.842E-01  0.245E-02  0.242E-03   60   5   3    7  13   5   14  14  18    0.000    0.000       0
     23  1    109 00:01:22.52   0.7560A  -2.3813  -1.5018   1.6676  0.073  0.76  -0.100E+01     46.  0.136E+01  0.844E-01  0.247E-02  0.239E-03   60   4   1   13  13   0   14  14  18    0.000    0.000       0
     24  1    110 00:01:23.28   0.7560A  -2.3814  -1.4993   1.6690  0.073  0.75  -0.100E+01     46.  0.136E+01  0.846E-01  0.249E-02  0.236E-03   60   4   0   13  13   0   15  14  18    0.000    0.000       0
     25  1    111 00:01:24.04   0.7560A  -2.3814  -1.4928   1.6874  0.073  0.75  -0.101E+01     46.  0.136E+01  0.848E-01  0.251E-02  0.233E-03   60   3  19   13  13   0   16  14  17    0.000    0.000       0
     26  1    112 00:01:24.79   0.7560A  -2.3815  -1.4843   1.7107  0.073  0.77  -0.101E+01     50.  0.135E+01  0.850E-01  0.252E-02  0.230E-03   60   3  18   13  13   0   16  14  17    0.000    0.000       0
     27  1    113 00:01:25.55   0.7560A  -2.3816  -1.4754   1.7311  0.073  0.77  -0.101E+01     52.  0.135E+01  0.851E-01  0.254E-02  0.227E-03   60   3  17   13  13   0   16  14  17    0.000    0.000       0
     28  1    114 00:01:26.30   0.7560A  -2.3817  -1.4669   1.7503  0.073  0.77  -0.101E+01     52.  0.135E+01  0.853E-01  0.255E-02  0.224E-03   60   2  16   13  13   0   17  14  17    0.000    0.000       0
     29  1    115 00:01:27.06   0.7560A  -2.3817  -1.4690   1.7722  0.073  0.77  -0.101E+01     52.  0.135E+01  0.854E-01  0.256E-02  0.221E-03   60   2  15   12  12   0   17  14  17    0.000    0.000       0
     30  1    116 00:01:27.82   0.7560A  -2.3818  -1.4712   1.7916  0.073  0.76  -0.102E+01     52.  0.135E+01  0.856E-01  0.256E-02  0.218E-03   60   2  14   12  12   0   17  14  17    0.000    0.000       0
     31  1    117 00:01:28.57   0.7560A  -2.3819  -1.4685   1.8139  0.073  0.76  -0.102E+01     52.  0.135E+01  0.857E-01  0.257E-02  0.215E-03   60   1  13   12  12   0   18  14  17    0.000    0.000       0
     32  1    118 00:01:29.33   0.7560A  -2.3820  -1.4625   1.8333  0.073  0.75  -0.102E+01     52.  0.135E+01  0.858E-01  0.257E-02  0.213E-03   60   1  12   12  12   0   18  14  17    0.000    0.000       0
     33  1    119 00:01:30.08   0.7560A  -2.3821  -1.4548   1.8501  0.073  0.75  -0.102E+01     52.  0.135E+01  0.859E-01  0.258E-02  0.211E-03   60   1  12   12  12   0   18  14  17    0.000    0.000       0
     34  1    120 00:01:30.84   0.7560A  -2.3821  -1.4461   1.8641  0.073  0.74  -0.102E+01     52.  0.135E+01  0.860E-01  0.258E-02  0.209E-03   60   1  11   12  12   0   18  14  17    0.000    0.000       0
     35  1    121 00:01:31.60   0.7560A  -2.3822  -1.4366   1.8753  0.073  0.74  -0.102E+01     52.  0.135E+01  0.860E-01  0.258E-02  0.207E-03   60   0  10   12  12   0   18  14  17    0.000    0.000       0
     36  1    122 00:01:32.35   0.7560A  -2.3823  -1.4260   1.8837  0.073  0.73  -0.104E+01     52.  0.135E+01  0.861E-01  0.259E-02  0.205E-03   60   0   9   12  12   0   18  14  17    0.000    0.000       0
     37  1    123 00:01:33.11   0.7560A  -2.3824  -1.4140   1.8943  0.073  0.73  -0.104E+01     52.  0.135E+01  0.862E-01  0.260E-02  0.204E-03   60   0   9   12  12   0   19  14  17    0.000    0.000       0
     38  1    124 00:01:33.86   0.7560A  -2.3824  -1.4007   1.9063  0.073  0.72  -0.103E+01     52.  0.135E+01  0.862E-01  0.260E-02  0.202E-03   60   0   8   15  13   0   19  14  17    0.000    0.000       0
     39  1    125 00:01:34.62   0.7560A  -2.3825  -1.4006   1.9161  0.073  0.72  -0.103E+01     52.  0.135E+01  0.862E-01  0.260E-02  0.201E-03   60  19   7   15  13   0   19  14  17    0.000    0.000       0
     40  1    126 00:01:35.38   0.7550A  -2.3826  -1.3994   1.9239  0.073  0.71  -0.103E+01     52.  0.134E+01  0.862E-01  0.261E-02  0.200E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
     41  1    127 00:01:36.13   0.7550A  -2.3827  -1.3969   1.9298  0.073  0.72  -0.102E+01     54.  0.134E+01  0.862E-01  0.261E-02  0.199E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
     42  1    128 00:01:36.89   0.7550A  -2.3827  -1.3930   1.9337  0.073  0.72  -0.102E+01     54.  0.134E+01  0.862E-01  0.262E-02  0.197E-03   60  19   6   15  13   0   19  14  17    0.000    0.000       0
     43  1    129 00:01:37.64   0.7550A  -2.3828  -1.3875   1.9351  0.073  0.72  -0.101E+01     54.  0.134E+01  0.861E-01  0.262E-02  0.196E-03   60  19   5   15  13   0   19  14  17    0.000    0.000       0
     44  1    130 00:01:38.40   0.7550A  -2.3829  -1.3798   1.9398  0.073  0.71  -0.101E+01     54.  0.134E+01  0.861E-01  0.263E-02  0.195E-03   60  18   4   15  13   0   20  14  16    0.000    0.000       0
     45  1    131 00:01:39.15   0.7550A  -2.3829  -1.3727   1.9581  0.073  0.71  -0.101E+01     54.  0.134E+01  0.860E-01  0.263E-02  0.194E-03   60  18   3   16  12   0   20  14  16    0.000    0.000       0
     46  1    132 00:01:39.91   0.7550A  -2.3830  -1.3741   1.9733  0.074  0.71  -0.100E+01     54.  0.134E+01  0.858E-01  0.263E-02  0.193E-03   60  18   3   16  12   0   20  14  16    0.000    0.000       0
     47  1    133 00:01:40.66   0.7550A  -2.3831  -1.3720   1.9969  0.074  0.71  -0.100E+01     54.  0.134E+01  0.857E-01  0.263E-02  0.192E-03   60  18   3   16  12   0   22  14  16    0.000    0.000       0
     48  1    134 00:01:41.42   0.7550A  -2.3831  -1.3642   2.0263  0.074  0.71  -0.999E+00     54.  0.134E+01  0.856E-01  0.263E-02  0.191E-03   60  18   2   16  12   0   22  14  16    0.000    0.000       0
     49  1    135 00:01:42.17   0.7550A  -2.3832  -1.3542   2.0499  0.074  0.71  -0.994E+00     54.  0.134E+01  0.854E-01  0.263E-02  0.190E-03   60  18   2   17  12   0   22  14  16    0.000    0.000       0
     50  1    136 00:01:42.93   0.7550A  -2.3833  -1.3461   2.0661  0.074  0.71  -0.989E+00     54.  0.134E+01  0.852E-01  0.262E-02  0.190E-03   60  17   0   17  12   0   22  14  16    0.000    0.000       0
     51  1    137 00:01:43.68   0.7550A  -2.3834  -1.3453   2.0742  0.074  0.71  -0.987E+00     54.  0.134E+01  0.850E-01  0.262E-02  0.189E-03   60  17   0   18  10  19   22  14  16    0.000    0.000       0
     52  1    138 00:01:44.44   0.7550A  -2.3834  -1.3633   2.0734  0.075  0.72  -0.989E+00     54.  0.134E+01  0.847E-01  0.262E-02  0.189E-03   60   5   4   18  10  19   22  14  16    0.000    0.000       0
     53  1    139 00:01:45.19   0.7550A  -2.3835  -1.3715   2.0642  0.075  0.73  -0.996E+00     56.  0.133E+01  0.845E-01  0.261E-02  0.188E-03   60   5   3   18  10  19   22  14  16    0.000    0.000       0
     54  1    140 00:01:45.95   0.7550A  -2.3836  -1.3685   2.0470  0.075  0.73  -0.100E+01     56.  0.133E+01  0.842E-01  0.261E-02  0.188E-03   60   5   3   18  10  19   22  14  16    0.000    0.000       0
     55  1    141 00:01:46.70   0.7550A  -2.3836  -1.3628   2.0231  0.075  0.73  -0.984E+00     56.  0.133E+01  0.839E-01  0.260E-02  0.188E-03   60   5   3   17   9  19   22  14  16    0.000    0.000       0
     56  1    142 00:01:47.46   0.7550A  -2.3837  -1.3481   1.9952  0.075  0.73  -0.974E+00     56.  0.133E+01  0.836E-01  0.260E-02  0.189E-03   60   4   1   17   9  19   22  14  16    0.000    0.000       0
     57  1    143 00:01:48.21   0.7550A  -2.3838  -1.3145   1.9654  0.075  0.73  -0.970E+00     56.  0.133E+01  0.833E-01  0.260E-02  0.190E-03   60   4   1   17   9  19   22  14  16    0.000    0.000       0
     58  1    144 00:01:48.97   0.7550A  -2.3839  -1.2749   1.9358  0.075  0.74  -0.968E+00     56.  0.133E+01  0.829E-01  0.259E-02  0.192E-03   60   4   1   17   8  18   22  14  16    0.000    0.000       0
     59  1    145 00:01:49.72   0.7550A  -2.3839  -1.2782   1.9151  0.075  0.74  -0.968E+00     56.  0.133E+01  0.826E-01  0.259E-02  0.194E-03   60   4   1   17   8  18   20  14  16    0.000    0.000       0
     60  1    146 00:01:50.48   0.7550A  -2.3840  -1.2657   1.8973  0.075  0.74  -0.970E+00     56.  0.133E+01  0.822E-01  0.258E-02  0.195E-03   60   4   0   17   8  18   20  14  16    0.000    0.000       0
     61  1    147 00:01:51.23   0.7550A  -2.3841  -1.2373   1.8782  0.075  0.74  -0.972E+00     56.  0.133E+01  0.818E-01  0.258E-02  0.197E-03   60   4   0   17   8  18   20  14  16    0.000    0.000       0
     62  1    148 00:01:51.99   0.7550A  -2.3842  -1.2352   1.8575  0.075  0.74  -0.976E+00     56.  0.133E+01  0.814E-01  0.256E-02  0.198E-03   60   4   0   18  14   0   20  14  16    0.000    0.000       0
     63  1    149 00:01:52.74   0.7550A  -2.3842  -1.2426   1.8352  0.075  0.74  -0.980E+00     56.  0.132E+01  0.810E-01  0.255E-02  0.199E-03   60   3  18   18  14   0   20  14  16    0.000    0.000       0
     64  1    150 00:01:53.50   0.7550A  -2.3843  -1.2463   1.8107  0.075  0.74  -0.984E+00     56.  0.132E+01  0.806E-01  0.254E-02  0.200E-03   60   3  18   18  14   0   20  14  16    0.000    0.000       0
     65  1    151 00:01:54.25   0.7550A  -2.3844  -1.2453   1.7839  0.074  0.75  -0.987E+00     58.  0.132E+01  0.802E-01  0.253E-02  0.200E-03   60   3  18   18  14   0   20  14  16    0.000    0.000       0
     66  1    152 00:01:55.01   0.7550A  -2.3845  -1.2494   1.7547  0.074  0.75  -0.992E+00     58.  0.132E+01  0.798E-01  0.251E-02  0.199E-03   60   3  18   17  14   0   20  14  16    0.000    0.000       0
     67  1    153 00:01:55.76   0.7550A  -2.3846  -1.2523   1.7235  0.074  0.75  -0.100E+01     58.  0.132E+01  0.794E-01  0.250E-02  0.199E-03   60   3  17   17  14   0   20  14  16    0.000    0.000       0
     68  1    154 00:01:56.52   0.7550A  -2.3846  -1.2479   1.6903  0.073  0.74  -0.105E+01     58.  0.132E+01  0.790E-01  0.250E-02  0.198E-03   60   3  17   15  12   4   20  14  16    0.000    0.000       0
     69  1    155 00:01:57.27   0.7550A  -2.3847  -1.2462   1.6556  0.073  0.74  -0.104E+01     58.  0.132E+01  0.787E-01  0.249E-02  0.197E-03   60   9  17   15  12   4   20  14  16    0.000    0.000       0
     70  1    156 00:01:58.03   0.7550A  -2.3848  -1.2408   1.6196  0.073  0.74  -0.104E+01     58.  0.132E+01  0.783E-01  0.248E-02  0.196E-03   60   9  17   14  10   4   20  14  16    0.000    0.000       0
     71  1    157 00:01:58.78   0.7550A  -2.3849  -1.2635   1.5819  0.072  0.74  -0.111E+01     58.  0.132E+01  0.780E-01  0.247E-02  0.196E-03   60   9  16   15  11   4   20  14  16    0.000    0.000       0
     72  1    158 00:01:59.54   0.7550A  -2.3849  -1.2827   1.5433  0.072  0.74  -0.102E+01     58.  0.132E+01  0.776E-01  0.247E-02  0.195E-03   60   9  16   14  10   4   22  14  16    0.000    0.000       0
     73  1    159 00:02:00.29   0.7550A  -2.3850  -1.2973   1.5082  0.072  0.74  -0.102E+01     58.  0.132E+01  0.773E-01  0.246E-02  0.195E-03   60   9  15   14  10   4   21  14  16    0.000    0.000       0
     74  1    160 00:02:01.05   0.7550A  -2.3851  -1.3035   1.5344  0.071  0.74  -0.101E+01     58.  0.131E+01  0.770E-01  0.246E-02  0.194E-03   60   8  14   14  10   4    6  14   2    0.000    0.000       0
     75  1    161 00:02:01.80   0.7550A  -2.3852  -1.2995   1.5661  0.071  0.75  -0.101E+01     60.  0.131E+01  0.767E-01  0.245E-02  0.194E-03   60   8  13   14  10   4    7  14   2    0.000    0.000       0
     76  1    162 00:02:02.56   0.7550A  -2.3853  -1.2959   1.6162  0.071  0.75  -0.101E+01     60.  0.131E+01  0.764E-01  0.245E-02  0.193E-03   60   8  13   14   9   4    7  14   2    0.000    0.000       0
     77  1    163 00:02:03.31   0.7550A  -2.3853  -1.3113   1.6575  0.070  0.75  -0.101E+01     60.  0.131E+01  0.761E-01  0.244E-02  0.192E-03   60   8  12   14   9   4    7  14   2    0.000    0.000       0
     78  1    164 00:02:04.07   0.7550A  -2.3854  -1.3115   1.6885  0.070  0.68  -0.102E+01     46.  0.131E+01  0.758E-01  0.242E-02  0.191E-03   60   7  11   14   9   4    7  14   2    0.000    0.000       0
     79  1    165 00:02:04.82   0.7550A  -2.3855  -1.3016   1.7097  0.070  0.69  -0.103E+01     46.  0.131E+01  0.756E-01  0.241E-02  0.191E-03   60   7  11   15  10   4    7  14   2    0.000    0.000       0
     80  1    166 00:02:05.58   0.7550A  -2.3856  -1.2824   1.7295  0.069  0.69  -0.107E+01     46.  0.131E+01  0.754E-01  0.240E-02  0.191E-03   60   7  10   15  10   4    8  14   2    0.000    0.000       0
     81  1    167 00:02:06.33   0.7550A  -2.3856  -1.2684   1.7661  0.069  0.70  -0.104E+01     46.  0.131E+01  0.751E-01  0.239E-02  0.191E-03   60   7  10    8  14   5    8  14   2    0.000    0.000       0
     82  1    168 00:02:07.09   0.7540A  -2.3857  -1.3015   1.7909  0.069  0.70  -0.102E+01     46.  0.131E+01  0.749E-01  0.238E-02  0.191E-03   60   7   9    8  14   5    8  14   2    0.000    0.000       0
     83  1    169 00:02:07.84   0.7540A  -2.3858  -1.3252   1.8063  0.068  0.71  -0.102E+01     46.  0.131E+01  0.747E-01  0.236E-02  0.192E-03   60   7   9    8  14   5    8  14   2    0.000    0.000       0
     84  1    170 00:02:08.59   0.7540A  -2.3858  -1.3400   1.8150  0.068  0.71  -0.101E+01     46.  0.131E+01  0.746E-01  0.234E-02  0.193E-03   60   6   8    8  14   5    8  14   2    0.000    0.000       0
     85  1    171 00:02:09.35   0.7540A  -2.3859  -1.3469   1.8187  0.068  0.71  -0.101E+01     46.  0.131E+01  0.744E-01  0.232E-02  0.194E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
     86  1    172 00:02:10.10   0.7540A  -2.3860  -1.3476   1.8180  0.068  0.72  -0.101E+01     46.  0.131E+01  0.742E-01  0.230E-02  0.196E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
     87  1    173 00:02:10.86   0.7540A  -2.3861  -1.3440   1.8130  0.068  0.72  -0.101E+01     46.  0.130E+01  0.741E-01  0.229E-02  0.197E-03   60   6   7    8  14   5    8  14   2    0.000    0.000       0
     88  1    174 00:02:11.61   0.7540A  -2.3861  -1.3381   1.8036  0.068  0.72  -0.101E+01     46.  0.130E+01  0.740E-01  0.227E-02  0.198E-03   60   6   6    8  14   5    8  14   2    0.000    0.000       0
     89  1    175 00:02:12.36   0.7540A  -2.3862  -1.3321   1.7997  0.068  0.72  -0.101E+01     46.  0.130E+01  0.738E-01  0.225E-02  0.199E-03   60   6   6    8  14   5   12  14  17    0.000    0.000       0
     90  1    176 00:02:13.12   0.7540A  -2.3863  -1.3283   1.8253  0.068  0.72  -0.102E+01     46.  0.130E+01  0.737E-01  0.224E-02  0.199E-03   60   5   5    8  14   5   12  14  17    0.000    0.000       0
     91  1    177 00:02:13.87   0.7540A  -2.3863  -1.3551   1.8597  0.068  0.73  -0.103E+01     46.  0.130E+01  0.736E-01  0.223E-02  0.199E-03   60   5   4   13  11   4   13  14  17    0.000    0.000       0
     92  1    178 00:02:14.63   0.7540A  -2.3864  -1.3705   1.8771  0.068  0.74  -0.103E+01     48.  0.130E+01  0.736E-01  0.221E-02  0.199E-03   60   5   4   13  11   4   13  14  17    0.000    0.000       0
     93  1    179 00:02:15.38   0.7540A  -2.3865  -1.3545   1.9080  0.068  0.74  -0.101E+01     48.  0.130E+01  0.735E-01  0.220E-02  0.199E-03   60   5   4   13  11   4   14  14  17    0.000    0.000       0
     94  1    180 00:02:16.13   0.7540A  -2.3865  -1.3690   1.9433  0.068  0.74  -0.101E+01     48.  0.130E+01  0.735E-01  0.219E-02  0.198E-03   60   5   4    8  14   5   13  14  16    0.000    0.000       0
     95  1    181 00:02:16.89   0.7540A  -2.3866  -1.3922   1.9860  0.068  0.74  -0.999E+00     48.  0.130E+01  0.734E-01  0.217E-02  0.197E-03   60   5   3    8  14   5   14  14  16    0.000    0.000       0
     96  1    182 00:02:17.64   0.7540A  -2.3867  -1.4187   2.0186  0.069  0.75  -0.993E+00     48.  0.130E+01  0.734E-01  0.216E-02  0.196E-03   60   5   3    7  14   5   15  14  16    0.000    0.000       0
     97  1    183 00:02:18.40   0.7540A  -2.3868  -1.4489   2.0535  0.069  0.75  -0.989E+00     48.  0.130E+01  0.734E-01  0.215E-02  0.195E-03   60  12   3    7  14   5   15  14  16    0.000    0.000       0
     98  1    184 00:02:19.15   0.7540A  -2.3868  -1.4716   2.0831  0.069  0.75  -0.986E+00     48.  0.130E+01  0.735E-01  0.215E-02  0.195E-03   60  11   2    7  14   5   16  14  16    0.000    0.000       0
     99  1    185 00:02:19.90   0.7540A  -2.3869  -1.4855   2.0896  0.069  0.75  -0.984E+00     48.  0.130E+01  0.735E-01  0.215E-02  0.195E-03   60  11   2    7  14   5   16  14  16    0.000    0.000       0
     100  1    186 00:02:20.66   0.7540A  -2.3870  -1.4913   2.1242  0.069  0.77  -0.983E+00     50.  0.130E+01  0.735E-01  0.215E-02  0.195E-03   60  11   1    7  14   5   17  14  16    0.000    0.000       0
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |