APRIL.FIG.R:april.analA<-function(APALL)
APRIL.FIG.R:april.fig<-function(AP)
APRIL.FIG.R:PLOT.MATAP<-function(ascd, tim=1, dt=1,  WIN=WIN, labs=LABS, notes=notes, sfact=1, LOG="", COL=col)
APRIL.FIG.R:dump.ALOOK<-function(mSW, STA="", FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, sfact=1 )
APRIL.FIG.R:dump.QLOOK<-function(mSW, STA="", FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, sfact=1 )
bosai.R:setJAPmap<-function()
bosai.R:japmap<-function(add=FALSE, STA=FALSE)
bosai.R:japmape<-function(add=FALSE, STA=FALSE)
bosai.R:xeq<-function(x, y, z, wd=10)
bosai.R:jcont<-function(x,y,z, d=0.7)
bosai.R:jimag<-function(x,y,z, d=0.7, col=terrain.colors(50))
bosai.R:jmap<-function()
bosai.R:jsta<-function()
bosai.R:SEESTA<-function(s)
bosai.R:get.jap.res<-function(file)
bosai.R:get.lq.res<-function(infile="residuals.P")
bosai.R:DO.jap.res<-function(files, dir='/home/lees/Site/Fuji/UW/ATEST')
bosai.R:map.jap1.res<-function(file, cmp=rainbow(100) )
bosai.R:plot.jap.res<-function(res, kind=1)
bosai.R:prep.ROSE<-function(infile)
bosai.R:##   warning:  this function sets global variables
bosai.R:STA.ROSE<-function()
bosai.R:STA1.ROSE<-function(i)
bosai.R:STA.ALL.ROSE<-function()
bosai.R:MAP.ROSE<-function()
bosai.R:MAPALL.ROSE<-function(WIN=win, ROSES=MROSE, rscale=1)
bosai.R:pageohist<-function(res)
bosai.R:histores<-function(res, disp=c(5,2), xlim=xlim)
bosai.R:acstats<-function(aclq, flag=TRUE, pch=pch, col=col)
bosai.R:getac<-function(fn)
bosai.R:compac<-function(aclq, acbos, lim=rep(TRUE, length(aclq$id)), zoom=list(x=c(0,1), y=c(0,1))  )
bosai.R:compac1<-function(aclq, acbos, lim=rep(TRUE, length(aclq$id)), zoom=list(x=c(0,1), y=c(0,1))  )
bosai.R:compacf<-function(aclq1, aclq0)
bosai.R:fix.delays<-function(pix="FQ1.pickfiles", fin="fuj1.del")
BUfilter.R:butfiltA<-function(a, ord=ord, fl=fl, fh=fh, deltat=dt, type=type, proto=pro)
BUfilter.R:butfilt<-function(a, fl, fh, deltat, type, proto)
BUfilter.R:FILT.spread<-function(x,y, dt, fl=fl, fh=fh, sfact=1, swin=NULL)
BUfilter.R:testbutter<-function(x,y, dt)
CAZI.R:plotBAR<-function(u, n)
CAZI.R:ValBAR<-function(u, n, v)
CAZI.R:LabelBAR<-function(u, n, labs)
CAZI.R:SortAZ<-function(FLS)
CAZI.R:doGetAZ<-function(FLS)
CAZI.R:AutoGetAmps<-function(FLS)
CAZI.R:AutoGetAZ<-function(FLS, ask=FALSE, PS=TRUE)
CAZI.R:Rsection<-function(LPF, COMP=1)
CAZI.R:getUWwin<-function(PF)
CAZI.R:slider.rot<-function(a1)
CAZI.R:img <-tkrplot(tt, function() srot(ascd, bb))
CAZI.R:f<-function(...) {
CAZI.R:doCAZI<-function(i, PF)
CAZI.R:setCOSOmap<-function(FN, eqs=NULL)
CAZI.R:COSOmap<-function(cmp, EQS=NULL, STA=NULL, LAB=NULL)
CAZI.R:CAZI<-function(ain, len, shift, p )
CAZI.R:TESTBAR<-function()
circ.R:rose.jml <- function(x, bins, main = "", prop = 1, pts = FALSE,
circ.R:Krose.jml <- function(x, wts=1, bins=36,rscale=1,  main = "", prop = 1, pts = FALSE,
circ.R:prose<-function(kros, ex=0, why=0, main = "", add=FALSE, prop = 1, perim=TRUE, style=1,
circ.R:SPLrose.jml <- function(x, wts=1, tags=NULL, bins=36,rscale=1,  main = "", prop = 1, pts = FALSE,
cluster.R:do.tree<-function(dx, method)
cluster.R:comp.tree4<-function(dx)
cluster.R:comp.tree<-function(dx, meth=c("single","average","median","centroid")   )
cluster.R:see.tree<-function(dend1)
cluster.R:plot.silhouette<-function (x, nmax.lab = 50, max.strlen = 5, main = NULL, sub = NULL,
Clust.R:clust.meth<-function(dis=d)
Clust.R:clust.k<-function(dis=d)
Default.R:plot.default<-function (x, y = NULL, type = "l", xlim = NULL, ylim = NULL,
DUMPSEIS.R:dump.get<-function(infofile, datfile, DIR=".", kname=1, kcomp=2)
DUMPSEIS.R:dump.plot<-function(S1)
DUMPSEIS.R:dump.Aplot<-function(S1, COMP=c(1:4), WIN=c(0,1), RMINST=FALSE, INTEG=FALSE, FILT=f, rot=FALSE, CINTEG=c(FALSE), LABS=LABS, sfact=1)
DUMPSEIS.R:get.dump.sta<-function(SIN, stpos, comppos, STA=stastromb)
FUNCTIONS.R:############  special functions for  R (splus)
FUNCTIONS.R:erf<-function(x)
General.R:ans<-function()
General.R:dnext<-function()
General.R:par1<-function()
General.R:one<-function()
General.R:no.dimnames<-function(a)
General.R:color.pr<-function(color="white")
General.R:plot.default<-function(x, y = NULL, type = "p", xlim = NULL, ylim = NULL,
General.R:letter.it<-function(a, corn)
General.R:label.it<-function(a, corn)
General.R:addtix<-function(side=3, pos=0,   tck=0.005, at=ttics, labels=FALSE, col=2, ...)
General.R:itoxyz<-function(i, nx, ny, nlay)
General.R:xyztoi<-function(ix, iy, iz, nx, ny, nlay)
General.R:jgrid<-function(x,y, col=col)
General.R:Fresh.Screens<-function(n)
General.R:Get.Screens<-function(n)
HEAT.EQN.R:############  functions for performing cross correlations in R (splus)
HEAT.EQN.R:erf<-function(x)
HEAT.EQN.R:heat.sol<-function(x, T0, k, t)
HEAT.EQN.R:DO.HALFSPACE<-function()
HEAT.EQN.R:DO.DYKE<-function(a=a, x=x, t=t, k=k, T0=T0, NDIM=TRUE)
INSTR.R:ReadSet.Instr<-function(file)
INSTR.R:PreSet.Instr<-function()
INSTR.R:poly<-function(x)
INSTR.R:polyval<-function(p,x)
INSTR.R:makefreq<-function(n,dt)
INSTR.R:FREQS<-function(b,a,w)
INSTR.R:response<-function(Kal,key,ff,tt=tt,plotkey=FALSE)
INSTR.R:  d=rep(0, length=n)                  ## % d is a delta function at time=0
INSTR.R:  DRET=ft(d,n,tstart,sintr);      ##% fourier transform of delta function
INSTR.R:        title(main='transfer function', xlab="frequency (Hz)", ylab='amplitude' )
INSTR.R:        title(main='transfer function', xlab="frequency (Hz)", ylab='phase' )
INSTR.R:        plot(RESP$t, resp, type='l', main='impulse response function', xlab="time (s)", ylab='amplitude')
INSTR.R:###   response      compute impulse response and transfer function
INSTR.R:### Convert poles and zeros of transfer function to polynomial coefficients, 
INSTR.R:### Calculate transfer function only if length(ff)>0
INSTR.R:### ff          = evaluate transfer function at the frequencies (Hz) in ff
INSTR.R:### transfer    = transfer function evaluated at frequencies (f)
INSTR.R:### f           = frequency vector corresponding to transfer function (Hz)
INSTR.R:ft<-function(g,n,tstart,dt)
INSTR.R:ift =function(G,n,tstart,dt)
INSTR.R:  ### function [g,f,t]=ift(G,n,tstart,dt);
INSTR.R:logspace<-function(d1, d2, n=n)
INSTR.R:LocalUnwrap<-function(p,cutoff=cutoff)
INSTR.R:deconinst<-function(data, sintr, KAL,key, Calibnew, waterlevel=1.e-8)
INSTR.R:finteg<-function(data, dt)
interp.R:kamplot<-function()
interp.R:ccw<-function(p0, p1, p2)
interp.R:intersect<-function( l1,  l2)
interp.R:inside<-function(A, POK )
interp.R:inpoly<-function(x,y, POK )
interp.R:jiden<-function(x, y, a, tol)
interp.R:jperim<-function( x, y, tol)
interp.R:jdist<-function(x, y)
interp.R:seeit<-function(j)
interp.R:quads<-function(x,y, mag)
interp.R:jfind<-function(KPT)
interp.R:jident<-function(KPT)
interp.R:jscal<-function(a, k)
interp.R:vscal<-function()
Inv.R:## these are functions for reading in data from a ray file  
Inv.R:get.ray.res<-function(file)
Inv.R:set.res<-function(a, sta)
Inv.R:jscale<-function(a, j)
Inv.R:plot.staparm<-function(lat, lon, nm, ista, g)
Inv.R:plot.stares<-function(res,  stas)
Inv.R:plot.res<-function(a)
Inv.R:get.stas<-function(file)
Inv.R:show.stas<-function(a)
Inv.R:tomo.col<-function(n)
iscale.R:image.SCALE<-function (z, col, x, y = NULL, size = NULL, digits = 2, labels = c("breaks",
iscale.R:HOZscale<-function(z, col, units="")
iscale.R:see.pal<-function(col)
KAM.R:readkammap<-function()
KAM.R:readkamSKS<-function()
KAM.R:pkammap<-function(VMAP)
KAM.R:plotSKS<-function()
KAM.R:plotSKSproj<-function( dfac=0.1)
kar.R:KAR.get<-function(A, DIR=".", j)
kar.R:KAR.it<-function(A, DIR=".", ZOOM=FALSE, PICK=1, Pause=0)
kar.R:vlocator<-function(COL=1, NUM=FALSE)
kar.R:plocator<-function(COL=1, NUM=FALSE, Yranges=NULL, style=0)
kar.R:KAR.view<-function(A, DIR=".", ZOOM=FALSE, PICK=1, Pause=0, XSTA=NULL)
kar.R:matchstas<-function(GFIL, STA, COMP)
kar.R:checkpick<-function(fin, dir="./")
kar.R:KAR.postview<-function(B, DIR=".", ZOOM=FALSE, PICK=1, Pause=0, XSTA=NULL, APLOT=FALSE, BPLOT=TRUE, COMP=4, s2nwin=2)
kar.R:Xpostview<-function(K1)
kar.R:KAR.picker<-function(jump=25, begin=1, flist=f98.250, dir='.', ZOOM=TRUE, pref="TES00", sta="1025" )
kar.R:KAR.calc<-function(B, DIR=".", rot=FALSE)
kar.R:UNPAK.calc<-function(LUV)
kar.R:KAR.Dchug<-function(A, B, CH=choose)
kar.R:KAR.chug<-function(A, B, DIR=".")
kar.R:KAR.look<-function(B, DIR=".")
kar.R:KAR.Getdat<-function(ifile, DIR)
kar.R:KAR.LIM<-function(FIG, WIN=c(0,1))
kar.R:KAR.4plot<-function(FIG, WIN=c(0,1), LABS=LABS)
kar.R:KAR.Gdat<-function(ifile, DIR, stloc=1, cmploc=2)
kar.R:KAR.Pdat<-function(GFIL)
kar.R:KAR.door<-function(B, DIR=".")
kar.R:KAR.Nplot<-function(FIG, COMP=c(1:4), WIN=c(0,1), LABS=LABS, COL=1, ZERO=FALSE)
kar.R:SETAIN<-function(A1, J)
kar.R:PMO.LOOK<-function(ain1, p)
kar.R:KAR.rose<-function(PM, let=1)
kar.R:KARAZI<-function(ain, len, shift, p )
kar.R:brune.fit<-function(amp, dt=1, f1=0.01, f2=15, nw=list(x=c(0,1)), PLOTA=FALSE, PLOTB=FALSE , tit="", PICK=FALSE)
kar.R:win.spec<-function(amp, dt, wins=list(x=c(0,1)) )
kar.R:KAR.MOM<-function(BFILES, DIR=".", VEC=VEC,  f1=0.01, f2=15, STA=c("1025"), COMP=4, PICK=FALSE, Pause=0 )
kar.R:SEE.MOM<-function(MOM, DIR=".", VEC=VEC)
LOC.R:## locator and other R-functions for interactive
LOC.R:vloc<-function()
MAP.R:MAPconstants<-function()
MAP.R:setEQSLOCS<-function(F)
MAP.R:setEQSLLZM<-function(FL, LLIM=NULL )
MAP.R:setMarkup<-function(LABS=NULL)
MAP.R:Markup<-function(MM=MM, cex=cex, ...)
MAP.R:MAPID<-function()
MAP.R:KMscale<-function(kmL=kmL, len=500, units="km")
MAP.R:VANE<-function(vane=vane)
MAP.R:setGenmap<-function(FN)
MAP.R:GENmap<-function(cmp, ADD=TRUE, ASP=TRUE, DIR=0)
MAP.R:GPROJm<-function(cmp, LAT=x, LON=y, pts=NULL, ADD=TRUE, WIN=NULL,
MAP.R:setPROJ<-function(type=1, LAT0=LAT, LON0=LON , LATS=NULL, LONS=NULL, DLAT=NULL, DLON=NULL)
MAP.R:GLOB.XY<-function(LAT, LON)
MAP.R:XY.GLOB<-function(x, y)
MAP.R:InRect<-function(x, y, X1, Y1, X2, Y2)
MAP.R:SELmap<-function(cmp,WIN=NULL)
MAP.R:PROJmap<-function(cmp, pts=NULL, ADD=TRUE, WIN=NULL, ASP=TRUE, NUM=FALSE,
MAP.R:AXESmap<-function(cmp, GRID=FALSE, WIN=NULL)
MAP.R:SUBmap<-function(cmp, tbox, pts=null)
MAP.R:merc.sphr.xy<-function(LON0, phi, lam)
MAP.R:merc.sphr.ll<-function(LON0, x, y)
MAP.R:utm.sphr.xy<-function(  phi,  lam)
MAP.R:utm.sphr.ll<-function( x, y)
MiscSeis.R:twostaZ<-function(dep=dep, sta1=list(d=40, r=120), sta2 = list(d=340,r=1060), vp = 3000.0, pts=pts)
MLAB.R:focmec<-function(az1, dip1, az2, dip2)
MLAB.R:dip.points<-function(az1, dip1, az2, dip2 )
MLAB.R:pole.points<-function(az1, dip1, az2, dip2 )
MLAB.R:TP.2.XYZ<-function(trend, dip)
MLAB.R:null.plane<-function(az1, dip1, az2, dip2)
MLAB.R:dip.point<-function(az1, dip1  )
MLAB.R:circle<- function()
MLAB.R:circ.tics<-function()
MLAB.R:fmtraj<-function(dip, az, iang)
myakejima.R:prep.cor<-function(MJALL2, MSTAS,MCOMPS )
myakejima.R:set.cor<-function(allcor, env=e1)
myakejima.R:seeV<-function(MYALL,STA="", N=4, fl=fl, fh=fh)
myakejima.R:seeMYSPEC<-function(MYALL,STA="", N=4, fl=fl, fh=fh)
myakejima.R:seeMYFILT<-function(MYALL,STA="", N=4, fl=fl, fh=fh)
myakejima.R:compfilt<-function(MYALL,STA="", comps)
myakejima.R:myak.Xtract<-function(mSW, STA="KAC", COMP="Vertical", BEF=1, AFT=5, FV=FV, PLT=FALSE )
NANDES.R:slabcont<-function(x,y,z, levs=25,  POL=POL)
NANDES.R:NANMAP <-function(PLOT=TRUE, WIN=GL)
NANDES.R:donan<-function(WIN=ECzoom)
NANDES.R:MAXCUT<-function(YX,YZ,thresh=5, dkm=100)
NANDES.R:NAMARKS<-function()
net.R:GreatDist<-function( LON1,  LAT1,  LON2,  LAT2)
net.R:AlongGreat<-function( LON1,  LAT1,  km1, ang )
net.R:focpoint<-function(az1, dip1, col=2, pch=5, lab="", UP=FALSE, PLOT=TRUE )
net.R:uppoint<-function(A)
net.R:upvec<-function(A, col=2)
net.R:fpoint<-function(az1, dip1, col=2, pch=5, lab="", UP=FALSE )
net.R:qpoint<-function(az, iang, col=2, pch=5, lab="", POS=4, UP=FALSE )
net.R:circ<-function()
net.R:HFOC<-function(H)
net.R:rectify.angle<-function(iang)
net.R:FixDip<-function(A)
net.R:DOT<-function(A1, A2)
net.R:LEN<-function(A1)
net.R:NORM<-function(A1)
net.R:ROT<-function(N, M)
net.R:CROSS<-function(A1, A2)
net.R:CROSS.DIP<-function(A1, A2)
net.R:FLIP<-function(A1)
net.R:XprodLL<-function(  deglon1, deglat1,   deglon2, deglat2  )
net.R:DANG<-function(A1,A2)
net.R:TOCART<-function(az, nadir)
net.R:TOSPHERE<-function(x, y, z)
net.R:tocart<-function(A)
net.R:tosphere<-function(A)
net.R:TOCART.DIP<-function(az, dip)
net.R:TOSPHERE.DIP<-function(x, y, z)
net.R:AXpoint<-function(UP=TRUE)
net.R:EApoint<-function()
net.R:pcirc<-function(gcol)
net.R:rotx<-function( deg )
net.R:REFLECT<-function(A)
net.R:Preflect<-function(az, dip )
net.R:CONVERTSDR<-function(strike, dip, rake)
net.R:net<-function(add)
net.R:getM<-function(FROM, TO)
net.R:prang<-function(a)
net.R:PRINTMC<-function(M)
net.R:MOD<-function(k, m)
net.R:fmod<-function(k, m)
net.R:faultplane<-function(az,  dip,  col = par("col"), PLOT=TRUE, UP=FALSE) 
net.R:POLYFOC<-function(H, PLOT=FALSE)
net.R:bang<-function(x1,y1,x2,y2)
net.R:polyfoc<-function(strike1, dip1, strike2, dip2, PLOT=FALSE)
net.R:lowplane<-function( az, dip, col = par("col"), UP=FALSE) 
net.R:GROTATE<-function()
net.R:JROTATE<-function(pfile)
net.R:MROTATE<-function(pfile)
net.R:AROTATE<-function(North, East, Down)
net.R:pplanes<-function(B1)
Nimage.R:imatestrgb<-function (mat3d,ngris=64,stretch=" ",ver=T)
Nimage.R:# apply(tstarr/256,c(1,2),function(z)do.call("rgb",as.list(z)))
Nimage.R:#used, but should be a fast function.
Nimage.R:rescale<-function(vector, oldmin = min(vector), oldmax = max(vector), newmin = 0,
Nimage.R:function(matriz, origen, imagen, directo = T)
Nimage.R:imagen<-function(x,col="bn",add=F)
p3d.R:see3comp_function(man.pick)
p3d.R:replot3<-function(rx)
p3d.R:replotSO<-function()
p3d.R:a3comp_function(ain)
Pal.R:tomo.col<-function(n)
Pal.R:shade.col<-function(n, acol=c(1,0,0) )
Pal.R:plot.pal<-function(col)
particle_motion.R:doHODO<-function(FLS, rot=FALSE, PS=TRUE)
particle_motion.R:autoHODO<-function(FLS, rot=FALSE, PS=TRUE)
particle_motion.R:testparm<-function(x, rot=FALSE, PS=FALSE)
particle_motion.R:HODO<-function(a1, rot=FALSE, PS=FALSE)
particle_motion.R:ptrace3<-function(ain, HODO=FALSE, rot=FALSE  )
particle_motion.R:plot3comp<-function(ain)
particle_motion.R:plot3co<-function(ain)
particle_motion.R:tracepick<-function( ascd, ex, labs=labs, ID="")
particle_motion.R:PMOtrace<-function(ascd, tim=1, labs=LAB, PS=FALSE, ID="")
particle_motion.R:ptrace<-function(ascd, tim=1, labs=LAB)
particle_motion.R:ptraceA<-function(ascd, tim=1, labs=LABS)
particle_motion.R:seltrace<-function(ascd)
particle_motion.R:DO.PMOT.ARR<-function(E, N)
particle_motion.R:pmotion<-function(ascd, X, Y , xlim=xx, ylim=yy, labs=labs, tit=TIT, PLPOINTS=FALSE, COL=0)
particle_motion.R:pmot3<-function(ascd, X, Y , xlim=xx, ylim=yy, labs=labs, tit=TIT)
particle_motion.R:ptraceB<-function(ascd, K, tim=1)
particle_motion.R:ptrace2<-function(ascd, K, tim=1)
particle_motion.R:slider.hodo<-function(abaz, ex)
particle_motion.R:img <-tkrplot(tt, hscale=1.5, vscale=1.5, function() ptrace2(abaz, bb, tim=ex))
particle_motion.R:f<-function(...) {
Part.Mo.R:wazi<-function(ain, len, shift, p )
Part.Mo.R:Cazi<-function(ain, len, shift, p )
Part.Mo.R:plot.t1t2<-function(ain)
Part.Mo.R:plot.ps<-function(ain)
Part.Mo.R:plot.medbars<-function(ain, wex, aaz)
Part.Mo.R:grotseis<-function(ang, flip=FALSE)
Part.Mo.R:read3comp<-function(len, shift, show.it, man.pick)
Part.Mo.R:POL.AZ<-function(Ofile)
radiation.R:printmc1<-function(M)
radiation.R:printmcrak<-function(RAK)
radiation.R:printmc0<-function(M)
radiation.R:seeoneM<-function(pf)
radiation.R:see1M<-function(pf)
radiation.R:Fmotion<-function(RAK)
radiation.R:seeMech<-function(pf, PLOT=TRUE)
radiation.R:seeALLM<-function(Mpf)
radiation.R:PPLOT<-function(pf)
radiation.R:get.ref.file<-function()
radiation.R:plotfoc<-function(RAK)
radiation.R:Simplefoc<-function(RAK)
radiation.R:Simplerfoc<-function(RAK)
radiation.R:SHFTNSCALE<-function(x,x1, minx, maxx, siz)
radiation.R:PLfoc<-function(x1, y1, siz,  GF, COL=gray(0.95)  )
radiation.R:Calcfoc<-function(RAK)
radiation.R:radP<-function( del, phiS, lam, ichi, phi)
radiation.R:radSV<-function( del, phiS, lam, ichi, phi)
radiation.R:radSH<-function( del, phiS, lam, ichi, phi)
radiation.R:Plot.Planes<-function(RAK)
radiation.R:Plot.UWPlanes<-function(M)
radiation.R:radiateP<-function(RAK, SCALE=FALSE, col=col)
radiation.R:radiateSV<-function(RAK, SCALE=FALSE, col=col)
radiation.R:radiateSH<-function(RAK, SCALE=FALSE, col=col)
radiation.R:imageP<-function(phiS, del , lam, SCALE=FALSE, UP=FALSE, col=NULL )
radiation.R:imageSV<-function( phiS, del, lam, SCALE=FALSE, UP=FALSE, col=NULL )
radiation.R:imageSH<-function( phiS, del, lam, SCALE=FALSE, UP=FALSE, col=NULL )
radiation.R:mc2cart<-function(az, dip)
radiation.R:GetRakeSense<-function(uaz, upl, vaz,  vpl , paz,  ppl,  taz, tpl)
radiation.R:MRake<-function(M)
radiation.R:SetMech<- function(RAK)
radiation.R:GetRake<- function( az1, dip1, az2, dip2,  dir)
radiation.R:plotfocPAPER<-function(RAK, SCALE=FALSE, col=rainbow(20) )
radiation.R:radiationRATIO<-function(RAK )
READ3.R:psplot_function(file)
READ3.R:psplot_function(file)
READ3.R:read3comp_function(len, shift, show.it, man.pick)
READ3.R:readsocomp_function()
READ3.R:liner_function()
READ3.R:   letter.it_function(a, corn)
READ3.R:FR_function(t)
READ3.R:gclc_function(phiorg, lamorg, phi, lam)
rotate.R:rotx <- function( deg )
rotate.R:roty <- function( deg )
rotate.R:rotz <- function( deg )
rotate.R:Motate<-function(x,y,z, a1, i)
rotate.R:Jotate<-function(x,y,z, a1, iseq=seq(from=0, to=360, by=10), sta=Dsta )
rotate.R:Xotate<-function(x,y,z, a1, iseq=seq(from=0, to=360, by=10), sta=Dsta )
scat3d.R:       cubedraw <- function(res3d, min = 0, max = 255, cex = 2)
scat3d.R:setHELmap<-function()
scat3d.R:helmap<-function(add=FALSE, STA=FALSE)
seis.R:getSEGY<-function(file)
seis.R:fillsegy<-function(SEGY.h)
seis.R:seisrot<-function(dat, sleeper=0.002)
seis.R:srot<-function(dat, i)
seis.R:pseis<-function(dat)
seis.R:pseis2<-function(ascd)
sig.R:lot<-function(x,y){
sig.R:jot<-function(x){
sig.R:zlot<-function(x, dt)
sig.R:next2<-function(x)
sig.R:which.outlier<-function(vec, m=1)
sig.R:RESCALE<-function(x, nx1, nx2, minx, maxx)
sig.R:addlines<-function(ex, y, OLDY, lty=2, AX=TRUE, LAB="", col=1)
sig.R:stats<-function(d)
sig.R:tufteboxplot<-function(istr, OR=1:length(istr), ALAB=FALSE, OL=FALSE)
sig.R:Get.S2N<-function(g , dt, t1, twin )
sig.R:ampsp<-function(x){
sig.R:Xspec<-function(x1, x2)
sig.R:GT<-function()
sig.R:envelope<-function(x)
sig.R:hilbert<-function(x)
sig.R:ra<-function(fn){
sig.R:MTMdisp<-function(a, f1=f1, f2=f2, len2=1024, PLOT=FALSE)
sig.R:spec<-function(a, f1=f1, f2=f2, len2=1024, PLOT=FALSE)
sig.R:specl<-function(a, f1=f1, f2=f2)
sig.R:add.brune0<-function(a,tstar,fc,Omega0, icol, f1=f1, f2=f2, plog=FALSE)
sig.R:add.brune<-function(a,tstar,fc,Omega0, gamma=2, icol=2, f1=f1, f2=f2, plog=FALSE)
sig.R:tsplot<-function(a)
sig.R:brune<-function(f,tstar,fc,Omega0)
sig.R:gettstar<-function(H, f,fc,Omega0)
sig.R:extrema<-function(x, y, smooth=TRUE , plot=FALSE)
sig.R:trapz<-function(y, dt)
sig.R:applytaper<-function(f, p=0.05)
sig.R:detrend<-function(x)
sig.R:mtapspec<-function(a, dt, klen=length(a), MTP=NULL)
sig.R:plot.mtm<-function(ma, PLOT=1, LOG='x')
sig.R:evolfft<-function(a, dt=0, Nfft=0,  Ns=0, Nov=0, fl=0, fh=10 )
sig.R:plotevol<-function(DEVOL, log=0,  fl=0, fh=10 , col=col, ylog=FALSE)
sourcer.R:                          command=function()cat("OW!\n"))
sourcer.R:f<-function(...) {
sourcer.R:img <-tkrplot(tt, function() plot(1:20,(1:20)^bb))
sourcer.R:f<-function(...) {
sourcer.R:img <-tkrplot(tt, function() srot(ascd, bb))
sourcer.R:f<-function(...) {
SPEC.dir.R:geign <- function(){
SPEC.dir.R:perplot <- function()
SPEC.dir.R:plcohere <- function()
SPEC.dir.R:cohspec <- function(blah, fil) {
SPEC.dir.R:cohere <- function(a, b)
SPEC.dir.R:spikeplot<-function(x,y)
SPEC.dir.R:jps <- function()
Specfit.R:brune.func<-function( freq,  omega0,  tstar0, fc, alpha, gamma)
Specfit.R:brune.score<-function(f, spec, jbrune, f1, f2 )
Specfit.R:brune.search1<-function(f, sp, omega0, tstar0 , fcorn,  alpha, gamma )
Specfit.R:brune.search<-function(infreq, inspec, f1, f2, omega0, fcorn, tstar0, gamma  )
Specfit.R:get.cornerA<-function( x, y)
Specfit.R:get.corner<-function( INfreq, INspec, dt, f1, f2, PLOT=FALSE)
stromb.R:dump.NLOOK<-function(mSW, COMP=c(1:4), FV=NULL, RMINST=RMINST, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, sfact=1, Pause=0 )
stromb.R:dump.LOOK<-function(mSW, STA="", FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, rot=FALSE, sfact=1, Pause=0 )
stromb.R:dump.CLOOK<-function(mSW, STA="", FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, sfact=1 )
stromb.R:dump.SPIN<-function(mSW, STA="", FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),  INTEG=FALSE, rot=FALSE, sfact=1 )
stromb.R:dump.Nplot<-function(FIG, COMP=c(1:4), WIN=c(0,1), RMINST=FALSE, INTEG=FALSE, FILT=f, rot=FALSE, CINTEG=c(FALSE), LABS=LABS, sfact=1)
stromb.R:ename<-function(NNN)
stromb.R:stromb.rot<-function(SIN, comps=c(2,3,4), vent=1, STA=stastromb, INTEG=FALSE, FILT=FALSE, rot=FALSE, sfact=1)
stromb.R:spin.rot<-function(SIN, comps=c(2,3,4),  angles=a, INTEG=FALSE, FILT=FALSE , sfact=1, PLOT=TRUE, LOOK=TRUE, FPARM=fparm, WIN=w)
stromb.R:dump.rot<-function(SIN, comps=c(2,3,4), ang=45, INTEG=FALSE, FILT=FALSE , PLOT=TRUE, sfact=1, FPARM=fparm, WIN=w)
stromb.R:ROTSTROM<-function(SALL, lens , CMP, ANGS)
stromb.R:SEESTROM<-function(SALL, lens , CMP, ANGS)
stromb.R:PLOT.ADDN<-function(n, pch=1, col=1)
stromb.R:PLOT.MATN<-function(ascd, tim=1, dt=1,  WIN=WIN, labs=LABS, notes=notes, sfact=1, LOG="", COL=col)
stromb.R:PLOT.LISTN<-function(ldat , N=N, tim=1,    WIN=WIN, labs=LABS, notes=notes, sfact=1,  COL=col)
stromb.R:PLOT.MAT3<-function(ascd, tim=1, labs=LABS, sfact=1)
stromb.R:PLOT.MAT.Aold<-function(ascd, tim=1, labs=LABS)
stromb.R:PLOTSTROMB<-function(NSTRMAP, lonrange=lonrange, latrange=latrange, STA=stastromb)
stromb.R:strmbcont<-function(x,y,z, levs=15,  POL=POL)
stromb.R:pstromb<-function(STRMAP=STRMAP, lonrange=lonrange, latrange=lonrange, STA=stastromb, tight=FALSE)
stromb.R:check.stromb<-function(SN, INTEG=FALSE, FILT=FALSE, rot=FALSE)
stromb.R:check.stromb<-function(SN, J, INTEG=FALSE, FILT=FALSE, rot=FALSE)
stromb.R:stromb.AIN<-function(A1, J, vent, STA, INTEG=FALSE, FILT=f)
stromb.R:check.response<-function(EJALL)
stromb.R:stromb.STACKER<-function(SELDA, tvert=30, prevert=7, tacoust=3, preacoust=1, WINDOWS=NULL)
stromb.R:etna.PIX<-function(mSW, STA=seq(from=2, to=7), FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),
stromb.R:dump.Nstagath<-function(mSW, STA="KAC", COMP="Vertical", NUM=5, BEF=1, AFT=5, FV=NULL)
stromb.R:dump.NPIX<-function(mSW, STA=seq(from=2, to=7), FV=NULL, RMINST=FALSE, CINTEG=c(FALSE, TRUE, TRUE, TRUE),
stromb.R:etna.Xtract<-function(mSW, STA=2, Wlen=30, FV=FV, sfact=1, Pause=0, PLT=FALSE )
stromb.R:etna.Xtract3<-function(mSW, STA=c(2,3,4),  ORD=ORD , Wlen=30, FV=FV, sfact=1, Pause=0, PLT=FALSE )
stromb.R:etna.Xtract4<-function(mSW, STA=c(2,3,4),  ORD=FALSE , Wlen=30, FV=FV, sfact=1, Pause=0, PLT=FALSE )
stromb.R:etna.test<-function(SMAT)
stromb.R:etna.clust1<-function(SMAT, FH=FH, dt=dt, KOR=KOR, KCLUST=temp, ORD=FALSE, DECON=FALSE, PRE=NULL, TIT="", GROUPS=3)
stromb.R:etna.exam<-function(amp, dt, flow = 1/50,  fhigh = 1/5)
stromb.R:etna.3comp<-function(DA1, DA2, DA3)
SUBS.R:libs<-function()
SUBS.R:sourcer<-function()
SUBS.R:recdate<-function(jday, hr, mi, sec)
SUBS.R:secdif<-function(jday, hr, mi, sec, jday2, hr2, mi2, sec2)
SUBS.R:jday<-function(yr, mo, day)
SUBS.R:plotK<-function(K)
SUBS.R:plotK3<-function(K)
SUBS.R:error.bar<-function(x, y, lo, hi, pch=sym, col=clr, add=FALSE)
SUBS.R:read.DUMP<-function(infof, dataf)
SUBS.R:hypot<-function(x1,y1,x2,y2)
SUBS.R:hypotL<-function(x)
SUBS.R:rotmat<-function(alph)
SUBS.R:rotseis<-function(k)
SUBS.R:rotseis2<-function(A, k)
SUBS.R:plotseis3<-function(A, amp)
SUBS.R:pwd<-function()
SUBS.R:find.file <- function(str){
SUBS.R:cormat<-function(n, h)
SUBS.R:slice<-function(x,y,r,a1, a2, col)
SUBS.R:flower<-function(x, y, z)
SUBS.R:greatAz<-function(lat1, lon1, lat2, lon2)
SUBS.R:GCLCFR<-function(t)
SUBS.R:gclc<-function(phiorg, lamorg, phi, lam)
SUBS.R:STA.dist<-function(  azis, zvals, asta)
SUBS.R:distaz<-function( olat, olon, tlat, tlon)
SUBS.R:reg<-function(x,y)
SUBS.R:VarRed<-function(X,Y, SEQ=SEQ)
SUBS.R:peaks<-function(series,span=3)
topo.R:CALCOL<-function()
topo.R:xytopper<-function(jtop)
topo.R:jtoppersp<-function(jtop)
topo.R:trans3d <- function(x,y,z, pmat) {
topo.R:makeimag<-function(x,y,z, d=0.7, col=terrain.colors(50), dx=100)
topo.R:PersPROJmap<-function(cmp, PMAT, pts=NULL, ADD=TRUE, WIN=NULL, ASP=TRUE, NUM=FALSE,
topo.R:PersAXESmap<-function(cmp, PMAT, GRID=FALSE, WIN=NULL)
UWFILES.R:printmc<-function(MC)
UWFILES.R:getmcard<-function(pf)
UWFILES.R:getpfile<-function(pf)
UWFILES.R:Grake<-function(MC)
UWFILES.R:GOALLM<-function(Mpf)
UWFILES.R:Get1Dvel<-function(infile)
UWFILES.R:Plot1Dvel<-function(v)
UWFILES.R:Comp1Dvel<-function(v, v2)
UWFILES.R:jlegend<-function (x, y, legend, fill, col = "black", lty, lwd, pch, angle = NULL,
UWFILES.R:    rect2 <- function(left, top, dx, dy, angle, ...) {
UWFILES.R:    segments2 <- function(x1, y1, dx, dy, ...) {
UWFILES.R:    points2 <- function(x, y, ...) {
UWFILES.R:    text2 <- function(x, y, ...) {
UWFILES.R:        catn <- function(...) do.call("cat", c(lapply(list(...),
volc.R:vdistF<-function(file)
volc.R:vdist1<-function(vkam)
volc.R:vproj<-function(vkam)
volc.R:vdist2<-function(M)
volc.R:vpicker<-function(ALLVOL)
XCOR.R:############  functions for performing cross correlations in R (splus)
XCOR.R:getlim<-function(a1, a2, dt)
XCOR.R:xchug<-function(a1,a2, lim1, dt,tskip)
XCOR.R:xchug1<-function(i4, s4, lim1, dt, tskip)
XCOR.R:xcrude<-function(a1, a2, dt, tskip)
XCOR.R:xcor2<-function(a1, a2, dt, PLOT=FALSE, LAG=100)
XCOR.R:xcor3<-function(a1, a2, dt, PLOT=FALSE)
XCOR.R:JXCOR<-function(dat, dt, PLOT=FALSE, PIC=FALSE)
XCOR.R:PLXCOR<-function(dat, dt, lags, v)
XMA.R:####  functions for reading in tomographic images
XMA.R:GXMA<-function(name)
XMA.R:GXMA3D<-function(name)
XMA.R:pltomo<-function(x,y,MOD,i)
XMA.R:pMOD<-function(x,y,D, MOD)
XMA.R:CMOD<-function(KMOD)
xmcolor.R:setXMCOL<-function()
xmcolor.R:SHOWPAL<-function(XMCOL)
ZOEP.R:#  function output= [angle, rmat, rra, ria, ang] =
ZOEP.R:newzoep<-function(icoef, vp1, vp2, vs1, vs2, rho1, rho2, incw)
ZOEP.R:plotzpr<-function( chtype , alpha1, alpha2, beta1, beta2, rho1 ,rho2, chincw, choutkind)
