Ignore:
Timestamp:
Dec 19, 2016 4:30:25 PM (4 years ago)
Author:
suehring
Message:

Bugfixes in initial salinity profile and generic topography definition in case of ocean simulations

File:
1 edited

Legend:

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

    r2038 r2088  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix in generic topography in case of ocean simulations
    2323!
    2424! Former revisions:
     
    574574          blx = NINT( building_length_x / dx )
    575575          bly = NINT( building_length_y / dy )
    576           bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
     576          IF ( .NOT. ocean )  THEN
     577             bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
     578          ELSE
     579             bh  = MINLOC( ABS( zw - zw(0) - building_height ), 1 ) - 1
     580          ENDIF
     581
    577582          IF ( ABS( zw(bh  ) - building_height ) == &
    578583               ABS( zw(bh+1) - building_height )    )  bh = bh + 1
     
    640645          ENDIF
    641646
    642           ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
     647          IF ( .NOT. ocean )  THEN
     648             ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
     649          ELSE
     650             ch  = MINLOC( ABS( zw - zw(0) - canyon_height ), 1 ) - 1
     651          ENDIF
     652
    643653          IF ( ABS( zw(ch  ) - canyon_height ) == &
    644654               ABS( zw(ch+1) - canyon_height )    )  ch = ch + 1
Note: See TracChangeset for help on using the changeset viewer.