Ignore:
Timestamp:
Mar 7, 2008 1:42:18 PM (16 years ago)
Author:
raasch
Message:

preliminary update for the turbulence recycling method

File:
1 edited

Legend:

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

    r110 r151  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix in volume flow control for non-cyclic boundary conditions
    77!
    88! Former revisions:
     
    110110                               / volume_flow_area(1)
    111111
    112        DO  j = nys, nyn
     112       DO  j = nys-1, nyn+1
    113113          DO  k = nzb_v_inner(j,i) + 1, nzt
    114114             u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
    115115          ENDDO
    116116       ENDDO
    117 
    118        CALL exchange_horiz( u )
    119117
    120118    ENDIF
     
    150148                               / volume_flow_area(2)
    151149
    152        DO  i = nxl, nxr
     150       DO  i = nxl-1, nxr+1
    153151          DO  k = nzb_v_inner(j,i) + 1, nzt
    154152             v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
    155153          ENDDO
    156154       ENDDO
    157 
    158        CALL exchange_horiz( v )
    159155
    160156    ENDIF
Note: See TracChangeset for help on using the changeset viewer.