Changeset 1691 for palm/trunk/SOURCE/lpm_advec.f90

Timestamp:

Oct 26, 2015 4:17:44 PM (9 years ago)

Author:

maronga

Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:

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

palm/trunk/SOURCE/lpm_advec.f90

@@ -14 +14 @@
 ! PALM. If not, see <http://www.gnu.org/licenses/>.
 !
-! Copyright 1997-2014 Leibniz Universitaet Hannover
+! Copyright 1997-2015 Leibniz Universitaet Hannover
 !--------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! ------------------
-!
+! Renamed prandtl_layer to constant_flux_layer.
 !
 ! Former revisions:
@@ -83 +83 @@
 
     USE control_parameters,                                                    &
-        ONLY:  atmos_ocean_sign, cloud_droplets, dt_3d, dt_3d_reached_l, dz,   &
-               g, kappa, molecular_viscosity, prandtl_layer, topography,       &
+        ONLY:  atmos_ocean_sign, cloud_droplets, constant_flux_layer, dt_3d,   &
+               dt_3d_reached_l, dz, g, kappa, molecular_viscosity, topography, &
                u_gtrans, v_gtrans, simulated_time
 
@@ -224 +224 @@
 !--       First, check if particle is located below first vertical grid level 
 !--       (Prandtl-layer height)
-          IF ( prandtl_layer  .AND.  particles(n)%z < z_p )  THEN
+          IF ( constant_flux_layer  .AND.  particles(n)%z < z_p )  THEN
 !
 !--          Resolved-scale horizontal particle velocity is zero below z0. 
@@ -294 +294 @@
        DO  n = start_index(nb), end_index(nb)
 
-          IF ( prandtl_layer  .AND.  particles(n)%z < z_p )  THEN
+          IF ( constant_flux_layer  .AND.  particles(n)%z < z_p )  THEN
 
              IF ( particles(n)%z < z0_av_global )  THEN