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

Timestamp:

Sep 23, 2009 9:40:33 AM (15 years ago)

Author:

raasch

Message:

in-situ AND potential density are calculated and used in the ocean version

File:

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

palm/trunk/SOURCE/init_ocean.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Bugfix: First calculation of hyp(0) changed
+! Bugfix: Initial profiles of hydrostatic pressure and density are calculated
+! iteratively. First calculation of hyp(0) changed.
 !
 ! Former revisions:
@@ -32 +33 @@
     INTEGER ::  k, n
 
-    REAL    ::  sa_l, pt_l, rho_l
+    REAL    ::  sa_l, pt_l
 
     REAL, DIMENSION(nzb:nzt+1) ::  rho_init
@@ -45 +46 @@
 !-- Calculate initial vertical profile of hydrostatic pressure (in Pa)
 !-- and the reference density (used later in buoyancy term)
+!-- First step: Calculate pressure using reference density
     hyp(nzt+1) = surface_pressure * 100.0
 
@@ -55 +57 @@
     hyp(0) = hyp(1) + rho_surface * g * dzu(1)
 
-    IF ( myid == 0 )  THEN
-       print*,'hydro pres using rho_surface'
-       DO  k = nzt+1, 0, -1
-          print*, 'k = ', k, ' hyp = ', hyp(k)
-       ENDDO
-       print*, ' '
-    ENDIF
-
+!
+!-- Second step: Iteratively calculate in situ density (based on presssure)
+!-- and pressure (based on in situ density)
     DO  n = 1, 5
 
@@ -85 +82 @@
        ENDDO
 
-       IF ( myid == 0 )  THEN
-          print*,'hydro pres / rho  n = ', n
-          DO  k = nzt+1, 0, -1
-             print*, 'k = ', k, ' hyp = ', hyp(k), ' rho = ', rho_init(k)
-          ENDDO
-          print*, ' '
-       ENDIF
-
     ENDDO
 
@@ -111 +100 @@
 
 !
-!-- Calculate the initial potential density, based on the initial
-!-- temperature and salinity profile
+!-- Calculate the 3d array of initial in situ and potential density,
+!-- based on the initial temperature and salinity profile
     CALL eqn_state_seawater