Ignore:
Timestamp:
Oct 19, 2007 3:47:46 PM (14 years ago)
Author:
raasch
Message:

further bugfixes for the ocean part

File:
1 edited

Legend:

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

    r110 r124  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: calculation of density flux in the ocean now starts from nzb+1
    77!
    88! Former revisions:
     
    382382!--                the bottom and top surface layer
    383383                   DO  j = nys, nyn
    384                       DO  k = nzb_s_inner(j,i), nzt
     384                      DO  k = nzb_s_inner(j,i)+1, nzt
    385385                         tend(k,j,i) = tend(k,j,i) +                    &
    386386                                       kh(k,j,i) * g / prho_reference * &
     
    419419!--                the bottom and top surface layer
    420420                   DO  j = nys, nyn
    421                       DO  k = nzb_s_inner(j,i), nzt
     421                      DO  k = nzb_s_inner(j,i)+1, nzt
    422422                         tend(k,j,i) = tend(k,j,i) -                &
    423423                                       kh(k,j,i) * g / rho(k,j,i) * &
     
    814814!--             So far in the ocean no special treatment of density flux in the
    815815!--             bottom and top surface layer
    816                 k= nzb_s_inner(j,i)
    817                 tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
    818                      ( rho(k+1,j,i) - rho(k,j,i) ) * dzu(k+1)
    819816                DO  k = nzb_s_inner(j,i)+1, nzt
    820817                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
     
    847844!--             So far in the ocean no special treatment of density flux in the
    848845!--             bottom and top surface layer
    849                 DO  k = nzb_s_inner(j,i), nzt
     846                DO  k = nzb_s_inner(j,i)+1, nzt
    850847                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / rho(k,j,i) * &
    851848                                      ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
Note: See TracChangeset for help on using the changeset viewer.