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

Timestamp:

Jul 13, 2009 1:14:25 PM (15 years ago)

Author:

maronga

Message:

xy cross section output of sensible and latent heatflux is now available

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -2365 +2365 @@
              unit = 'psu'
 
-          CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*' )
+          CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'qsws*', 'shf*', 'z0*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
                 message_string = 'illegal value for data_output: "' // &
@@ -2392 +2392 @@
                 CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
              ENDIF
-
-
-             IF ( TRIM( var ) == 'u*'   )  unit = 'm/s'
-             IF ( TRIM( var ) == 't*'   )  unit = 'K'
-             IF ( TRIM( var ) == 'lwp*' )  unit = 'kg/kg*m'
-             IF ( TRIM( var ) == 'pra*' )  unit = 'mm'
-             IF ( TRIM( var ) == 'prr*' )  unit = 'mm/s'
-             IF ( TRIM( var ) == 'z0*'  )  unit = 'm'
+             IF ( TRIM( var ) == 'qsws*'  .AND.  .NOT. humidity )  THEN
+                message_string = 'output of "' // TRIM( var ) // '" requi' // &
+                                 'res humidity = .TRUE.'
+                CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
+             ENDIF
+
+             IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/kg*m'
+             IF ( TRIM( var ) == 'pra*'   )  unit = 'mm'
+             IF ( TRIM( var ) == 'prr*'   )  unit = 'mm/s'
+             IF ( TRIM( var ) == 'qsws*'  )  unit = 'kgm/kgs'
+             IF ( TRIM( var ) == 'shf*'   )  unit = 'K*m/s'
+             IF ( TRIM( var ) == 't*'     )  unit = 'K'
+             IF ( TRIM( var ) == 'u*'     )  unit = 'm/s'
+             IF ( TRIM( var ) == 'z0*'    )  unit = 'm'
+
 
           CASE ( 'p', 'pt', 'u', 'v', 'w' )