Ignore:
Timestamp:
Mar 27, 2012 3:10:33 PM (12 years ago)
Author:
gryschka
Message:

Bugifx in init_grid in case of ocean runs

File:
1 edited

Legend:

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

    r863 r864  
    44! Current revisions:
    55! -----------------
     6!
     7! In case of ocean and dirichlet bc for u and v at the ground
     8! dzu_mg and ddzu_pres where not defined correctly for k=1.
    69!
    710! Former revisions:
     
    246249       ALLOCATE( ddzu_pres(1:nzt+1) )
    247250       ddzu_pres = ddzu
    248        IF( .NOT. ocean )  ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
     251       ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
    249252    ENDIF   
    250253
     
    263266       dzu_mg(:,maximum_grid_level) = dzu
    264267!       
    265 !--    Next line to ensure an equally spaced grid. For ocean runs this is not
    266 !--    necessary,
    267 !--    because zu(0) does not changed so far. Also this would cause errors
    268 !--    if a vertical stretching for ocean runs is used.
    269        IF ( .NOT. ocean ) dzu_mg(1,maximum_grid_level) = dzu(2)
     268!--    Next line to ensure an equally spaced grid.
     269       dzu_mg(1,maximum_grid_level) = dzu(2)
    270270
    271271       dzw_mg(:,maximum_grid_level) = dzw
Note: See TracChangeset for help on using the changeset viewer.