/test/plot_lnP.R

https://bitbucket.org/NeilMiller/mixeos · R · 28 lines · 25 code · 3 blank · 0 comment · 0 complexity · 2a93b85859b11314d3df0a9ebbc1d94f MD5 · raw file

  1. filename = "test_lnP_lnRho.dat"
  2. lnPvect = read.table(filename)
  3. plotfile = "plot_lnP_lnRho.ps"
  4. postscript(plotfile)
  5. plot(lnPvect[,1],lnPvect[,3], xlab="log10(Rho)", ylab="lnP",type="l")
  6. plot(lnPvect[,1],lnPvect[,4], xlab="log10(Rho)", ylab="dlnP_dlnRho",type="l")
  7. lines(lnPvect[,1],lnPvect[,5], col="blue")
  8. lines(lnPvect[,1],lnPvect[,6], col="red")
  9. plot(lnPvect[,1],lnPvect[,7], xlab="log10(Rho)", ylab="d2lnP_dlnRho2", type="l")
  10. lines(lnPvect[,1],lnPvect[,8], col="red")
  11. dev.off()
  12. filename = "test_lnP_lnT.dat"
  13. lnPvect = read.table(filename)
  14. plotfile = "plot_lnP_lnT.ps"
  15. postscript(plotfile)
  16. plot(lnPvect[,1],lnPvect[,3], xlab="log10(T)", ylab="lnP",type="l")
  17. plot(lnPvect[,1],lnPvect[,4], xlab="log10(T)", ylab="dlnP_dlnT",type="l")
  18. lines(lnPvect[,1],lnPvect[,5], col="blue")
  19. lines(lnPvect[,1],lnPvect[,6], col="red")
  20. plot(lnPvect[,1],lnPvect[,7], xlab="log10(T)", ylab="d2lnP_dlnT2", type="l")
  21. lines(lnPvect[,1],lnPvect[,8], col="red")
  22. plot(lnPvect[,1],lnPvect[,9], xlab="log10(T)", ylab="d2lnP_dlnRhodlnT", type="l")
  23. lines(lnPvect[,1],lnPvect[,10], col="blue")
  24. lines(lnPvect[,1],lnPvect[,11], col="red")
  25. dev.off()