Changeset 3231 for palm/trunk

Timestamp:

Sep 6, 2018 11:53:51 AM (6 years ago)

Author:

schwenkel

Message:

Bugfix for microphysic scheme = saturation_adjustment

File:

• palm/trunk/SOURCE/diagnostic_quantities_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/diagnostic_quantities_mod.f90

@@ -25 +25 @@
 ! ----------------- 
 ! $Id$
+! Bugfix for microphysic scheme = 'saturation_adjustment'
+! 
+! 3026 2018-05-22 10:30:53Z schwenkel
 ! Changed the name specific humidity to mixing ratio, since we are computing
 ! mixing ratios.
@@ -72 +75 @@
 ! Description:
 ! ------------
-!> Computation of the diagnostic supersaturation sat, actual temperature t_l
-!< and saturation water vapor mixing ratio q_
+!> Computation of the diagnostic supersaturation sat, actual liquid water 
+!< temperature t_l and saturation water vapor mixing ratio q_s
 !------------------------------------------------------------------------------!
     SUBROUTINE supersaturation ( i,j,k )
@@ -84 +87 @@
 
        USE control_parameters,                                                 &
-           ONLY:  microphysics_kessler                                       
+           ONLY:  microphysics_kessler, microphysics_sat_adjust                                       
 
        IMPLICIT NONE
@@ -107 +110 @@
 !--    (see: Cuijpers + Duynkerke, 1993, JAS, 23)
        q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
-
 !
 !--    Supersaturation:
-!--    Not in case of microphysics_kessler since qr is unallocated
-       IF ( .NOT. microphysics_kessler ) THEN
+!--    Not in case of microphysics_kessler or microphysics_sat_adjust 
+!--    since qr is unallocated
+       IF ( .NOT. microphysics_kessler  .AND.                                  & 
+            .NOT.  microphysics_sat_adjust )                                THEN
           sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
        ENDIF