Changeset 2839

Timestamp:

Feb 27, 2018 9:49:06 AM (7 years ago)

Author:

schwenkel

Message:

Bugfix for Kessler microphysics

File:

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

palm/trunk/SOURCE/diagnostic_quantities_mod.f90

@@ -25 +25 @@
 ! -----------------! 
 ! $Id$
+! Bugfix for Kessler microphysics
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -74 +77 @@
 
        USE cloud_parameters,                                                   &
-           ONLY:  l_d_cp, l_d_r, t_d_pt                                              
+           ONLY:  l_d_cp, l_d_r, t_d_pt       
+
+       USE control_parameters,                                                 &
+           ONLY:  microphysics_kessler                                       
 
        IMPLICIT NONE
@@ -97 +103 @@
 !--    (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:
-       sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
+!--    Not in case of microphysics_kessler since qr is unallocated
+       IF ( .NOT. microphysics_kessler ) THEN
+          sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
+       ENDIF
 
     END SUBROUTINE supersaturation