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

Timestamp:

Oct 26, 2020 10:23:38 AM (4 years ago)

Author:

schwenkel

Message:

implement relative humidity as diagnostic output quantity

File:

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

palm/trunk/SOURCE/dynamics_mod.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Implement relative humidity as diagnostic output quantity
+!
+! 4731 2020-10-07 13:25:11Z schwenkel
 ! Move exchange_horiz from time_integration to modules
 !
@@ -82 +85 @@
                u, u_1, u_2, u_init, u_p, u_m_l, u_m_n, u_m_r, u_m_s,                               &
                v, v_1, v_2, v_p, v_init, v_m_l, v_m_n, v_m_r, v_m_s,                               &
-               w, w_1, w_2, w_p, w_m_l, w_m_n, w_m_r, w_m_s
+               w, w_1, w_2, w_p, w_m_l, w_m_n, w_m_r, w_m_s,                                       &
+               zu
 
     USE basic_constants_and_equations_mod,                                                         &
@@ -518 +522 @@
     REAL(wp)     ::  q_s !< saturation mixing ratio
     REAL(wp)     ::  t_l !< actual temperature
+    REAL(wp)     ::  rh     = 9999999.9_wp !< relative humidity
+    REAL(wp)     ::  rh_min = 9999999.9_wp !< max relative humidity
+    REAL(wp)     ::  height = 9999999.9_wp !< height of supersaturated regions
+    REAL(wp)     ::  min_height = 9999999.9_wp !< height of supersaturated regions
 
 !
@@ -537 +545 @@
                 q_s = rd_d_rv * e_s / ( hyp(k) - e_s )
 
-                IF ( q(k,j,i) > 1.02_wp * q_s )  realistic_q = .FALSE.
+                IF ( q(k,j,i) > 1.02_wp * q_s )  THEN
+                   realistic_q = .FALSE.
+                   rh = q(k,j,i) / q_s * 100.0_wp
+                   height = zu(k)
+                ENDIF
              ENDDO
           ENDDO
@@ -546 +558 @@
 #if defined( __parallel )
        CALL MPI_ALLREDUCE( MPI_IN_PLACE, realistic_q, 1, MPI_LOGICAL, MPI_LAND, comm2d, ierr)
+       CALL MPI_ALLREDUCE( rh, rh_min, 1, MPI_REAL, MPI_MIN, comm2d, ierr )
+       CALL MPI_ALLREDUCE( height, min_height, 1, MPI_REAL, MPI_MIN, comm2d, ierr )
 #endif
 
        IF ( .NOT. realistic_q  .AND.                                                               &
             TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
-          message_string = 'The initial mixing ratio exceeds the saturation mixing ratio.'
+            WRITE( message_string, * ) 'The initial mixing ratio exceeds the saturation mixing' // &
+               'ratio, with rh =', rh_min, '% in a height of', min_height, 'm for the first time'
           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.'
+          WRITE( message_string, * ) 'The initial mixing ratio exceeds the saturation mixing' //   &
+              'ratio, with rh =', rh_min, '% in a height of', min_height, 'm for the first time'
           CALL message( 'dynamic_init_checks', 'PA0697', 0, 1, 0, 6, 0 )
        ENDIF