Ignore:
Timestamp:
Mar 30, 2011 9:31:40 AM (10 years ago)
Author:
raasch
Message:

formatting adjustments

File:
1 edited

Legend:

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

    r708 r709  
    44! Current revisions:
    55! -----------------
    6 !
     6! formatting adjustments
    77!
    88! Former revisions:
     
    8989   
    9090!
    91 !   Computation of the array bounds.
     91!-- Calculation of horizontal array bounds including ghost layers
    9292    nxlg = nxl - nbgp
    9393    nxrg = nxr + nbgp
    9494    nysg = nys - nbgp
    9595    nyng = nyn + nbgp
     96
    9697!
    9798!-- Allocate grid arrays
     
    196197   
    197198!   
    198 !-- In case of FFT method or SOR swap inverse grid lenght ddzu to ddzu_fft and
    199 !-- modify the lowest entry. This is necessary for atmosphere runs, because
    200 !-- zu(0) and so ddzu(1) changed. Accompanied with this modified arrays
    201 !-- pressure solver uses wrong values and this causes kinks in the profiles
    202 !-- of turbulent quantities. 
    203     IF ( psolver /= 'multigrid' ) THEN
     199!-- The FFT- SOR-pressure solvers assume grid spacings of a staggered grid
     200!-- everywhere. For the actual grid, the grid spacing at the lowest level
     201!-- is only dz/2, but should be dz. Therefore, an additional array
     202!-- containing with appropriate grid information is created for these
     203!-- solvers.
     204    IF ( psolver /= 'multigrid' )  THEN
    204205       ALLOCATE( ddzu_pres(1:nzt+1) )
    205206       ddzu_pres = ddzu
    206        IF( .NOT. ocean ) ddzu_pres(1) = ddzu_pres(2)
     207       IF( .NOT. ocean )  ddzu_pres(1) = ddzu_pres(2)  ! change for lowest level
    207208    ENDIF   
    208209
     
    221222       dzu_mg(:,maximum_grid_level) = dzu
    222223!       
    223 !--    To ensure a equally spaced grid. For ocean runs this is not necessary,
     224!--    Next line to ensure an equally spaced grid. For ocean runs this is not
     225!--    necessary,
    224226!--    because zu(0) does not changed so far. Also this would cause errors
    225227!--    if a vertical stretching for ocean runs is used.
    226228       IF ( .NOT. ocean ) dzu_mg(1,maximum_grid_level) = dzu(2)
     229
    227230       dzw_mg(:,maximum_grid_level) = dzw
    228231       nzt_l = nzt
     
    11321135
    11331136!
    1134 !-- Need to set lateral boundary conditions for l_wall
    1135 
     1137!-- Set lateral boundary conditions for l_wall
    11361138    CALL exchange_horiz( l_wall, nbgp )
    11371139
Note: See TracChangeset for help on using the changeset viewer.