MAKE_DEPEND=yes
include ${SEPINC}/SEP.top
#include FIGMakefilegrey
include FIGMakefile

LATOPTS  = style=html
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
#LATOPTS=class=geophysics
SIGDIR = /mnt/sepdata/1a/sigsbee_multiple
OLEVEL=2
#UF90FLAGS =-mp -check bounds
UF90FLAGS =-parallel -openmp -openmp_report2  -O3 -axW -vec_report3
MAINLIBS=-lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -L${LIBDIR}  -lsepgeef90
UF90LIBS=${MPILIBS} ${MAINLIBS} -lsepparf90 -lseppar_mpi -lsepgeef90 -lsep2df90 -static-libcxa -lguide#-Bstatic
UF77LIBS = ${BEILIB}

COLOR=yes

RESULTSER = join16 join24 join50 join00 join15 join30 join40 Ijoin16 Ijoin24 Ijoin50 Ijoin00 Ijoin15 Ijoin30 Ijoin40 Astk Cstk Stk jDoff01 jDoff02 jDoff03 Ojoin00 Ojoin15 Ojoin16 Ojoin24 Ojoin50 Ojoin30 Ojoin40 Sigsbee

BIN=./Bin
PAR=./Par
RESDIR=~claudio/Research/Ill_ADCIG/Fig/New
FIG=~claudio/Research/Ill_ADCIG/Fig/New
DATA=./Data
SRC=./Src

#UF90LIBS=-lsep2df90 -lsepparf90 -lseppar -lsupersetf90 -lsuperset -lsepgeef90
#UF90FLAGS=-O1

########
# Data

# LtL_Sis_az1_ax1_wave_marine_old_off.H
# mig_Sis_Parallel_pad_2048_off.H


ODCIG_lim.H:
	Transp reshape=2,4,5 plane=13 < ~valencia/research/sep124/mig_Sis_Parallel_pad_800.H |Window3d squeeze=n min1=7975 min4=32000 n4=70 > $@

ADCIG_orig.H:
	<ODCIG_lim.H Pad end1=100 beg2=43 end2=44 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 >junk.H; <junk.H Window squeeze=n min2=-20 > $@
	rm -f junk.H

ILL_off.H:
#	cp ~valencia/research/sep124/LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.H .
#	sftp gomorrah; get /scr4/alejandro/LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.H@;
	<LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.H Window n1=382 squeeze=n |Transp reshape=1,2,4 plane=23 > $@

ADCIG_mult.H:
	Math file1=ODCIG_lim.H file2=ILL_off.H exp="file1*file2" > junk.H;<junk.H Pad end1=100 beg2=43 end2=44 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y \
	eps=30 >junk1.H; <junk1.H Window squeeze=n min2=-20 > $@
	rm -f junk*.H

ADCIG_const.H:
	Math file1=ODCIG_lim.H exp="10*file1" > junk.H;<junk.H Pad end1=100 beg2=43 end2=44 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y \
	eps=30 >junk1.H; <junk1.H Window squeeze=n min2=-20 > $@
	rm -f junk*.H

Wgh_ang_eps001.H:
	 Math file1=ADCIG_mult.H file2=ADCIG_orig.H exp="file1*file2/(file2*file2+0.001)" > $@

Wgh_angcte_eps001.H:
	 Math file1=ADCIG_const.H file2=ADCIG_orig.H exp="file1*file2/(file2*file2+0.001)" > $@

Wgh_ang_eps005.H:
	 Math file1=ADCIG_mult.H file2=ADCIG_orig.H exp="file1*file2/(file2*file2+0.005)" > $@

Wgh_ang_eps0.H:
	 Math file1=ADCIG_mult.H file2=ADCIG_orig.H exp="file1/file2" > $@

Wgh_ang_eps01.H:
	 Math file1=ADCIG_mult.H file2=ADCIG_orig.H exp="file1*file2/(file2*file2+0.01)" > $@

Wgh_ang_eps001_med11.H:
	Transp plane=34 < Wgh_ang_eps001.H | Median wind1=11 wind2=1  > $@

Wgh_ang_eps005_med11.H:
	Transp plane=34 < Wgh_ang_eps005.H | Median wind1=11 wind2=1  > $@

Wgh_ang_eps005_med31.H:
	Transp plane=34 < Wgh_ang_eps005.H | Median wind1=31 wind2=1  > $@

Wgh_ang_eps0_med51.H:
	Transp plane=34 < Wgh_ang_eps0.H | Median wind1=51 wind2=1  > $@

Wgh_ang_eps01_med11.H:
	Transp plane=34 < Wgh_ang_eps01.H | Median wind1=11 wind2=1  > $@

Wgh_ang_eps005_med2107.H:
	Transp plane=34 < Wgh_ang_eps005.H | Transp plane=23 | Median wind1=21 wind2=7  > $@

Wgh_angcte_eps001_med2107.H:
	Transp plane=34 < Wgh_angcte_eps001.H | Transp plane=23 | Median wind1=21 wind2=7  > $@

ADCIG_corr_eps005_01.H:
	Transp plane=34 < ADCIG_orig.H > junk.H ; Math file1=junk.H file2=Wgh_ang_eps005_med11.H exp="file1*file2/(file2*file2+0.1)" > $@

ADCIG_corr_eps005_2107_01.H:
	Transp plane=34 < ADCIG_orig.H |Transp plane=23 > junk.H ; Math file1=junk.H file2=Wgh_ang_eps005_med2107.H exp="file1*file2/(file2*file2+0.1)" > $@

ADCIG_corr_eps005_2107_001.H:
	Transp plane=34 < ADCIG_orig.H |Transp plane=23 > junk.H ; Math file1=junk.H file2=Wgh_ang_eps005_med2107.H exp="file1*file2/(file2*file2+0.01)" > $@

CorrStk_eps005_2107_01.H:
	<ADCIG_corr_eps005_2107_01.H Transp plane=23 | Stack > $@


##########################################################
#
# Computations using the correct diagonal of the Hessian
#

001-ODCIG.H:
	Window min1=7975 min4=32000 n4=50 j2=3 < mig_Sis_Parallel_pad_2048_off.H > $@

002-ILLOFF.H:
	Window n2=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag.H | Transp reshape=1,3 plane=12 > $@

003-ILLANG_prev.H:
	Transp plane=34 < 002-ILLOFF.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 > junk1.H; <junk1.H Window max1=16000 > $@
	Rm junk1.H

004-ILLANG.H:
	Math file1=001-ODCIG.H file2=002-ILLOFF.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=21 wind2=7 | Transp plane=23 > $@
	Rm j*H


004-ILLANG1.H:
	Math file1=001-ODCIG.H file2=002-ILLOFF.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > $@
	Rm j*H

004-ILLANG1a.H:
	Math file1=001-ODCIG.H file2=002-ILLOFF.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Math file1=junk1.H file2=junk2.H file3=junk3.H exp="file1*file2/file3" > junk4.H
	Transp plane=23 < junk4.H | Median wind1=51 wind2=5 | Transp plane=23 > $@

004-ILLANG2.H:
	Math file1=001-ODCIG.H exp="file1*10." > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > $@
	Rm j*H

005-ADCIG.H:
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > $@

001-ODCIG_inv.H:
	<005-ADCIG.H Pad beg2=43 end2=44 |Transp plane=34 | /net/koko/bob/SEP_BOB/bin/OFF2ANG inverse=y hmax=1200 verb=y eps=1 maxmem=1000 > $@

005-ADCIG_eps01.H:
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=1 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > $@

005-ADCIG_diag.H:
	Math file1=001-ODCIG.H file2=002-ILLOFF.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > $@
	Rm funk1.H j1.H

006-ADCIG_corr.H:
	Math file1=005-ADCIG.H file2=004-ILLANG1.H exp="file1*file2/(file2*file2+0.01)" > $@

diagonal_off:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=0 n6=1 f6=0 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=1 n6=1 f6=1 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=2 n6=1 f6=2 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=3 n6=1 f6=3 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=4 n6=1 f6=4 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=5 n6=1 f6=5 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=6 n6=1 f6=6 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=7 n6=1 f6=7 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=8 n6=1 f6=8 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=9 n6=1 f6=9 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=10 n6=1 f6=10 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=11 n6=1 f6=11 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=12 n6=1 f6=12 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=13 n6=1 f6=13 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=14 n6=1 f6=14 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=15 n6=1 f6=15 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=16 n6=1 f6=16 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=17 n6=1 f6=17 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=18 n6=1 f6=18 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=19 n6=1 f6=19 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=20 n6=1 f6=20 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=21 n6=1 f6=21 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=22 n6=1 f6=22 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=23 n6=1 f6=23 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=24 n6=1 f6=24 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=25 n6=1 f6=25 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=26 n6=1 f6=26 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=27 n6=1 f6=27 > n27.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=28 n6=1 f6=28 > n28.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=29 n6=1 f6=29 > n29.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=30 n6=1 f6=30 > n30.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=31 n6=1 f6=31 > n31.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f4=32 n6=1 f6=32 > n32.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      n30.H n31.H n32.H axis=3 >LtL_Sis_az1_ax1_wave_marine_old_off_diag.H
	echo o3=-1200. >> LtL_Sis_az1_ax1_wave_marine_old_off_diag.H
#	Rm n*.H


diagonal01:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=1 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=2 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=3 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=4 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=5 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=6 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=7 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=8 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=9 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=10 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=11 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=12 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=13 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=14 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=15 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=16 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=17 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=18 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=19 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=20 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=21 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=22 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=23 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=24 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=25 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=26 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=27 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=28 > n27.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=28 n6=1 f4=29 > n28.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=29 n6=1 f4=30 > n29.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=30 n6=1 f4=31 > n30.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=31 n6=1 f4=32 > n31.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      n30.H n31.H axis=3 > junk1.H
	echo o3=-1125. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=1 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag01.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF01.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG01.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag01.H > 002-DIAGOFF01.H
	Rm j*H
	Rm n*.H junk1.H

diagonal02:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=2 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=3 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=4 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=5 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=6 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=7 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=8 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=9 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=10 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=11 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=12 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=13 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=14 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=15 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=16 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=17 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=18 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=19 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=20 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=21 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=22 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=23 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=24 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=25 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=26 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=27 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=28 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=29 > n27.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=28 n6=1 f4=30 > n28.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=29 n6=1 f4=31 > n29.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=30 n6=1 f4=32 > n30.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      n30.H axis=3 > junk1.H
	echo o3=-1125. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=1 end2=1 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag02.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag02.H > 002-DIAGOFF02.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF02.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG02.H
	Rm j*H
	Rm n*.H junk1.H

diagonal03:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=3 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=4 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=5 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=6 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=7 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=8 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=9 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=10 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=11 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=12 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=13 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=14 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=15 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=16 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=17 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=18 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=19 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=20 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=21 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=22 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=23 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=24 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=25 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=26 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=27 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=28 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=29 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=30 > n27.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=28 n6=1 f4=31 > n28.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=29 n6=1 f4=32 > n29.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      axis=3 > junk1.H
	echo o3=-1050. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=2 end2=1 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag03.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag03.H > 002-DIAGOFF03.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF03.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG03.H
	Rm j*H
	Rm n*.H junk1.H

diagonal04:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=4 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=5 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=6 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=7 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=8 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=9 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=10 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=11 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=12 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=13 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=14 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=15 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=16 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=17 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=18 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=19 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=20 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=21 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=22 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=23 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=24 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=25 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=26 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=27 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=28 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=29 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=30 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=31 > n27.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=28 n6=1 f4=32 > n28.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H \
	      axis=3 > junk1.H
	echo o3=-1050. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=2 end2=2 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag04.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag04.H > 002-DIAGOFF04.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF04.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG04.H
	Rm j*H
	Rm n*.H junk1.H

diagonal05:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=5 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=6 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=7 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=8 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=9 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=10 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=11 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=12 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=13 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=14 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=15 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=16 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=17 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=18 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=19 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=20 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=21 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=22> n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=23 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=24 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=25 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=26 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=27 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=28 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=29 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=30 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=31 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=32 > n27.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H \
	      axis=3 > junk1.H
	echo o3=-975. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=3 end2=2 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag05.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag05.H > 002-DIAGOFF05.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF05.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG05.H
	Rm j*H
	Rm n*.H junk1.H

diagonal05a:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=5 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=6 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=7 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=8 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=9 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=10 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=11 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=12 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=13 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=14 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=15 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=16 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=17 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=18 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=19 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=20 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=21 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=22> n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=23 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=24 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=25 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=26 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=27 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=28 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=29 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=30 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=31 > n26.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=27 n6=1 f4=32 > n27.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H \
	      axis=3 > junk1.H
	echo o3=-975. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=3 end2=2 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag05.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag05.H > 002-DIAGOFF05.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF05.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk2.H exp="file1*file1" | Transp plane=23 | Median wind1=7 wind2=3 | Transp plane=23 > junk3.H
	Math file1=junk1.H file2=junk2.H file3=junk3.H exp="file1*file2/file3" > junk4.H
	Transp plane=23 < junk4.H | Median wind1=21 wind2=11 | Transp plane=23 > 004-DIAGANG05a.H
#	Rm j*H
#	Rm n*.H junk1.H

diagonal06:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=6 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=7 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=8 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=9 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=10 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=11 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=12 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=13 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=14 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=15 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=16 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=17 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=18 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=19 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=20 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=21 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=22 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=23 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=24 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=25 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=26 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=27 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=28 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=29 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=30 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=31 > n25.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=26 n6=1 f4=32 > n26.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H \
	      axis=3 > junk1.H
	echo o3=-975. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=3 end2=3 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag06.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag06.H > 002-DIAGOFF06.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF06.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG06.H
	Rm j*H
	Rm n*.H junk1.H

diagonal07:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=7 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=8 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=9 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=10 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=11 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=12 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=13 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=14 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=15 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=16 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=17 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=18 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=19 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=20 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=21 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=22 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=23 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=24 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=25 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=26 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=27 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=28 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=29 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=30 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=31 > n24.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=25 n6=1 f4=32 > n25.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H \
	      axis=3 > junk1.H
	echo o3=-900. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=4 end2=3 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag07.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag07.H > 002-DIAGOFF07.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF07.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG07.H
	Rm j*H
	Rm n*.H junk1.H

diagonal08:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=8 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=9 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=10 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=11 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=12 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=13 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=14 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=15 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=16 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=17 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=18 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=19 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=20 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=21 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=22 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=23 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=24 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=25 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=26 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=27 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=28 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=29 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=30 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=31 > n23.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=24 n6=1 f4=32 > n24.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H \
	      axis=3 > junk1.H
	echo o3=-900. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=4 end2=4 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag08.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag08.H > 002-DIAGOFF08.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF08.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG08.H
	Rm j*H
	Rm n*.H junk1.H

diagonal09:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=9 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=10 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=11 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=12 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=13 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=14 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=15 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=16 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=17 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=18 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=19 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=20 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=21 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=22 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=23 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=24 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=25 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=26 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=27 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=28 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=29 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=30 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=31 > n22.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=23 n6=1 f4=32 > n23.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H \
	      axis=3 > junk1.H
	echo o3=-825. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=5 end2=4 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag09.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag09.H > 002-DIAGOFF09.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF09.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG09.H
	Rm j*H
	Rm n*.H junk1.H

diagonal10:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=10 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=11 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=12 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=13 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=14 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=15 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=16 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=17 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=18 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=19 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=20 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=21 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=22 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=23 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=24 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=25 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=26 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=27 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=28 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=29 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=30 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=31 > n21.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=22 n6=1 f4=32 > n22.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H \
	      axis=3 > junk1.H
	echo o3=-825. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=5 end2=5 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag10.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag10.H > 002-DIAGOFF10.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF10.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG10.H
	Rm j*H
	Rm n*.H junk1.H

diagonal11:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=11 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=12 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=13 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=14 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=15 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=16 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=17 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=18 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=19 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=20 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=21 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=22 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=23 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=24 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=25 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=26 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=27 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=28 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=29 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=30 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=31 > n20.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=21 n6=1 f4=32 > n21.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H \
	      axis=3 > junk1.H
	echo o3=-750. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=6 end2=5 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag11.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag11.H > 002-DIAGOFF11.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF11.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG11.H
	Rm j*H
	Rm n*.H junk1.H

diagonal12:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=12 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=13 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=14 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=15 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=16 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=17 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=18 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=19 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=20 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=21 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=22 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=23 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=24 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=25 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=26 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=27 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=28 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=29 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=20 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=31 > n19.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=20 n6=1 f4=32 > n20.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H \
	      axis=3 > junk1.H
	echo o3=-750. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=6 end2=6 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag12.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag12.H > 002-DIAGOFF12.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF12.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG12.H
	Rm j*H
	Rm n*.H junk1.H

diagonal13:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=13 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=14 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=15 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=16 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=17 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=18 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=19 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=20 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=21 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=22 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=23 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=24 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=25 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=26 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=27 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=28 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=29 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=30 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=31 > n18.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=19 n6=1 f4=32 > n19.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	      axis=3 > junk1.H
	echo o3=-675. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=7 end2=6 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag13.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag13.H > 002-DIAGOFF13.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF13.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG13.H
	Rm j*H
	Rm n*.H junk1.H

diagonal14:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=14 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=15 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=16 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=17 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=18 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=19 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=20 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=21 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=22 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=23 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=24 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=25 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=26 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=27 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=28 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=29 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=30 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=31 > n17.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=18 n6=1 f4=32 > n18.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H \
	      axis=3 > junk1.H
	echo o3=-675. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=7 end2=7 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag14.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag14.H > 002-DIAGOFF14.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF14.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG14.H
	Rm j*H
	Rm n*.H junk1.H

diagonal15:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=15 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=16 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=17 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=18 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=19 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=20 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=21 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=22 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=23 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=24 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=25 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=26 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=27 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=28 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=29 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=30 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=31 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=32 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > junk1.H
	echo o3=-600. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=8 end2=7 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag15.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag15.H > 002-DIAGOFF15.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF15.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG15.H
	Rm j*H
	Rm n*.H junk1.H


diagonal15a:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=15 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=16 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=17 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=18 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=19 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=20 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=21 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=22 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=23 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=24 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=25 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=26 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=27 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=28 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=29 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=30 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=31 > n16.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=17 n6=1 f4=32 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > junk1.H
	echo o3=-600. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=8 end2=7 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag15.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag15.H > 002-DIAGOFF15.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF15.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Math file1=junk1.H file2=junk2.H file3=junk3.H exp="file1*file2/file3" > junk4.H
	Transp plane=23 < junk4.H | Median wind1=11 wind2=5 | Transp plane=23 > 004-DIAGANG15a.H
#	Rm j*H
#	Rm n*.H junk1.H

diagonal16:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=16 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=17 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=18 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=19 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=20 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=21 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=22 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=23 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=24 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=25 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=26 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=27 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=28 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=29 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=30 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=31 > n15.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=16 n6=1 f4=32 > n16.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H \
	      axis=3 > junk1.H
	echo o3=-600. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=8 end2=8 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag16.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag16.H > 002-DIAGOFF16.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF16.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG16.H
	Rm j*H
	Rm n*.H junk1.H

diagonal17:
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=0 n6=1 f4=17 > n0.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=1 n6=1 f4=18 > n1.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=2 n6=1 f4=19 > n2.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=3 n6=1 f4=20 > n3.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=4 n6=1 f4=21 > n4.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=5 n6=1 f4=22 > n5.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=6 n6=1 f4=23 > n6.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=7 n6=1 f4=24 > n7.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=8 n6=1 f4=25 > n8.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=9 n6=1 f4=26 > n9.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=10 n6=1 f4=27 > n10.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=11 n6=1 f4=28 > n11.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=12 n6=1 f4=29 > n12.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=13 n6=1 f4=30 > n13.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=14 n6=1 f4=31 > n14.H
	Window3d <LtL_Sis_az1_ax1_wave_marine_old_off.H n4=1 f6=15 n6=1 f4=32 > n15.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H \
	      axis=3 > junk1.H
	echo o3=-525. >> junk1.H
	Transp reshape=1,3,5 plane=12 <junk1.H |Pad beg2=9 end2=8 extend=1 >LtL_Sis_az1_ax1_wave_marine_old_off_diag17.H
	Window n1=382 < LtL_Sis_az1_ax1_wave_marine_old_off_diag17.H > 002-DIAGOFF17.H
	Math file1=001-ODCIG.H file2=002-DIAGOFF17.H exp="file1*file2" > funk1.H
	Transp plane=34 < funk1.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j1.H ; <j1.H Window max1=16000 > junk1.H
	Transp plane=34 < 001-ODCIG.H | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 maxmem=1000 >j2.H ; <j2.H Window max1=16000 > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1*file2/(file2*file2+0.005)" > junk3.H
	Transp plane=23 < junk3.H | Median wind1=31 wind2=11 | Transp plane=23 > 004-DIAGANG17.H
#	Rm j*H
#	Rm n*.H junk1.H


Hessiang:
	Cat 004-ILLANG1.H 004-DIAGANG01.H 004-DIAGANG02.H 004-DIAGANG03.H 004-DIAGANG04.H 004-DIAGANG05.H 004-DIAGANG06.H 004-DIAGANG07.H 004-DIAGANG08.H 004-DIAGANG09.H 004-DIAGANG10.H 004-DIAGANG11.H 004-DIAGANG12.H 004-DIAGANG13.H 004-DIAGANG14.H 004-DIAGANG15.H 004-DIAGANG16.H 004-DIAGANG17.H > 007-HESSIANG.H

Hess01:
	Window3d <004-ILLANG1.H   n3=1 f3=0 > n0.H
	Window3d <004-DIAGANG01.H n3=1 f3=0 > n1.H
	Window3d <004-DIAGANG02.H n3=1 f3=0 > n2.H
	Window3d <004-DIAGANG03.H n3=1 f3=0 > n3.H
	Window3d <004-DIAGANG04.H n3=1 f3=0 > n4.H
	Window3d <004-DIAGANG05.H n3=1 f3=0 > n5.H
	Window3d <004-DIAGANG06.H n3=1 f3=0 > n6.H
	Window3d <004-DIAGANG07.H n3=1 f3=0 > n7.H
	Window3d <004-DIAGANG08.H n3=1 f3=0 > n8.H
	Window3d <004-DIAGANG09.H n3=1 f3=0 > n9.H
	Window3d <004-DIAGANG10.H n3=1 f3=0 > n10.H
	Window3d <004-DIAGANG11.H n3=1 f3=0 > n11.H
	Window3d <004-DIAGANG12.H n3=1 f3=0 > n12.H
	Window3d <004-DIAGANG13.H n3=1 f3=0 > n13.H
	Window3d <004-DIAGANG14.H n3=1 f3=0 > n14.H
	Window3d <004-DIAGANG15.H n3=1 f3=0 > n15.H
	Window3d <004-DIAGANG16.H n3=1 f3=0 > n16.H
	Window3d <004-DIAGANG17.H n3=1 f3=0 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > 007-HESSIAN01.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> 007-HESSIAN01.H
	Rm n*.H

Hess16:
	Window3d <004-ILLANG1.H   n3=1 f3=17 > n0.H
	Window3d <004-DIAGANG01.H n3=1 f3=17 > n1.H
	Window3d <004-DIAGANG02.H n3=1 f3=17 > n2.H
	Window3d <004-DIAGANG03.H n3=1 f3=17 > n3.H
	Window3d <004-DIAGANG04.H n3=1 f3=17 > n4.H
	Window3d <004-DIAGANG05.H n3=1 f3=17 > n5.H
	Window3d <004-DIAGANG06.H n3=1 f3=17 > n6.H
	Window3d <004-DIAGANG07.H n3=1 f3=17 > n7.H
	Window3d <004-DIAGANG08.H n3=1 f3=17 > n8.H
	Window3d <004-DIAGANG09.H n3=1 f3=17 > n9.H
	Window3d <004-DIAGANG10.H n3=1 f3=17 > n10.H
	Window3d <004-DIAGANG11.H n3=1 f3=17 > n11.H
	Window3d <004-DIAGANG12.H n3=1 f3=17 > n12.H
	Window3d <004-DIAGANG13.H n3=1 f3=17 > n13.H
	Window3d <004-DIAGANG14.H n3=1 f3=17 > n14.H
	Window3d <004-DIAGANG15.H n3=1 f3=17 > n15.H
	Window3d <004-DIAGANG16.H n3=1 f3=17 > n16.H
	Window3d <004-DIAGANG17.H n3=1 f3=17 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > 007-HESSIAN16.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> 007-HESSIAN16.H
	Rm n*.H

Hess24:
	Window3d <004-ILLANG1.H   n3=1 f3=25 > n0.H
	Window3d <004-DIAGANG01.H n3=1 f3=25 > n1.H
	Window3d <004-DIAGANG02.H n3=1 f3=25 > n2.H
	Window3d <004-DIAGANG03.H n3=1 f3=25 > n3.H
	Window3d <004-DIAGANG04.H n3=1 f3=25 > n4.H
	Window3d <004-DIAGANG05.H n3=1 f3=25 > n5.H
	Window3d <004-DIAGANG06.H n3=1 f3=25 > n6.H
	Window3d <004-DIAGANG07.H n3=1 f3=25 > n7.H
	Window3d <004-DIAGANG08.H n3=1 f3=25 > n8.H
	Window3d <004-DIAGANG09.H n3=1 f3=25 > n9.H
	Window3d <004-DIAGANG10.H n3=1 f3=25 > n10.H
	Window3d <004-DIAGANG11.H n3=1 f3=25 > n11.H
	Window3d <004-DIAGANG12.H n3=1 f3=25 > n12.H
	Window3d <004-DIAGANG13.H n3=1 f3=25 > n13.H
	Window3d <004-DIAGANG14.H n3=1 f3=25 > n14.H
	Window3d <004-DIAGANG15.H n3=1 f3=25 > n15.H
	Window3d <004-DIAGANG16.H n3=1 f3=25 > n16.H
	Window3d <004-DIAGANG17.H n3=1 f3=25 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > 007-HESSIAN24.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> 007-HESSIAN24.H
	Rm n*.H

Hess40:
	Window3d <004-ILLANG1.H   n3=1 f3=41 > n0.H
	Window3d <004-DIAGANG01.H n3=1 f3=41 > n1.H
	Window3d <004-DIAGANG02.H n3=1 f3=41 > n2.H
	Window3d <004-DIAGANG03.H n3=1 f3=41 > n3.H
	Window3d <004-DIAGANG04.H n3=1 f3=41 > n4.H
	Window3d <004-DIAGANG05.H n3=1 f3=41 > n5.H
	Window3d <004-DIAGANG06.H n3=1 f3=41 > n6.H
	Window3d <004-DIAGANG07.H n3=1 f3=41 > n7.H
	Window3d <004-DIAGANG08.H n3=1 f3=41 > n8.H
	Window3d <004-DIAGANG09.H n3=1 f3=41 > n9.H
	Window3d <004-DIAGANG10.H n3=1 f3=41 > n10.H
	Window3d <004-DIAGANG11.H n3=1 f3=41 > n11.H
	Window3d <004-DIAGANG12.H n3=1 f3=41 > n12.H
	Window3d <004-DIAGANG13.H n3=1 f3=41 > n13.H
	Window3d <004-DIAGANG14.H n3=1 f3=41 > n14.H
	Window3d <004-DIAGANG15.H n3=1 f3=41 > n15.H
	Window3d <004-DIAGANG16.H n3=1 f3=41 > n16.H
	Window3d <004-DIAGANG17.H n3=1 f3=41 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > 007-HESSIAN40.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> 007-HESSIAN40.H
	Rm n*.H

Hess48:
	Window3d <004-ILLANG1.H   n3=1 f3=49 > n0.H
	Window3d <004-DIAGANG01.H n3=1 f3=49 > n1.H
	Window3d <004-DIAGANG02.H n3=1 f3=49 > n2.H
	Window3d <004-DIAGANG03.H n3=1 f3=49 > n3.H
	Window3d <004-DIAGANG04.H n3=1 f3=49 > n4.H
	Window3d <004-DIAGANG05.H n3=1 f3=49 > n5.H
	Window3d <004-DIAGANG06.H n3=1 f3=49 > n6.H
	Window3d <004-DIAGANG07.H n3=1 f3=49 > n7.H
	Window3d <004-DIAGANG08.H n3=1 f3=49 > n8.H
	Window3d <004-DIAGANG09.H n3=1 f3=49 > n9.H
	Window3d <004-DIAGANG10.H n3=1 f3=49 > n10.H
	Window3d <004-DIAGANG11.H n3=1 f3=49 > n11.H
	Window3d <004-DIAGANG12.H n3=1 f3=49 > n12.H
	Window3d <004-DIAGANG13.H n3=1 f3=49 > n13.H
	Window3d <004-DIAGANG14.H n3=1 f3=49 > n14.H
	Window3d <004-DIAGANG15.H n3=1 f3=49 > n15.H
	Window3d <004-DIAGANG16.H n3=1 f3=49 > n16.H
	Window3d <004-DIAGANG17.H n3=1 f3=49 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > 007-HESSIAN48.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> 007-HESSIAN48.H
	Rm n*.H

HessCA30:
	Transp plane=23  <004-ILLANG1.H   | Window3d n3=1 f3=46 > n0.H
	Transp plane=23  <004-DIAGANG01.H | Window3d n3=1 f3=46 > n1.H
	Transp plane=23  <004-DIAGANG02.H | Window3d n3=1 f3=46 > n2.H
	Transp plane=23  <004-DIAGANG03.H | Window3d n3=1 f3=46 > n3.H
	Transp plane=23  <004-DIAGANG04.H | Window3d n3=1 f3=46 > n4.H
	Transp plane=23  <004-DIAGANG05.H | Window3d n3=1 f3=46 > n5.H
	Transp plane=23  <004-DIAGANG06.H | Window3d n3=1 f3=46 > n6.H
	Transp plane=23  <004-DIAGANG07.H | Window3d n3=1 f3=46 > n7.H
	Transp plane=23  <004-DIAGANG08.H | Window3d n3=1 f3=46 > n8.H
	Transp plane=23  <004-DIAGANG09.H | Window3d n3=1 f3=46 > n9.H
	Transp plane=23  <004-DIAGANG10.H | Window3d n3=1 f3=46 > n10.H
	Transp plane=23  <004-DIAGANG11.H | Window3d n3=1 f3=46 > n11.H
	Transp plane=23  <004-DIAGANG12.H | Window3d n3=1 f3=46 > n12.H
	Transp plane=23  <004-DIAGANG13.H | Window3d n3=1 f3=46 > n13.H
	Transp plane=23  <004-DIAGANG14.H | Window3d n3=1 f3=46 > n14.H
	Transp plane=23  <004-DIAGANG15.H | Window3d n3=1 f3=46 > n15.H
	Transp plane=23  <004-DIAGANG16.H | Window3d n3=1 f3=46 > n16.H
	Transp plane=23  <004-DIAGANG17.H | Window3d n3=1 f3=46 > n17.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H \
	      axis=3 > junk.H
	echo n3=18 o3=1 d3=1 label3=diagonal >> junk.H
	Transp plane=23 < junk.H > 007-HESSIAN_CA30.H
	Rm n*.H junk.H

diagonal_ang:
	Window3d <008-Matrix.H n4=1 f4=0 n6=1 f6=0 > n0.H
	Window3d <008-Matrix.H n4=1 f4=1 n6=1 f6=1 > n1.H
	Window3d <008-Matrix.H n4=1 f4=2 n6=1 f6=2 > n2.H
	Window3d <008-Matrix.H n4=1 f4=3 n6=1 f6=3 > n3.H
	Window3d <008-Matrix.H n4=1 f4=4 n6=1 f6=4 > n4.H
	Window3d <008-Matrix.H n4=1 f4=5 n6=1 f6=5 > n5.H
	Window3d <008-Matrix.H n4=1 f4=6 n6=1 f6=6 > n6.H
	Window3d <008-Matrix.H n4=1 f4=7 n6=1 f6=7 > n7.H
	Window3d <008-Matrix.H n4=1 f4=8 n6=1 f6=8 > n8.H
	Window3d <008-Matrix.H n4=1 f4=9 n6=1 f6=9 > n9.H
	Window3d <008-Matrix.H n4=1 f4=10 n6=1 f6=10 > n10.H
	Window3d <008-Matrix.H n4=1 f4=11 n6=1 f6=11 > n11.H
	Window3d <008-Matrix.H n4=1 f4=12 n6=1 f6=12 > n12.H
	Window3d <008-Matrix.H n4=1 f4=13 n6=1 f6=13 > n13.H
	Window3d <008-Matrix.H n4=1 f4=14 n6=1 f6=14 > n14.H
	Window3d <008-Matrix.H n4=1 f4=15 n6=1 f6=15 > n15.H
	Window3d <008-Matrix.H n4=1 f4=16 n6=1 f6=16 > n16.H
	Window3d <008-Matrix.H n4=1 f4=17 n6=1 f6=17 > n17.H
	Window3d <008-Matrix.H n4=1 f4=18 n6=1 f6=18 > n18.H
	Window3d <008-Matrix.H n4=1 f4=19 n6=1 f6=19 > n19.H
	Window3d <008-Matrix.H n4=1 f4=20 n6=1 f6=20 > n20.H
	Window3d <008-Matrix.H n4=1 f4=21 n6=1 f6=21 > n21.H
	Window3d <008-Matrix.H n4=1 f4=22 n6=1 f6=22 > n22.H
	Window3d <008-Matrix.H n4=1 f4=23 n6=1 f6=23 > n23.H
	Window3d <008-Matrix.H n4=1 f4=24 n6=1 f6=24 > n24.H
	Window3d <008-Matrix.H n4=1 f4=25 n6=1 f6=25 > n25.H
	Window3d <008-Matrix.H n4=1 f4=26 n6=1 f6=26 > n26.H
	Window3d <008-Matrix.H n4=1 f4=27 n6=1 f6=27 > n27.H
	Window3d <008-Matrix.H n4=1 f4=28 n6=1 f6=28 > n28.H
	Window3d <008-Matrix.H n4=1 f4=29 n6=1 f6=29 > n29.H
	Window3d <008-Matrix.H n4=1 f4=30 n6=1 f6=30 > n30.H
	Window3d <008-Matrix.H n4=1 f4=31 n6=1 f6=31 > n31.H
	Window3d <008-Matrix.H n4=1 f4=32 n6=1 f6=32 > n32.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
	  n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      n30.H n31.H n32.H axis=3 > xdiag.H
	echo o3=-1200. >> xdiag.H

teste1:
	Attr <004-ILLANG1.H |awk '/rms/ {print $$3}' > lixo
	Attr <004-DIAGANG01.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG02.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG03.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG04.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG05.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG06.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG07.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG08.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG09.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG10.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG11.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG12.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG13.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG14.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG15.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG16.H |awk '/rms/ {print $$3}' >> lixo
	Attr <004-DIAGANG17.H |awk '/rms/ {print $$3}' >> lixo
	

teste2:
	Math file1=005-ADCIG.H file2=004-ILLANG1.H exp="file1*file2/(file2*file2+0.07)*0.735906*0.735906" > n0.H
	Math file1=005-ADCIG.H file2=004-DIAGANG01.H exp="file1*file2/(file2*file2+0.03)*0.371056*0.371056" > n1.H
	Math file1=005-ADCIG.H file2=004-DIAGANG02.H exp="file1*file2/(file2*file2+0.01)*0.12876*0.12876" > n2.H
	Math file1=005-ADCIG.H file2=004-DIAGANG03.H exp="file1*file2/(file2*file2+0.008)*0.0822569*0.0822569" > n3.H
	Math file1=005-ADCIG.H file2=004-DIAGANG04.H exp="file1*file2/(file2*file2+0.005)*0.0527149*0.0527149" > n4.H
	Math file1=005-ADCIG.H file2=004-DIAGANG05.H exp="file1*file2/(file2*file2+0.003)*0.0356188*0.0356188" > n5.H
	Math file1=005-ADCIG.H file2=004-DIAGANG06.H exp="file1*file2/(file2*file2+0.003)*0.0303835*0.0303835" > n6.H
	Math file1=005-ADCIG.H file2=004-DIAGANG07.H exp="file1*file2/(file2*file2+0.002)*0.0256147*0.0256147" > n7.H
	Math file1=005-ADCIG.H file2=004-DIAGANG08.H exp="file1*file2/(file2*file2+0.002)*0.0237133*0.0237133" > n8.H
	Math file1=005-ADCIG.H file2=004-DIAGANG09.H exp="file1*file2/(file2*file2+0.002)*0.0199597*0.0199597" > n9.H
	Math file1=005-ADCIG.H file2=004-DIAGANG10.H exp="file1*file2/(file2*file2+0.002)*0.0188286*0.0188286" > n10.H
	Math file1=005-ADCIG.H file2=004-DIAGANG11.H exp="file1*file2/(file2*file2+0.0015)*0.0155852*0.0155852" > n11.H
	Math file1=005-ADCIG.H file2=004-DIAGANG12.H exp="file1*file2/(file2*file2+0.0017)*0.0170111*0.0170111" > n12.H
	Math file1=005-ADCIG.H file2=004-DIAGANG13.H exp="file1*file2/(file2*file2+0.001)*0.0124157*0.0124157" > n13.H
	Math file1=005-ADCIG.H file2=004-DIAGANG14.H exp="file1*file2/(file2*file2+0.001)*0.012704*0.012704" > n14.H
	Math file1=005-ADCIG.H file2=004-DIAGANG15.H exp="file1*file2/(file2*file2+0.001)*0.0105007*0.0105007" > n15.H
	Math file1=005-ADCIG.H file2=004-DIAGANG16.H exp="file1*file2/(file2*file2+0.001)*0.010937*0.010937" > n16.H
	Math file1=005-ADCIG.H file2=004-DIAGANG17.H exp="file1*file2/(file2*file2+0.0008)*0.00868882*0.00868882" > n17.H
	Math file1=n0.H file2=n1.H file3=n2.H file4=n3.H file5=n4.H file6=n5.H file7=n6.H file8=n7.H file9=n8.H file10=n9.H file11=n10.H file12=n11.H file13=n12.H \
	file14=n13.H file15=n14.H file16=n15.H file17=n16.H file18=n17.H exp="file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15+file16+file17+file18" \
	> 011-Arr.H
	Stack <001-Arr.H |Window min1=10000|Grey|Tube&

teste3:
	Math file1=005-ADCIG.H file2=004-ILLANG1.H exp="file1*file2/(file2*file2+0.07)*0.735906" > n0.H
	Math file1=005-ADCIG.H file2=004-DIAGANG01.H exp="file1*file2/(file2*file2+0.03)*0.371056" > n1.H
	Math file1=005-ADCIG.H file2=004-DIAGANG02.H exp="file1*file2/(file2*file2+0.01)*0.12876" > n2.H
	Math file1=005-ADCIG.H file2=004-DIAGANG03.H exp="file1*file2/(file2*file2+0.008)*0.0822569" > n3.H
	Math file1=005-ADCIG.H file2=004-DIAGANG04.H exp="file1*file2/(file2*file2+0.005)*0.0527149" > n4.H
	Math file1=005-ADCIG.H file2=004-DIAGANG05.H exp="file1*file2/(file2*file2+0.003)*0.0356188" > n5.H
	Math file1=005-ADCIG.H file2=004-DIAGANG06.H exp="file1*file2/(file2*file2+0.003)*0.0303835" > n6.H
	Math file1=005-ADCIG.H file2=004-DIAGANG07.H exp="file1*file2/(file2*file2+0.002)*0.0256147" > n7.H
	Math file1=005-ADCIG.H file2=004-DIAGANG08.H exp="file1*file2/(file2*file2+0.002)*0.0237133" > n8.H
	Math file1=005-ADCIG.H file2=004-DIAGANG09.H exp="file1*file2/(file2*file2+0.002)*0.0199597" > n9.H
	Math file1=005-ADCIG.H file2=004-DIAGANG10.H exp="file1*file2/(file2*file2+0.002)*0.0188286" > n10.H
	Math file1=005-ADCIG.H file2=004-DIAGANG11.H exp="file1*file2/(file2*file2+0.0015)*0.0155852" > n11.H
	Math file1=005-ADCIG.H file2=004-DIAGANG12.H exp="file1*file2/(file2*file2+0.0017)*0.0170111" > n12.H
	Math file1=005-ADCIG.H file2=004-DIAGANG13.H exp="file1*file2/(file2*file2+0.001)*0.0124157" > n13.H
	Math file1=005-ADCIG.H file2=004-DIAGANG14.H exp="file1*file2/(file2*file2+0.001)*0.012704" > n14.H
	Math file1=005-ADCIG.H file2=004-DIAGANG15.H exp="file1*file2/(file2*file2+0.001)*0.0105007" > n15.H
	Math file1=005-ADCIG.H file2=004-DIAGANG16.H exp="file1*file2/(file2*file2+0.001)*0.010937" > n16.H
	Math file1=005-ADCIG.H file2=004-DIAGANG17.H exp="file1*file2/(file2*file2+0.0008)*0.00868882" > n17.H
	Math file1=n0.H file2=n1.H file3=n2.H file4=n3.H file5=n4.H file6=n5.H file7=n6.H file8=n7.H file9=n8.H file10=n9.H file11=n10.H file12=n11.H file13=n12.H file14=n13.H file15=n14.H file16=n15.H file17=n16.H file18=n17.H exp="file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15+file16+file17+file18" > 011-Arr.H

teste4:
	Math file1=005-ADCIG.H file2=004-ILLANG1.H exp="file1*file2/(file2*file2+0.007)" > n0.H
	Math file1=005-ADCIG.H file2=004-DIAGANG01.H exp="file1*file2/(file2*file2+0.003)" > n1.H
	Math file1=005-ADCIG.H file2=004-DIAGANG02.H exp="file1*file2/(file2*file2+0.001)" > n2.H
	Math file1=005-ADCIG.H file2=004-DIAGANG03.H exp="file1*file2/(file2*file2+0.0008)" > n3.H
	Math file1=005-ADCIG.H file2=004-DIAGANG04.H exp="file1*file2/(file2*file2+0.0005)" > n4.H
	Math file1=005-ADCIG.H file2=004-DIAGANG05.H exp="file1*file2/(file2*file2+0.0003)" > n5.H
	Math file1=005-ADCIG.H file2=004-DIAGANG06.H exp="file1*file2/(file2*file2+0.0003)" > n6.H
	Math file1=005-ADCIG.H file2=004-DIAGANG07.H exp="file1*file2/(file2*file2+0.0002)" > n7.H
	Math file1=005-ADCIG.H file2=004-DIAGANG08.H exp="file1*file2/(file2*file2+0.0002)" > n8.H
	Math file1=005-ADCIG.H file2=004-DIAGANG09.H exp="file1*file2/(file2*file2+0.0002)" > n9.H
	Math file1=005-ADCIG.H file2=004-DIAGANG10.H exp="file1*file2/(file2*file2+0.0002)" > n10.H
	Math file1=005-ADCIG.H file2=004-DIAGANG11.H exp="file1*file2/(file2*file2+0.00015)" > n11.H
	Math file1=005-ADCIG.H file2=004-DIAGANG12.H exp="file1*file2/(file2*file2+0.00017)" > n12.H
	Math file1=005-ADCIG.H file2=004-DIAGANG13.H exp="file1*file2/(file2*file2+0.0001)" > n13.H
	Math file1=005-ADCIG.H file2=004-DIAGANG14.H exp="file1*file2/(file2*file2+0.0001)" > n14.H
	Math file1=005-ADCIG.H file2=004-DIAGANG15.H exp="file1*file2/(file2*file2+0.0001)" > n15.H
	Math file1=005-ADCIG.H file2=004-DIAGANG16.H exp="file1*file2/(file2*file2+0.0001)" > n16.H
	Math file1=005-ADCIG.H file2=004-DIAGANG17.H exp="file1*file2/(file2*file2+0.00008)" > n17.H
	Math file1=n0.H file2=n1.H file3=n2.H file4=n3.H file5=n4.H file6=n5.H file7=n6.H file8=n7.H file9=n8.H file10=n9.H file11=n10.H file12=n11.H file13=n12.H file14=n13.H file15=n14.H file16=n15.H file17=n16.H file18=n17.H exp="file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15+file16+file17+file18" > 011-Arr.H
	Stack <011-Arr.H |Window min1=10000|Grey|Tube&

teste5:
	Math file1=004-ILLANG1.H file2=004-DIAGANG01.H file3=004-DIAGANG02.H file4=004-DIAGANG03.H file5=004-DIAGANG04.H file6=004-DIAGANG05.H file7=004-DIAGANG06.H file8=004-DIAGANG07.H file9=004-DIAGANG08.H file10=004-DIAGANG09.H file11=004-DIAGANG10.H file12=004-DIAGANG11.H file13=004-DIAGANG12.H file14=004-DIAGANG13.H file15=004-DIAGANG14.H file16=004-DIAGANG15.H file17=004-DIAGANG16.H file18=004-DIAGANG17.H exp="file1+file2+file3+file4+file5+file6+file7+file8+file9+file10+file11+file12+file13+file14+file15+file16+file17+file18" > 004-SumDiags.H
	Math file1=005-ADCIG.H file2=004-SumDiags.H exp="file1*file2/(file2*file2+0.007)" > 011-Arr.H
	Stack <011-Arr.H |Window min1=10000|Grey|Tube&
######################################################
#
# Inverse problem
#

008-Matrix.H:
	./Bin/built_hessian.x < 007-HESSIANG.H n1=322 o1=7975 d1=25 n2=120 o2=-60 d2=1 n3=900 o3=32000 d3=75 n4=18 d4=1 o4=-17 > $@

008a-Matrix.H:
	./Bin/build_hessian.x < 007-HESSIANG.H n1=322 o1=7975 d1=25 n2=120 o2=-60 d2=1 n3=900 o3=32000 d3=75 n4=18 d4=1 o4=-17 > $@

005-InADCIG.H:
	Transp reshape=1,3,5 plane=12 <005-ADCIG.H | Transp plane=24 > $@

009-Inversion1.H:
#	< 008-Matrix.H ~valencia/research/NoStaDecon_save/Bin/LINUX86_64/test_symmetry.x m=pp3.H >$@
	< 008-Matrix.H /homes/sep/alejandro/research/NoStaDecon_save/Bin/LINUX86_64/test_symmetry.x m=005-InADCIG.H >$@

# to run in gomorrah
009-Inversion.H:
	< 008a-Matrix.H ~valencia/research/NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=005-InADCIG.H niter=20 eps=.1 \ 
	alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov.H wmov=wmov.H > $@

010-convol.H:
	Spike n1=1 n2=1 n3=322 n4=50 n5=1 o1=0 o2=0 o3=7975 o4=32000 o5=1 d1=1 d2=1 d3=25 d4=75 d5=1 k1=1 k3=160 k4=30 mag=10 nsp=1 > zspk.H
	Window3d < 008a-Matrix.H n4=1 f4=17 n6=1 f6=61 squeeze=n > zdiag00_00deg.H 
	< zdiag00_00deg.H ~valencia/research/NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=zspk.H >$@
#	<008a-Matrix.H ~valencia/research/NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=zspk.H >$@

diag1.H:
	Transp plane=34 <ADCIG_orig.H |Transp plane=23 >junk1.H
	Transp plane=34 <ADCIG_mult.H |Transp plane=23 >junk2.H
#	./Bin/Hessang.x <junk1.H adcig2=junk2.H niter=20 eps=10. > $@

diagtest.H:
	./Bin/Hessang.x adcig1=test.H adcig2=test10.H n1=1500 o1=0 d1=0.008 n2=695 o2=-26025 d2=75 n3=1 o3=25925 d3=150 niter=500 eps=.00001 > $@

diagtest1.H:
	./Bin/Hessang1.x adcig1=test.H adcig2=test10.H n1=1500 o1=0 d1=0.008 n2=695 o2=-26025 d2=75 n3=1 o3=25925 d3=150 niter=5000 eps=0.001> $@

diagtest2.H:
#	Math file1=005-ADCIG_diag.H file2=005-ADCIG.H exp="file1/file2" > junk0.H
	Window min1=10000 max1=16000 min2=-15 n3=1 <005-ADCIG.H >jtest1.H
	Window min1=10000 max1=16000 min2=-15 n3=1 <005-ADCIG_diag.H >jtest2.H
	./Bin/Hessang1.x adcig1=jtest1.H adcig2=jtest2.H n1=241 o1=10000 d1=25 n2=75 o2=-15 d2=1 n3=1 o3=32000 d3=75 niter=60 eps=0.001> $@
#	./Bin/Hessang1.x adcig1=test.H adcig2=test10.H n1=1500 o1=0 d1=0.008 n2=695 o2=-26025 d2=75 n3=1 o3=25925 d3=150 niter=20 eps=0.001> $@

diag.H:
	./Bin/Hessang.x adcig1=j1.H adcig2=j2.H n1=251 o1=4 d1=0.008 n2=267 o2=-9975 d2=75 n3=1 o3=25925 d3=150 niter=2000 eps=0.0 > $@

clean:
	rm -f *.o *.fix.* *.mod *.x .make.dependencies.LINUX86_64

include ${SEPINC}/SEP.bottom