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

Timestamp:

Oct 12, 2016 4:42:37 PM (8 years ago)

Author:

kanani

Message:

changes related to urban surface model and output of ssws

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added missing CASE, error message and unit for ssws*,
+! increased number of possible output quantities to 500.
 ! 
 ! Former revisions:
@@ -2834 +2835 @@
     IF ( data_output_user(1) /= ' ' )  THEN
        i = 1
-       DO  WHILE ( data_output(i) /= ' '  .AND.  i <= 100 )
+       DO  WHILE ( data_output(i) /= ' '  .AND.  i <= 500 )
           i = i + 1
        ENDDO
        j = 1
-       DO  WHILE ( data_output_user(j) /= ' '  .AND.  j <= 100 )
-          IF ( i > 100 )  THEN
+       DO  WHILE ( data_output_user(j) /= ' '  .AND.  j <= 500 )
+          IF ( i > 500 )  THEN
              message_string = 'number of output quantitities given by data' // &
-                '_output and data_output_user exceeds the limit of 100'
+                '_output and data_output_user exceeds the limit of 500'
              CALL message( 'check_parameters', 'PA0102', 1, 2, 0, 6, 0 )
           ENDIF
@@ -3018 +3019 @@
              CONTINUE
 
-          CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'shf*', 't*', &
+          CASE ( 'lwp*', 'ol*', 'pra*', 'prr*', 'qsws*', 'shf*', 'ssws*', 't*', &
                  'u*', 'z0*', 'z0h*', 'z0q*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
@@ -3054 +3055 @@
                 CALL message( 'check_parameters', 'PA0322', 1, 2, 0, 6, 0 )
              ENDIF
+             IF ( TRIM( var ) == 'ssws*'  .AND.  .NOT.  passive_scalar )  THEN
+                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                                 'res passive_scalar = .TRUE.'
+                CALL message( 'check_parameters', 'PA0361', 1, 2, 0, 6, 0 )
+             ENDIF             
 
              IF ( TRIM( var ) == 'lwp*'   )  unit = 'kg/m2'
@@ -3061 +3067 @@
              IF ( TRIM( var ) == 'qsws*'  )  unit = 'kgm/kgs'
              IF ( TRIM( var ) == 'shf*'   )  unit = 'K*m/s'
+             IF ( TRIM( var ) == 'ssws*'  )  unit = 'kg/m2*s'      
              IF ( TRIM( var ) == 't*'     )  unit = 'K'
              IF ( TRIM( var ) == 'u*'     )  unit = 'm/s'