Changeset 761 for palm/trunk

Timestamp:

Oct 5, 2011 5:58:52 PM (13 years ago)

Author:

suehring

Message:

Bugfix: Avoid divisions by zero in case of using a 'neumann' bc for the pressure at the top of the model domain.

File:

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

palm/trunk/SOURCE/poisfft.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: Avoid divisions by zero in case of using a 'neumann' bc for the
+! pressure at the top of the model domain.
+!
 !
 ! Former revisions:
@@ -1545 +1547 @@
 !--       Backward substitution
           DO  i = 0, nx
-             ar(i,nz) = ar1(i,nz-1) / tri(4,i,nz-1)
+             ar(i,nz) = ar1(i,nz-1) / ( tri(4,i,nz-1) + 1.0E-20 )
           ENDDO
           DO  k = nz-2, 0, -1