Changeset 1427

Timestamp:

Jul 7, 2014 2:04:59 PM (10 years ago)

Author:

maronga

Message:

bugfix in surface coupler

File:

• palm/trunk/SOURCE/surface_coupler.f90 (modified) (5 diffs)

palm/trunk/SOURCE/surface_coupler.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: value of l_v corrected.
 ! 
 ! Former revisions:
@@ -72 +72 @@
                total_2d_o, tswst, u, usws, uswst, v, vsws, vswst
 
+    USE cloud_parameters,                                                      &
+        ONLY:  cp, l_v
+
     USE control_parameters,                                                    &
         ONLY:  coupling_mode, coupling_mode_remote, coupling_topology,         &
@@ -92 +95 @@
     REAL(wp)    ::  time_since_reference_point_rem        !: 
     REAL(wp)    ::  total_2d(-nbgp:ny+nbgp,-nbgp:nx+nbgp) !: 
+
+    REAL(wp)    ::  cpw = 4218.0_wp !: heat capacity of water at constant pressure
 
 #if defined( __parallel )
@@ -399 +404 @@
 !--       * latent heat of vaporization in m2/s2, or 540 cal/g, or 40.65 kJ/mol
 !--       /(rho_atm(=1.0)*c_p)
-          tswst = tswst + qswst_remote * 2.2626108E6_wp / 1005.0_wp
+          tswst = tswst + qswst_remote * l_v / cp
 !
 !--        ...and convert it to a salinity flux at the sea surface (top)
 !--       following Steinhorn (1991), JPO 21, pp. 1681-1683:
 !--       S'w' = -S * evaporation / ( rho_water * ( 1 - S ) )
-          saswst = -1.0_wp * sa(nzt,:,:) * 0.001 * qswst_remote /  &
-                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) * 0.001 ) )
+          saswst = -1.0_wp * sa(nzt,:,:) * 0.001_wp * qswst_remote /  &
+                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) * 0.001_wp ) )
        ENDIF
 
@@ -412 +417 @@
 !--    (constants are the specific heat capacities for air and water)
 !--    now tswst is the ocean top heat flux
-       tswst = tswst / rho(nzt,:,:) * 1005.0_wp / 4218.0_wp
+       tswst = tswst / rho(nzt,:,:) * cp / cpw
 
 !