Changeset 2502 for palm/trunk


Ignore:
Timestamp:
Sep 26, 2017 3:51:59 PM (7 years ago)
Author:
gronemeier
Message:

Bugfix: corrected reference state for rho_ocean

File:
1 edited

Legend:

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

    r2195 r2502  
    2121! -----------------
    2222!
    23 !
     23! 
    2424! Former revisions:
    2525! ------------------
    2626! $Id$
     27! Bugfix: use equation of state for seawater to calculate rho_ocean_init at nzt
     28!
     29! 2195 2017-03-23 08:15:17Z raasch
    2730!
    2831! 2194 2017-03-23 08:03:04Z raasch
     
    126129       rho_reference = rho_surface * 0.5_wp * dzu(nzt+1)
    127130
    128        DO  k = nzt-1, 0, -1
     131       DO  k = nzt, 0, -1
    129132
    130133          sa_l = 0.5_wp * ( sa_init(k) + sa_init(k+1) )
     
    139142       rho_reference = rho_reference / ( zw(nzt) - zu(nzb) )
    140143
    141        DO  k = nzt-1, 0, -1
    142           hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_ocean_init(k) + rho_ocean_init(k+1 ) ) * &
    143                               dzu(k+1)
     144       DO  k = nzt, 0, -1
     145          hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_ocean_init(k)                 &
     146                                           + rho_ocean_init(k+1) ) * dzu(k+1)
    144147       ENDDO
    145148
Note: See TracChangeset for help on using the changeset viewer.