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

Timestamp:

Oct 6, 2011 9:32:09 AM (13 years ago)

Author:

suehring

Message:

comment added

File:

• palm/trunk/SOURCE/poisfft.f90 (modified) (1 diff)

palm/trunk/SOURCE/poisfft.f90

@@ -1548 +1548 @@
 !
 !--       Backward substitution
+!--       Note, the add of 1.0E-20 in the denominator is due to avoid divisions
+!--       by zero appearing if the pressure bc is set to neumann at the top of
+!--       the model domain.
           DO  i = 0, nx
              ar(i,nz) = ar1(i,nz-1) / ( tri(4,i,nz-1) + 1.0E-20 )