Changeset 1322 for palm/trunk/SOURCE/lpm_droplet_collision.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_droplet_collision.f90 (modified) (11 diffs)

palm/trunk/SOURCE/lpm_droplet_collision.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! REAL constants defined as wp_kind
 !
 ! Former revisions:
@@ -207 +207 @@
 !--                Determine dissipation class index of this gridbox
                    IF ( wang_kernel )  THEN
-                      eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
+                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
                                     dissipation_classes ) + 1
                       epsilon = diss(k,j,i)
@@ -213 +213 @@
                       epsilon = 0.0
                    ENDIF
-                   IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
+                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
                       eclass = 0   ! Hall kernel is used
                    ELSE
@@ -270 +270 @@
 !--                   Change of the current droplet radius
                       rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
-                                       sum3 )**0.33333333333333
+                                       sum3 )**0.33333333333333_wp
 
                       IF ( weight(n-is+1) < 0.0 )  THEN
@@ -347 +347 @@
 !--                   Change of the current droplet radius
                       rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
-                                       sum3 )**0.33333333333333
+                                       sum3 )**0.33333333333333_wp
 
                       IF ( weight(n-is+1) < 0.0 )  THEN
@@ -534 +534 @@
 !
 !--                   Change in radius due to collision
-                      delta_r = effective_coll_efficiency / 3.0             &
+                      delta_r = effective_coll_efficiency / 3.0_wp          &
                                 * pi * sl_r3 * ddx * ddy / dz               &
                                 * SQRT( ( u_int - particles(n)%speed_x )**2 &
@@ -551 +551 @@
                       IF ( delta_v >= sl_r3  .AND.  sl_r3 > 0.0 )  THEN
 
-                         delta_r = ( ( sl_r3/particles(n)%weight_factor )   &
-                                     + particles(n)%radius**3 )**( 1./3. )  &
+                         delta_r = ( ( sl_r3/particles(n)%weight_factor )           &
+                                     + particles(n)%radius**3 )**( 1.0_wp/3.0_wp )  &
                                      - particles(n)%radius
 
@@ -612 +612 @@
 
                    IF ( wang_kernel )  THEN
-                      eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
+                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
                                     dissipation_classes ) + 1
                       epsilon = diss(k,j,i)
@@ -618 +618 @@
                       epsilon = 0.0
                    ENDIF
-                   IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
+                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
                       eclass = 0   ! Hall kernel is used
                    ELSE
@@ -631 +631 @@
 
                    particles(psi)%radius = ( particles(psi)%radius**3 / &
-                                             sum3 )**0.33333333333333
+                                             sum3 )**0.33333333333333_wp
                    particles(psi)%weight_factor = particles(psi)%weight_factor &
                                                   * sum3
@@ -652 +652 @@
 
                    particles(psi)%radius = ( particles(psi)%radius**3 / &
-                                             sum3 )**0.33333333333333
+                                             sum3 )**0.33333333333333_wp
                    particles(psi)%weight_factor = particles(psi)%weight_factor &
                                                   * sum3