Changeset 76 for palm/trunk/SOURCE/poisfft.f90

Timestamp:

Mar 29, 2007 12:58:32 AM (17 years ago)

Author:

raasch

Message:

changes for Neumann p conditions both at bottom and top

File:

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

palm/trunk/SOURCE/poisfft.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Tridiagonal coefficients adjusted for Neumann boundary conditions both at
+! the bottom and the top.
 !
 ! Former revisions:
@@ -393 +394 @@
 !------------------------------------------------------------------------------!
 
+          USE control_parameters
+
           IMPLICIT NONE
 
@@ -443 +446 @@
              ENDDO
           ENDDO
+
+!
+!--       Indices i=0, j=0 correspond to horizontally averaged pressure.
+!--       The respective values of ar should be zero at all k-levels if
+!--       acceleration of horizontally averaged vertical velocity is zero.
+          IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1 )  THEN
+             IF ( j == 0  .AND.  nxl_z == 0 )  THEN
+#if defined( __parallel )
+                DO  k = 1, nz
+                   ar(nxl_z,j,k) = 0.0
+                ENDDO
+#else
+                DO  k = 1, nz
+                   ar(k,j,nxl_z) = 0.0
+                ENDDO
+#endif
+             ENDIF
+          ENDIF
 
        END SUBROUTINE substi
@@ -1466 +1487 @@
           IMPLICIT NONE
 
-          INTEGER ::  i, j, k
+          INTEGER ::  i, k
 
           REAL, DIMENSION(0:nx,nz)       ::  ar
@@ -1495 +1516 @@
           ENDDO
 
+!
+!--       Indices i=0, j=0 correspond to horizontally averaged pressure.
+!--       The respective values of ar should be zero at all k-levels if
+!--       acceleration of horizontally averaged vertical velocity is zero.
+          IF ( ibc_p_b == 1  .AND.  ibc_p_t == 1 )  THEN
+             IF ( j == 0 )  THEN
+                DO  k = 1, nz
+                   ar(0,k) = 0.0
+                ENDDO
+             ENDIF
+          ENDIF
+
        END SUBROUTINE substi_1dd