Ignore:
Timestamp:
Oct 16, 2017 5:12:01 PM (4 years ago)
Author:
boeske
Message:

enable simulations with complex terrain

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_grid.f90

    r2478 r2550  
    2525! -----------------
    2626! $Id$
     27! Set lateral boundary conditions for topography on all three ghost layers
     28!
     29! 2478 2017-09-18 13:37:24Z suehring
    2730! Bugfix, correct flag for use_top
    2831!
     
    10411044         
    10421045          IF ( .NOT. bc_ns_cyc )  THEN
    1043              IF ( nys == 0  )  nzb_local(-1,:)   = nzb_local(0,:)
    1044              IF ( nyn == ny )  nzb_local(ny+1,:) = nzb_local(ny,:)
     1046             IF ( nys == 0  )  THEN
     1047                nzb_local(-1,:)   = nzb_local(0,:)
     1048                nzb_local(-2,:)   = nzb_local(0,:)
     1049                nzb_local(-3,:)   = nzb_local(0,:)
     1050             ELSEIF ( nyn == ny )  THEN
     1051                nzb_local(ny+1,:) = nzb_local(ny,:)
     1052                nzb_local(ny+2,:) = nzb_local(ny,:)
     1053                nzb_local(ny+3,:) = nzb_local(ny,:)
     1054             ENDIF
    10451055          ENDIF
    10461056
    10471057          IF ( .NOT. bc_lr_cyc )  THEN
    1048              IF ( nxl == 0  )  nzb_local(:,-1)   = nzb_local(:,0)
    1049              IF ( nxr == nx )  nzb_local(:,nx+1) = nzb_local(:,nx)         
     1058             IF ( nxl == 0  )  THEN
     1059                nzb_local(:,-1)   = nzb_local(:,0)
     1060                nzb_local(:,-2)   = nzb_local(:,0)
     1061                nzb_local(:,-3)   = nzb_local(:,0)
     1062             ELSEIF ( nxr == nx )  THEN
     1063                nzb_local(:,nx+1) = nzb_local(:,nx)
     1064                nzb_local(:,nx+2) = nzb_local(:,nx)         
     1065                nzb_local(:,nx+3) = nzb_local(:,nx)         
     1066             ENDIF         
    10501067          ENDIF
    10511068
Note: See TracChangeset for help on using the changeset viewer.