Changeset 2502 for palm/trunk/SOURCE/init_ocean.f90

Timestamp:

Sep 26, 2017 3:51:59 PM (7 years ago)

Author:

gronemeier

Message:

Bugfix: corrected reference state for rho_ocean

File:

• palm/trunk/SOURCE/init_ocean.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_ocean.f90

@@ -21 +21 @@
 ! -----------------
 ! 
-!
+! 
 ! Former revisions:
 ! ------------------
 ! $Id$
+! Bugfix: use equation of state for seawater to calculate rho_ocean_init at nzt
+!
+! 2195 2017-03-23 08:15:17Z raasch
 !
 ! 2194 2017-03-23 08:03:04Z raasch
@@ -126 +129 @@
        rho_reference = rho_surface * 0.5_wp * dzu(nzt+1)
 
-       DO  k = nzt-1, 0, -1
+       DO  k = nzt, 0, -1
 
           sa_l = 0.5_wp * ( sa_init(k) + sa_init(k+1) )
@@ -139 +142 @@
        rho_reference = rho_reference / ( zw(nzt) - zu(nzb) )
 
-       DO  k = nzt-1, 0, -1
-          hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_ocean_init(k) + rho_ocean_init(k+1 ) ) * &
-                              dzu(k+1)
+       DO  k = nzt, 0, -1
+          hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_ocean_init(k)                 &
+                                           + rho_ocean_init(k+1) ) * dzu(k+1)
        ENDDO