/test/plot_lnP.R
R | 28 lines | 25 code | 3 blank | 0 comment | 0 complexity | 2a93b85859b11314d3df0a9ebbc1d94f MD5 | raw file
- filename = "test_lnP_lnRho.dat"
- lnPvect = read.table(filename)
- plotfile = "plot_lnP_lnRho.ps"
- postscript(plotfile)
- plot(lnPvect[,1],lnPvect[,3], xlab="log10(Rho)", ylab="lnP",type="l")
- plot(lnPvect[,1],lnPvect[,4], xlab="log10(Rho)", ylab="dlnP_dlnRho",type="l")
- lines(lnPvect[,1],lnPvect[,5], col="blue")
- lines(lnPvect[,1],lnPvect[,6], col="red")
- plot(lnPvect[,1],lnPvect[,7], xlab="log10(Rho)", ylab="d2lnP_dlnRho2", type="l")
- lines(lnPvect[,1],lnPvect[,8], col="red")
- dev.off()
- filename = "test_lnP_lnT.dat"
- lnPvect = read.table(filename)
- plotfile = "plot_lnP_lnT.ps"
- postscript(plotfile)
- plot(lnPvect[,1],lnPvect[,3], xlab="log10(T)", ylab="lnP",type="l")
- plot(lnPvect[,1],lnPvect[,4], xlab="log10(T)", ylab="dlnP_dlnT",type="l")
- lines(lnPvect[,1],lnPvect[,5], col="blue")
- lines(lnPvect[,1],lnPvect[,6], col="red")
- plot(lnPvect[,1],lnPvect[,7], xlab="log10(T)", ylab="d2lnP_dlnT2", type="l")
- lines(lnPvect[,1],lnPvect[,8], col="red")
- plot(lnPvect[,1],lnPvect[,9], xlab="log10(T)", ylab="d2lnP_dlnRhodlnT", type="l")
- lines(lnPvect[,1],lnPvect[,10], col="blue")
- lines(lnPvect[,1],lnPvect[,11], col="red")
- dev.off()