Changeset 2088 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Dec 19, 2016 4:30:25 PM (7 years ago)

Author:

suehring

Message:

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

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix in generic topography in case of ocean simulations
 ! 
 ! Former revisions:
@@ -574 +574 @@
           blx = NINT( building_length_x / dx )
           bly = NINT( building_length_y / dy )
-          bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
+          IF ( .NOT. ocean )  THEN
+             bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
+          ELSE
+             bh  = MINLOC( ABS( zw - zw(0) - building_height ), 1 ) - 1
+          ENDIF
+
           IF ( ABS( zw(bh  ) - building_height ) == &
                ABS( zw(bh+1) - building_height )    )  bh = bh + 1
@@ -640 +645 @@
           ENDIF
 
-          ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
+          IF ( .NOT. ocean )  THEN
+             ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
+          ELSE
+             ch  = MINLOC( ABS( zw - zw(0) - canyon_height ), 1 ) - 1
+          ENDIF
+
           IF ( ABS( zw(ch  ) - canyon_height ) == &
                ABS( zw(ch+1) - canyon_height )    )  ch = ch + 1