Changeset 3638 for palm/trunk/SOURCE/chemistry_model_mod.f90

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:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Added missing conversion factor fr2ppm for qvap
+! 
+! 
+! 3637 2018-12-20 01:51:36Z knoop
 ! Implementation of the PALM module interface
 ! 
@@ -1942 +1946 @@
     REAL(wp)                         ::  conv                                !< conversion factor
     REAL(wp), PARAMETER              ::  ppm2fr  = 1.0e-6_wp                 !< Conversion factor ppm to fraction
+    REAL(wp), PARAMETER              ::  fr2ppm  = 1.0e6_wp                  !< Conversion factor fraction to ppm
 !    REAL(wp), PARAMETER              ::  xm_air  = 28.96_wp                  !< Mole mass of dry air
 !    REAL(wp), PARAMETER              ::  xm_h2o  = 18.01528_wp               !< Mole mass of water vapor
@@ -1971 +1976 @@
        IF ( humidity ) THEN
           IF ( bulk_cloud_model )  THEN
-             tmp_qvap(:) = ( q(nzb+1:nzt,j,i) - ql(nzb+1:nzt,j,i) ) * xm_air/xm_h2o * tmp_fact(:)
+             tmp_qvap(:) = ( q(nzb+1:nzt,j,i) - ql(nzb+1:nzt,j,i) ) *  &
+                             xm_air/xm_h2o * fr2ppm * tmp_fact(:)
           ELSE
-             tmp_qvap(:) = q(nzb+1:nzt,j,i) * xm_air/xm_h2o * tmp_fact(:)
+             tmp_qvap(:) = q(nzb+1:nzt,j,i) * xm_air/xm_h2o * fr2ppm * tmp_fact(:)
           ENDIF
        ELSE
-          tmp_qvap(:) = 0.01 * tmp_fact(:)                          !< Constant value for q if water vapor is not computed
+          tmp_qvap(:) = 0.01 * xm_air/xm_h2o * fr2ppm * tmp_fact(:)          !< Constant value for q if water vapor is not computed
        ENDIF