Ignore:
Timestamp:
Aug 16, 2007 2:30:26 PM (14 years ago)
Author:
raasch
Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:
1 edited

Legend:

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

    r98 r106  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: wrong sign removed from the buoyancy production term in the case
     7! use_reference = .T.,
     8! u_0 and v_0 are calculated for nxr+1, nyn+1 also (otherwise these values are
     9! not available in case of non-cyclic boundary conditions)
    710!
    811! Former revisions:
     
    386389                   DO  j = nys, nyn
    387390                      DO  k = nzb_diff_s_inner(j,i), nzt_diff
    388                          tend(k,j,i) = tend(k,j,i) +                  &
     391                         tend(k,j,i) = tend(k,j,i) -                  &
    389392                                       kh(k,j,i) * g / pt_reference * &
    390393                                       ( pt(k+1,j,i) - pt(k-1,j,i) ) * dd2zu(k)
     
    976979!--       (otherwise the timestep may be significantly reduced by large
    977980!--       surface winds).
     981!--       Upper bounds are nxr+1 and nyn+1 because otherwise these values are
     982!--       not available in case of non-cyclic boundary conditions.
    978983!--       WARNING: the exact analytical solution would require the determination
    979984!--                of the eddy diffusivity by km = u* * kappa * zp / phi_m.
    980985          !$OMP PARALLEL DO PRIVATE( ku, kv )
    981           DO  i = nxl, nxr
    982              DO  j = nys, nyn
     986          DO  i = nxl, nxr+1
     987             DO  j = nys, nyn+1
    983988
    984989                ku = nzb_u_inner(j,i)+1
Note: See TracChangeset for help on using the changeset viewer.