Ignore:
Timestamp:
Apr 25, 2019 12:33:20 PM (2 years ago)
Author:
kanani
Message:

Bugfixes and clean-up for output quantity theta_2m*

File:
1 edited

Legend:

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

    r3885 r3933  
    2525! -----------------
    2626! $Id$
     27! Alphabetical resorting in CASE, condense settings for theta_2m* into one IF clause
     28!
     29! 3885 2019-04-11 11:29:34Z kanani
    2730! Changes related to global restructuring of location messages and introduction
    2831! of additional debug messages
     
    30773080             CONTINUE
    30783081
    3079           CASE ( 'ghf*', 'lwp*', 'ol*', 'qsws*', 'r_a*', 'theta_2m*',          &
    3080                  'shf*', 'ssws*', 't*', 'tsurf*', 'us*', 'z0*', 'z0h*', 'z0q*' )
     3082          CASE ( 'ghf*', 'lwp*', 'ol*', 'qsws*', 'r_a*',                       &
     3083                 'shf*', 'ssws*', 't*', 'theta_2m*', 'tsurf*', 'us*',          &
     3084                 'z0*', 'z0h*', 'z0q*' )
    30813085             IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
    30823086                message_string = 'illegal value for data_output: "' //         &
     
    31173121                CALL message( 'check_parameters', 'PA0361', 1, 2, 0, 6, 0 )
    31183122             ENDIF
    3119 
    3120            
    31213123!
    31223124!--          Activate calculation of 2m temperature if output is requested
    31233125             IF ( TRIM( var ) == 'theta_2m*' )  THEN
    31243126                do_output_at_2m = .TRUE.
    3125              ENDIF             
    3126 
     3127                unit = 'K'
     3128             ENDIF
    31273129
    31283130             IF ( TRIM( var ) == 'ghf*'   )  unit = 'W/m2'
    31293131             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/m2'
    31303132             IF ( TRIM( var ) == 'ol*'    )  unit = 'm'
    3131              IF ( TRIM( var ) == 'theta_2m*' )  unit = 'K'           
    31323133             IF ( TRIM( var ) == 'qsws*'  )  unit = 'kgm/kgs'
    3133              IF ( TRIM( var ) == 'r_a*'   )  unit = 's/m'     
     3134             IF ( TRIM( var ) == 'r_a*'   )  unit = 's/m'
    31343135             IF ( TRIM( var ) == 'shf*'   )  unit = 'K*m/s'
    31353136             IF ( TRIM( var ) == 'ssws*'  )  unit = 'kg/m2*s'
    31363137             IF ( TRIM( var ) == 't*'     )  unit = 'K'
    3137              IF ( TRIM( var ) == 'tsurf*' )  unit = 'K' 
     3138             IF ( TRIM( var ) == 'tsurf*' )  unit = 'K'
    31383139             IF ( TRIM( var ) == 'us*'    )  unit = 'm/s'
    31393140             IF ( TRIM( var ) == 'z0*'    )  unit = 'm'
Note: See TracChangeset for help on using the changeset viewer.