Ignore:
Timestamp:
Mar 6, 2007 12:31:10 PM (17 years ago)
Author:
raasch
Message:

one change missed in previous commit

File:
1 edited

Legend:

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

    r48 r49  
    370370!
    371371!-- Consistency checks and index array initialization are only required for
    372 !-- non-flat topography
     372!-- non-flat topography, also the initialization of topography heigth arrays
     373!-- zu_s_inner and zw_w_inner
    373374    IF ( TRIM( topography ) /= 'flat' )  THEN
    374375
     
    532533       CALL exchange_horiz_2d_int( nzb_w_outer )
    533534       CALL exchange_horiz_2d_int( nzb_s_outer )
     535
     536!
     537!--    Allocate and set the arrays containing the topography height
     538       IF ( myid == 0 )  THEN
     539
     540          ALLOCATE( zu_s_inner(0:nx+1,0:ny+1), zw_w_inner(0:nx+1,0:ny+1) )
     541
     542          DO  i = 0, nx + 1
     543             DO  j = 0, ny + 1
     544                zu_s_inner(i,j) = zu(nzb_local(j,i))
     545                zw_w_inner(i,j) = zw(nzb_local(j,i))
     546             ENDDO
     547          ENDDO
     548         
     549       ENDIF
    534550
    535551    ENDIF
Note: See TracChangeset for help on using the changeset viewer.