source("/home/lees/Progs/R_stuff/SUBS.R")
source("/home/lees/Progs/R_stuff/Part.Mo.R")

source("/home/lees/Progs/R_stuff/XCOR.R")
source("/home/lees/Progs/R_stuff/p3d.R")

source("/home/lees/Progs/R_stuff/CAZI.R")
source("/home/lees/Progs/R_stuff/interp.R")
source("/home/lees/Progs/R_stuff/iscale.R")
source("/home/lees/Progs/R_stuff/net.R")
source("/home/lees/Progs/R_stuff/p3d.R")
source("/home/lees/Progs/R_stuff/Part.Mo.R")
source("/home/lees/Progs/R_stuff/radiation.R")
source("/home/lees/Progs/R_stuff/READ3.R")
source("/home/lees/Progs/R_stuff/SUBS.R")
source("/home/lees/Progs/R_stuff/XCOR.R")


cmp = setCOSOmap()



library(mva)
library(cluster)
library(ts)



wazi(K, 100, 20, 15)
plotK3(K)

locator(1)

open two windows dev (2,3)



run 

w = wazi(K, 100, 20, 15)

y = locator(2)
mean(w$aaz[w$aex<y$x[2] & w$aex>y$x[1]])

par(mfrow=c(3, 1))
k = length(fs$f)/2
for(i in 1:(k))
{
   j = 2*(i)-1;
   print(c(j, fs$f[j], fs$f[j+1]))

K = read.DUMP( fs$f[j+1] ,  fs$f[j])
plotK3(K)
locator(1)




}

coso.sta$rot_rep(0,length(coso.sta$lon))

ANGS = rep(0,k)
X11()

for(i in 1:(k))
{
   j = 2*(i)-1;

K = read.DUMP( fs$f[j+1] ,  fs$f[j])
w = wazi(K, 100, 20, 15)
y = locator(2)
ANGS[i] = mean(w$aaz[w$aex<y$x[2] & w$aex>y$x[1]])
VARS[i] = var(w$aaz[w$aex<y$x[2] & w$aex>y$x[1]])

}

---------------------------------------

fs.etu = scan(file="etu.files", what=list(f=""))

ketu  = length(fs$f)/2
ANGS.etu = rep(0,ketu)
VARS.etu= rep(0,ketu)
#  X11()

for(i in 1:(ketu))
{
   j = 2*(i)-1;

K = read.DUMP( fs.etu$f[j+1] ,  fs.etu$f[j])
w = wazi(K, 100, 20, 15)
y = locator(2)
ANGS.etu[i] = mean(w$aaz[w$aex<y$x[2] & w$aex>y$x[1]])
VARS.etu[i] = var(w$aaz[w$aex<y$x[2] & w$aex>y$x[1]])

}

--------------------------------------
####  see /home/lees/Progs/R_stuff/XCOR.R
###  source("/home/lees/Progs/R_stuff/XCOR.R")

lim1 =  getlim(i4, s4, dt)
xchug(i4, s4, lim1, dt, tskip)

xchug1(i4, s4, lim1, dt, tskip)


D3_read.DUMP('SEIS_INFO_0003.dat', 'SEIS_DATA_0003.dat')
s3 = D3$data[,2]
i3 = D3$data[,1]
xchug1(i3, s3, lim1, dt, tskip)
############
D5_read.DUMP('SEIS_INFO_0005.dat', 'SEIS_DATA_0005.dat')
s5 = D5$data[,2]
i5 = D5$data[,1]

xchug1(i5, s5, lim1, dt, tskip)
############
############
D6_read.DUMP('SEIS_INFO_0006.dat', 'SEIS_DATA_0006.dat')
s6 = D6$data[,2]
i6 = D6$data[,1]

xchug1(i6, s6, lim1, dt, tskip)

############
D7_read.DUMP('SEIS_INFO_0007.dat', 'SEIS_DATA_0007.dat')
s7 = D7$data[,2]
i7 = D7$data[,1]

xchug1(i7, s7, lim1, dt, tskip)
############
D8_read.DUMP('SEIS_INFO_0008.dat', 'SEIS_DATA_0008.dat')
s8 = D8$data[,2]
i8 = D8$data[,1]

xchug1(i8, s8, lim1, dt, tskip)

############
D9_read.DUMP('SEIS_INFO_0009.dat', 'SEIS_DATA_0009.dat')
s9 = D9$data[,2]
i9 = D9$data[,1]

xchug1(i9, s9, lim1, dt, tskip)

############
D10_read.DUMP('SEIS_INFO_0010.dat', 'SEIS_DATA_0010.dat')
s10 = D10$data[,2]
i10 = D10$data[,1]

xchug1(i10, s10, lim1, dt, tskip)





xcor_function(a1, a2, dt)
{
b1 = ts(a1, deltat=dt)
b2 = ts(a2, deltat=dt)
ccf(b1, b2, lag.max = round(length(a1)/4) , type ="correlation",  plot = TRUE )
}

dlt = 0.008


pls1_function(a1, dt)
{
par(mfrow=c(1,1))
n1 = length(a1)
plot(dt*(1:n1), a1, type = 'l')
}


plsi_function(a1, dt)
{
par(mfrow=c(1,1))
n1 = length(a1)
plot((1:n1), a1, type = 'l')
}

plsi(s4,dlt)

plsi(i4,dlt)

lim = round(xlim())

plot(i4[lim[1]:lim[2]], type = 'l')

plot(s4[lim[1]:lim[2]], type = 'l')

xcor2(i4[lim[1]:lim[2]], s4[lim[1]:lim[2]], dlt)


xcor3(i4[lim[1]:lim[2]], s4, dlt)

plsi(i4[lim[1]:lim[2]],dlt)



par(mfrow = c(1,1))

par(mfrow = c(2,1))

xcs = xcrude(chug, s4, 0.008, 0.1)

xci = xcrude(chug, i4, 0.008, 0.1)

max(xci[!is.na(xci)])

xtims[xci[!is.na(xci)]==max(xci[!is.na(xci)])]

xtims[xcs[!is.na(xcs)]==max(xcs[!is.na(xcs)])]

xtims = xt[!is.na(xci)]

xt[!is.na(xcs)]


Fxci_xci!='NA'
xt[Fxci]




plot(xt, xcs, type = 'p')

plot(xt,xci, type = 'p')



T2pf = scan(file="/home/beer/lees/Coso/Scat_Data/S4T2.pfiles", what="")
T2Wf_paste(substring(T2pf, 1, nchar(T2pf)-1), "W", sep="")


T1pf = scan(file="/home/beer/lees/Coso/Scat_Data/S4T1.pfiles", what="")
T1Wf_paste(substring(T1pf, 1, nchar(T1pf)-1), "W", sep="")




xlim_function()
{
y = locator(2)
return(y$x)
}




ccf(i4, s4, lag.max = NULL, type ="correlation",  plot = TRUE )



.First <- function() 
{
cat("\n   Welcome to R!\n\n")
source("/home/lees/Progs/R_stuff/SUBS.R")
source("/home/lees/Progs/R_stuff/Part.Mo.R")
source("/home/lees/Progs/R_stuff/p3d.R")
source("/home/lees/Progs/R_stuff/CAZI.R")

library(mva)
library(cluster)
library(ts)


}
#############################

graphics.off()

#  source("/home/lees/Progs/R_stuff/CAZI.R")


# labs = c("Quit", "Next", "P", "T1", "S", "T2")
TESTBAR()

#  t1out = doGetAZ(T1pf)
#  t2out = doGetAZ(T2pf)

 


#  Z = doCAZI(j)   LabelBAR(u, 5, labs)  print(ValBAR(u, 5,locator(1)) )

#   doGetAZ() 

#     

getUW_AM.prl 19*p > AM.cards


source("/home/lees/Progs/R_stuff/CAZI.R")

source("/home/lees/Progs/R_stuff/SUBS.R")

eqs = setEQSLOCS("EQ.LOCS")
stas = setEQSLOCS("/home/lees/Site/Coso/coso_sta")


stas = setEQSLOCS("/home/lees/Progs/snaps/Tables/ActiveStas/coso_station.LLZ")



cmp = setCOSOmap()

COSOmap(cmp, EQS=eqs, STA=stas)

doGetAZ(T2pf)

cmdM = paste(" ", "grep ^M",  "199410100803p" , "> MCARD")

M = system(cmdM )

system("cat MCARD")

M = scan(file="MCARD", list(id="", F="", az1=0, d1=0, G="", az2=0, d2=0, 
U="", uaz=0, ud=0, V="", vaz=0, vd=0, P="", paz=0, pd =0, T="", taz=0, td=0))

M$az1

-----------------------

cmdM = paste(" ", "",  "199410100803p" , "> AMCARD")







doCAZI(9, T2pf)


sta.org=list(lat=stas$lat[stas$name=="S4"], lon=stas$lon[stas$name=="S4"])

S = gclc(sta.org$lat, sta.org$lon, stas$lat , stas$lon)

S = gclc( sta.org$lon, sta.org$lat, stas$lon,  stas$lat )



plot(stas$lon, stas$lat, type='n')
text(stas$lon, stas$lat,labels=stas$name)

dev.set(dev.next())
plot(S$x, S$y,type='n' )
text( S$x, S$y ,labels=stas$name)

u = par("usr")

A = gclc( sta.org$lon, sta.org$lat,  c(u[1],u[2]) , c(u[3], u[4])  )


Rasp = (A$y[2]-A$y[1])/(A$x[2]-A$x[1])
Rasp = (A$x[2]-A$x[1])/(A$y[2]-A$y[1])

plot(stas$lon, stas$lat, asp=Rasp)

text(stas$lon, stas$lat,labels=stas$name)


help(plot.window)


for( i in 1:length(Mpf)) {
RAK = see1M(Mpf[i])
locator(1)
}
