Changeset 95 for palm/trunk/SOURCE/init_cloud_physics.f90

Timestamp:

Jun 2, 2007 4:48:38 PM (17 years ago)

Author:

raasch

Message:

further preliminary uncomplete changes for ocean version

File:

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

palm/trunk/SOURCE/init_cloud_physics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! hydro_press renamed hyp
 !
 ! Former revisions:
@@ -35 +35 @@
     REAL    ::  t_surface
 
-    ALLOCATE( hydro_press(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
+    ALLOCATE( hyp(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
 
 !
@@ -54 +54 @@
 !-- pt / t : ratio of potential and actual temperature (pt_d_t)
 !-- t / pt : ratio of actual and potential temperature (t_d_pt)
-!-- p_0(z) : vertical profile of the hydrostatic pressure (hydro_press)
+!-- p_0(z) : vertical profile of the hydrostatic pressure (hyp)
     t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
     DO  k = nzb, nzt+1
-       hydro_press(k) = surface_pressure * 100.0 * &
-                        ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
-       pt_d_t(k)      = ( 100000.0 / hydro_press(k) )**0.286
-       t_d_pt(k)      = 1.0 / pt_d_t(k)       
+       hyp(k)    = surface_pressure * 100.0 * &
+                   ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
+       pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
+       t_d_pt(k) = 1.0 / pt_d_t(k)       
     ENDDO