Changeset 1322 for palm/trunk/SOURCE/lpm_advec.f90

Timestamp:

Mar 20, 2014 4:38:49 PM (10 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

• palm/trunk/SOURCE/lpm_advec.f90 (modified) (10 diffs)

palm/trunk/SOURCE/lpm_advec.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
+! REAL constants defined as wp_kind
 !
 ! Former revisions:
@@ -186 +187 @@
 !
 !--          Determine the sublayer. Further used as index.
-             height_p     = ( particles(n)%z - z0_av_global )           &
-                                  * REAL( number_of_sublayers )         &
+             height_p     = ( particles(n)%z - z0_av_global )            &
+                                  * REAL( number_of_sublayers, KIND=wp ) &
                                   * d_z_p_z0 
 !
@@ -210 +211 @@
              us_int = 0.5 * ( us(j,i) + us(j,i-1) )  
 
-             u_int  = -usws(j,i) / ( us_int * kappa + 1E-10 )           & 
+             u_int  = -usws(j,i) / ( us_int * kappa + 1E-10_wp )           & 
                       * log_z_z0_int 
 
@@ -277 +278 @@
               us_int = 0.5 * ( us(j,i) + us(j-1,i) )    
 
-              v_int  = -vsws(j,i) / ( us_int * kappa + 1E-10 )             &
+              v_int  = -vsws(j,i) / ( us_int * kappa + 1E-10_wp )             &
                        * log_z_z0_int 
 
@@ -1093 +1094 @@
           ENDIF
 
-          vv_int = ( 1.0 / 3.0 ) * ( aa + bb + cc )
-
-          fs_int = ( 2.0 / 3.0 ) * e_mean_int / &
-                      ( vv_int + ( 2.0 / 3.0 ) * e_mean_int )
+          vv_int = ( 1.0_wp / 3.0_wp ) * ( aa + bb + cc )
+
+          fs_int = ( 2.0_wp / 3.0_wp ) * e_mean_int / &
+                      ( vv_int + ( 2.0_wp / 3.0_wp ) * e_mean_int )
 
 !
@@ -1126 +1127 @@
 !--          from becoming unrealistically large.
              particles(n)%rvar1 = SQRT( 2.0 * sgs_wfu_part * e_int ) * &
-                                        ( random_gauss( iran_part, 5.0 ) - 1.0 )
+                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
              particles(n)%rvar2 = SQRT( 2.0 * sgs_wfv_part * e_int ) * &
-                                        ( random_gauss( iran_part, 5.0 ) - 1.0 )
+                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
              particles(n)%rvar3 = SQRT( 2.0 * sgs_wfw_part * e_int ) * &
-                                        ( random_gauss( iran_part, 5.0 ) - 1.0 )
+                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
 
           ELSE
@@ -1168 +1169 @@
                                   ) * dt_particle / 2.0          +             &
                                   SQRT( fs_int * c_0 * diss_int ) *            &
-                                  ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
+                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
                                   SQRT( dt_particle )
 
@@ -1177 +1178 @@
                                     particles(n)%rvar2 /                       &
                                     ( 2.0 * sgs_wfv_part * e_int ) + de_dy_int &
-                                  ) * dt_particle / 2.0          +             &
+                                  ) * dt_particle / 2.0_wp          +             &
                                   SQRT( fs_int * c_0 * diss_int ) *            &
-                                  ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
+                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
                                   SQRT( dt_particle )
 
@@ -1188 +1189 @@
                                     particles(n)%rvar3 /                       &
                                     ( 2.0 * sgs_wfw_part * e_int ) + de_dz_int &
-                                  ) * dt_particle / 2.0          +             &
+                                  ) * dt_particle / 2.0_wp         +             &
                                   SQRT( fs_int * c_0 * diss_int ) *            &
-                                  ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
+                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
                                   SQRT( dt_particle )
 
@@ -1252 +1253 @@
                                 ( v_int - particles(n)%speed_y )**2 +   &
                                 ( w_int - particles(n)%speed_z )**2 ) / &
-                                         molecular_viscosity )**0.687   &
+                                         molecular_viscosity )**0.687_wp   &
                         )
              exp_term = EXP( -exp_arg * dt_particle )