Changeset 107 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Aug 17, 2007 1:54:45 PM (17 years ago)

Author:

raasch

Message:

further bugfix for non-cyclic BCs

File:

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

palm/trunk/SOURCE/boundary_conds.f90

@@ -285 +285 @@
                                        ( u(k,-1,i) - u(k,0,i) ) * ddy
 
-             v_p(k,-1,i) = v(k,-1,i) - dt_3d * tsc(2) * c_v(k,i) * &
+             v_p(k,0,i)  = v(k,0,i)  - dt_3d * tsc(2) * c_v(k,i) * &
                                        ( v(k,0,i) - v(k,1,i) ) * ddy
 
@@ -433 +433 @@
                 IF ( denom /= 0.0 )  THEN
                    c_u(k,j) = -c_max * ( u(k,j,1) - u_m_l(k,j,1) ) / denom
-                   IF ( c_u(k,j) > 0.0 )  THEN
+                   IF ( c_u(k,j) < 0.0 )  THEN
                       c_u(k,j) = 0.0
-                   ELSEIF ( c_u(k,j) < -c_max )  THEN
-                      c_u(k,j) = -c_max
-                   ENDIF
-                ELSE
-                   c_u(k,j) = -c_max
+                   ELSEIF ( c_u(k,j) > c_max )  THEN
+                      c_u(k,j) = c_max
+                   ENDIF
+                ELSE
+                   c_u(k,j) = c_max
                 ENDIF