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

Timestamp:

Oct 26, 2015 4:17:44 PM (8 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/init_grid.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_grid.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:
@@ -144 +144 @@
                building_length_y, building_wall_left, building_wall_south,     &
                canyon_height, canyon_wall_left, canyon_wall_south,             &
-               canyon_width_x, canyon_width_y, coupling_char, dp_level_ind_b,  &
-               dz, dz_max, dz_stretch_factor, dz_stretch_level,                &
-               dz_stretch_level_index, ibc_uv_b, io_blocks, io_group,          &
-               inflow_l, inflow_n, inflow_r, inflow_s, masking_method,         &
-               maximum_grid_level, message_string, momentum_advec, ocean,      &
-               outflow_l, outflow_n, outflow_r, outflow_s, prandtl_layer,      &
-               psolver, scalar_advec, topography, topography_grid_convention,  &
-               use_surface_fluxes, use_top_fluxes, wall_adjustment_factor 
+               canyon_width_x, canyon_width_y, constant_flux_layer,            &
+               coupling_char, dp_level_ind_b, dz, dz_max, dz_stretch_factor,   &
+               dz_stretch_level, dz_stretch_level_index, ibc_uv_b, io_blocks,  &
+               io_group, inflow_l, inflow_n, inflow_r, inflow_s,               &
+               masking_method, maximum_grid_level, message_string,             &
+               momentum_advec, ocean, outflow_l, outflow_n, outflow_r,         &
+               outflow_s, psolver, scalar_advec, topography,                   &
+               topography_grid_convention, use_surface_fluxes, use_top_fluxes, &
+               wall_adjustment_factor 
         
     USE grid_variables,                                                        &
@@ -477 +478 @@
 !-- Define vertical gridpoint from (or to) which on the usual finite difference
 !-- form (which does not use surface fluxes) is applied 
-    IF ( prandtl_layer  .OR.  use_surface_fluxes )  THEN
+    IF ( constant_flux_layer  .OR.  use_surface_fluxes )  THEN
        nzb_diff = nzb + 2
     ELSE
@@ -955 +956 @@
 !-- the usual finite difference form (which does not use surface fluxes) is
 !-- applied 
-    IF ( prandtl_layer  .OR.  use_surface_fluxes )  THEN
+    IF ( constant_flux_layer  .OR.  use_surface_fluxes )  THEN
        nzb_diff_u         = nzb_u_inner + 2
        nzb_diff_v         = nzb_v_inner + 2