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

Timestamp:

Sep 24, 2018 6:28:20 AM (6 years ago)

Author:

raasch

Message:

bugfix for r2502: wrong calculation of hyp(nzt) in second iterative step

File:

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

palm/trunk/SOURCE/init_ocean.f90

@@ -25 +25 @@
 ! ------------------
 ! $Id$
+! bugfix for r2502: wrong calculation of hyp(nzt) in second iterative step
+! 
+! 2846 2018-03-01 08:48:47Z raasch
 ! FORTRAN bugfix for r2845
 ! 
@@ -132 +135 @@
 !-- First step: Calculate pressure using reference density
     hyp(nzt+1) = surface_pressure * 100.0_wp
-
-    hyp(nzt)      = hyp(nzt+1) + rho_surface * g * 0.5_wp * dzu(nzt+1)
-    rho_ocean_init(nzt) = rho_surface
+    hyp(nzt)   = hyp(nzt+1) + rho_surface * g * 0.5_wp * dzu(nzt+1)
+    rho_ocean_init(nzt)   = rho_surface
+    rho_ocean_init(nzt+1) = rho_surface  ! only required for output
 
     DO  k = nzt-1, 1, -1
@@ -161 +164 @@
        rho_reference = rho_reference / ( zw(nzt) - zu(nzb) )
 
-       DO  k = nzt, 0, -1
+       hyp(nzt) = hyp(nzt+1) + rho_surface * g * 0.5_wp * dzu(nzt+1)
+       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)