Changeset 875

Timestamp:

Apr 2, 2012 3:35:15 PM (12 years ago)

Author:

gryschka

Message:

Bugfix in case of dirichlet inflow bc at the right or north boundary

File:

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

palm/trunk/SOURCE/boundary_conds.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+!
+! Bugfix in case of dirichlet inflow bc at the right or north boundary
 ! 
 !  
@@ -190 +192 @@
 
 !
-!--    Lateral boundary conditions at the inflow. Quasi Neumann conditions
-!--    are needed for the wall normal velocity in order to ensure zero
-!--    divergence. Dirichlet conditions are used for all other quantities.
+!--    In case of inflow at the south boundary the boundary for v is at nys
+!--    and in case of inflow at the left boundary the boundary for u is at nxl.
+!--    Since in prognostic_equations (cache optimized version) these levels are
+!--    handled as a prognostic level, boundary values have to be restored here.
        IF ( inflow_s )  THEN
           v_p(:,nys,:) = v_p(:,nys-1,:)
-       ELSEIF ( inflow_n ) THEN
-          v_p(:,nyn,:) = v_p(:,nyn+1,:)
        ELSEIF ( inflow_l ) THEN
           u_p(:,:,nxl) = u_p(:,:,nxl-1)
-       ELSEIF ( inflow_r ) THEN
-          u_p(:,:,nxr) = u_p(:,:,nxr+1)
        ENDIF