/test/plot_lnP.R
https://bitbucket.org/NeilMiller/mixeos · 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()