Index: palm/trunk/SOURCE/init_grid.f90
===================================================================
--- palm/trunk/SOURCE/init_grid.f90	(revision 4168)
+++ palm/trunk/SOURCE/init_grid.f90	(revision 4180)
@@ -67,316 +67,4 @@
 ! Comment added
 ! 
-! 3294 2018-10-01 02:37:10Z raasch
-! ocean renamed ocean_mode
-! 
-! 3241 2018-09-12 15:02:00Z raasch
-! unused variables removed
-! 
-! 3200 2018-08-17 14:46:36Z suehring
-! Bugfix, missing pre-processor directive
-! 
-! 3183 2018-07-27 14:25:55Z suehring
-! Rename variables in mesoscale-offline nesting mode
-! 
-! 3182 2018-07-27 13:36:03Z suehring
-! Bugfix in referencing buildings on orography top
-! 
-! 3139 2018-07-17 11:30:10Z Giersch
-! Bugfix in case of restarts and grid stretching
-! 
-! 3115 2018-07-10 12:49:26Z suehring
-! Referencing of buildings onto top of terrain - special treatment for bridges. 
-! 
-! 3103 2018-07-04 17:30:52Z suehring
-! Reference lowest terrain height to zero level
-! 
-! 3068 2018-06-12 14:49:41Z Giersch
-! New warning message concerning grid stretching has been introduced
-! 
-! 3066 2018-06-12 08:55:55Z Giersch
-! Bugfix in IF statement before error message
-! 
-! 3065 2018-06-12 07:03:02Z Giersch
-! New vertical stretching mechanism introduced
-! 
-! 3051 2018-05-30 17:43:55Z suehring
-! Minor bugfix concerning mapping 3D buildings on top of terrain
-! 
-! 3045 2018-05-28 07:55:41Z Giersch
-! Error messages revised
-! 
-! 3045 2018-05-28 07:55:41Z Giersch
-! Error messages revised
-! 
-! 2968 2018-04-13 11:52:24Z suehring
-! Bugfix in initialization in case of elevated model surface. Introduce
-! index for minimum topography-top. 
-! 
-! 2955 2018-04-09 15:14:01Z suehring
-! Improve topography filter routine and add ghost-point exchange for building 
-! ID and building type. 
-! 
-! 2927 2018-03-23 15:13:00Z suehring
-! Bugfix, setting boundary conditions for topography index array.
-! 
-! 2918 2018-03-21 15:52:14Z gronemeier
-! Moved init_mixing_length to turbulence_closure_mod.f90
-! 
-! 2897 2018-03-15 11:47:16Z suehring
-! Relax restrictions for topography input, terrain and building heights can be
-! input separately and are not mandatory any more. 
-! 
-! 2893 2018-03-14 16:20:52Z suehring
-! Revise informative message concerning filtered topography (1 grid-point 
-! holes). 
-! 
-! 2892 2018-03-14 15:06:29Z suehring
-! Bugfix, uninitialized array in case of single_building.
-! 
-! 2867 2018-03-09 09:40:23Z suehring
-! Revise mapping of 3D buildings onto onto orography.
-! 
-! 2823 2018-02-20 15:31:45Z Giersch
-! Set boundary conditions for 3D topography in case of non-cyclic boundary 
-! conditions
-! 
-! 2796 2018-02-08 12:25:39Z suehring
-! Bugfix in 3D building initialization
-! 
-! 2747 2018-01-15 12:44:17Z suehring
-! Bugfix, topography height is rounded to the nearest discrete grid level
-! 
-! 2718 2018-01-02 08:49:38Z maronga
-! Corrected "Former revisions" section
-! 
-! 2701 2017-12-15 15:40:50Z suehring
-! Changes from last commit documented
-! 
-! 2698 2017-12-14 18:46:24Z suehring
-! Bugfix in get_topography_top_index
-!
-! 2696 2017-12-14 17:12:51Z kanani
-! Change in file header (GPL part)
-! Revised topography input
-! Set nzb_max not for the entire nest domain, only for boundary PEs
-! Re-organize routine, split-up into several subroutines
-! Modularize poismg_noopt
-! Remove setting bit 26, 27, 28 in wall_flags_0, indicating former '_outer'
-! arrays (not required any more).  
-! Bugfix in generic tunnel setup (MS)
-! 
-! 2550 2017-10-16 17:12:01Z boeske
-! Set lateral boundary conditions for topography on all three ghost layers
-! 
-! 2478 2017-09-18 13:37:24Z suehring
-! Bugfix, correct flag for use_top
-! 
-! 2365 2017-08-21 14:59:59Z kanani
-! Vertical nesting implemented (SadiqHuq)
-! 
-! 2319 2017-07-20 17:33:17Z suehring
-! Remove print statements
-! 
-! 2318 2017-07-20 17:27:44Z suehring
-! Get topography top index via Function call 
-! 
-! 2317 2017-07-20 17:27:19Z suehring
-! Bugfixes in reading 3D topography from file
-! 
-! 2274 2017-06-09 13:27:48Z Giersch
-! Changed error messages
-! 
-! 2233 2017-05-30 18:08:54Z suehring
-!
-! 2232 2017-05-30 17:47:52Z suehring
-! - Adjustments according to new topography representation
-! - Bugfix: Move determination of nzb_max behind topography modification in 
-!   cell-edge case 
-! - Get rid off global arrays required for topography output
-! - Enable topography input via netcdf 
-! - Generic tunnel set-up added
-! 
-! 2200 2017-04-11 11:37:51Z suehring
-! monotonic_adjustment removed
-! 
-! 2169 2017-03-06 18:16:35Z suehring
-! Bugfix, move setting for topography grid convention to init_grid, else, if no
-! value is set, the simulation may abort in case of restarts
-! 
-! 2128 2017-01-23 15:00:03Z suehring
-! Bugfix in setting topography from file in case of ocean simulations
-! 
-! 2088 2016-12-19 16:30:25Z suehring
-! Bugfix in generic topography in case of ocean simulations
-! 
-! 2037 2016-10-26 11:15:40Z knoop
-! Anelastic approximation implemented
-! 
-! 2021 2016-10-07 14:08:57Z suehring
-! Bugfix: setting Neumann boundary conditions for topography required for 
-! topography flags in multigrid_noopt solver
-! 
-! 2000 2016-08-20 18:09:15Z knoop
-! Forced header and separation lines into 80 columns
-! 
-! 1994 2016-08-15 09:52:21Z suehring
-! Bugfix in definition of generic topography
-! 
-! 1982 2016-08-01 11:04:48Z suehring
-! Bugfix concering consistency check for topography
-! 
-! 1968 2016-07-18 12:01:49Z suehring
-! Changed: PE-wise reading of topography file in order to avoid global definition
-! of arrays nzb_local and nzb_tmp. Thereby, topography definition for single
-! buildings and street canyons has changed, as well as flag setting for 
-! multigrid scheme.
-!
-! Bugfix in checking l_grid anisotropy. 
-! Simplify initial computation of lwall and vertical_influence, i.e. remove 
-! nzb_s_inner as it is still zero at this point.
-! 
-! 1942 2016-06-14 12:18:18Z suehring
-! Topography filter implemented to fill holes resolved by only one grid point.
-! Initialization of flags for ws-scheme moved to advec_ws.  
-! 
-! 1931 2016-06-10 12:06:59Z suehring
-! Rename multigrid into multigrid_noopt and multigrid_fast into multigrid
-!
-! 1910 2016-05-26 06:49:46Z raasch
-! Bugfix: if topography is read from file, Neumann conditions are used for the
-! nzb_local array (instead of cyclic conditions) in case that non-cyclic
-! boundary conditions are switched on for the run
-!
-! 1902 2016-05-09 11:18:56Z suehring
-! Set topography flags for multigrid solver only (not for multigrid_fast)
-!
-! 1886 2016-04-21 11:20:47Z suehring
-! Bugfix: setting advection flags near walls
-! reformulated index values for nzb_v_inner
-! variable discriptions added in declaration block
-!
-! 1845 2016-04-08 08:29:13Z raasch
-! nzb_2d removed
-! 
-! 1804 2016-04-05 16:30:18Z maronga
-! Removed code for parameter file check (__check)
-!
-! 1779 2016-03-03 08:01:28Z raasch
-! coupling_char is trimmed at every place it occurs, because it can have
-! different length now
-!
-! 1762 2016-02-25 12:31:13Z hellstea
-! Introduction of nested domain feature
-!
-! 1743 2016-01-13 10:23:51Z raasch
-! Bugfix for calculation of nzb_s_outer and nzb_u_outer at north boundary of
-! total domain
-!
-! 1691 2015-10-26 16:17:44Z maronga
-! Renamed prandtl_layer to constant_flux_layer.
-! 
-! 1682 2015-10-07 23:56:08Z knoop
-! Code annotations made doxygen readable 
-! 
-! 1677 2015-10-02 13:25:23Z boeske
-! Bugfix: Ghost points are included in wall_flags_0 and wall_flags_00
-! 
-! 1675 2015-10-02 08:28:59Z gronemeier
-! Bugfix: Definition of topography grid levels
-! 
-! 1660 2015-09-21 08:15:16Z gronemeier
-! Bugfix: Definition of topography grid levels if vertical grid stretching
-!         starts below the maximum topography height.
-!
-! 1580 2015-04-10 13:43:49Z suehring
-! Bugfix: setting flags for 5th order scheme near buildings
-!
-! 1575 2015-03-27 09:56:27Z raasch
-! adjustments for psolver-queries
-!
-! 1557 2015-03-05 16:43:04Z suehring
-! Adjustment for monotoinic limiter
-!
-! 1418 2014-06-06 13:05:08Z fricke
-! Bugfix: Change if-condition for stretched grid in the ocean, with the old
-!          condition and a negative value for dz_stretch_level the condition
-!          was always true for the whole model domain
-!
-! 1409 2014-05-23 12:11:32Z suehring
-! Bugfix: set wall_flags_0 at inflow and outflow boundary also for i <= nxlu 
-! j <= nysv
-! 
-! 1353 2014-04-08 15:21:23Z heinze
-! REAL constants provided with KIND-attribute 
-!
-! 1322 2014-03-20 16:38:49Z raasch
-! REAL constants defined as wp-kind
-!
-! 1320 2014-03-20 08:40:49Z raasch
-! ONLY-attribute added to USE-statements,
-! kind-parameters added to all INTEGER and REAL declaration statements, 
-! kinds are defined in new module kinds, 
-! revision history before 2012 removed,
-! comment fields (!:) to be used for variable explanations added to
-! all variable declaration statements 
-! 
-! 1221 2013-09-10 08:59:13Z raasch
-! wall_flags_00 introduced to hold bits 32-63,
-! additional 3D-flag arrays for replacing the 2D-index array nzb_s_inner in
-! loops optimized for openACC (pres + flow_statistics)
-!
-! 1092 2013-02-02 11:24:22Z raasch
-! unused variables removed
-!
-! 1069 2012-11-28 16:18:43Z maronga
-! bugfix: added coupling_char to TOPOGRAPHY_DATA to allow topography in the
-!         ocean model in case of coupled runs
-!
-! 1036 2012-10-22 13:43:42Z raasch
-! code put under GPL (PALM 3.9)
-!
-! 1015 2012-09-27 09:23:24Z raasch
-! lower index for calculating wall_flags_0 set to nzb_w_inner instead of
-! nzb_w_inner+1
-!
-! 996 2012-09-07 10:41:47Z raasch
-! little reformatting
-!
-! 978 2012-08-09 08:28:32Z fricke
-! Bugfix: nzb_max is set to nzt at non-cyclic lateral boundaries
-! Bugfix: Set wall_flags_0 for inflow boundary
-!
-! 927 2012-06-06 19:15:04Z raasch
-! Wall flags are not set for multigrid method in case of masking method
-!
-! 864 2012-03-27 15:10:33Z gryschka 
-! In case of ocean and Dirichlet bottom bc for u and v dzu_mg and ddzu_pres
-! were not correctly defined for k=1.
-!
-! 861 2012-03-26 14:18:34Z suehring
-! Set wall_flags_0. The array is needed for degradation in ws-scheme near walls,
-! inflow and outflow boundaries as well as near the bottom and the top of the
-! model domain.!
-! Initialization of nzb_s_inner and nzb_w_inner.
-! gls has to be at least nbgp to do not exceed the array bounds of nzb_local
-! while setting wall_flags_0
-!
-! 843 2012-02-29 15:16:21Z gryschka
-! In case of ocean and dirichlet bc for u and v at the bottom
-! the first u-level ist defined at same height as the first w-level
-!
-! 818 2012-02-08 16:11:23Z maronga
-! Bugfix: topo_height is only required if topography is used. It is thus now
-! allocated in the topography branch
-!
-! 809 2012-01-30 13:32:58Z maronga
-! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
-!
-! 807 2012-01-25 11:53:51Z maronga
-! New cpp directive "__check" implemented which is used by check_namelist_files
-!
-! Revision 1.1  1997/08/11 06:17:45  raasch
-! Initial revision (Testversion)
-!
 !
 ! Description:
