Changeset 1742 for palm

Timestamp:

Jan 13, 2016 9:50:06 AM (9 years ago)

Author:

raasch

Message:

bugfix for outflow Neumann boundary conditions at bottom and top

File:

• palm/trunk/SOURCE/boundary_conds.f90 (modified) (3 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! bugfix for outflow Neumann boundary conditions at bottom and top
 !
 ! Former revisions:
@@ -521 +521 @@
              v_p(nzt+1,0,:)  = v_init(nzt+1)
           ELSE
-             u_p(nzt+1,-1,:) = u(nzt,-1,:)
-             v_p(nzt+1,0,:)  = v(nzt,0,:)
+             u_p(nzt+1,-1,:) = u_p(nzt,-1,:)
+             v_p(nzt+1,0,:)  = v_p(nzt,0,:)
           ENDIF
           w_p(nzt:nzt+1,-1,:) = 0.0_wp
@@ -941 +941 @@
              v_p(nzt+1,:,nx+1) = v_p(nzt,:,nx+1)
           ENDIF
-          w(nzt:nzt+1,:,nx+1) = 0.0_wp
+          w_p(nzt:nzt+1,:,nx+1) = 0.0_wp
 
        ENDIF