Index: palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
===================================================================
--- palm/trunk/SOURCE/surface_layer_fluxes_mod.f90	(revision 3987)
+++ palm/trunk/SOURCE/surface_layer_fluxes_mod.f90	(revision 4180)
@@ -54,200 +54,4 @@
 ! OpenACC port for SPEC
 ! 
-! 3597 2018-12-04 08:40:18Z maronga
-! Added routine for calculating near surface air potential temperature (moved
-! from urban_surface_mod)
-! 
-! 3547 2018-11-21 13:21:24Z suehring
-! variables documented
-! 
-! 3361 2018-10-16 20:39:37Z knoop
-! Modularization of all bulk cloud physics code components
-! 
-! 3271 2018-09-24 08:20:34Z suehring
-! Comment revised
-! 
-! 3157 2018-07-19 21:08:49Z maronga
-! Added local free convection velocity scale w_lfc in calculation of uvw_abs.
-! This can be switche on/off by the user via the flag namelist parameter 
-! use_free_convection_scaling
-! 
-! 3152 2018-07-19 13:26:52Z suehring
-! q_surface is now part of surface structure
-! 
-! 3149 2018-07-19 05:48:09Z maronga
-! Correct accidental last commit
-! 
-! 3148 2018-07-19 05:45:25Z maronga
-! Major bugfix in calculation of Obukhov length
-! 
-! 3146 2018-07-18 22:36:19Z maronga
-! move phi_m from turbulence_closure_mod
-! 
-! 3045 2018-05-28 07:55:41Z Giersch
-! Error message revised
-! 
-! 2766 2018-01-22 17:17:47Z kanani
-! Removed preprocessor directive __chem
-! 
-! 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)
-! - Implementation of chemistry module (FK)
-! - Added calculation of pt1 and qv1 for all surface types. Added calculation of
-!   pt_surface for default-type surfaces (BM)
-! - Add flag to disable computation of qsws in case of urban surface (MS)
-! 
-! 2547 2017-10-16 12:41:56Z schwenkel
-! extended by cloud_droplets option
-! 
-! 2321 2017-07-24 15:57:07Z schwenkel
-! Bugfix: Correct index in lookup table for Obukhov length 
-! 
-! 2299 2017-06-29 10:14:38Z suehring
-! Adjusted for allow separate spinups of LSM and atmosphere code
-! 
-! 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). 
-! 
-! 2281 2017-06-13 11:34:50Z suehring
-! Clean-up unnecessary index access to surface type
-! 
-! 2233 2017-05-30 18:08:54Z suehring
-!
-! 2232 2017-05-30 17:47:52Z suehring
-! Adjustments to new surface concept
-! OpenMP bugfix
-! 
-! 2118 2017-01-17 16:38:49Z raasch
-! OpenACC directives and related code removed
-! 
-! 2091 2016-12-21 16:38:18Z suehring
-! Bugfix in calculation of vsws ( incorrect linear interpolation of us ) 
-! 
-! 2076 2016-12-02 13:54:20Z raasch
-! further openmp bugfix for lookup method
-! 
-! 2073 2016-11-30 14:34:05Z raasch
-! openmp bugfix for lookup method
-! 
-! 2037 2016-10-26 11:15:40Z knoop
-! Anelastic approximation implemented
-! 
-! 2011 2016-09-19 17:29:57Z kanani
-! Flag urban_surface is now defined in module control_parameters.
-! 
-! 2007 2016-08-24 15:47:17Z kanani
-! Account for urban surface model in computation of vertical kinematic heatflux
-! 
-! 2000 2016-08-20 18:09:15Z knoop
-! Forced header and separation lines into 80 columns
-! 
-! 1992 2016-08-12 15:14:59Z suehring
-! Minor bug, declaration of look-up index as INTEGER
-! 
-! 1960 2016-07-12 16:34:24Z suehring
-! Treat humidity and passive scalar separately
-! 
-! 1929 2016-06-09 16:25:25Z suehring
-! Bugfix: avoid segmentation fault in case one grid point is horizontally 
-! completely surrounded by topography
-! 
-! 1920 2016-05-30 10:50:15Z suehring
-! Avoid segmentation fault (see change in 1915) by different initialization of 
-! us instead of adding a very small number in the denominator
-!
-! 1915 2016-05-27 11:05:02Z suehring
-! Bugfix: avoid segmentation fault in case of most_method = 'circular' at first
-! timestep
-!
-! 1850 2016-04-08 13:29:27Z maronga
-! Module renamed
-!
-! 
-! 1822 2016-04-07 07:49:42Z hoffmann
-! icloud_scheme replaced by microphysics_*
-!
-! 1788 2016-03-10 11:01:04Z maronga
-! Added parameter z0q which replaces z0h in the similarity functions for
-! humidity.
-! Syntax layout improved.
-! 
-! 1757 2016-02-22 15:49:32Z maronga
-! Minor fixes.
-!
-! 1749 2016-02-09 12:19:56Z raasch
-! further OpenACC adjustments
-!
-! 1747 2016-02-08 12:25:53Z raasch
-! adjustments for OpenACC usage
-!
-! 1709 2015-11-04 14:47:01Z maronga
-! Bugfix: division by zero could occur when calculating rib at low wind speeds
-! Bugfix: calculation of uv_total for neutral = .T., initial value for ol for
-! neutral = .T.
-! 
-! 1705 2015-11-02 14:28:56Z maronga
-! Typo removed
-!
-! 1697 2015-10-28 17:14:10Z raasch
-! FORTRAN and OpenMP errors removed
-!
-! 1696 2015-10-27 10:03:34Z maronga
-! Modularized and completely re-written version of prandtl_fluxes.f90. In the
-! course of the re-writing two additional methods have been implemented. See
-! updated description.
-!
-! 1551 2015-03-03 14:18:16Z maronga
-! Removed land surface model part. The surface fluxes are now always calculated
-! within prandtl_fluxes, based on the given surface temperature/humidity (which 
-! is either provided by the land surface model, by large scale forcing data, or
-! directly prescribed by the user.
-! 
-! 1496 2014-12-02 17:25:50Z maronga
-! Adapted for land surface model
-! 
-! 1494 2014-11-21 17:14:03Z maronga
-! Bugfixes: qs is now calculated before calculation of Rif. calculation of
-! buoyancy flux in Rif corrected (added missing humidity term), allow use of 
-! topography for coupled runs (not tested)
-! 
-! 1361 2014-04-16 15:17:48Z hoffmann
-! Bugfix: calculation of turbulent fluxes of rain water content (qrsws) and rain 
-! drop concentration (nrsws) added
-! 
-! 1340 2014-03-25 19:45:13Z kanani
-! 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
-!
-! 1276 2014-01-15 13:40:41Z heinze
-! Use LSF_DATA also in case of Dirichlet bottom boundary condition for scalars
-!
-! 1257 2013-11-08 15:18:40Z raasch
-! openACC "kernels do" replaced by "kernels loop", "loop independent" added
-!
-! 1036 2012-10-22 13:43:42Z raasch
-! code put under GPL (PALM 3.9)
-!
-! 1015 2012-09-27 09:23:24Z raasch
-! OpenACC statements added
-!
-! 978 2012-08-09 08:28:32Z fricke
-! roughness length for scalar quantities z0h added
-!
-! Revision 1.1  1998/01/23 10:06:06  raasch
-! Initial revision
-!
 !
 ! Description:
