next up previous

17.3.3 Model and Figure Generation

The Reef velocity model can be generated using the reef.f program and the model.dat file. The user must specify the grid size, grid spacing, the length of time for the model run and the time step. A model name is also required.

The output can be viewed with a postscript viewer such as ghostview. A Unix shell program reef_reef.sh is provided, and it generates a reef_reef.ps file. The first process of the shell is to convert the unformatted output file from reef.f into an SU format file using suaddhead. The shell uses SU program supsimage to convert the SU file to postscript.

The seismic image file with the time data is generated using the seismic.f program. The output file is processed by reef_seismic.sh, which is a shell program that is similar to the reef_reef.sh shell program. It converts the program output into an SU file and then into a .ps file.

The reverse time imaging which converts the time data back into depth data is generated using the revtime.f program. The output file is processed by the reef_revtime.sh shell similarly to reef_reef.sh. The depth image data file is converted from a Fortran unformatted file to an SU file and then to a .ps file.

The Salt Intrusion model and the Dipping layer models have Unix scripts and programs identical to those for the Reef model. The names of the files used for all three models are shown in Table 2.

Table 2: Fortran Codes and Shell Scripts Used in Seismic Modeling and Imaging of Data.

A comment about the shell programs described above: they all require parameters to create output data images. Any changes to model.dat will probably require changes in the shell programs.