Index: palm/trunk/SOURCE/advec_ws.f90
===================================================================
--- palm/trunk/SOURCE/advec_ws.f90	(revision 4110)
+++ palm/trunk/SOURCE/advec_ws.f90	(revision 4180)
@@ -77,218 +77,4 @@
 ! OpenACC port for SPEC
 ! 
-! 3589 2018-11-30 15:09:51Z suehring
-! Move the control parameter "salsa" from salsa_mod to control_parameters 
-! (M. Kurppa)
-! 
-! 3582 2018-11-29 19:16:36Z suehring
-! - Computation of vertical fluxes separated from computation of horizontal 
-!   fluxes. Loops are splitted in order to avoid indirect indexing and allow 
-!   for better vectorization. 
-! - Accelerate code by remove type-conversions of ibits
-! - replace pointer definition in scalar routine by simple explicit definition
-!   of the passed array
-! 
-! 3547 2018-11-21 13:21:24Z suehring
-! variables documented
-! 
-! 3467 2018-10-30 19:05:21Z suehring
-! Implementation of a new aerosol module salsa.
-! 
-! 3302 2018-10-03 02:39:40Z raasch
-! Stokes drift velocity added in scalar advection
-! 
-! 3298 2018-10-02 12:21:11Z kanani
-! Add formatted note from ketelsen about statistics for chemistry
-!
-! 3294 2018-10-01 02:37:10Z raasch
-! ocean renamed ocean_mode
-! 
-! 3274 2018-09-24 15:42:55Z knoop
-! Modularization of all bulk cloud physics code components
-! 
-! 3241 2018-09-12 15:02:00Z raasch
-! unused variables removed
-! 
-! 3183 2018-07-27 14:25:55Z suehring
-! Rename variables for boundary-condition flags and for offline nesting mode 
-! 
-! 3182 2018-07-27 13:36:03Z suehring
-! Bugfix in calculation of left-sided fluxes for u-component in OpenMP case.
-! 
-! 2731 2018-01-09 17:44:02Z suehring
-! Enable loop vectorization by splitting the k-loop
-!
-! 2718 2018-01-02 08:49:38Z maronga
-! Corrected "Former revisions" section
-! 
-! 2696 2017-12-14 17:12:51Z kanani
-! Change in file header (GPL part)
-! Implement advection for TKE-dissipation in case of RANS-TKE-e closure (TG)
-! Allocate advc_flags_1/2 within ws_init_flags instead of init_grid
-! Change argument list for exchange_horiz_2d_int (MS)
-! 
-! 2329 2017-08-03 14:24:56Z knoop
-! Bugfix concerning density in divergence correction close to buildings
-! 
-! 2292 2017-06-20 09:51:42Z schwenkel
-! Implementation of new microphysic scheme: cloud_scheme = 'morrison' 
-! includes two more prognostic equations for cloud drop concentration (nc)  
-! and cloud water content (qc). 
-! 
-! 2233 2017-05-30 18:08:54Z suehring
-! 
-! 2232 2017-05-30 17:47:52Z suehring
-! Rename wall_flags_0 and wall_flags_00 into advc_flags_m and advc_flags_1, 
-! respectively.
-! Set advc_flags_1/2 on basis of wall_flags_0/00 instead of nzb_s/u/v/w_inner.
-! Setting advc_flags_1/2 also for downward-facing walls
-! 
-! 2200 2017-04-11 11:37:51Z suehring
-! monotonic_adjustment removed
-! 
-! 2118 2017-01-17 16:38:49Z raasch
-! OpenACC version of subroutines removed
-! 
-! 2037 2016-10-26 11:15:40Z knoop
-! Anelastic approximation implemented
-! 
-! 2000 2016-08-20 18:09:15Z knoop
-! Forced header and separation lines into 80 columns
-! 
-! 1996 2016-08-18 11:42:29Z suehring
-! Bugfix concerning calculation of turbulent of turbulent fluxes
-! 
-! 1960 2016-07-12 16:34:24Z suehring
-! Separate humidity and passive scalar
-! 
-! 1942 2016-06-14 12:18:18Z suehring
-! Initialization of flags for ws-scheme moved from init_grid.
-! 
-! 1873 2016-04-18 14:50:06Z maronga
-! Module renamed (removed _mod)
-! 
-! 
-! 1850 2016-04-08 13:29:27Z maronga
-! Module renamed
-!
-! 
-! 1822 2016-04-07 07:49:42Z hoffmann
-! icloud_scheme removed, microphysics_seifert added
-!
-! 1682 2015-10-07 23:56:08Z knoop
-! Code annotations made doxygen readable 
-! 
-! 1630 2015-08-26 16:57:23Z suehring
-! 
-! 
-! 1629 2015-08-26 16:56:11Z suehring
-! Bugfix concerning wall_flags at left and south PE boundaries
-!
-! 1581 2015-04-10 13:45:59Z suehring
-! 
-! 
-! 1580 2015-04-10 13:43:49Z suehring
-! Bugfix: statistical evaluation of scalar fluxes in case of monotonic limiter
-!
-! 1567 2015-03-10 17:57:55Z suehring
-! Bugfixes in monotonic limiter.
-!
-! 2015-03-09 13:10:37Z heinze
-! Bugfix: REAL constants provided with KIND-attribute in call of 
-! intrinsic functions like MAX and MIN 
-! 
-! 1557 2015-03-05 16:43:04Z suehring
-! Enable monotone advection for scalars using monotonic limiter
-!
-! 1374 2014-04-25 12:55:07Z raasch
-! missing variables added to ONLY list
-! 
-! 1361 2014-04-16 15:17:48Z hoffmann
-! accelerator and vector version for qr and nr added
-! 
-! 1353 2014-04-08 15:21:23Z heinze
-! REAL constants provided with KIND-attribute,
-! module kinds added
-! some formatting adjustments
-!
-! 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,
-! old module precision_kind is removed,
-! revision history before 2012 removed,
-! comment fields (!:) to be used for variable explanations added to
-! all variable declaration statements
-!
-! 1257 2013-11-08 15:18:40Z raasch
-! accelerator loop directives removed
-!
-! 1221 2013-09-10 08:59:13Z raasch
-! wall_flags_00 introduced, which holds bits 32-...
-!
-! 1128 2013-04-12 06:19:32Z raasch
-! loop index bounds in accelerator version replaced by i_left, i_right, j_south,
-! j_north
-!
-! 1115 2013-03-26 18:16:16Z hoffmann
-! calculation of qr and nr is restricted to precipitation
-!
-! 1053 2012-11-13 17:11:03Z hoffmann
-! necessary expansions according to the two new prognostic equations (nr, qr) 
-! of the two-moment cloud physics scheme:
-! +flux_l_*, flux_s_*, diss_l_*, diss_s_*, sums_ws*s_ws_l 
-!
-! 1036 2012-10-22 13:43:42Z raasch
-! code put under GPL (PALM 3.9)
-!
-! 1027 2012-10-15 17:18:39Z suehring
-! Bugfix in calculation indices k_mm, k_pp in accelerator version
-!
-! 1019 2012-09-28 06:46:45Z raasch
-! small change in comment lines
-!
-! 1015 2012-09-27 09:23:24Z raasch
-! accelerator versions (*_acc) added
-!
-! 1010 2012-09-20 07:59:54Z raasch
-! cpp switch __nopointer added for pointer free version
-!
-! 888 2012-04-20 15:03:46Z suehring
-! Number of REAL( IBITS() calls with identical arguments is reduced.
-!
-! 862 2012-03-26 14:21:38Z suehring
-! ws-scheme also work with topography in combination with vector version.
-! ws-scheme also work with outflow boundaries in combination with
-! vector version.
-! Degradation of the applied order of scheme is now steered by multiplying with
-! Integer advc_flags_1. 2nd order scheme, WS3 and WS5 are calculated on each
-! grid point and mulitplied with the appropriate flag.
-! 2nd order numerical dissipation term changed. Now the appropriate 2nd order
-! term derived according to the 4th and 6th order terms is applied. It turns
-! out that diss_2nd does not provide sufficient dissipation near walls.
-! Therefore, the function diss_2nd is removed.
-! Near walls a divergence correction is necessary to overcome numerical
-! instabilities due to too less divergence reduction of the velocity field.
-! boundary_flags and logicals steering the degradation are removed.
-! Empty SUBROUTINE local_diss removed.
-! Further formatting adjustments.
-!
-! 801 2012-01-10 17:30:36Z suehring
-! Bugfix concerning OpenMP parallelization. Summation of sums_wsus_ws_l,
-! sums_wsvs_ws_l, sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wspts_ws_l,
-! sums_wsqs_ws_l, sums_wssas_ws_l is now thread-safe by adding an additional
-! dimension.
-!
-! Initial revision
-!
-! 411 2009-12-11 12:31:43 Z suehring
-!
-!
-!
-! @author Matthias Suehring
-!
 !
 ! Description:
