Changeset 4347 for palm/trunk/SOURCE

Timestamp:

Dec 18, 2019 1:18:33 PM (4 years ago)

Author:

suehring

Message:

Check for realistically initial values of mixing ratio implemented. Mixing ratio must not exceed its saturation value, else surface fluxes in the land-surface scheme may become unrealistic. Test cases updated

Location:

palm/trunk/SOURCE

Files:

• Makefile (modified) (2 diffs)
• dynamics_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/Makefile

@@ -25 +25 @@
 # -----------------
 # $Id$
+# add dependency to basic_constants_and_equations_mod for dynamics_mod
+# 
+# 4331 2019-12-10 18:25:02Z suehring
 # Move diagnostic surface output to diagnostic_output_quantities
 # 
@@ -553 +556 @@
         surface_mod.o
 dynamics_mod.o: \
+        basic_constants_and_equations_mod.o \
         mod_kinds.o \
         surface_mod.o \

palm/trunk/SOURCE/dynamics_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implement post-initialization check for realistically inital values of mixing ratio
+! 
+! 4281 2019-10-29 15:15:39Z schwenkel
 ! Moved boundary conditions in dynamics module
 ! 
@@ -44 +47 @@
         ONLY:  c_u, c_u_m, c_u_m_l, c_v, c_v_m, c_v_m_l, c_w, c_w_m, c_w_m_l,  &
                dzu, &
+               exner, &
+               hyp, &
                pt, pt_1, pt_2, pt_init, pt_p, &
                q, q_1, q_2, q_p, &
@@ -50 +55 @@
                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
+
+    USE basic_constants_and_equations_mod,                                                         &
+        ONLY:  magnus,                                                                             &
+               rd_d_rv
 
     USE control_parameters, &
@@ -74 +83 @@
                intermediate_timestep_count, &
                length, &
+               message_string, &
                nesting_offline, &
                nudging, &
@@ -455 +465 @@
  SUBROUTINE dynamics_init_checks
 
+    INTEGER(iwp) ::  i !< loop index in x-direction
+    INTEGER(iwp) ::  j !< loop index in y-direction
+    INTEGER(iwp) ::  k !< loop index in z-direction
+
+    LOGICAL      ::  realistic_q = .TRUE. !< flag indicating realistic mixing ratios
+
+    REAL(wp)     ::  e_s !< saturation water vapor pressure
+    REAL(wp)     ::  q_s !< saturation mixing ratio
+    REAL(wp)     ::  t_l !< actual temperature
+
+!
+!-- Check for realistic initial mixing ratio. This must be in a realistic phyiscial range and must 
+!-- 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.
+    IF ( humidity  .AND.  .NOT. neutral )  THEN
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             DO  k = nzb+1, nzt
+!
+!--             Calculate actual temperature, water vapor saturation pressure, and based on this
+!--             the saturation mixing ratio.
+                t_l = exner(k) * pt(k,j,i)
+                e_s = magnus( t_l )
+                q_s = rd_d_rv * e_s / ( hyp(k) - e_s )
+
+                IF ( q(k,j,i) > 1.02_wp * q_s )  realistic_q = .FALSE. 
+             ENDDO
+          ENDDO
+       ENDDO
+!
+!--    Since the check is performed locally, merge the logical flag from all mpi ranks, 
+!--    in order to do not print the error message multiple times.
+#if defined( __parallel )
+       CALL MPI_ALLREDUCE( MPI_IN_PLACE, realistic_q, 1, MPI_LOGICAL, MPI_LAND, comm2d, ierr)
+#endif
+
+       IF ( .NOT. realistic_q )  THEN
+          message_string = 'The initial mixing ratio exceeds the saturation mixing ratio.'
+          CALL message( 'dynamic_init_checks', 'PA0697', 2, 2, 0, 6, 0 )
+       ENDIF
+    ENDIF
 
  END SUBROUTINE dynamics_init_checks