Changeset 864

Timestamp:

Mar 27, 2012 3:10:33 PM (13 years ago)

Author:

gryschka

Message:

Bugifx in init_grid in case of ocean runs

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+!
+! In case of ocean and dirichlet bc for u and v at the ground
+! dzu_mg and ddzu_pres where not defined correctly for k=1.
 !
 ! Former revisions:
@@ -246 +249 @@
        ALLOCATE( ddzu_pres(1:nzt+1) )
        ddzu_pres = ddzu
-       IF( .NOT. ocean )  ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
+       ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
     ENDIF   
 
@@ -263 +266 @@
        dzu_mg(:,maximum_grid_level) = dzu
 !       
-!--    Next line to ensure an equally spaced grid. For ocean runs this is not
-!--    necessary, 
-!--    because zu(0) does not changed so far. Also this would cause errors 
-!--    if a vertical stretching for ocean runs is used.
-       IF ( .NOT. ocean ) dzu_mg(1,maximum_grid_level) = dzu(2)
+!--    Next line to ensure an equally spaced grid. 
+       dzu_mg(1,maximum_grid_level) = dzu(2)
 
        dzw_mg(:,maximum_grid_level) = dzw