Ignore:
Timestamp:
Jan 11, 2018 12:01:27 PM (6 years ago)
Author:
suehring
Message:

Output of resistance also urban-type surfaces

File:
1 edited

Legend:

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

    r2718 r2735  
    2525! -----------------
    2626! $Id$
     27! output of r_a moved from land-surface to consider also urban-type surfaces
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    32823285             CONTINUE
    32833286
    3284           CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'shf*', 'ssws*', 't*', &
    3285                  'u*', 'z0*', 'z0h*', 'z0q*' )
     3287          CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'r_a*', 'shf*',      &
     3288                 'ssws*', 't*', 'u*', 'z0*', 'z0h*', 'z0q*' )
    32863289             IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
    32873290                message_string = 'illegal value for data_output: "' //         &
     
    33183321                CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
    33193322             ENDIF
     3323             IF ( TRIM( var ) == 'r_a*'  .AND.  .NOT.  land_surface  .AND.     &
     3324                                                .NOT.  urban_surface )         &         
     3325             THEN
     3326                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     3327                                 'res land_surface = .TRUE. or ' //            &
     3328                                 'urban_surface = .TRUE.'
     3329                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
     3330             ENDIF
    33203331             IF ( TRIM( var ) == 'ssws*'  .AND.  .NOT.  passive_scalar )  THEN
    33213332                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     
    33253336
    33263337             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/m2'
    3327              IF ( TRIM( var ) == 'ol*'   )   unit = 'm'
     3338             IF ( TRIM( var ) == 'ol*'    )  unit = 'm'
    33283339             IF ( TRIM( var ) == 'pra*'   )  unit = 'mm'
    33293340             IF ( TRIM( var ) == 'prr*'   )  unit = 'mm/s'
    33303341             IF ( TRIM( var ) == 'qsws*'  )  unit = 'kgm/kgs'
     3342             IF ( TRIM( var ) == 'r_a*'   )  unit = 's/m'     
    33313343             IF ( TRIM( var ) == 'shf*'   )  unit = 'K*m/s'
    33323344             IF ( TRIM( var ) == 'ssws*'  )  unit = 'kg/m2*s'
Note: See TracChangeset for help on using the changeset viewer.