Changeset 1320 for palm/trunk/SOURCE/diffusivities.f90

Timestamp:

Mar 20, 2014 8:40:49 AM (10 years ago)

Author:

raasch

Message:

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

File:

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

palm/trunk/SOURCE/diffusivities.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! 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 
 !
 ! Former revisions:
@@ -36 +42 @@
 ! adjustment of mixing length to the Prandtl mixing length at first grid point
 ! above ground removed
-!
-! 667 2010-12-23 12:06:00Z suehring/gryschka
-! nxl-1, nxr+1, nys-1, nyn+1 replaced by nxlg, nxrg, nysg, nyng
-!
-! 137 2007-11-28 08:50:10Z letzel
-! Bugfix for summation of sums_l_l for flow_statistics
-! Vertical scalar profiles now based on nzb_s_inner and ngp_2dh_s_inner.
-!
-! 97 2007-06-21 08:23:15Z raasch
-! Adjustment of mixing length calculation for the ocean version.
-! This is also a bugfix, because the height above the topography is now
-! used instead of the height above level k=0.
-! theta renamed var, dpt_dz renamed dvar_dz, +new argument var_reference
-! use_pt_reference renamed use_reference
-!
-! 57 2007-03-09 12:05:41Z raasch
-! Reference temperature pt_reference can be used in buoyancy term
-!
-! RCS Log replace by Id keyword, revision history cleaned up
-!
-! Revision 1.24  2006/04/26 12:16:26  raasch
-! OpenMP optimization (+sums_l_l_t), sqrt_e must be private
 !
 ! Revision 1.1  1997/09/19 07:41:10  raasch
@@ -69 +53 @@
 !------------------------------------------------------------------------------!
 
-    USE arrays_3d
-    USE control_parameters
-    USE grid_variables
-    USE indices
+    USE arrays_3d,                                                             &
+        ONLY:  dd2zu, e, kh, km, l_grid, l_wall
+        
+    USE control_parameters,                                                    &
+        ONLY:  atmos_ocean_sign, e_min, g, outflow_l, outflow_n, outflow_r,    &
+                outflow_s, use_single_reference_value, wall_adjustment
+                
+    USE indices,                                                               &
+        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb_s_inner, nzb, nzt
+    USE kinds
+    
     USE pegrid
-    USE statistics
+    
+    USE statistics,                                                            &
+        ONLY :  rmask, statistic_regions, sums_l_l
 
     IMPLICIT NONE
 
-    INTEGER ::  i, j, k, omp_get_thread_num, sr, tn
-
-    REAL    ::  dvar_dz, l, ll, l_stable, sqrt_e, var_reference
-
-    REAL    ::  var(nzb:nzt+1,nysg:nyng,nxlg:nxrg)
+    INTEGER(iwp) ::  i                   !:
+    INTEGER(iwp) ::  j                   !:
+    INTEGER(iwp) ::  k                   !:
+    INTEGER(iwp) ::  omp_get_thread_num  !:
+    INTEGER(iwp) ::  sr                  !:
+    INTEGER(iwp) ::  tn                  !:
+
+    REAL(wp)     ::  dvar_dz             !:
+    REAL(wp)     ::  l                   !:
+    REAL(wp)     ::  ll                  !:
+    REAL(wp)     ::  l_stable            !:
+    REAL(wp)     ::  sqrt_e              !:
+    REAL(wp)     ::  var_reference       !:
+
+    REAL(wp)     ::  var(nzb:nzt+1,nysg:nyng,nxlg:nxrg)  !:
 
 
@@ -136 +139 @@
                 IF ( dvar_dz > 0.0 ) THEN
                    IF ( use_single_reference_value )  THEN
-                      l_stable = 0.76 * sqrt_e / &
+                      l_stable = 0.76 * sqrt_e /                               &
                                  SQRT( g / var_reference * dvar_dz ) + 1E-5
                    ELSE
-                      l_stable = 0.76 * sqrt_e / &
+                      l_stable = 0.76 * sqrt_e /                               &
                                  SQRT( g / var(k,j,i) * dvar_dz ) + 1E-5
                    ENDIF