Changeset 1342 for palm/trunk/SOURCE/pres.f90

Timestamp:

Mar 26, 2014 5:04:47 PM (10 years ago)

Author:

kanani

Message:

REAL constants defined as wp-kind

File:

• palm/trunk/SOURCE/pres.f90 (modified) (12 diffs)

palm/trunk/SOURCE/pres.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! REAL constants defined as wp-kind
 !
 ! Former revisions:
@@ -135 +135 @@
 
 
-    ddt_3d = 1.0 / dt_3d
-    d_weight_pres = 1.0 / weight_pres(intermediate_timestep_count)
+    ddt_3d = 1.0_wp / dt_3d
+    d_weight_pres = 1.0_wp / weight_pres(intermediate_timestep_count)
 
 !
@@ -180 +180 @@
     IF ( conserve_volume_flow  .AND.  ( outflow_l .OR. outflow_r ) )  THEN
 
-       volume_flow(1)   = 0.0
-       volume_flow_l(1) = 0.0
+       volume_flow(1)   = 0.0_wp
+       volume_flow_l(1) = 0.0_wp
 
        IF ( outflow_l )  THEN
@@ -219 +219 @@
     IF ( conserve_volume_flow  .AND.  ( outflow_n .OR. outflow_s ) )  THEN
 
-       volume_flow(2)   = 0.0
-       volume_flow_l(2) = 0.0
+       volume_flow(2)   = 0.0_wp
+       volume_flow_l(2) = 0.0_wp
 
        IF ( outflow_s )  THEN
@@ -257 +257 @@
 !-- Remove mean vertical velocity
     IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1 )  THEN
-       IF ( simulated_time > 0.0 )  THEN ! otherwise nzb_w_inner not yet known
-          w_l = 0.0;  w_l_l = 0.0
+       IF ( simulated_time > 0.0_wp )  THEN ! otherwise nzb_w_inner not yet known
+          w_l = 0.0_wp;  w_l_l = 0.0_wp
           DO  i = nxl, nxr
              DO  j = nys, nyn
@@ -295 +295 @@
           DO  j = nys-1, nyn+1
              DO  k = nzb, nzt+1
-                d(k,j,i) = 0.0
+                d(k,j,i) = 0.0_wp
              ENDDO
           ENDDO
@@ -305 +305 @@
           DO  j = nys, nyn
              DO  k = nzb+1, nzt
-                d(k,j,i) = 0.0
+                d(k,j,i) = 0.0_wp
              ENDDO
           ENDDO
@@ -312 +312 @@
     ENDIF
 
-    localsum  = 0.0
-    threadsum = 0.0
+    localsum  = 0.0_wp
+    threadsum = 0.0_wp
 
 #if defined( __ibm )
@@ -432 +432 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                tend(nzb_s_inner(j,i),j,i) = 0.0
+                tend(nzb_s_inner(j,i),j,i) = 0.0_wp
              ENDDO
           ENDDO
@@ -460 +460 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                tend(nzt+1,j,i) = 0.0
+                tend(nzt+1,j,i) = 0.0_wp
              ENDDO
           ENDDO
@@ -577 +577 @@
 !-- pressure just computed
     IF ( conserve_volume_flow  .AND.  ( bc_lr_cyc .OR. bc_ns_cyc ) )  THEN
-       volume_flow_l(1) = 0.0
-       volume_flow_l(2) = 0.0
+       volume_flow_l(1) = 0.0_wp
+       volume_flow_l(2) = 0.0_wp
     ENDIF
 
@@ -704 +704 @@
 !-- a possible PE-sum is computed in flow_statistics 
     CALL cpu_log( log_point_s(1), 'divergence', 'start' )
-    sums_divnew_l = 0.0
+    sums_divnew_l = 0.0_wp
 
 !
 !-- d must be reset to zero because it can contain nonzero values below the
 !-- topography
-    IF ( topography /= 'flat' )  d = 0.0
-
-    localsum  = 0.0
-    threadsum = 0.0
+    IF ( topography /= 'flat' )  d = 0.0_wp
+
+    localsum  = 0.0_wp
+    threadsum = 0.0_wp
 
     !$OMP PARALLEL PRIVATE (i,j,k) FIRSTPRIVATE(threadsum) REDUCTION(+:localsum)