MAKE_DEPEND=yes
include ${SEPINC}/SEP.top
include FIGMakefile1
include MigMakefile
include SeminarMake

LATOPTS  = style=html
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
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 Istk

BIN=./Bin
PAR=./Par
RESDIR=./Fig/New
FIG=./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-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 > $@

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

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

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 |Window min1=10000 min2=-15 > 007-HESSIANG.H
	< 005-ADCIG.H Window min1=10000 min2=-15 > 007-GADCIG.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

######################################################
#
# Inverse problem
#
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