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

Timestamp:

Aug 25, 2009 8:06:27 AM (15 years ago)

Author:

raasch

Message:

speed optomizations +bugfix in init_ocean

File:

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

palm/trunk/SOURCE/init_ocean.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: First calculation of hyp(0) changed
 !
 ! Former revisions:
@@ -50 +50 @@
     rho_init(nzt) = rho_surface
 
-    DO  k = nzt-1, 0, -1
+    DO  k = nzt-1, 1, -1
        hyp(k) = hyp(k+1) + rho_surface * g * dzu(k)
     ENDDO
+    hyp(0) = hyp(1) + rho_surface * g * dzu(1)
 
     IF ( myid == 0 )  THEN