Ignore:
Timestamp:
May 9, 2014 2:03:54 PM (7 years ago)
Author:
knoop
Message:

Parallel random number generator added (preliminary version).

File:
1 edited

Legend:

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

    r1354 r1400  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Parallel random number generator added
    2323!
    2424! Former revisions:
     
    7272    USE random_function_mod,                                                   &
    7373        ONLY: random_function
     74       
     75    USE random_generator_parallel,                                             &
     76        ONLY:  random_number_parallel, random_seed_parallel, random_dummy,     &
     77               id_random_array, seq_random_array
    7478
    7579    IMPLICIT NONE
     
    112116                ENDIF
    113117             ENDDO
     118          ENDDO
     119       ENDDO
     120    ELSEIF ( random_generator == 'random-parallel' )  THEN
     121       DO  i = dist_nxl(dist_range), dist_nxr(dist_range)
     122          DO  j = dist_nys(dist_range), dist_nyn(dist_range)
     123             CALL random_seed_parallel( put=seq_random_array(:, j, i) )
     124             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
     125                CALL random_number_parallel( random_dummy )
     126                randomnumber = 3.0_wp * disturbance_amplitude *                &
     127                               ( random_dummy - 0.5_wp )
     128                IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.         &
     129                     nyn >= j )                                                &
     130                THEN
     131                   dist1(k,j,i) = randomnumber
     132                ENDIF
     133             ENDDO
     134!              IF (myid == 0) PRINT*, nxl, i, nxr, i, nys, j, nyn, j
     135             CALL random_seed_parallel( get=seq_random_array(:, j, i) )
    114136          ENDDO
    115137       ENDDO
Note: See TracChangeset for help on using the changeset viewer.