Ignore:
Timestamp:
Mar 23, 2018 3:13:00 PM (3 years ago)
Author:
suehring
Message:

Bugfix in setting boundary condition for topography index array in case of non-cyclic boundary conditions and Piascek-Williams advection scheme is used.

File:
1 edited

Legend:

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

    r2918 r2927  
    2525! -----------------
    2626! $Id$
     27! Bugfix, setting boundary conditions for topography index array.
     28!
     29! 2918 2018-03-21 15:52:14Z gronemeier
    2730! Moved init_mixing_length to turbulence_closure_mod.f90
    2831!
     
    655658    IF ( .NOT. bc_ns_cyc )  THEN
    656659       IF ( nys == 0  )  THEN
    657           nzb_local(-1,:)   = nzb_local(0,:)
     660          DO  i = 1, nbgp 
     661             nzb_local(nys-i,:)   = nzb_local(nys,:)
     662          ENDDO
    658663       ELSEIF ( nyn == ny )  THEN
    659           nzb_local(ny+1,:) = nzb_local(ny,:)
     664          DO  i = 1, nbgp 
     665             nzb_local(ny+i,:) = nzb_local(ny,:)
     666          ENDDO
    660667       ENDIF
    661668    ENDIF
     
    663670    IF ( .NOT. bc_lr_cyc )  THEN
    664671       IF ( nxl == 0  )  THEN
    665           nzb_local(:,-1)   = nzb_local(:,0)
    666           nzb_local(:,-2)   = nzb_local(:,0)
    667           nzb_local(:,-3)   = nzb_local(:,0)
     672          DO  i = 1, nbgp 
     673             nzb_local(:,nxl-i)   = nzb_local(:,nxl)
     674          ENDDO
    668675       ELSEIF ( nxr == nx )  THEN
    669           nzb_local(:,nx+1) = nzb_local(:,nx)
    670           nzb_local(:,nx+2) = nzb_local(:,nx)         
    671           nzb_local(:,nx+3) = nzb_local(:,nx)         
     676          DO  i = 1, nbgp 
     677             nzb_local(:,nx+i) = nzb_local(:,nx)
     678          ENDDO
    672679       ENDIF         
    673680    ENDIF
Note: See TracChangeset for help on using the changeset viewer.