The Santa Barbara Channel Experiment (SBCX)

 SAIC Software

• First Order Quality Control Check Program Time Check:

This program reads a FFP file that is in CASDE format and writes the start
time of the file to a disk based log file. It is written in C and requires
the "casde.h" include file. One can use this program in conjuction with the 2 Unix scripts:
"dumphd.cmd" and "timeit.cmd" to batch process an entire FFP tape (17 files). (See below)

• Second Order Quality Control Check Program

It comes in 2 flavors: C and Fortran 90.
C Version of the Program:Block Check (C):
Fortran 90 Version of the Program:Block Check (Fortran 90):
This program reads a FFP file that is in CASDE format and writes the start
time as well as the total number of Acoustic and Non-Acoustic sensor IDs that
were recorded within the file. The results are written to a disk based log file.
These programs require the following
include files: "casde.h" and "bbnlook.inc" respectively. (See below)

• Third Order Quality Control Check Program:
  Pre-Processing/Diagnostics (with byteswapping) Diagnostics (bbnlook.f90)

This program can perform several diagnostic functions on the FFP data files.
The functions include:
1. Time series plots of LF and HF data from a single VLA
2. FFT spectral analysis plots of Phone v Frequency for LF and HF data within a single VLA
3. Spectrogram plots of LF and HF data from a single phone within a single VLA
4. Text or Binary files of the spectra from each phone over time from a single VLA
5. Peak spectral values (files) for LF data (used for shape estimation from tilt sensors)
6. Average spectral values (files) for LF data (used for shape estimation from tilt sensors)
7. K-Omega surface (plot) for LF data from a single array
8. Vertical Angle v Frequency surface (plot) for LF data from a single array
9. Continuity check (plot) of LF reference tone (used to confirm good vs dropped data)
10. Ping detections (plot and file) from HF data in a single VLA (used for shape estimation from pingers)
11. Compass data logger (file) for all VLAs (used to support the AEL derived from the tilt sensors)
12. Averaged Complex FFT Spectral Samples (file) for (used to create covariance matrices by a separate program)

For this program to work, FFP tapes are first dumped from their DLT 7000
media to disk. Each file requires about 1.79 GB of space. This represents
about 6.25 minutes of data. (See "dumpit.cmd" below).

The output plot files can be subsequently viewed with a public domain
Windows 95/NT graphics tool called DPLOT95 .
The plots can also be viewed with MatLab:
MatLab XY Graph viewerMatLab XY Graph viewer
MatLab 2D Surface viewerMatLab 2D Surface viewer
This program comes in 1 flavor: Fortran 90, and it requires
the include file: "bbnlook.inc" and a set of subroutines in the file "looksubs.f90"
A sample input parameter file is in a file "bbnlook.in". (See below).

• Fourth Order Quality Control Check Program

Non-Acoustic Data (NAD) Check Program:NAD Check:
This program reads a FFP file that is in CASDE format and writes the summary
information as well as the NAD blocks that it encounters inside the file.
The results are written to the screen or can be re-directed to an output disk based log file.
One can use this program to quickly evaluate the integrity of the compass data strings.
Some where between FFP tape #29 and FFP tape #34, the compass data strings
were fixed, such that no buffer corruption occurred. Prior to FFP tape #26,
no compass data was available on tape at all.
The program requires the following file:
"casde.h" (See below)

• A data compression program is available in 2 flavors (C and Fortran 90):

The "lfonly" program will strip out only the LF data and NAD data
(if it is contained within the LF frame) and write out a new file that
is about one-half as big as the original file.
There are 4 parameters used to drive the program:
1. Input file name
2. Output file name
3. Number of frames to skip (Max = 675309)
4. Number of frames to process (Max = 675309)

Data Compression Program (C Version LF data output only): Data Compression Program (LF data output only) (lfonly.c)) :
Data Compression Program (Fortran Version LF data output only): Data Compression Program (LF data output only) (lfonly.f90)) :
Data Compression Program Parameter File (Fortran only): Data Compression Program Parameter File (Fortran only) (lfonly.in)) :
Data Compression Program Linking File (Fortran only): Data Compression Program Linking File (Fortran only) (lfonly.lnk)) :

• Support files for the programs described above are located at:

CASDE header file (C language): CASDE C header file (casde.h) :
Pre-Processing/Diagnostics include file (Fortran 90): Pre-Processing/Diagnostics F90 include file (bbnlook.inc) :
Pre-Processing/Diagnostic F90 Subroutines (with byteswapping) Diagnostic F90 Subroutines (looksubs.f90)
Lahey Fortran specific routines Lahey Fortran specific Subroutines (lahey.f90)
SGI specific routines SGI specific Subroutines (sgi.f90)
Pre-Processing/Diagnostics Linking file (Fortran 90): Pre-Processing/Diagnostics F90 Linking file (bbnlook.lnk) :
This program comes in 1 flavor: Fortran 90, and it requires an input Parameter file: bbnlook.in :
MatLab XY Graph viewerMatLab XY Graph viewer
MatLab 2D Surface viewerMatLab 2D Surface viewer
Unix DLT tape dumping script: Unix Tape Dump Script (dumpit.cmd) :
Unix FFP file header dumping script: Unix FFP File Header Dump Script (dumphd.cmd) :
Unix Time Check script: Unix Time Check script (timeit.cmd) :
Unix DLT tape Dubbing script: Unix Tape dubbing script (dubit.cmd) :
DLT Tape Copying Program for SUN: Tape Copying Program for SUN (copyit.c)) :
DLT Tape Copying Subroutines for SUN: Tape Copying Subroutines for SUN (sun_ops.c)) :
DLT Tape Copying Header file for SUN: Tape Copying Header file for SUN (sun_ops.h)) :
Solaris SCSI Config File: Solaris SCSI Config File (st.conf) :

How to Run these programs

• Directions for running these programs can be found in the following files:

Unix DLT tape dumping script: Unix Tape Dump Script (dumpit.cmd) :
Unix FFP file header dumping script: Unix FFP File Header Dump Script (dumphd.cmd) :
Unix Time Check script: Unix Time Check script (timeit.cmd) :
Unix Time Check program: Unix Time Check program (bbntime.txt) :
Unix Data Gap program: Unix Data Gap program (bbngap.txt) :
Sample Unix Session: Sample Unix Session (qc.txt) :
Unix NAD Check program: Unix Time NAD check program (bbnnad.txt) :
Unix DLT tape Dubbing script: Unix Tape dubbing script (dubit.cmd) :
F90 Diagnostic Program: Diagnostics (bbnlook.txt) :
DLT Tape Copying Program for SUN: Tape Copying Program for SUN (copyit.txt)) :

AEL

To download the Slope Sensor Based AEL Files go to: Slope Sensor Based AEL Files

• A Fortran 90 program is available to convert from UTM coordinates to Lat/Lon coordinates (and vice versa):

1. Main Demonstration Program: Main Demo Program (utmlltst.f90):
2. Conversion Subroutines: Conversion Subroutines (utm2ll.f90):
2. Linking File: Conversion Subroutines (utmlltst.lnk):

• A Fortran 90 program is available to Average the Noisy Compass Data:

1. Program to Average Compass data: Program to Average Compass data (avgnad.f90):

The output of this program is needed by the MatLab AEL program (see the .avc file below).

• Some MatLab programs are available to determine the AEL of 5 FFP arrays:

The only method available currently uses the "slope sensors" on each
FFP array (there are 30 of them within each array) to derive the AEL.
There are 6 types of files needed to complete the processing:
1. Main Program: Main MatLab Program (slopeav.m):
2. Compass Data Reader: Compass Data Reader (avghdt.m):
3. Slope Sensor Calibration File: Slope Sensor Calibration File (ffa.cal):
4. Slope Sensor Serial Number File: Slope Sensor Serial Number File (ffahydsn.dat):
5. Average Compass Data File: Average Compass Data File (f344.avc):
6a. Peak Frequencies File: Peak Frequencies Data File (f3441lff.mpb):
6b. Peak Frequencies File: Peak Frequencies Data File (f3442lff.mpb):
6c. Peak Frequencies File: Peak Frequencies Data File (f3443lff.mpb):
6d. Peak Frequencies File: Peak Frequencies Data File (f3444lff.mpb):
6e. Peak Frequencies File: Peak Frequencies Data File (f3445lff.mpb):

The output of each run are 4 JPEG files and 1 ASCII file that contains the (xyz) coordinates of each phone of each array.
The JPEG files can be viewed with a WWW browser. The ASCII file can be edited.
The files are as follows:
1. X v Z Plot: X v Z Plot (south_sl.jpg)
2. Y v Z Plot: X v Z Plot (east_sl.jpg)
3. Map View Plot: Map View Plot (map_sl.jpg)
4. Water Temperature v Z Plot: Water Temperature v Z Plot (wtemp_sl.jpg)
5. X,Y,Z (Tape 34,File 4): X,Y,Z (Tape 34, File 4)

• Array Shape Animations:

1. X2 in Map View: X2 in Map View
2. X2 from a Southerly View: X2 Southerly View

• Time History of FFP VLA motion:

Map View of X2 Time Period: X2 Map Views
Map View of X9 Time Period: X9 Map Views

BATHYMETRY

• Bathymetry Extraction Program for the SBCX area:

There is a Fortran 90 program that can extract bathymetric values
from a database file for the SBCX area.

There are 5 files needed to complete this processing:

1. Main Program: Main F90 Program (getbath.f90):
2. Geodesy subroutines: Geodesy subroutines (geodesy.f90):
3. Geodesy Include File: Geodesy Include File (geodesy.inc):
4. Parameter File: Parameter File (getbath.in):
5. SBCX Bathymetry Data File: SBCX Bathymetry Data File (ASCII) (sbcxbath.dos):

SBCX Bathymetry in Colorfill Map View
SBCX Bathymetry in Contour Map View
SBCX Bathymetry (Zoomed Version) in Contour Map View
SBCX Bathymetry in 3D Perspective View

Last updated: 000703  Comments/Questions: pmd@alum.mit.edu