Index: palm/trunk/SOURCE/turbulence_closure_mod.f90
===================================================================
--- palm/trunk/SOURCE/turbulence_closure_mod.f90	(revision 3359)
+++ palm/trunk/SOURCE/turbulence_closure_mod.f90	(revision 3361)
@@ -170,5 +170,5 @@
 
     USE basic_constants_and_equations_mod,                                     &
-        ONLY:  g, kappa, lv_d_cp, lv_d_rd
+        ONLY:  g, kappa, lv_d_cp, lv_d_rd, rd_d_rv
 
     USE control_parameters,                                                    &
@@ -3006,6 +3006,6 @@
                          k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                  &
                                        ( q(k,j,i) - ql(k,j,i) ) *              &
-                              ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
-                              ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
+                              ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /        &
+                              ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *         &
                               ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                          k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3066,6 +3066,6 @@
                                k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *            &
                                           ( q(k,j,i) - ql(k,j,i) ) *           &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
                                  ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                                k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3103,6 +3103,6 @@
                             k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
                                           ( q(k,j,i) - ql(k,j,i) ) *           &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
                                  ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                             k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3139,6 +3139,6 @@
                             k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
                                           ( q(k,j,i) - ql(k,j,i) ) *           &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
-                                 ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
+                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
                                  ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                             k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3178,6 +3178,6 @@
                             k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
                                        ( q(k,j,i) - ql(k,j,i) ) *              &
-                              ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
-                              ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
+                              ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /        &
+                              ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *         &
                               ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                             k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3794,6 +3794,6 @@
                    k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                        &
                                  ( q(k,j,i) - ql(k,j,i) ) *                    &
-                        ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
-                        ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
+                        ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /              &
+                        ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *               &
                         ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                    k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3851,6 +3851,6 @@
                         k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                   &
                                    ( q(k,j,i) - ql(k,j,i) ) *                  &
-                          ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /           &
-                          ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *            &
+                          ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /            &
+                          ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *             &
                           ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                         k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3888,6 +3888,6 @@
                        k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
                                      ( q(k,j,i) - ql(k,j,i) ) *                &
-                            ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
-                            ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
+                            ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /          &
+                            ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *           &
                             ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                        k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3924,6 +3924,6 @@
                        k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
                                      ( q(k,j,i) - ql(k,j,i) ) *                &
-                            ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
-                            ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
+                            ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /          &
+                            ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *           &
                             ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                        k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
@@ -3964,6 +3964,6 @@
                       k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                     &
                                  ( q(k,j,i) - ql(k,j,i) ) *                    &
-                        ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
-                        ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
+                        ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /              &
+                        ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *               &
                         ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
                       k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
