Extracting Values from the T and S <BR>Database Files
Execute the shell script TS_extract
and answer the prompts for the longitudes and latitudes bounding the region of
interest.
This script:
- first creates an updated version of a fortran source file
newextract_ts.f
- then compiles and runs newextract_ts, which
reads the temperature and salinity
database files and extracts T and S data for this geographic region,
writing the output to file tands
- then creates an updated version of a fortran source file
newplotform.f
- finally compiles and runs newplotform,
which
reads tands and writes the (x,y,z)-format files
fort.21, ..., fort.25 (longitude-latitude-temperature values)
and
fort.31, ..., fort.35 (longitude-latitude-salinity)
for input to the GMT plotting routines.
(At present, these fort.* files contain the values of the T and S
fields only
at the ocean surface, but a comment in newplotform.f
indicates the change needed to plot the fields at any of the other depths).