Home


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

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

LSM output of r_a and r_s added

File:
1 edited

Legend:

Unmodified
Added
Removed

  • 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'
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |