Changeset 1418 for palm/trunk/SOURCE/surface_coupler.f90

Timestamp:

Jun 6, 2014 1:05:08 PM (10 years ago)

Author:

fricke

Message:

Bugfixes concerning grid stretching for the ocean and calculation of the salinity flux in routine surface_coupler

File:

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

palm/trunk/SOURCE/surface_coupler.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: For caluclation of the salinity flux at the sea surface, 
+!          the given value for salinity must be in percent and not in psu
 ! 
 ! Former revisions:
@@ -400 +401 @@
 !--       following Steinhorn (1991), JPO 21, pp. 1681-1683:
 !--       S'w' = -S * evaporation / ( rho_water * ( 1 - S ) )
-          saswst = -1.0_wp * sa(nzt,:,:) * qswst_remote /  &
-                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) ) )
+          saswst = -1.0_wp * sa(nzt,:,:) * 0.001 * qswst_remote /  &
+                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) * 0.001 ) )
        ENDIF