Changeset 4577 for palm/trunk/SOURCE/biometeorology_mod.f90

Timestamp:

Jun 25, 2020 9:53:58 AM (4 years ago)

Author:

raasch

Message:

further re-formatting to follow the PALM coding standard

File:

• palm/trunk/SOURCE/biometeorology_mod.f90 (modified) (14 diffs)

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! further re-formatting concerning Fortran parameter variables
+!
+! 4540 2020-05-18 15:23:29Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -187 +190 @@
 !
 !-- Declare all global variables within the module (alphabetical order)
+    REAL(wp), PARAMETER ::  bio_fill_value = -9999.0_wp  !< set module fill value, replace by global fill value as soon as available
+    REAL(wp), PARAMETER ::  human_absorb = 0.7_wp  !< SW absorbtivity of a human body (Fanger 1972)
+    REAL(wp), PARAMETER ::  human_emiss = 0.97_wp  !< LW emissivity of a human body after (Fanger 1972)
+
     INTEGER(iwp) ::  bio_cell_level     !< cell level biom calculates for
 
@@ -209 +216 @@
 
     REAL(wp)    ::  bio_output_height  !< height output is calculated in m
-
-    REAL(wp), PARAMETER ::  bio_fill_value = -9999.0_wp  !< set module fill value, replace by global fill value as soon as available
-    REAL(wp), PARAMETER ::  human_absorb = 0.7_wp  !< SW absorbtivity of a human body (Fanger 1972)
-    REAL(wp), PARAMETER ::  human_emiss = 0.97_wp  !< LW emissivity of a human body after (Fanger 1972)
 
     REAL(wp), DIMENSION(:), ALLOCATABLE   ::  mrt_av_grid   !< time average mean
@@ -2270 +2273 @@
 !
 !-- Parameters for standard "Klima-Michel"
-    REAL(wp), PARAMETER :: actlev = 134.6862_wp  !< Workload by activity per standardized surface 
-                                                 !< (A_Du)
+    REAL(wp), PARAMETER :: actlev = 134.6862_wp  !< Workload by activity per standardized surface (A_Du)
     REAL(wp), PARAMETER :: eta = 0.0_wp          !< Mechanical work efficiency for walking on flat
                                                  !< ground (compare to Fanger (1972) pp 24f)
 !
 !-- Type of program variables
+    REAL(wp), PARAMETER :: eps = 0.0005  !< Accuracy in clothing insulation (clo) for evaluation the root of Fanger's PMV (pmva=0)
+
     INTEGER(iwp) :: ncount      !< running index
     INTEGER(iwp) :: nerr_cold   !< error number (cold conditions)
@@ -2281 +2285 @@
 
     LOGICAL :: sultrieness
-
-    REAL(wp), PARAMETER :: eps = 0.0005  !< Accuracy in clothing insulation (clo) for evaluation the
-                                         !< root of Fanger's PMV (pmva=0)
 
     REAL(wp) ::  clon           !< clo for neutral conditions   (clo)
@@ -2504 +2505 @@
 !-- Type of program variables
     INTEGER(iwp), PARAMETER  ::  max_iteration = 15_iwp       !< max number of iterations
+
+    REAL(wp),     PARAMETER  ::  guess_0       = -1.11e30_wp  !< initial guess
+
     INTEGER(iwp) ::  j       !< running index
-
-    REAL(wp),     PARAMETER  ::  guess_0       = -1.11e30_wp  !< initial guess
 
     REAL(wp) ::  clo_lower   !< lower limit of clothing insulation      (clo)
@@ -2788 +2790 @@
     REAL(wp) ::  dtmrt        !< difference mean radiation to air temperature
     REAL(wp) ::  pa           !< vapor pressure (hPa) with hard bounds
-    REAL(wp) ::  pa_p50       !< ratio actual water vapour pressure to that of relative humidity of 
+    REAL(wp) ::  pa_p50       !< ratio actual water vapour pressure to that of relative humidity of
                               !< 50 %
     REAL(wp) ::  pmv          !< temp storage og predicted mean vote
@@ -2973 +2975 @@
 !
 !-- Additional output variables of argument list:
-    REAL(wp), INTENT ( OUT ) ::  dperctm    !< Mean deviation perct (classical gt) to gt* (rational 
+    REAL(wp), INTENT ( OUT ) ::  dperctm    !< Mean deviation perct (classical gt) to gt* (rational
                                             !< gt calculated based on Gagge's rational PMV*)
     REAL(wp), INTENT ( OUT ) ::  dperctstd  !< dperctm plus its standard deviation times a factor
@@ -3039 +3041 @@
 !
 !-- Type of output argument
-    INTEGER(iwp), INTENT ( OUT ) ::  nerr !< Error indicator: 0 = o.k., +1 = denominator for 
+    INTEGER(iwp), INTENT ( OUT ) ::  nerr !< Error indicator: 0 = o.k., +1 = denominator for
                                           !< intersection = 0
 
-    REAL(wp),     INTENT ( OUT ) ::  dpmv_cold_res    !< Increment to adjust pmva according to the 
+    REAL(wp),     INTENT ( OUT ) ::  dpmv_cold_res    !< Increment to adjust pmva according to the
                                                       !< results of Gagge's 2 node model depending on the input
 !
@@ -3344 +3346 @@
 !
 !-- Internal variables
+    REAL(wp), PARAMETER :: eps = 0.0005_wp
+    REAL(wp), PARAMETER :: eta = 0.0_wp
+
     INTEGER(iwp) ::  ncount
     INTEGER(iwp) ::  nerr_cold
@@ -3349 +3354 @@
 
     LOGICAL ::  sultrieness
-
-    REAL(wp), PARAMETER :: eps = 0.0005_wp
-    REAL(wp), PARAMETER :: eta = 0.0_wp
 
 !    REAL(wp) ::  acti
@@ -3639 +3641 @@
 !
 !-- Internal variables
+    REAL(wp), PARAMETER  ::  time_equil = 7200.0_wp
+
     INTEGER(iwp) :: i         !< running index
     INTEGER(iwp) ::  niter    !< Running index
-
-    REAL(wp), PARAMETER  ::  time_equil = 7200.0_wp
 
     REAL(wp) ::  adjustrate        !< Max storage adjustment rate
@@ -3661 +3663 @@
     REAL(wp) ::  ws                !< wind speed                             (m/s)
     REAL(wp) ::  z1                !< Empiric factor for the adaption of the heat
-                                   !< ballance equation to the psycho-physical scale 
+                                   !< ballance equation to the psycho-physical scale
                                    !< (Equ. 40 in FANGER)
     REAL(wp) ::  z2                !< Water vapour diffution through the skin
@@ -4087 +4089 @@
           IF ( eswdif <= 0.0_wp )  esw = eswpot     !< Limit is evaporation
           IF ( eswdif > 0.0_wp )   esw = eswphy     !< Limit is sweat production
-          IF ( esw  > 0.0_wp )     esw = 0.0_wp     !< Sweat can't be evaporated, no more cooling 
+          IF ( esw  > 0.0_wp )     esw = 0.0_wp     !< Sweat can't be evaporated, no more cooling
                                                     !< effect
 !