#!/bin/sh

#####
#	this is a test of the Trillium SN 0036
#	Cal pulses are put through individual sensors in the UVW mode
#	and also in the XYZ mode: This permits a test of the uvwtril option
#	introduced to gsac on March 30, 2009
#####
# The 0036 trillium was jumpered into UVW mode and
# 0.78v positive steps were applied to the 13.6KOhm cal input.
# from 17:39 to 1745 a pulse was applied to all components.
# i.e. all cal enables were enabled.
# 
# from 18:14 to 18:20 a pulse was applied with the U cal enabled
# data on the HNE channel
# from 18:34 to 18:40 a pulse was applied with the V cal enabled.
# data on the HNN channel.
# from 19:04 to 19:10 a pulse was applied with the W cal enabled.
# data on the HNZ  channel.
# About 19:  the seis was switched to XYZ mode.
# from 19:24 to 19:29 a pulse was applied with the W cal enabled
# 
# from 20:04 to 20:10 a pulse was applied with the V cal enabled.
# 
# from 20:24 to 20:30 a pulse was applied with the U cal enabled.
# data on ew-slu1 under the directory /cnmsn/FVM/  currently the active file. 
#######
mkdir TEST
gsac << EOF
# 18 14 59 U cal in UVW mode
cut cal 2009 04 01 18 14 49 0000 cal 2009 04 01 18 19 10 000
r Sac/*
cd TEST
w append .UVW.Ucal
mv FVMHNE.NM..S.UVW.Ucal FVMHNU.NM..S.UVW.Ucal
mv FVMHNN.NM..S.UVW.Ucal FVMHNV.NM..S.UVW.Ucal
mv FVMHNZ.NM..S.UVW.Ucal FVMHNW.NM..S.UVW.Ucal
cd ..

# 18 35 00 V cal in UVW mode
cut cal 2009 04 01 18 34 50 0000 cal 2009 04 01 18 39 10 000
r Sac/*
cd TEST
w append .UVW.Vcal
mv FVMHNE.NM..S.UVW.Vcal FVMHNU.NM..S.UVW.Vcal
mv FVMHNN.NM..S.UVW.Vcal FVMHNV.NM..S.UVW.Vcal
mv FVMHNZ.NM..S.UVW.Vcal FVMHNW.NM..S.UVW.Vcal
cd ..


# 19 14 05 W cal in UVW mode
cut cal 2009 04 01 19 04 50 0000 cal 2009 04 01 19 09 01 000
r Sac/*
cd TEST
w append .UVW.Wcal
mv FVMHNE.NM..S.UVW.Wcal FVMHNU.NM..S.UVW.Wcal
mv FVMHNN.NM..S.UVW.Wcal FVMHNV.NM..S.UVW.Wcal
mv FVMHNZ.NM..S.UVW.Wcal FVMHNW.NM..S.UVW.Wcal
cd ..

# 19 25 00 W cal in XYZ mode
cut cal 2009 04 01 19 24 50 0000 cal 2009 04 01 19 29 03 000
r Sac/*
cd TEST
w append .XYZ.Wcal
rotate3 to uvwtril
w
cut off
r FVMHNU FVMHNV FVMHNW
w append .xyz.Wcal
cd ..

# 20 04 59 V cal in XYZ mode
cut cal 2009 04 01 20 04 49 0000 cal 2009 04 01 20 09 01 000
r Sac/*
cd TEST
w append .XYZ.Vcal
rotate3 to uvwtril
w
cut off
r FVMHNU FVMHNV FVMHNW
w append .xyz.Vcal
cd ..

# 20 24 51 U cal in XYZ mode
cut cal 2009 04 01 20 24 41 000 cal 2009 04 01 20 29 02 000
r Sac/*
cd TEST
w append .XYZ.Ucal
rotate3 to uvwtril
w
cut off
r FVMHNU FVMHNV FVMHNW
w append .xyz.Ucal
cd ..

#####
#	now do the graphics
#####
cd TEST
fileid name
r *UVW.Ucal *xyz.Ucal
ylim all
TITLE ON L Top S S Text "U sensor calibration: (top) UVW recordings, (bot) UVW from xyz recording"
bg plt
p relative
r *UVW.Vcal *xyz.Vcal
TITLE ON L Top S S Text "V sensor calibration: (top) UVW recordings, (bot) UVW from xyz recording"
p relative
r *UVW.Wcal *xyz.Wcal
TITLE ON L Top S S Text 'W sensor calibration: (top) UVW recordings, (bot) UVW from xyz recording'
p relative

quit
EOF