Changeset 3231 for palm/trunk/SOURCE


Ignore:
Timestamp:
Sep 6, 2018 11:53:51 AM (4 years ago)
Author:
schwenkel
Message:

Bugfix for microphysic scheme = saturation_adjustment

File:
1 edited

Legend:

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

    r3026 r3231  
    2525! -----------------
    2626! $Id$
     27! Bugfix for microphysic scheme = 'saturation_adjustment'
     28!
     29! 3026 2018-05-22 10:30:53Z schwenkel
    2730! Changed the name specific humidity to mixing ratio, since we are computing
    2831! mixing ratios.
     
    7275! Description:
    7376! ------------
    74 !> Computation of the diagnostic supersaturation sat, actual temperature t_l
    75 !< and saturation water vapor mixing ratio q_
     77!> Computation of the diagnostic supersaturation sat, actual liquid water
     78!< temperature t_l and saturation water vapor mixing ratio q_s
    7679!------------------------------------------------------------------------------!
    7780    SUBROUTINE supersaturation ( i,j,k )
     
    8487
    8588       USE control_parameters,                                                 &
    86            ONLY:  microphysics_kessler                                       
     89           ONLY:  microphysics_kessler, microphysics_sat_adjust                                       
    8790
    8891       IMPLICIT NONE
     
    107110!--    (see: Cuijpers + Duynkerke, 1993, JAS, 23)
    108111       q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
    109 
    110112!
    111113!--    Supersaturation:
    112 !--    Not in case of microphysics_kessler since qr is unallocated
    113        IF ( .NOT. microphysics_kessler ) THEN
     114!--    Not in case of microphysics_kessler or microphysics_sat_adjust
     115!--    since qr is unallocated
     116       IF ( .NOT. microphysics_kessler  .AND.                                  &
     117            .NOT.  microphysics_sat_adjust )                                THEN
    114118          sat   = ( q(k,j,i) - qr(k,j,i) - qc(k,j,i) ) / q_s - 1.0_wp
    115119       ENDIF
Note: See TracChangeset for help on using the changeset viewer.