C
C           Demonstrates the use of make_seed with MPI          
C'make_sprng_seed' is called to produce a new seed each time the program
C is run. The same seed is produced on each process. 
C

       program seedf_simple_mpi
       implicit none

#define SIMPLE_SPRNG	
#define USE_MPI 1
#include <mpif.h>
#include "sprng_f.h"

       integer seed, i, myid, ierror, junk
       SPRNG_POINTER junkPtr
       real*8 rn
C---
       integer gtype
C---                                                                           
      call MPI_INIT(ierror)
      call MPI_COMM_RANK(MPI_COMM_WORLD, myid, ierror)

       seed = make_sprng_seed() !produce a new seed each time program is run
C--- node 0 is reading in a generator type
       if (myid .eq. 0) then
C--- You may not need this. This is for generator type explanation.
#include "genf_types_menu.h"
         print *,'Type in an integer for the generator type:  '
         read *, gtype
       endif
       call MPI_BCAST(gtype,1, MPI_INTEGER,0,MPI_COMM_WORLD,ierror)
C---                                                                          

C Seed should be the same on all processes
       write(6, 911) myid, seed
 911   format("Process", i2, ": seed = ", i16)
    
       junkPtr = init_sprng(gtype,seed,SPRNG_DEFAULT)
       write(6, 922) myid
 922   format("Process", i2, 
     &        ": Print information about stream:")
       junk = print_sprng()

       do 100 i = 1, 3
          rn = sprng()
          write(6, 933) myid, i, rn
 100   continue
 933   format("process", i2, ", random number", i2, ": ", f8.6)

       call MPI_FINALIZE(ierror)

       end