Changeset 1400 for palm/trunk


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

Parallel random number generator added (preliminary version).

Location:
palm/trunk/SOURCE
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r1381 r1400  
    2020# Current revisions:
    2121# ------------------
    22 #
     22# Added new module random_generator_parallel
    2323#
    2424# Former revisions:
     
    206206        poismg.f90 prandtl_fluxes.f90 pres.f90 print_1d.f90 \
    207207        production_e.f90 prognostic_equations.f90 random_function.f90 \
    208         random_gauss.f90 read_3d_binary.f90 read_var_list.f90 run_control.f90 \
     208        random_gauss.f90 random_generator_parallel.f90 read_3d_binary.f90 \
     209        read_var_list.f90 run_control.f90 \
    209210        set_slicer_attributes_dvrp.f90 singleton.f90 sor.f90 \
    210211        subsidence.f90 sum_up_3d_data.f90 \
     
    297298diffusion_w.o: modules.o mod_kinds.o wall_fluxes.o
    298299diffusivities.o: modules.o mod_kinds.o
    299 disturb_field.o: modules.o cpulog.o mod_kinds.o random_function.o
     300disturb_field.o: modules.o cpulog.o mod_kinds.o random_function.o random_generator_parallel.o
    300301disturb_heatflux.o: modules.o cpulog.o mod_kinds.o
    301302eqn_state_seawater.o: modules.o mod_kinds.o
     
    310311init_1d_model.o: modules.o mod_kinds.o
    311312init_3d_model.o: modules.o cpulog.o mod_kinds.o random_function.o advec_ws.o \
    312         ls_forcing.o lpm_init.o
     313        ls_forcing.o lpm_init.o random_generator_parallel.o
    313314init_advec.o: modules.o mod_kinds.o
    314315init_cloud_physics.o: modules.o mod_kinds.o
     
    348349lpm_init.o: modules.o lpm_collision_kernels.o mod_kinds.o \
    349350        random_function.o mod_particle_attributes.o lpm_exchange_horiz.o \
    350         lpm_pack_arrays.o
     351        lpm_pack_arrays.o random_generator_parallel.o
    351352lpm_init_sgs_tke.o: modules.o mod_kinds.o mod_particle_attributes.o
    352353lpm_pack_arrays.o: modules.o cpulog.o mod_kinds.o mod_particle_attributes.o
     
    354355        lpm_pack_arrays.o
    355356lpm_release_set.o: modules.o mod_kinds.o random_function.o \
    356         mod_particle_attributes.o lpm_init.o
     357        mod_particle_attributes.o lpm_init.o  random_generator_parallel.o
    357358lpm_set_attributes.o: modules.o cpulog.o mod_kinds.o mod_particle_attributes.o
    358359lpm_write_exchange_statistics.o: modules.o mod_kinds.o mod_particle_attributes.o
     
    383384        nudging.o plant_canopy_model.o production_e.o subsidence.o user_actions.o
    384385random_function.o: mod_kinds.o
    385 random_gauss.o: mod_kinds.o random_function.o
    386 read_3d_binary.o: modules.o cpulog.o mod_kinds.o random_function.o
     386random_gauss.o: mod_kinds.o random_function.o random_generator_parallel.o
     387random_generator_parallel.o: mod_kinds.o
     388read_3d_binary.o: modules.o cpulog.o mod_kinds.o random_function.o random_generator_parallel.o
    387389read_var_list.o: modules.o mod_kinds.o
    388390run_control.o: modules.o cpulog.o mod_kinds.o
     
    430432user_statistics.o: modules.o mod_kinds.o user_module.o
    431433wall_fluxes.o: modules.o mod_kinds.o
    432 write_3d_binary.o: modules.o cpulog.o mod_kinds.o random_function.o
     434write_3d_binary.o: modules.o cpulog.o mod_kinds.o random_function.o random_generator_parallel.o
    433435write_var_list.o: modules.o mod_kinds.o
  • palm/trunk/SOURCE/check_parameters.f90

    r1385 r1400  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Check random generator extended by option random-parallel
    2323!
    2424! Former revisions:
     
    34853485!
    34863486!-- Check random generator
    3487     IF ( random_generator /= 'system-specific'  .AND. &
    3488          random_generator /= 'numerical-recipes' )  THEN
     3487    IF ( (random_generator /= 'system-specific'     .AND.                      &
     3488          random_generator /= 'random-parallel'   ) .AND.                      &
     3489          random_generator /= 'numerical-recipes' )  THEN
    34893490       message_string = 'unknown random generator: random_generator = "' // &
    34903491                        TRIM( random_generator ) // '"'
  • 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
  • palm/trunk/SOURCE/init_3d_model.f90

    r1385 r1400  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Parallel random number generator added
    2323!
    2424! Former revisions:
     
    198198   
    199199    USE random_function_mod
     200   
     201    USE random_generator_parallel,                                             &
     202        ONLY:  random_number_parallel, random_seed_parallel, random_dummy,     &
     203               id_random_array, seq_random_array
    200204   
    201205    USE statistics,                                                            &
     
    481485
    482486!
     487!-- ID-array and state-space-array for the parallel random number generator
     488    IF ( random_generator == 'random-parallel' )  THEN
     489       ALLOCATE ( seq_random_array(5,nysg:nyng,nxlg:nxrg) )
     490       ALLOCATE ( id_random_array(0:ny,0:nx) )
     491       seq_random_array = 0
     492       id_random_array  = 0
     493    ENDIF
     494   
     495!
    483496!-- 3D-arrays for the leaf area density and the canopy drag coefficient
    484497    IF ( plant_canopy ) THEN
     
    913926
    914927!
     928!--    Initialize the random number generators (from numerical recipes)
     929       CALL random_function_ini
     930       IF ( random_generator == 'random-parallel' )  THEN
     931          CALL random_function_ini
     932          !- asigning an ID to every vertical gridpoint column
     933          DO j=0,ny
     934             DO i=0,nx
     935                id_random_array(j,i) = 1000000*i + j
     936             END DO
     937          ENDDO
     938          !- executing the initial ran_parallel
     939          DO j = nysg, nyng
     940             DO i = nxlg, nxrg
     941                CALL random_seed_parallel (random_sequence=id_random_array(j, i))
     942                CALL random_number_parallel (random_dummy)
     943                CALL random_seed_parallel (get=seq_random_array(:, j, i))
     944             END DO
     945          ENDDO
     946       ENDIF
     947
     948!
    915949!--    Initialize fluxes at bottom surface
    916950       IF ( use_surface_fluxes )  THEN
     
    10981132          q(nzb,:,:) = q(nzb,:,:) + q_surface_initial_change
    10991133       ENDIF
    1100 !
    1101 !--    Initialize the random number generator (from numerical recipes)
    1102        CALL random_function_ini
    11031134
    11041135!
  • palm/trunk/SOURCE/read_3d_binary.f90

    r1359 r1400  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! reading of arrays for random_generator_parallel added
    2323!
    2424! Former revisions:
     
    9595    USE random_function_mod,                                                   &
    9696        ONLY:  random_iv, random_iy
     97
     98    USE random_generator_parallel,                                             &
     99        ONLY:  id_random_array, seq_random_array
    97100
    98101    USE statistics,                                                            &
     
    681684                   IF ( k == 1 )  READ ( 13 )  random_iv
    682685                   IF ( k == 1 )  READ ( 13 )  random_iy
     686                   
     687                CASE ( 'seq_random_array' )  ! still unresolved issue
     688                   IF ( k == 1 )  READ ( 13 )  id_random_array
     689                   IF ( k == 1 )  READ ( 13 )  seq_random_array
    683690
    684691                CASE ( 'rho_av' )
  • palm/trunk/SOURCE/write_3d_binary.f90

    r1360 r1400  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! writing of arrays for random_generator_parallel added
    2323!
    2424! Former revisions:
     
    9191    USE random_function_mod,                                                   &
    9292        ONLY:  random_iv, random_iy
     93
     94    USE random_generator_parallel,                                             &
     95        ONLY:  id_random_array, seq_random_array
    9396       
    9497    USE statistics,                                                            &
     
    223226    WRITE ( 14 )  'random_iv           ';  WRITE ( 14 )  random_iv
    224227                                           WRITE ( 14 )  random_iy
     228    IF ( ALLOCATED( seq_random_array ) )  THEN
     229    WRITE ( 14 )  'seq_random_array    ';  WRITE ( 14 )  id_random_array
     230                                           WRITE ( 14 )  seq_random_array
     231    ENDIF
    225232    WRITE ( 14 )  'rif                 ';  WRITE ( 14 )  rif
    226233    IF ( topography /= 'flat' )  THEN
Note: See TracChangeset for help on using the changeset viewer.