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

Timestamp:

Jun 4, 2007 8:07:41 AM (17 years ago)

Author:

raasch

Message:

more preliminary uncomplete changes for ocean version

File:

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

palm/trunk/SOURCE/init_ocean.f90

@@ -8 +8 @@
 ! Former revisions:
 ! ------------------
-! $Id:$
+! $Id$
 !
 ! Initial revision (raasch 01/06/07)
@@ -37 +37 @@
 !
 !-- Calculate initial vertical profile of hydrostatic pressure (in Pa)
+!-- and the reference density (used later in buoyancy term)
     hyp(nzt+1) = surface_pressure * 100.0
 
     hyp(nzt)   = hyp(nzt+1) + rho_surface * g * 0.5 * dzu(nzt+1)
+    rho_ref    = rho_surface * 0.5 * dzu(nzt+1)
 
     DO  k = nzt-1, 0, -1
@@ -46 +48 @@
        pt_l = 0.5 * ( pt_init(k) + pt_init(k+1) )
 
-       rho_l = eqn_state_seawater_func( hyp(k+1), pt_l, sa_l )
+       rho_l   = eqn_state_seawater_func( hyp(k+1), pt_l, sa_l )
 
-       hyp(k) = hyp(k+1) + rho_l * g * dzu(k+1)
+       hyp(k)  = hyp(k+1) + rho_l * g * dzu(k+1)
+       rho_ref = rho_ref + rho_l * dzu(k+1)
 
     ENDDO
 
+    rho_ref = rho_ref / ( zw(nzt) - zu(nzb) )
+    print*, '*** rho_ref = ', rho_ref
+
 
  END SUBROUTINE init_ocean