as-500/hw21/.Rhistory

fe = read.table("shock.fe")
fe
plot(fe)
plot(fe[1,])
fe[1.]
fe[1,]
fe[1,]
fe[1,1]
fe[1,2]
fe[1,,]
fe[2,]
fe[0,]
fe[1,]
plot,x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28),y=fe[1,])
plot(x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28),y=fe[1,])
1..28
1...28
?seq
seq(1..28)
seq(1,28)
plot(seq(1,28),fe[1,])
plot(seq(1,27),fe[1,2..28])
plot(seq(1,27),fe[1,2..27])
plot(seq(1,27),fe[1,2...27])
plot(seq(1,27),fe[1,])
plot(seq(1,28),fe[1,])
plot(seq(1,28),fe[1,2])
plot(seq(1,28),fe[1,0])
plot(seq(1,28),fe[1,])
plot(seq(1,28),fe[1,>1])
plot(seq(1,28),subset(fe[1,],>1)
?subset
fe[c(seq(2,28)]
fe[c(seq(2,28))]
fe[c(seq(1,28))]
fe[c(seq(2,28)]
plot(seq(1,28),fe[1,c(seq(2,28))])
plot(seq(2,28),fe[1,c(seq(2,28))])
plot(seq(2,28),fe[1,c(seq(2,28))],type="b")
pdf("fe_stages.pdf")
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction")
dv.off()
dev.off()
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction")
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(seq(1,28)))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,1))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,1,28))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,28))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,27))
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,27))
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,27))
fe
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,0))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(1,28,0))
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,0))
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
?for
?
for
?for
for
axis(1,"yo")
?axis
axis(1,at=1,"yo"0
axis(1,at=1,"yo)
axis(1,at=1,"yo")
axis(1,at=1,label="yo")
fe[1,1]
fe[0,1]
fe[0,2]
fe[1,1]
fe[1,0]
fe[1,]
fe[1,1]
fe[2,1]
fe[2,2]
fe[2,3]
fe[1,3]
fe[0,3]
fe[0,5]
fe[1,5]
fe[1,1]
name(fe[1,1])
names(fe[1,1])
names(fe[1,2])
names(fe[1,3])
fe
fe[1,]
fe[1,1]
fe[1,2]
fe[1,0]
fe[1,1]
fe[1,2]
fe
title
title(fe[1,1])
fe[1]
fe[2]
fe[3]
names(fe[3])
names(fe[4])
names(fe[5])
names(fe[6])
names(fe[1])
names(fe[2])
fe
names(fe[3])
names(fe[4])
axis(1,at=2,label="Fe")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
axis(1,at=2,label="Fe")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages.")
axis(1,at=2,label="Fe")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds")
subtitle("v=10,000 km/s, T = 2e9 K")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K")
?title
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(3,28),axis(1,at=coord,label=paste("FE+",coord-2)))
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(3,28)){axis(1,at=coord,label=paste("FE+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(3,28,2)){axis(1,at=coord,label=paste("FE+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste("FE+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste("FE-",coord-2))}
?paste
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE-",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE^+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE",^+,coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE",^"+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE","^+",coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE",parse("^+"),coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=paste0("FE",parse(text="^+"),coord-2))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("FE","^+",coord-2)))
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("FE","^+",coord-2))}
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("FE","^+",coord-2)))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("FE","^+",coord-2)))}
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
#title("Ionization Fraction for Iron Stages, AGN Shock Winds\rv=10,000 km/s, T = 2e9 K",main="")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
dev.off()
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
dev.off()
oxy= read.table("shock.oxy",header=T)
oxy
pdf("oxy_ion.pdf")
oxy= read.table("shock.oxy",header=T)
plot(seq(2,12),oxy[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Oxygen Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="O")
for(coord in seq(4,12,2)){axis(1,at=coord,label=parse(text=paste0("O","^+",coord-2)))}
dev.off()
oxy= read.table("shock.oxy",header=T)
plot(seq(2,12),oxy[1,c(seq(2,12))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
xy[1,c(seq(2,12))]
oxy[1,c(seq(2,12))]
pdf("oxy_ion.pdf")
oxy= read.table("shock.oxy",header=T)
plot(seq(2,12),oxy[1,c(seq(2,12))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Oxygen Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="O")
for(coord in seq(4,12,2)){axis(1,at=coord,label=parse(text=paste0("O","^+",coord-2)))}
dev.off()
pdf("oxy_ion.pdf")
oxy= read.table("shock.oxy",header=T)
plot(seq(2,10),oxy[1,c(seq(2,10))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Oxygen Stages, AGN Shock Winds (v=10,000 km/s, T = 2e9 K)")
axis(1,at=2,label="O")
for(coord in seq(4,10,2)){axis(1,at=coord,label=parse(text=paste0("O","^+",coord-2)))}
dev.off()
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Iron Stages, AGN Shock Winds (v=10,000 km/s)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
dev.off()
pdf("oxy_ion.pdf")
oxy= read.table("shock.oxy",header=T)
plot(seq(2,10),oxy[1,c(seq(2,10))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Oxygen Stages, AGN Shock Winds (v=10,000 km/s)")
axis(1,at=2,label="O")
for(coord in seq(4,10,2)){axis(1,at=coord,label=parse(text=paste0("O","^+",coord-2)))}
dev.off()
pdf("fe_ion.pdf")
fe = read.table("shock.fe",header=T)
plot(seq(2,28),fe[1,c(seq(2,28))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for Fe Stages, AGN Shock Winds (v=10,000 km/s)")
axis(1,at=2,label="Fe")
for(coord in seq(4,28,2)){axis(1,at=coord,label=parse(text=paste0("Fe","^+",coord-2)))}
dev.off()
pdf("oxy_ion.pdf")
oxy= read.table("shock.oxy",header=T)
plot(seq(2,10),oxy[1,c(seq(2,10))],type="b",xlab="Ionization Stage",ylab="Ionization Fraction",xaxp=c(0,0,1))
title("Ionization Fraction for O Stages, AGN Shock Winds (v=10,000 km/s)")
axis(1,at=2,label="O")
for(coord in seq(4,10,2)){axis(1,at=coord,label=parse(text=paste0("O","^+",coord-2)))}
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con")
plot(cont[,1],cont[,7])
cont
pdf("cont.pdf")
cont = read.table("shock.con")
plot(cont[,1],cont[,7])
dev.off()
?plot
?Read
?read
?read
?read.table
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,7])
dev.off()
cont
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4])
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="xy")
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x")
dev.off()
plot(cont)
plot(cont)
plot(cont)
plot(cont[,1],cont[,4],log="x")
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab="Intensity (ergs/cm/cm/s)")
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab="Flux (ergs/cm/cm/s)")
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab="Flux (ergs/cm^2/s)")
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text="Flux (ergs/cm^2/s)"))
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux,(ergs/cm^2/s)")))
dev.off()
pdf("cont.pdf")
cont = read.table("shock.con",fill=T)
plot(cont[,1],cont[,4],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
dev.off()
cont[,1]
cont[,4]
cont[,3]
cont[,2]
cont = read.table("cont.preplot")
plot(cont)
plot(cont)
cont
cont[1]
cont[2]
plot(cont[1],cont[2])
cont[1,]
cont[,1]
cont[,2]
pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
dev.off()
pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],type="l",log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
dev.off()
pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],type="l",log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
title("Emitted Continuum, AGN Shock Winds (v=10,000 km/s, T=2e9 K) ")
dev.off()pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],type="l",log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
title("Emitted Continuum, AGN Shock Winds (v=10,000 km/s, T=2e9 K) ")
dev.off()
pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],type="l",log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
title("Emitted Continuum, AGN Shock Winds (v=10,000 km/s, T=2e9 K) ")
dev.off()
pdf("cont.pdf")
cont = read.table("cont.preplot",fill=T)
plot(cont[,1],cont[,2],type="l",log="x",xlim=c(1e-5,1e4),xlab="Wavelength (microns)",ylab=parse(text=paste("Flux","(ergs/cm^2/s)")))
title("Emitted Continuum, AGN Shock Winds (v=10,000 km/s, T=2e9 K) ")
dev.off()
exit
quit()