Changeset 49 for palm

Timestamp:

Mar 6, 2007 12:31:10 PM (18 years ago)

Author:

raasch

Message:

one change missed in previous commit

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -370 +370 @@
 !
 !-- Consistency checks and index array initialization are only required for
-!-- non-flat topography
+!-- non-flat topography, also the initialization of topography heigth arrays
+!-- zu_s_inner and zw_w_inner
     IF ( TRIM( topography ) /= 'flat' )  THEN
 
@@ -532 +533 @@
        CALL exchange_horiz_2d_int( nzb_w_outer )
        CALL exchange_horiz_2d_int( nzb_s_outer )
+
+!
+!--    Allocate and set the arrays containing the topography height
+       IF ( myid == 0 )  THEN
+
+          ALLOCATE( zu_s_inner(0:nx+1,0:ny+1), zw_w_inner(0:nx+1,0:ny+1) )
+
+          DO  i = 0, nx + 1
+             DO  j = 0, ny + 1
+                zu_s_inner(i,j) = zu(nzb_local(j,i))
+                zw_w_inner(i,j) = zw(nzb_local(j,i))
+             ENDDO
+          ENDDO
+          
+       ENDIF
 
     ENDIF