Changeset 3638


Ignore:
Timestamp:
Dec 20, 2018 1:18:23 PM (6 years ago)
Author:
forkel
Message:

chemistry_model_mod: Added missing conversion factor fr2ppm for qvap

File:
1 edited

Legend:

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

    r3637 r3638  
    2727! -----------------
    2828! $Id$
     29! Added missing conversion factor fr2ppm for qvap
     30!
     31!
     32! 3637 2018-12-20 01:51:36Z knoop
    2933! Implementation of the PALM module interface
    3034!
     
    19421946    REAL(wp)                         ::  conv                                !< conversion factor
    19431947    REAL(wp), PARAMETER              ::  ppm2fr  = 1.0e-6_wp                 !< Conversion factor ppm to fraction
     1948    REAL(wp), PARAMETER              ::  fr2ppm  = 1.0e6_wp                  !< Conversion factor fraction to ppm
    19441949!    REAL(wp), PARAMETER              ::  xm_air  = 28.96_wp                  !< Mole mass of dry air
    19451950!    REAL(wp), PARAMETER              ::  xm_h2o  = 18.01528_wp               !< Mole mass of water vapor
     
    19711976       IF ( humidity ) THEN
    19721977          IF ( bulk_cloud_model )  THEN
    1973              tmp_qvap(:) = ( q(nzb+1:nzt,j,i) - ql(nzb+1:nzt,j,i) ) * xm_air/xm_h2o * tmp_fact(:)
     1978             tmp_qvap(:) = ( q(nzb+1:nzt,j,i) - ql(nzb+1:nzt,j,i) ) *  &
     1979                             xm_air/xm_h2o * fr2ppm * tmp_fact(:)
    19741980          ELSE
    1975              tmp_qvap(:) = q(nzb+1:nzt,j,i) * xm_air/xm_h2o * tmp_fact(:)
     1981             tmp_qvap(:) = q(nzb+1:nzt,j,i) * xm_air/xm_h2o * fr2ppm * tmp_fact(:)
    19761982          ENDIF
    19771983       ELSE
    1978           tmp_qvap(:) = 0.01 * tmp_fact(:)                          !< Constant value for q if water vapor is not computed
     1984          tmp_qvap(:) = 0.01 * xm_air/xm_h2o * fr2ppm * tmp_fact(:)          !< Constant value for q if water vapor is not computed
    19791985       ENDIF
    19801986
Note: See TracChangeset for help on using the changeset viewer.