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

Timestamp:

Mar 30, 2011 9:31:40 AM (13 years ago)

Author:

raasch

Message:

formatting adjustments

File:

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

palm/trunk/SOURCE/pres.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! formatting adjustments
 !
 ! Former revisions:
@@ -104 +104 @@
 
     ddt_3d = 1.0 / dt_3d
-    d_weight_pres = 1. / weight_pres(intermediate_timestep_count)
+    d_weight_pres = 1.0 / weight_pres(intermediate_timestep_count)
 
 !
@@ -146 +146 @@
 !
 !-- Left/right
-
-    IF ( conserve_volume_flow  .AND.  ( outflow_l  .OR. outflow_r ) )  THEN
+    IF ( conserve_volume_flow  .AND.  ( outflow_l .OR. outflow_r ) )  THEN
 
        volume_flow(1)   = 0.0
@@ -161 +160 @@
 !
 !--       Sum up the volume flow through the south/north boundary
-          DO  k = nzb_2d(j,i) + 1, nzt
+          DO  k = nzb_2d(j,i)+1, nzt
              volume_flow_l(1) = volume_flow_l(1) + u(k,j,i) * dzw(k)
           ENDDO
@@ -173 +172 @@
        volume_flow = volume_flow_l  
 #endif
-       volume_flow_offset(1) = ( volume_flow_initial(1) - volume_flow(1) )    &
+       volume_flow_offset(1) = ( volume_flow_initial(1) - volume_flow(1) ) &
                                / volume_flow_area(1)
 
        DO  j = nysg, nyng
-          DO  k = nzb_2d(j,i) + 1, nzt
+          DO  k = nzb_2d(j,i)+1, nzt
              u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
           ENDDO
@@ -186 +185 @@
 !
 !-- South/north
-    IF ( conserve_volume_flow  .AND.  ( outflow_n  .OR. outflow_s ) )  THEN
+    IF ( conserve_volume_flow  .AND.  ( outflow_n .OR. outflow_s ) )  THEN
 
        volume_flow(2)   = 0.0
@@ -200 +199 @@
 !
 !--       Sum up the volume flow through the south/north boundary
-          DO  k = nzb_2d(j,i) + 1, nzt
+          DO  k = nzb_2d(j,i)+1, nzt
              volume_flow_l(2) = volume_flow_l(2) + v(k,j,i) * dzw(k)
           ENDDO
@@ -216 +215 @@
 
        DO  i = nxlg, nxrg
-          DO  k = nzb_v_inner(j,i) + 1, nzt
+          DO  k = nzb_v_inner(j,i)+1, nzt
              v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
           ENDDO
@@ -226 +225 @@
 !-- 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 is not yet known
+       IF ( simulated_time > 0.0 )  THEN ! otherwise nzb_w_inner not yet known
           w_l = 0.0;  w_l_l = 0.0
           DO  i = nxl, nxr
@@ -237 +236 @@
 #if defined( __parallel )   
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-          CALL MPI_ALLREDUCE( w_l_l(1), w_l(1), nzt, MPI_REAL, MPI_SUM, comm2d, &
-                              ierr )
+          CALL MPI_ALLREDUCE( w_l_l(1), w_l(1), nzt, MPI_REAL, MPI_SUM, &
+                              comm2d, ierr )
 #else
           w_l = w_l_l  
@@ -574 +573 @@
 !-- Correction of the provisional velocities with the current perturbation 
 !-- pressure just computed
-    IF ( conserve_volume_flow  .AND.  ( bc_lr_cyc  .OR.  bc_ns_cyc ) )  THEN
+    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