Creating a Grid Program using MPI.

Timothy H. Kaiser, Ph.D.
tkaiser@sdsc.edu
www.sdsc.edu/~tkaiser
San Diego Supercomputer Center
 
 
This page has links to the slides and source code for my talk given January 99 at NAVO on creating a grid program using MPI.  See the notes below for program changes.

If you are not familiar with MPI you may want to look at http://www.sdsc.edu/~tkaiser/my_mpi.html 

If you are not familiar with Fortran 90  you may want to look at http://www.sdsc.edu/~tkaiser/f90.html
 
  • talk1.pdf:  Stommel model. (202554 bytes )
  • talk2.pdf:  Stommel Model in Fortran 90. (96345 bytes )
  • talk3.pdf:  Stommel Model in MPI with 1d and 2d decompisition. (316496 bytes )
  • my_stommel.tar: Contains the source code, the pdf files of the talks and the notes shown below. (720896 bytes )
  • The plot generated by Luke's plotit.c and my parallel stommel program, st_4.f

    •  
     
    Notes included in my_stommel.tar:
     I made a few changes to the programs: 
    (1) I changed the forcing function to force=-alpha*sin(y*a6)). 
    (2) I changed the output format slightly for stommel.f, st_3.f, 
        and st_4.f to be able to be read by the plotting program. 
    (3) new input values 

    golden % cat st.in 
    200 200 
    2000000 2000000 
    1.0e-9 2.25e-11 3.0e-6 
    75000 

    contents: 
     

    talk1.pdf, talk2.pdf, talk3.pdf  - new versions of the stommel/mpi 
             talks i gave 
     
    stommel.f -  serial stommel program 

    st_1.f - 1d decomposition mpi version of stommel program 

    st_2.f - 2d decomposition mpi version 

    st_3.f - same as st_2.f except all data is sent back to 
             mpi_master and it write the output file.  this 
             version produces the exact same results as stommel.f 
             but it prints the residual only every 1000 steps 
             I added slides to talk3 to discuss how the output 
             works using mpi_gather and mpi_gatherv 
     
    st_4.f - same as st_3.f but the forcing function is stored 
             in an array so it is faster 
     
    plotit.c - slight modification of Luke's plotting programs 
             to use the program run st_4d to produce out3d_000. 
             run plotter as shown below.... 
     

    golden % plotit 
     This program will convert a Plot3D Whole IEEE Unformatted Function 
     file into a pair of files suitable for use with GNUplot 
     

     Please Enter the Plot3D format input file name: out3d_000 

     Dimensions of input file: NI:   202 NJ:   202 NK:     1 
    0   5.08452e+07 

     You should now be able to execute the command: 
                    "gnuplot gnuplot_commandfile" 
     to get a 2D XY contour plot and postscript printable output. 
     

     Normal Termination. 
    golden % 
    golden % 
    golden % gnuplot gnuplot_commandfile 

    generates a file plotfile.ps which can be seen which ghostview. 
     
     
     
     
    Last Revised: January 29, 1999
    tkaiser@sdsc.edu