MACRO gengen
par(mfrow=c(2,2))
?gen(ntp.isx)
?gen(ntp.ncr)
?gen(ntp.umd)
END
MACRO gengen3
par(mfrow=c(2,1))
?gen2(min.cal,2,.2)
?gen2(min.eas,2,.2)
?gen2(min.eur,8,2)
END
MACRO gengen1
par(mfrow=c(2,2))
?gen(ntp.ncry)
?gen(ntp.pscy)
?gen(ntp.uiuy)
?gen(ntp.umdy)
END
MACRO gengen2
par(mfrow=c(2,2))
?gen(ntp.frdz)
?gen(ntp.isxz)
?gen(ntp.ncrz)
?gen(ntp.umdz)
END
MACRO gengen4
par(mfrow=c(2,2))
?gen(ntp.frdw)
?gen(ntp.isxw)
?gen(ntp.ncrw)
?gen(ntp.umdw)
END
MACRO gen(f)
par(cex=1.6,mex=1.6)
?series(f,3)
?wedge(f)
?series(f,4)
?rank(f,3,4)
END
MACRO gen2(f,g,h)
?series2(f,3,g)
?series2(f,4,h)
END
MACRO ntp(f,g)
g<-matrix(read("f"),ncol=4,byrow=TRUE)
g<-g[8:nrow(g),]
END
MACRO ntpgen
?ntp(isx.flz,ntp.isx)
?ntp(ncr.flz,ntp.ncr)
?ntp(umd.flz,ntp.umd)
?ntp(cal.flz,ntp.cal)
?ntp(eas.flz,ntp.eas)
?ntp(eur.flz,ntp.eur)
?ntp(ncry.flz,ntp.ncry)
?ntp(pscy.flz,ntp.pscy)
?ntp(uiuy.flz,ntp.uiuy)
?ntp(umdy.flz,ntp.umdy)
?ntp(frdz.flz,ntp.frdz)
?ntp(isxz.flz,ntp.isxz)
?ntp(ncrz.flz,ntp.ncrz)
?ntp(umdz.flz,ntp.umdz)
END
MACRO min(f,g)
y2 <- f[,3]+f[,4]
y4 <- y2+f[,5]+f[,6]
y8 <- y4+f[,7]+f[,8]+f[,9]+f[,10]
g <- cbind(f[,1],f[,2],f[,3],y2/2,y4/4,y8/8)
END
MACRO mingen
?min(rank.cal,mean.cal)
?min(rank.eas,mean.eas)
?min(rank.eur,mean.eur)
?min(rank.isx,mean.isx)
?min(rank.ncr,mean.ncr)
?min(rank.umd,mean.umd)
END
MACRO mingen1
?min(rank.ncry,mean.ncry)
?min(rank.pscy,mean.pscy)
?min(rank.uiuy,mean.uiuy)
?min(rank.umdy,mean.umdy)
END
MACRO med(f,g)
x <- f[,3:10]
y3 <- x[,order(x[,1:3])[2]]
y5 <- x[,order(x[,1:5])[3]]
y9 <- x[,order(x[,1:9])[5]]
g <- cbind(f[,1],f[,2],f[,3],y3,y5,y9)
END
MACRO rank(f,g,h)
x1 <- sort(abs(f[,g])) ; x1 <- x1/1000
y1 <- 100*(len(x1)+1-rank(x1))/len(x1)
x1[x1<.001] <- .001 ; x1[x1>9.99] <- NA
plot(x1,y1,type="l",log="xy",xlab="Offset (sec)",ylab="Prob (%)",main="f",
xlim=c(.001,10),ylim=c(.01,100))
x1 <- sort(abs(f[,h])) ; x1 <- x1/1000
y1 <- 100*(len(x1)+1-rank(x1))/len(x1)
x1[x1<.001] <- .001 ; x1[x1>9.99] <- NA
lines(x1,y1)
END
MACRO rankgen
par(mfrow=c(2,2))
?rank(min.cal,3,6,8)
?rank(min.eas,3,6,8)
?rank(min.eur,3,6,8)
?rank(min.isx,3,6,8)
?rank(min.ncr,3,6,8)
?rank(min.umd,3,6,8)
END
MACRO rankgen1
par(mfrow=c(2,2))
?rank(min.ncry,3,6,8)
?rank(min.pscy,3,6,8)
?rank(min.uiuy,3,6,8)
?rank(min.umdy,3,6,8)
END
MACRO resp(f)
x <- f[,1] ; y <- f[,4]
plot(x,y,type="l",xlab="Time (min)",ylab="Response (%)",main="f")
y <- f[,5] ; lines(x,y)
!y <- f[,4]
!lines(x,y)
!y <- f[,5]
!lines(x,y)
END
MACRO resp1(f)
x <- f[,1] ; y <- f[,2]
plot(x,y,type="l",xlab="Time (min)",ylab="Response (%)",main="f")
!y <- f[,4]
!lines(x,y)
!y <- f[,5]*1000/60
!plot(x,y,type="l",xlab="Time (min)",ylab="Drift (ppm)",main="f")
y <- f[,3]*1000/60
lines(x,y)
END
MACRO respgen
par(mfrow=c(2,2))
!?resp(filt.1)
!?resp(filt.3)
?resp1(filt.4)
END
MACRO delay(f)
x <- f[,1] ; y <- f[,2]
y <- abs(y/1000) ; y[y>9.99] <- NA
plot(x,y,type="l",log="y",xlab="Time (UT hours)",ylab="Delay (sec)",main="f",
ylim=c(.1,9.99))
END
MACRO delgen
par(mfrow=c(2,2))
?delay(rank.cal)
?delay(rank.eas)
?delay(rank.eur)
?delay(rank.isx)
?delay(rank.ncr)
?delay(rank.umd)
END
MACRO delgen1
par(mfrow=c(2,2))
?delay(rank.ncry)
?delay(rank.pscy)
?delay(rank.uiuy)
?delay(rank.umdy)
END
MACRO series(f,g)
x <- f[,1]
y <- abs(f[,g]/1000) ; y[y<.001] <- .001 ; y[y>9.99] <- NA
plot(x,y,type="l",log="y",xlab="Time (UT hours)",ylab="Offset (sec)",main="f",
ylim=c(.001,9.99))
END
MACRO series2(f,g,h)
x <- f[,1]
y <- f[,g]/1000 ; y[y< -h] <- NA ; y[y>h] <- NA
plot(x,y,type="l",xlab="Time (UT hours)",ylab="Offset (sec)",main="f",
ylim=c(-h,h))
END
MACRO sergen
par(mfrow=c(2,2))
?series(rank.cal,3)
?series(rank.eas,3)
?series(rank.eur,3)
?series(rank.isx,3)
?series(rank.ncr,3)
?series(rank.umd,3)
END
MACRO sergen2
par(mfrow=c(2,2))
?series(ntp.cal,3) ; ?series(ntp.cal,4)
?series(ntp.eas,3) ; ?series(ntp.eas,4)
?series(ntp.eur,3) ; ?series(ntp.eur,4)
?series(rank.isx,3) ; ?series(ntp.isx,3)
?series(rank.umd,3) ; ?series(ntp.umd,3)
?series(rank.ncr,3) ; ?series(ntp.ncr,3)
END
MACRO sergen1
par(mfrow=c(2,2))
?series(rank.ncry,3)
?series(rank.pscy,3)
?series(rank.uiuy,3)
?series(rank.umdy,3)
END
MACRO wedge(f)
x <- f[,2]/1000 ; x[x>2] <- NA
y <- f[,3]/1000 ; y[abs(y)>1] <- NA
plot(x,y,pch=",",xlab="Delay (sec)",ylab="Offset (sec)",main="f",
xlim=c(0,2),ylim=c(-1,1))
END
MACRO wedgegen
par(mfrow=c(2,2))
?wedge(rank.cal)
?wedge(rank.eas)
?wedge(rank.eur)
?wedge(rank.isx)
?wedge(rank.ncr)
?wedge(rank.umd)
END
MACRO wedgegen1
par(mfrow=c(2,2))
?wedge(rank.ncry)
?wedge(rank.pscy)
?wedge(rank.uiuy)
?wedge(rank.umdy)
END
                                                                                                                                                                                                                                                                                                                                                                                 