PageRenderTime 28ms CodeModel.GetById 22ms RepoModel.GetById 1ms app.codeStats 0ms

/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()