Changeset 2742 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Jan 12, 2018 2:59:47 PM (6 years ago)

Author:

suehring

Message:

Output of surface temperature tsurf* at urban- and natural-type surfaces

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable output of surface temperature
+! 
+! 2735 2018-01-11 12:01:27Z suehring
 ! output of r_a moved from land-surface to consider also urban-type surfaces
 ! 
@@ -3286 +3289 @@
 
           CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'r_a*', 'shf*',       &
-                 'ssws*', 't*', 'u*', 'z0*', 'z0h*', 'z0q*' )
+                 'ssws*', 't*', 'tsurf*', 'u*', 'z0*', 'z0h*', 'z0q*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
                 message_string = 'illegal value for data_output: "' //         &
@@ -3321 +3324 @@
                 CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
              ENDIF
-             IF ( TRIM( var ) == 'r_a*'  .AND.  .NOT.  land_surface  .AND.     &
-                                                .NOT.  urban_surface )         &         
+             IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'tsurf*' )      &
+                 .AND.  .NOT.  land_surface  .AND.  .NOT.  urban_surface )     &         
              THEN
                 message_string = 'output of "' // TRIM( var ) // '" requi' //  &
@@ -3344 +3347 @@
              IF ( TRIM( var ) == 'ssws*'  )  unit = 'kg/m2*s'
              IF ( TRIM( var ) == 't*'     )  unit = 'K'
+             IF ( TRIM( var ) == 'tsurf*' )  unit = 'K'  
              IF ( TRIM( var ) == 'u*'     )  unit = 'm/s'
              IF ( TRIM( var ) == 'z0*'    )  unit = 'm'