Changeset 1691 for palm/trunk/SOURCE/wall_fluxes.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/wall_fluxes.f90 (modified) (8 diffs)

palm/trunk/SOURCE/wall_fluxes.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 rif_min and rif_max with zeta_min and zeta_max, respectively.
 ! 
 ! Former revisions:
@@ -69 +69 @@
 !> it gives slightly different results from the ij-version for some unknown
 !> reason.
+!>
+!> @todo Rename rif to zeta throughout the routine
 !------------------------------------------------------------------------------!
  MODULE wall_fluxes_mod
@@ -107 +109 @@
        
        USE control_parameters,                                                 &
-           ONLY:  g, kappa, rif_max, rif_min
+           ONLY:  g, kappa, zeta_max, zeta_min
        
        USE grid_variables,                                                     &
@@ -234 +236 @@
 !--                shear stresses and very small momentum fluxes (both are
 !--                generally unrealistic).
-                   IF ( rifs < rif_min )  rifs = rif_min
-                   IF ( rifs > rif_max )  rifs = rif_max
+                   IF ( rifs < zeta_min )  rifs = zeta_min
+                   IF ( rifs > zeta_max )  rifs = zeta_max
 
 !
@@ -291 +293 @@
        
        USE control_parameters,                                                 &
-           ONLY:  g, kappa, rif_max, rif_min
+           ONLY:  g, kappa, zeta_max, zeta_min
        
        USE grid_variables,                                                     &
@@ -428 +430 @@
 !--                shear stresses and very small momentum fluxes (both are
 !--                generally unrealistic).
-                   IF ( rifs < rif_min )  rifs = rif_min
-                   IF ( rifs > rif_max )  rifs = rif_max
+                   IF ( rifs < zeta_min )  rifs = zeta_min
+                   IF ( rifs > zeta_max )  rifs = zeta_max
 
 !
@@ -485 +487 @@
        
        USE control_parameters,                                                 &
-           ONLY:  g, kappa, rif_max, rif_min
+           ONLY:  g, kappa, zeta_max, zeta_min
        
        USE grid_variables,                                                     &
@@ -596 +598 @@
 !--       consequence would result in very large shear stresses and very
 !--       small momentum fluxes (both are generally unrealistic).
-          IF ( rifs < rif_min )  rifs = rif_min
-          IF ( rifs > rif_max )  rifs = rif_max
+          IF ( rifs < zeta_min )  rifs = zeta_min
+          IF ( rifs > zeta_max )  rifs = zeta_max
 
 !