Changeset 2550 for palm/trunk/SOURCE/header.f90

Timestamp:

Oct 16, 2017 5:12:01 PM (7 years ago)

Author:

boeske

Message:

enable simulations with complex terrain

File:

• palm/trunk/SOURCE/header.f90 (modified) (4 diffs)

palm/trunk/SOURCE/header.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added output for complex terrain simulations
+! 
+! 2544 2017-10-13 18:09:32Z maronga
 ! Moved initial day of year and time to inipar. 
 ! 
@@ -403 +406 @@
 
     USE surface_mod,                                                           &
-        ONLY:  surf_def_h
+        ONLY:  surf_def_h, get_topography_top_index
 
     USE synthetic_turbulence_generator_mod,                                    &
@@ -972 +975 @@
        ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' )  THEN
           WRITE ( io, 279 )
+       ENDIF
+    ENDIF
+
+!-- Complex terrain
+    IF ( complex_terrain )  THEN
+       WRITE( io, 280 ) 
+       IF ( turbulent_inflow )  THEN
+          WRITE( io, 281 )  zu( get_topography_top_index( 0, 0, 's' ) )
+       ENDIF
+       IF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
+          WRITE( io, 282 )
        ENDIF
     ENDIF
@@ -2051 +2065 @@
 279 FORMAT (' Topography grid definition convention:'/ &
             ' cell center (scalar grid points)' /)
+280 FORMAT (' Complex terrain simulation is activated.')
+281 FORMAT ('    --> Mean inflow profiles are adjusted.' / &
+            '    --> Elevation of inflow boundary: ', F7.1, ' m' )
+282 FORMAT ('    --> Initial data from 3D-precursor run is shifted' / &
+            '        vertically depending on local surface height.')
 300 FORMAT (//' Boundary conditions:'/ &
              ' -------------------'// &