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

palm/trunk/SOURCE/calc_precipitation.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,
+! revision history before 2012 removed,
+! comment fields (!:) to be used for variable explanations added to
+! all variable declaration statements
 !
 ! Former revisions:
@@ -28 +33 @@
 ! 1036 2012-10-22 13:43:42Z raasch
 ! code put under GPL (PALM 3.9)
-!
-! 403 2009-10-22 13:57:16Z franke
-! Bugfix in calculation of precipitation_rate(j,i)
-!
-! 73 2007-03-20 08:33:14Z raasch
-! Precipitation rate and amount are calculated/stored,
-! + module control_parameters
-!
-! 19 2007-02-23 04:53:48Z raasch
-! Calculation extended for gridpoint nzt
-!
-! RCS Log replace by Id keyword, revision history cleaned up
-!
-! Revision 1.5  2004/01/30 10:15:57  raasch
-! Scalar lower k index nzb replaced by 2d-array nzb_2d
 !
 ! Revision 1.1  2000/04/13 14:45:22  schroeter
@@ -71 +61 @@
     SUBROUTINE calc_precipitation
 
-       USE arrays_3d
-       USE cloud_parameters
-       USE constants
-       USE control_parameters
-       USE indices
+       USE arrays_3d,                                                          &
+           ONLY:  dzw, ql, tend
+
+       USE cloud_parameters,                                                   &
+           ONLY:  precipitation_amount, precipitation_rate, prec_time_const,   &
+                  ql_crit
+
+       USE control_parameters,                                                 &
+           ONLY:  dt_do2d_xy, dt_3d,                                           &
+                  intermediate_timestep_count, intermediate_timestep_count_max,&
+                  precipitation_amount_interval, time_do2d_xy
+
+       USE indices,                                                            &
+           ONLY:  nxl, nxr, nyn, nys, nzb_2d, nzt
+
+       USE kinds
+
 
        IMPLICIT NONE
 
-       INTEGER ::  i, j, k
-       REAL    ::  dqdt_precip
+       INTEGER(iwp) ::  i !:
+       INTEGER(iwp) ::  j !:
+       INTEGER(iwp) ::  k !:
+       
+       REAL(wp)    ::  dqdt_precip !:
 
        precipitation_rate = 0.0
@@ -96 +101 @@
 !
 !--             Precipitation rate in kg / m**2 / s (= mm/s)
-                precipitation_rate(j,i) = precipitation_rate(j,i) + &
+                precipitation_rate(j,i) = precipitation_rate(j,i) +            &
                                           dqdt_precip * dzw(k)
 
@@ -102 +107 @@
 !
 !--          Sum up the precipitation amount, unit kg / m**2 (= mm)
-             IF ( intermediate_timestep_count ==         &
-                  intermediate_timestep_count_max  .AND. &
+             IF ( intermediate_timestep_count ==                               &
+                  intermediate_timestep_count_max  .AND.                       &
                   ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&
              THEN
-                precipitation_amount(j,i) = precipitation_amount(j,i) + &
+                precipitation_amount(j,i) = precipitation_amount(j,i) +        &
                                             precipitation_rate(j,i) * dt_3d
              ENDIF
@@ -120 +125 @@
     SUBROUTINE calc_precipitation_ij( i, j )
 
-       USE arrays_3d
-       USE cloud_parameters
-       USE constants
-       USE control_parameters
-       USE indices
+       USE arrays_3d,                                                          &
+           ONLY:  dzw, ql, tend
+
+       USE cloud_parameters,                                                   &
+           ONLY:  precipitation_amount, precipitation_rate, prec_time_const,   &
+                  ql_crit
+
+       USE control_parameters,                                                 &
+           ONLY:  dt_do2d_xy, dt_3d,                                           &
+                  intermediate_timestep_count, intermediate_timestep_count_max,&
+                  precipitation_amount_interval, time_do2d_xy
+
+       USE indices,                                                            &
+           ONLY:  nzb_2d, nzt
+
+       USE kinds
+
 
        IMPLICIT NONE
 
-       INTEGER ::  i, j, k
-       REAL    ::  dqdt_precip
+       INTEGER(iwp) ::  i !:
+       INTEGER(iwp) ::  j !:
+       INTEGER(iwp) ::  k !:
+       
+       REAL(wp)    ::  dqdt_precip !:       
 
        precipitation_rate(j,i) = 0.0
@@ -147 +167 @@
 !
 !--       Precipitation rate in kg / m**2 / s (= mm/s)
-          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip * &
+          precipitation_rate(j,i) = precipitation_rate(j,i) + dqdt_precip *    &
                                                               dzw(k)
 
@@ -157 +177 @@
             .AND. ( dt_do2d_xy-time_do2d_xy ) < precipitation_amount_interval )&
        THEN
-          precipitation_amount(j,i) = precipitation_amount(j,i) + &
+          precipitation_amount(j,i) = precipitation_amount(j,i) +              &
                                       precipitation_rate(j,i) * dt_3d
        ENDIF