Ignore:
Timestamp:
Jan 13, 2010 3:10:53 PM (11 years ago)
Author:
franke
Message:

collision of cloud droplets has changed in advec_particles
bugfixes in advec_particles and collision_efficiency
change in disturb_field to make runs reproducible on HLRN

File:
1 edited

Legend:

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

    r77 r420  
    44! Actual revisions:
    55! -----------------
    6 !
     6! A loop has been splitted to make runs reproducible on HLRN systems
    77!
    88! Former revisions:
     
    9898!
    9999!-- Exchange of ghost points for the random perturbation
     100
    100101    CALL exchange_horiz( dist1 )
    101102
     
    104105!-- Neighboured grid points in all three directions are used for the
    105106!-- filter operation.
    106     DO  i = nxl, nxr
    107        DO  j = nys, nyn
     107!-- Loop has been splitted to make runs reproducible on HLRN systems using
     108!-- compiler option -O3
     109     DO  i = nxl, nxr
     110        DO  j = nys, nyn
    108111          DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
    109              dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1) + dist1(k,j-1,i) &
    110                             + dist1(k,j+1,i) + dist1(k+1,j,i) + dist1(k-1,j,i) &
     112             dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1) &
     113                            + dist1(k,j+1,i) + dist1(k+1,j,i) &
     114                            ) / 12.0
     115          ENDDO
     116          DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
     117              dist2(k,j,i) = dist2(k,j,i) + ( dist1(k,j-1,i) + dist1(k-1,j,i)  &
    111118                            + 6.0 * dist1(k,j,i)                               &
    112119                            ) / 12.0
    113120          ENDDO
    114        ENDDO
    115     ENDDO
     121        ENDDO
     122     ENDDO
    116123
    117124!
     
    119126!-- Afterwards, filter operation and exchange of ghost points are repeated.
    120127    CALL exchange_horiz( dist2 )
     128
    121129    DO  i = nxl, nxr
    122130       DO  j = nys, nyn
     
    129137       ENDDO
    130138    ENDDO
     139
    131140    CALL exchange_horiz( dist1 )
    132141
Note: See TracChangeset for help on using the changeset viewer.