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

Timestamp:

Mar 19, 2007 8:20:46 AM (17 years ago)

Author:

raasch

Message:

preliminary changes for precipitation output

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -6 +6 @@
 ! "by_user" allowed as initializing action, -data_output_ts,
 ! leapfrog with non-flat topography not allowed any more, loop_optimization
-! and pt_reference are checked
+! and pt_reference are checked,
+! output of precipitation amount/rate and roughnes length + check
+! possible negative humidities are avoided in initial profile
 !
 ! Former revisions:
@@ -351 +353 @@
 
     IF ( cloud_physics  .AND.  .NOT. moisture )  THEN
-       IF ( myid == 0 )  PRINT*, '+++ check_parameters: moisture =', &
-                                 moisture, ' is not allowed with ',  &
-                                 'cloud_physics=', cloud_physics
+       IF ( myid == 0 )  PRINT*, '+++ check_parameters: cloud_physics =', &
+                                 cloud_physics, ' is not allowed with ',  &
+                                 'moisture =', moisture
+       CALL local_stop
+    ENDIF
+
+    IF ( precipitation  .AND.  .NOT.  cloud_physics )  THEN
+       IF ( myid == 0 )  PRINT*, '+++ check_parameters: precipitation =', &
+                                 precipitation, ' is not allowed with ',  &
+                                 'cloud_physics =', cloud_physics
        CALL local_stop
     ENDIF
@@ -570 +579 @@
                 q_init(k) = q_init(k-1)
              ENDIF
+!
+!--          Avoid negative humidities
+             IF ( q_init(k) < 0.0 )  THEN
+                q_init(k) = 0.0
+             ENDIF
           ENDDO
 
@@ -1163 +1177 @@
 
 !
+!-- Set the default value for the integration interval of precipitation amount
+    IF ( precipitation )  THEN
+       IF ( precipitation_amount_interval == 9999999.9 )  THEN
+          precipitation_amount_interval = dt_do2d_xy
+       ELSE
+          IF ( precipitation_amount_interval > dt_do2d_xy )  THEN
+             IF ( myid == 0 )  PRINT*, '+++ check_parameters: ',              &
+                                       'precipitation_amount_interval =',     &
+                                        precipitation_amount_interval,        &
+                                       ' must not be larger than dt_do2d_xy', &
+                                       ' = ', dt_do2d_xy   
+       CALL local_stop
+          ENDIF
+       ENDIF
+    ENDIF
+
+!
 !-- Determine the number of output profiles and check whether they are
 !-- permissible
@@ -1840 +1871 @@
              unit = 'conc'
 
-          CASE ( 'u*', 't*', 'lwp*' )
+          CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
                 IF ( myid == 0 )  THEN
@@ -1856 +1887 @@
                 CALL local_stop
              ENDIF
+             IF ( TRIM( var ) == 'pra*'  .AND.  .NOT. precipitation )  THEN
+                IF ( myid == 0 )  THEN
+                   PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
+                                '" requires precipitation = .TRUE.'
+                ENDIF
+                CALL local_stop
+             ENDIF
+             IF ( TRIM( var ) == 'pra*'  .AND.  j == 1 )  THEN
+                IF ( myid == 0 )  THEN
+                   PRINT*, '+++ check_parameters: temporal averaging of ', &
+                           ' precipitation amount "', TRIM( var ),         &
+                           '" not possible' 
+                ENDIF
+                CALL local_stop
+             ENDIF
+             IF ( TRIM( var ) == 'prr*'  .AND.  .NOT. precipitation )  THEN
+                IF ( myid == 0 )  THEN
+                   PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
+                                '" requires precipitation = .TRUE.'
+                ENDIF
+                CALL local_stop
+             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'
 
           CASE ( 'p', 'pt', 'u', 'v', 'w' )