Changeset 4358 for palm/trunk/SOURCE/dynamics_mod.f90

Timestamp:

Dec 30, 2019 1:03:09 PM (5 years ago)

Author:

suehring

Message:

Bugfix, time coordinate is relative to origin_time rather than to 00:00:00 UTC; Refine post-initialization check for realistically inital values of mixing ratio. Give an error message for faulty initial values, but only a warning in a restart run.

File:

• palm/trunk/SOURCE/dynamics_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/dynamics_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Refine post-initialization check for realistically inital values of mixing ratio. Give an error
+! message for faulty initial values, but only a warning in a restart run. 
+! 
+! 4347 2019-12-18 13:18:33Z suehring
 ! Implement post-initialization check for realistically inital values of mixing ratio
 ! 
@@ -479 +483 @@
 !-- not exceed the saturation mixing ratio by more than 2 percent. Please note, the check is 
 !-- performed for each grid point (not just for a vertical profile), in order to cover also 
-!-- three-dimensional initialization.
+!-- three-dimensional initialization. Note, this check gives an error only for the initial run not
+!-- for a restart run. In case there are no cloud physics considered, the mixing ratio can exceed
+!-- the saturation moisture. This case a warning is given. 
     IF ( humidity  .AND.  .NOT. neutral )  THEN
        DO  i = nxl, nxr
@@ -502 +508 @@
 #endif
 
-       IF ( .NOT. realistic_q )  THEN
+       IF ( .NOT. realistic_q  .AND.                                                               &
+            TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
           message_string = 'The initial mixing ratio exceeds the saturation mixing ratio.'
           CALL message( 'dynamic_init_checks', 'PA0697', 2, 2, 0, 6, 0 )
+       ELSEIF ( .NOT. realistic_q  .AND.                                                           &
+                TRIM( initializing_actions ) == 'read_restart_data' )  THEN
+          message_string = 'The mixing ratio exceeds the saturation mixing ratio.'
+          CALL message( 'dynamic_init_checks', 'PA0697', 0, 1, 0, 6, 0 )
        ENDIF
     ENDIF