Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (8 years ago)
Author:
raasch
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/lpm_droplet_collision.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp_kind
    2323!
    2424! Former revisions:
     
    207207!--                Determine dissipation class index of this gridbox
    208208                   IF ( wang_kernel )  THEN
    209                       eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
     209                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
    210210                                    dissipation_classes ) + 1
    211211                      epsilon = diss(k,j,i)
     
    213213                      epsilon = 0.0
    214214                   ENDIF
    215                    IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
     215                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
    216216                      eclass = 0   ! Hall kernel is used
    217217                   ELSE
     
    270270!--                   Change of the current droplet radius
    271271                      rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
    272                                        sum3 )**0.33333333333333
     272                                       sum3 )**0.33333333333333_wp
    273273
    274274                      IF ( weight(n-is+1) < 0.0 )  THEN
     
    347347!--                   Change of the current droplet radius
    348348                      rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
    349                                        sum3 )**0.33333333333333
     349                                       sum3 )**0.33333333333333_wp
    350350
    351351                      IF ( weight(n-is+1) < 0.0 )  THEN
     
    534534!
    535535!--                   Change in radius due to collision
    536                       delta_r = effective_coll_efficiency / 3.0             &
     536                      delta_r = effective_coll_efficiency / 3.0_wp          &
    537537                                * pi * sl_r3 * ddx * ddy / dz               &
    538538                                * SQRT( ( u_int - particles(n)%speed_x )**2 &
     
    551551                      IF ( delta_v >= sl_r3  .AND.  sl_r3 > 0.0 )  THEN
    552552
    553                          delta_r = ( ( sl_r3/particles(n)%weight_factor )   &
    554                                      + particles(n)%radius**3 )**( 1./3. )  &
     553                         delta_r = ( ( sl_r3/particles(n)%weight_factor )           &
     554                                     + particles(n)%radius**3 )**( 1.0_wp/3.0_wp )  &
    555555                                     - particles(n)%radius
    556556
     
    612612
    613613                   IF ( wang_kernel )  THEN
    614                       eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
     614                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
    615615                                    dissipation_classes ) + 1
    616616                      epsilon = diss(k,j,i)
     
    618618                      epsilon = 0.0
    619619                   ENDIF
    620                    IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
     620                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
    621621                      eclass = 0   ! Hall kernel is used
    622622                   ELSE
     
    631631
    632632                   particles(psi)%radius = ( particles(psi)%radius**3 / &
    633                                              sum3 )**0.33333333333333
     633                                             sum3 )**0.33333333333333_wp
    634634                   particles(psi)%weight_factor = particles(psi)%weight_factor &
    635635                                                  * sum3
     
    652652
    653653                   particles(psi)%radius = ( particles(psi)%radius**3 / &
    654                                              sum3 )**0.33333333333333
     654                                             sum3 )**0.33333333333333_wp
    655655                   particles(psi)%weight_factor = particles(psi)%weight_factor &
    656656                                                  * sum3
Note: See TracChangeset for help on using the changeset viewer.