Changeset 2839 for palm/trunk


Ignore:
Timestamp:
Feb 27, 2018 9:49:06 AM (4 years ago)
Author:
schwenkel
Message:

Bugfix for Kessler microphysics

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/diagnostic_quantities_mod.f90

    r2718 r2839  
    2525! -----------------!
    2626! $Id$
     27! Bugfix for Kessler microphysics
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    7477
    7578       USE cloud_parameters,                                                   &
    76            ONLY:  l_d_cp, l_d_r, t_d_pt                                             
     79           ONLY:  l_d_cp, l_d_r, t_d_pt       
     80
     81       USE control_parameters,                                                 &
     82           ONLY:  microphysics_kessler                                       
    7783
    7884       IMPLICIT NONE
     
    97103!--    (see: Cuijpers + Duynkerke, 1993, JAS, 23)
    98104       q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
     105
    99106!
    100107!--    Supersaturation:
    101        sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
     108!--    Not in case of microphysics_kessler since qr is unallocated
     109       IF ( .NOT. microphysics_kessler ) THEN
     110          sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
     111       ENDIF
    102112
    103113    END SUBROUTINE supersaturation
Note: See TracChangeset for help on using the changeset viewer.