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

Timestamp:

Apr 8, 2014 3:21:23 PM (10 years ago)

Author:

heinze

Message:

REAL constants provided with KIND-attribute

File:

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

palm/trunk/SOURCE/surface_coupler.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -206 +207 @@
 !
 !--       Send heat flux at bottom surface to the ocean
-          total_2d_a = 0.0
-          total_2d   = 0.0
+          total_2d_a = 0.0_wp
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = shf(nys:nyn,nxl:nxr)
 
@@ -216 +217 @@
 !--       Send humidity flux at bottom surface to the ocean
           IF ( humidity )  THEN
-             total_2d_a = 0.0
-             total_2d   = 0.0
+             total_2d_a = 0.0_wp
+             total_2d   = 0.0_wp
              total_2d(nys:nyn,nxl:nxr) = qsws(nys:nyn,nxl:nxr)
 
@@ -236 +237 @@
 !
 !--       Send momentum flux (u) at bottom surface to the ocean
-          total_2d_a = 0.0 
-          total_2d   = 0.0
+          total_2d_a = 0.0_wp 
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = usws(nys:nyn,nxl:nxr)
           CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, &
@@ -244 +245 @@
 !
 !--       Send momentum flux (v) at bottom surface to the ocean
-          total_2d_a = 0.0
-          total_2d   = 0.0
+          total_2d_a = 0.0_wp
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = vsws(nys:nyn,nxl:nxr)
           CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, &
@@ -338 +339 @@
 !
 !--       Send surface temperature to atmosphere
-          total_2d_o = 0.0
-          total_2d   = 0.0
+          total_2d_o = 0.0_wp
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = pt(nzt,nys:nyn,nxl:nxr)
 
@@ -367 +368 @@
 !
 !--       Send u to atmosphere
-          total_2d_o = 0.0  
-          total_2d   = 0.0
+          total_2d_o = 0.0_wp  
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = u(nzt,nys:nyn,nxl:nxr)
           CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, &
@@ -375 +376 @@
 !
 !--       Send v to atmosphere
-          total_2d_o = 0.0
-          total_2d   = 0.0
+          total_2d_o = 0.0_wp
+          total_2d   = 0.0_wp
           total_2d(nys:nyn,nxl:nxr) = v(nzt,nys:nyn,nxl:nxr)
           CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, &
@@ -396 +397 @@
 !--       following Steinhorn (1991), JPO 21, pp. 1681-1683:
 !--       S'w' = -S * evaporation / ( rho_water * ( 1 - S ) )
-          saswst = -1.0 * sa(nzt,:,:) * qswst_remote /  &
-                    ( rho(nzt,:,:) * ( 1.0 - sa(nzt,:,:) ) )
+          saswst = -1.0_wp * sa(nzt,:,:) * qswst_remote /  &
+                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) ) )
        ENDIF
 
@@ -474 +475 @@
        dny2 = 2 * ( dny / 2 )
 
-       total_2d_a = 0.0
+       total_2d_a = 0.0_wp
 !
 !--    Interpolation from ocean-grid-layer to atmosphere-grid-layer