Ignore:
Timestamp:
Jun 2, 2007 4:48:38 PM (14 years ago)
Author:
raasch
Message:

further preliminary uncomplete changes for ocean version

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_cloud_physics.f90

    r4 r95  
    44! Actual revisions:
    55! -----------------
    6 !
     6! hydro_press renamed hyp
    77!
    88! Former revisions:
     
    3535    REAL    ::  t_surface
    3636
    37     ALLOCATE( hydro_press(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
     37    ALLOCATE( hyp(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
    3838
    3939!
     
    5454!-- pt / t : ratio of potential and actual temperature (pt_d_t)
    5555!-- t / pt : ratio of actual and potential temperature (t_d_pt)
    56 !-- p_0(z) : vertical profile of the hydrostatic pressure (hydro_press)
     56!-- p_0(z) : vertical profile of the hydrostatic pressure (hyp)
    5757    t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
    5858    DO  k = nzb, nzt+1
    59        hydro_press(k) = surface_pressure * 100.0 * &
    60                         ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
    61        pt_d_t(k)      = ( 100000.0 / hydro_press(k) )**0.286
    62        t_d_pt(k)      = 1.0 / pt_d_t(k)       
     59       hyp(k)    = surface_pressure * 100.0 * &
     60                   ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
     61       pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
     62       t_d_pt(k) = 1.0 / pt_d_t(k)       
    6363    ENDDO
    6464
Note: See TracChangeset for help on using the changeset viewer.