Changeset 1555


Ignore:
Timestamp:
Mar 4, 2015 5:44:27 PM (6 years ago)
Author:
maronga
Message:

LSM output of r_a and r_s added

Location:
palm/trunk/SOURCE
Files:
6 edited

Legend:

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

    r1552 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    8181        ONLY:  c_liq_av, c_soil_av, c_veg_av, ghf_eb_av, lai_av, m_liq_eb_av,  &
    8282               m_soil_av, nzb_soil, nzt_soil, qsws_eb_av, qsws_liq_eb_av,      &
    83                qsws_soil_eb_av, qsws_veg_eb_av, shf_eb_av, t_soil_av
     83               qsws_soil_eb_av, qsws_veg_eb_av, r_a, r_a_av, r_s, r_s_av,      &
     84               shf_eb_av, t_soil_av
    8485
    8586    USE radiation_model_mod,                                                   &
     
    359360             ENDDO
    360361
     362         CASE ( 'r_a*' )
     363             DO  i = nxlg, nxrg
     364                DO  j = nysg, nyng
     365                   r_a_av(j,i) = r_a_av(j,i) / REAL( average_count_3d, KIND=wp )
     366                ENDDO
     367             ENDDO
     368
     369         CASE ( 'r_s*' )
     370             DO  i = nxlg, nxrg
     371                DO  j = nysg, nyng
     372                   r_s_av(j,i) = r_s_av(j,i) / REAL( average_count_3d, KIND=wp )
     373                ENDDO
     374             ENDDO
     375
    361376          CASE ( 'rho' )
    362377             DO  i = nxlg, nxrg
  • palm/trunk/SOURCE/check_parameters.f90

    r1554 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s. Renumbering of LSM PA-messages.
    2323!
    2424! Former revisions:
     
    990990 
    991991          IF ( min_canopy_resistance == 9999999.9_wp)  THEN
    992              message_string = 'veg_type = 0 (user_defined)'//                  &
     992             message_string = 'veg_type = 0 (user defined)'//                  &
    993993                              'requires setting of min_canopy_resistance'//    &
    994994                              '/= 9999999.9'
    995              CALL message( 'check_parameters', 'PA0415', 1, 2, 0, 6, 0 )
     995             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    996996          ENDIF
    997997
     
    10001000                              'requires setting of leaf_area_index'//          &
    10011001                              '/= 9999999.9'
    1002              CALL message( 'check_parameters', 'PA0416', 1, 2, 0, 6, 0 )
     1002             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10031003          ENDIF
    10041004
     
    10071007                              'requires setting of vegetation_coverage'//      &
    10081008                              '/= 9999999.9'
    1009              CALL message( 'check_parameters', 'PA0417', 1, 2, 0, 6, 0 )
     1009             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10101010          ENDIF
    10111011
     
    10141014                              'requires setting of'//                          &
    10151015                              'canopy_resistance_coefficient /= 9999999.9'
    1016              CALL message( 'check_parameters', 'PA0418', 1, 2, 0, 6, 0 )
     1016             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10171017          ENDIF
    10181018
     
    10211021                              'requires setting of lambda_surface_stable'//    &
    10221022                              '/= 9999999.9'
    1023              CALL message( 'check_parameters', 'PA0419', 1, 2, 0, 6, 0 )
     1023             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10241024          ENDIF
    10251025
     
    10281028                              'requires setting of lambda_surface_unstable'//  &
    10291029                              '/= 9999999.9'
    1030              CALL message( 'check_parameters', 'PA0420', 1, 2, 0, 6, 0 )
     1030             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10311031          ENDIF
    10321032
     
    10351035                              'requires setting of f_shortwave_incoming'//     &
    10361036                              '/= 9999999.9'
    1037              CALL message( 'check_parameters', 'PA0421', 1, 2, 0, 6, 0 )
     1037             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10381038          ENDIF
    10391039
     
    10421042                              'requires setting of z0_eb'//                   &
    10431043                              '/= 9999999.9'
    1044              CALL message( 'check_parameters', 'PA0422', 1, 2, 0, 6, 0 )
     1044             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10451045          ENDIF
    10461046
     
    10491049                              'requires setting of z0h_eb'//                  &
    10501050                              '/= 9999999.9'
    1051              CALL message( 'check_parameters', 'PA0423', 1, 2, 0, 6, 0 )
     1051             CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    10521052          ENDIF
    10531053
     
    10611061                              'requires setting of alpha_vangenuchten'//       &
    10621062                              '/= 9999999.9'
    1063              CALL message( 'check_parameters', 'PA0422', 1, 2, 0, 6, 0 )
     1063             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10641064          ENDIF
    10651065
     
    10681068                              'requires setting of l_vangenuchten'//           &
    10691069                              '/= 9999999.9'
    1070              CALL message( 'check_parameters', 'PA0423', 1, 2, 0, 6, 0 )
     1070             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10711071          ENDIF
    10721072
     
    10751075                              'requires setting of n_vangenuchten'//           &
    10761076                              '/= 9999999.9'
    1077              CALL message( 'check_parameters', 'PA0424', 1, 2, 0, 6, 0 )
     1077             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10781078          ENDIF
    10791079
     
    10821082                              'requires setting of hydraulic_conductivity'//   &
    10831083                              '/= 9999999.9'
    1084              CALL message( 'check_parameters', 'PA0425', 1, 2, 0, 6, 0 )
     1084             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10851085          ENDIF
    10861086
     
    10891089                              'requires setting of saturation_moisture'//      &
    10901090                              '/= 9999999.9'
    1091              CALL message( 'check_parameters', 'PA0426', 1, 2, 0, 6, 0 )
     1091             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10921092          ENDIF
    10931093
     
    10961096                              'requires setting of field_capacity'//           &
    10971097                              '/= 9999999.9'
    1098              CALL message( 'check_parameters', 'PA0427', 1, 2, 0, 6, 0 )
     1098             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    10991099          ENDIF
    11001100
     
    11031103                              'requires setting of wilting_point'//            &
    11041104                              '/= 9999999.9'
    1105              CALL message( 'check_parameters', 'PA0428', 1, 2, 0, 6, 0 )
     1105             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    11061106          ENDIF
    11071107
     
    11101110                              'requires setting of residual_moisture'//        &
    11111111                              '/= 9999999.9'
    1112              CALL message( 'check_parameters', 'PA0429', 1, 2, 0, 6, 0 )
     1112             CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
    11131113          ENDIF
    11141114
     
    11181118          message_string = 'lsm requires '//                                   &
    11191119                           'radiation = .T.'
    1120           CALL message( 'check_parameters', 'PA0402', 1, 2, 0, 6, 0 )
     1120          CALL message( 'check_parameters', 'PA0400', 1, 2, 0, 6, 0 )
    11211121       ENDIF
    11221122
     
    11331133          message_string = 'unknown radiation_scheme = '//                     &
    11341134                           TRIM( radiation_scheme )
    1135           CALL message( 'check_parameters', 'PA0430', 1, 2, 0, 6, 0 )
     1135          CALL message( 'check_parameters', 'PA0405', 1, 2, 0, 6, 0 )
    11361136       ENDIF
    11371137    ENDIF
     
    30053005                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    30063006                                 'lemented for land_surface = .FALSE.'
    3007                 CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
     3007                CALL message( 'check_parameters', 'PA0402', 1, 2, 0, 6, 0 )
    30083008             ELSE
    30093009                dopr_index(i) = 89
     
    30223022                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    30233023                                 'lemented for land_surface = .FALSE.'
    3024                 CALL message( 'check_parameters', 'PA0403', 1, 2, 0, 6, 0 )
     3024                CALL message( 'check_parameters', 'PA0402', 1, 2, 0, 6, 0 )
    30253025             ELSE
    30263026                dopr_index(i) = 91
     
    32683268          CASE ( 'c_liq*', 'c_soil*', 'c_veg*', 'ghf_eb*', 'lai*', 'lwp*',     &
    32693269                 'm_liq_eb*', 'pra*', 'prr*', 'qsws*', 'qsws_eb*',             &
    3270                  'qsws_liq_eb*', 'qsws_soil_eb*', 'qsws_veg_eb*',              &
    3271                  'rad_net*', 'rad_sw_in*', 'shf*', 'shf_eb*', 't*', 'u*',      &
    3272                  'z0*', 'z0h*' )
     3270                 'qsws_liq_eb*', 'qsws_soil_eb*', 'qsws_veg_eb*', 'rad_net*',  &
     3271                 'rad_sw_in*', 'r_a*', 'r_s*', 'shf*', 'shf_eb*', 't*',        &
     3272                 'u*', 'z0*', 'z0h*' )
    32733273             IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
    32743274                message_string = 'illegal value for data_output: "' //         &
     
    32803280                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    32813281                                 'res land_surface = .TRUE.'
    3282                 CALL message( 'check_parameters', 'PA0411', 1, 2, 0, 6, 0 )
     3282                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    32833283             ENDIF
    32843284             IF ( TRIM( var ) == 'c_soil*'  .AND.  .NOT. land_surface )  THEN
    32853285                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    32863286                                 'res land_surface = .TRUE.'
    3287                 CALL message( 'check_parameters', 'PA0412', 1, 2, 0, 6, 0 )
     3287                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    32883288             ENDIF
    32893289             IF ( TRIM( var ) == 'c_veg*'  .AND.  .NOT. land_surface )  THEN
    32903290                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    32913291                                 'res land_surface = .TRUE.'
    3292                 CALL message( 'check_parameters', 'PA0413', 1, 2, 0, 6, 0 )
     3292                CALL message( 'check_parameters', 'PA0401', 1, 2, 0, 6, 0 )
    32933293             ENDIF
    32943294             IF ( TRIM( var ) == 'ghf_eb*'  .AND.  .NOT. land_surface )  THEN
    32953295                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    32963296                                 'res land_surface = .TRUE.'
    3297                 CALL message( 'check_parameters', 'PA0405', 1, 2, 0, 6, 0 )
     3297                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    32983298             ENDIF
    32993299             IF ( TRIM( var ) == 'lai*'  .AND.  .NOT. land_surface )  THEN
    33003300                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33013301                                 'res land_surface = .TRUE.'
    3302                 CALL message( 'check_parameters', 'PA0414', 1, 2, 0, 6, 0 )
     3302                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33033303             ENDIF
    33043304             IF ( TRIM( var ) == 'lwp*'  .AND.  .NOT. cloud_physics )  THEN
     
    33103310                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33113311                                 'res land_surface = .TRUE.'
    3312                 CALL message( 'check_parameters', 'PA0406', 1, 2, 0, 6, 0 )
     3312                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33133313             ENDIF
    33143314             IF ( TRIM( var ) == 'pra*'  .AND.  .NOT. precipitation )  THEN
     
    33353335                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33363336                                 'res land_surface = .TRUE.'
    3337                 CALL message( 'check_parameters', 'PA0407', 1, 2, 0, 6, 0 )
     3337                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33383338             ENDIF
    33393339             IF ( TRIM( var ) == 'qsws_liq_eb*'  .AND.  .NOT. land_surface )  &
     
    33413341                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33423342                                 'res land_surface = .TRUE.'
    3343                 CALL message( 'check_parameters', 'PA0408', 1, 2, 0, 6, 0 )
     3343                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33443344             ENDIF
    33453345             IF ( TRIM( var ) == 'qsws_soil_eb*'  .AND.  .NOT. land_surface ) &
     
    33473347                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33483348                                 'res land_surface = .TRUE.'
    3349                 CALL message( 'check_parameters', 'PA0409', 1, 2, 0, 6, 0 )
     3349                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33503350             ENDIF
    33513351             IF ( TRIM( var ) == 'qsws_veg_eb*'  .AND.  .NOT. land_surface )  &
     
    33533353                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    33543354                                 'res land_surface = .TRUE.'
    3355                 CALL message( 'check_parameters', 'PA0410', 1, 2, 0, 6, 0 )
     3355                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
     3356             ENDIF
     3357             IF ( TRIM( var ) == 'r_a*'  .AND.  .NOT. land_surface ) &
     3358             THEN
     3359                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     3360                                 'res land_surface = .TRUE.'
     3361                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
     3362             ENDIF
     3363             IF ( TRIM( var ) == 'r_s*'  .AND.  .NOT. land_surface ) &
     3364             THEN
     3365                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     3366                                 'res land_surface = .TRUE.'
     3367                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    33563368             ENDIF
    33573369
     
    33703382             IF ( TRIM( var ) == 'qsws_veg_eb*'  ) unit = 'W/m2'
    33713383             IF ( TRIM( var ) == 'rad_net*')       unit = 'W/m2'     
    3372              IF ( TRIM( var ) == 'rad_sw_in*')     unit = 'W/m2'   
     3384             IF ( TRIM( var ) == 'rad_sw_in*')     unit = 'W/m2'
     3385             IF ( TRIM( var ) == 'r_a*')     unit = 's/m'     
     3386             IF ( TRIM( var ) == 'r_s*')     unit = 's/m'
    33733387             IF ( TRIM( var ) == 'shf*'   )  unit = 'K*m/s'
    33743388             IF ( TRIM( var ) == 'shf_eb*')  unit = 'W/m2'
  • palm/trunk/SOURCE/data_output_2d.f90

    r1552 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    139139               m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
    140140               qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
    141                qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
    142                t_soil_av, zs
     141               qsws_veg_eb, qsws_veg_eb_av, r_a, r_a_av, r_s, r_s_av, shf_eb,  &
     142               shf_eb_av, t_soil, t_soil_av, zs
    143143   
    144144    USE netcdf_control
     
    934934                   to_be_resorted => rho_av
    935935                ENDIF
     936
     937             CASE ( 'r_a*_xy' )        ! 2d-array
     938                IF ( av == 0 )  THEN
     939                   DO  i = nxlg, nxrg
     940                      DO  j = nysg, nyng
     941                         local_pf(i,j,nzb+1) = r_a(j,i)
     942                      ENDDO
     943                   ENDDO
     944                ELSE
     945                   DO  i = nxlg, nxrg
     946                      DO  j = nysg, nyng
     947                         local_pf(i,j,nzb+1) = r_a_av(j,i)
     948                      ENDDO
     949                   ENDDO
     950                ENDIF
     951                resorted = .TRUE.
     952                two_d = .TRUE.
     953                level_z(nzb+1) = zu(nzb+1)
     954
     955             CASE ( 'r_s*_xy' )        ! 2d-array
     956                IF ( av == 0 )  THEN
     957                   DO  i = nxlg, nxrg
     958                      DO  j = nysg, nyng
     959                         local_pf(i,j,nzb+1) = r_s(j,i)
     960                      ENDDO
     961                   ENDDO
     962                ELSE
     963                   DO  i = nxlg, nxrg
     964                      DO  j = nysg, nyng
     965                         local_pf(i,j,nzb+1) = r_s_av(j,i)
     966                      ENDDO
     967                   ENDDO
     968                ENDIF
     969                resorted = .TRUE.
     970                two_d = .TRUE.
     971                level_z(nzb+1) = zu(nzb+1)
    936972
    937973             CASE ( 's_xy', 's_xz', 's_yz' )
  • palm/trunk/SOURCE/flow_statistics.f90

    r1552 r1555  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! Added output of r_a and r_s.
    2424!
    2525! Former revisions:
     
    160160    USE land_surface_model_mod,                                                &
    161161        ONLY:   dots_soil, ghf_eb, land_surface, m_soil, nzb_soil, nzt_soil,   &
    162                 qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, shf_eb,       &
    163                 t_soil
     162                qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, r_a, r_s,     &
     163                shf_eb, t_soil
    164164
    165165    USE pegrid
     
    717717
    718718             IF ( land_surface )  THEN
    719                 sums_l(nzb,93,tn) = sums_l(nzb,93,tn) + ghf_eb(j,i)
    720                 sums_l(nzb,94,tn) = sums_l(nzb,94,tn) + shf_eb(j,i)
    721                 sums_l(nzb,95,tn) = sums_l(nzb,95,tn) + qsws_eb(j,i)
    722                 sums_l(nzb,96,tn) = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
    723                 sums_l(nzb,97,tn) = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
    724                 sums_l(nzb,98,tn) = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
     719                sums_l(nzb,93,tn)  = sums_l(nzb,93,tn) + ghf_eb(j,i)
     720                sums_l(nzb,94,tn)  = sums_l(nzb,94,tn) + shf_eb(j,i)
     721                sums_l(nzb,95,tn)  = sums_l(nzb,95,tn) + qsws_eb(j,i)
     722                sums_l(nzb,96,tn)  = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
     723                sums_l(nzb,97,tn)  = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
     724                sums_l(nzb,98,tn)  = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
     725                sums_l(nzb,99,tn)  = sums_l(nzb,99,tn) + r_a(j,i)
     726                sums_l(nzb,100,tn) = sums_l(nzb,100,tn)+ r_s(j,i)
    725727             ENDIF
    726728
    727729             IF ( radiation )  THEN
    728                 sums_l(nzb,99,tn)  = sums_l(nzb,99,tn)  + rad_net(j,i)
    729                 sums_l(nzb,100,tn) = sums_l(nzb,100,tn) + rad_sw_in(j,i)
     730                sums_l(nzb,101,tn)  = sums_l(nzb,101,tn)  + rad_net(j,i)
     731                sums_l(nzb,102,tn) = sums_l(nzb,102,tn) + rad_sw_in(j,i)
    730732             ENDIF
    731733
     
    12931295          hom(:,1,91,sr) = sums(:,91)              ! m_soil
    12941296                                                   ! 92 is initial m_soil profile
    1295           hom(:,1,93,sr) = sums(:,93)              ! ghf_eb
    1296           hom(:,1,94,sr) = sums(:,94)              ! shf_eb
    1297           hom(:,1,95,sr) = sums(:,95)              ! qsws_eb
    1298           hom(:,1,96,sr) = sums(:,96)              ! qsws_liq_eb
    1299           hom(:,1,97,sr) = sums(:,97)              ! qsws_soil_eb
    1300           hom(:,1,98,sr) = sums(:,98)              ! qsws_veg_eb
     1297          hom(:,1,93,sr)  = sums(:,93)             ! ghf_eb
     1298          hom(:,1,94,sr)  = sums(:,94)             ! shf_eb
     1299          hom(:,1,95,sr)  = sums(:,95)             ! qsws_eb
     1300          hom(:,1,96,sr)  = sums(:,96)             ! qsws_liq_eb
     1301          hom(:,1,97,sr)  = sums(:,97)             ! qsws_soil_eb
     1302          hom(:,1,98,sr)  = sums(:,98)             ! qsws_veg_eb
     1303          hom(:,1,99,sr)  = sums(:,99)             ! r_a
     1304          hom(:,1,100,sr) = sums(:,100)            ! r_s
     1305
    13011306       ENDIF
    13021307
    13031308       IF ( radiation )  THEN
    1304           hom(:,1,99 ,sr) = sums(:,99)             ! rad_net
    1305           hom(:,1,100,sr) = sums(:,100)            ! rad_sw_in
     1309          hom(:,1,101 ,sr) = sums(:,101)            ! rad_net
     1310          hom(:,1,102,sr)  = sums(:,102)            ! rad_sw_in
    13061311       ENDIF
    13071312
     
    14601465          ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
    14611466          ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
     1467          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
     1468          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
    14621469       ENDIF
    14631470!
     
    32633270          ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
    32643271          ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
     3272          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
     3273          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
    32653274       ENDIF
    32663275!
  • palm/trunk/SOURCE/land_surface_model.f90

    r1554 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    236236              qsws_veg_eb_av,   & !: average of qsws_veg_eb
    237237              r_a,              & !: aerodynamic resistance
     238              r_a_av,           & !: avergae of r_a
    238239              r_canopy,         & !: canopy resistance
    239240              r_soil,           & !: soil resitance
    240241              r_soil_min,       & !: minimum soil resistance
    241               r_s,              & !: total surface resistance (combination of r_soil and r_canopy)         
     242              r_s,              & !: total surface resistance (combination of r_soil and r_canopy)
     243              r_s_av,           & !: avergae of r_s
    242244              r_canopy_min,     & !: minimum canopy (stomatal) resistance
    243245              shf_eb,           & !: surface flux of sensible heat
     
    474476           lai, lai_av, qsws_eb, qsws_eb_av, qsws_liq_eb, qsws_liq_eb_av,      &
    475477           qsws_soil_eb, qsws_soil_eb_av, qsws_veg_eb, qsws_veg_eb_av,         &
    476            shf_eb, shf_eb_av
     478           r_a, r_a_av, r_s, r_s_av, shf_eb, shf_eb_av
    477479
    478480
     
    888890       dots_unit(dots_num+1:dots_num+6) = "W/m2"
    889891
     892       dots_label(dots_num+7) = "r_a"
     893       dots_label(dots_num+8) = "r_s"
     894       dots_unit(dots_num+7:dots_num+8) = "s/m"
     895
    890896       dots_soil = dots_num + 1
    891        dots_num  = dots_num + 6
     897       dots_num  = dots_num + 8
    892898
    893899
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r1552 r1555  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Added output of r_a and r_s
    2323!
    2424! Former revisions:
     
    102102               m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
    103103               qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
    104                qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
    105                t_soil_av
     104               qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, r_a, r_a_av,    &
     105               r_s, r_s_av, t_soil, t_soil_av
    106106
    107107    USE particle_attributes,                                                   &
     
    334334                ENDIF
    335335                rho_av = 0.0_wp
     336
     337             CASE ( 'r_a*' )
     338                IF ( .NOT. ALLOCATED( r_a_av ) )  THEN
     339                   ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
     340                ENDIF
     341                r_a_av = 0.0_wp
     342
     343             CASE ( 'r_s*' )
     344                IF ( .NOT. ALLOCATED( r_s_av ) )  THEN
     345                   ALLOCATE( r_s_av(nysg:nyng,nxlg:nxrg) )
     346                ENDIF
     347                r_s_av = 0.0_wp
    336348
    337349             CASE ( 's' )
     
    724736             ENDDO
    725737
     738          CASE ( 'r_a*' )
     739             DO  i = nxlg, nxrg
     740                DO  j = nysg, nyng
     741                   r_a_av(j,i) = r_a_av(j,i) + r_a(j,i)
     742                ENDDO
     743             ENDDO
     744
     745          CASE ( 'r_s*' )
     746             DO  i = nxlg, nxrg
     747                DO  j = nysg, nyng
     748                   r_s_av(j,i) = r_s_av(j,i) + r_s(j,i)
     749                ENDDO
     750             ENDDO
     751
    726752          CASE ( 'rho' )
    727753             DO  i = nxlg, nxrg
Note: See TracChangeset for help on using the changeset viewer.