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. 

  2. lnPvect = read.table(filename)
    3. 

  3. 

  4. plotfile = "plot_lnP_lnRho.ps"
    5. 

  5. postscript(plotfile)
    6. 

  6. plot(lnPvect[,1],lnPvect[,3], xlab="log10(Rho)", ylab="lnP",type="l")
    7. 

  7. plot(lnPvect[,1],lnPvect[,4], xlab="log10(Rho)", ylab="dlnP_dlnRho",type="l")
    8. 

  8. lines(lnPvect[,1],lnPvect[,5], col="blue")
    9. 

  9. lines(lnPvect[,1],lnPvect[,6], col="red")
    10. 

  10. plot(lnPvect[,1],lnPvect[,7], xlab="log10(Rho)", ylab="d2lnP_dlnRho2", type="l")
    11. 

  11. lines(lnPvect[,1],lnPvect[,8], col="red")
    12. 

  12. dev.off()
    13. 

  13. 

  14. filename = "test_lnP_lnT.dat"
    15. 

  15. lnPvect = read.table(filename)
    16. 

  16. 

  17. plotfile = "plot_lnP_lnT.ps"
    18. 

  18. postscript(plotfile)
    19. 

  19. plot(lnPvect[,1],lnPvect[,3], xlab="log10(T)", ylab="lnP",type="l")
    20. 

  20. plot(lnPvect[,1],lnPvect[,4], xlab="log10(T)", ylab="dlnP_dlnT",type="l")
    21. 

  21. lines(lnPvect[,1],lnPvect[,5], col="blue")
    22. 

  22. lines(lnPvect[,1],lnPvect[,6], col="red")
    23. 

  23. plot(lnPvect[,1],lnPvect[,7], xlab="log10(T)", ylab="d2lnP_dlnT2", type="l")
    24. 

  24. lines(lnPvect[,1],lnPvect[,8], col="red")
    25. 

  25. plot(lnPvect[,1],lnPvect[,9], xlab="log10(T)", ylab="d2lnP_dlnRhodlnT", type="l")
    26. 

  26. lines(lnPvect[,1],lnPvect[,10], col="blue")
    27. 

  27. lines(lnPvect[,1],lnPvect[,11], col="red")
    28. 

  28. dev.off()
    29. 

  29.