Ignore:
Timestamp:
Sep 23, 2009 9:40:33 AM (12 years ago)
Author:
raasch
Message:

in-situ AND potential density are calculated and used in the ocean version

File:
1 edited

Legend:

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

    r364 r388  
    44! Actual revisions:
    55! -----------------
     6! Bugfix: wrong sign in buoyancy production of ocean part in case of not using
     7!         the reference density (only in 3D routine production_e)
    68! Bugfix to avoid zero division by km_neutral
    79!
     
    450452                   DO  j = nys, nyn
    451453                      DO  k = nzb_s_inner(j,i)+1, nzt
    452                          tend(k,j,i) = tend(k,j,i) +                    &
    453                                        kh(k,j,i) * g / prho_reference * &
     454                         tend(k,j,i) = tend(k,j,i) +                   &
     455                                       kh(k,j,i) * g / rho_reference * &
    454456                                       ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k)
    455457                      ENDDO
     
    487489                   DO  j = nys, nyn
    488490                      DO  k = nzb_s_inner(j,i)+1, nzt
    489                          tend(k,j,i) = tend(k,j,i) -                &
     491                         tend(k,j,i) = tend(k,j,i) +                &
    490492                                       kh(k,j,i) * g / rho(k,j,i) * &
    491493                                       ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k)
     
    946948!--             bottom and top surface layer
    947949                DO  k = nzb_s_inner(j,i)+1, nzt
    948                    tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
     950                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / rho_reference * &
    949951                                      ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
    950952                ENDDO
Note: See TracChangeset for help on using the changeset viewer.