Ignore:
Timestamp:
Mar 22, 2007 9:54:05 AM (15 years ago)
Author:
raasch
Message:

preliminary update for changes concerning non-cyclic boundary conditions

File:
1 edited

Legend:

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

    r60 r75  
    3737! Call for all grid points
    3838!------------------------------------------------------------------------------!
    39     SUBROUTINE wall_fluxes( wall_flux, a, b, c1, c2, ixp, jyp, nzb_uvw_inner, &
     39    SUBROUTINE wall_fluxes( wall_flux, a, b, c1, c2, nzb_uvw_inner, &
    4040                            nzb_uvw_outer, wall )
    4141
     
    4848       IMPLICIT NONE
    4949
    50        INTEGER ::  i, ixp, j, jyp, k, wall_index
     50       INTEGER ::  i, j, k, wall_index
    5151
    5252       INTEGER, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1) ::  nzb_uvw_inner, &
     
    5555       REAL ::  pts, pt_i, rifs, u_i, v_i, us_wall, vel_total, ws, wspts
    5656
    57        REAL, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1)           ::  wall
    58        REAL, DIMENSION(nzb:nzt+1,nys:nyn+jyp,nxl:nxr+ixp) ::  wall_flux
     57       REAL, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1)   ::  wall
     58       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  wall_flux
    5959
    6060
     
    6363       wall_index = NINT( a+ 2*b + 3*c1 + 4*c2 )
    6464
    65        DO  i = nxl, nxr+ixp
    66           DO  j = nys, nyn+jyp
     65       DO  i = nxl, nxr
     66          DO  j = nys, nyn
    6767
    6868             IF ( wall(j,i) /= 0.0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.