Ignore:
Timestamp:
Aug 25, 2009 8:06:27 AM (12 years ago)
Author:
raasch
Message:

speed optomizations +bugfix in init_ocean

File:
1 edited

Legend:

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

    r336 r366  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: First calculation of hyp(0) changed
    77!
    88! Former revisions:
     
    5050    rho_init(nzt) = rho_surface
    5151
    52     DO  k = nzt-1, 0, -1
     52    DO  k = nzt-1, 1, -1
    5353       hyp(k) = hyp(k+1) + rho_surface * g * dzu(k)
    5454    ENDDO
     55    hyp(0) = hyp(1) + rho_surface * g * dzu(1)
    5556
    5657    IF ( myid == 0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.