Nothing
### This function is a component of astrochron: An R Package for Astrochronology
### Copyright (C) 2018 Stephen R. Meyers
###
###########################################################################
### getLaskar: download Laskar et al. (2004; 2011) astronomical solutions.
### (SRM: January 22, 2015; January 4, 2016; February 7, 2017;
### June 19, 2018)
###
###########################################################################
getLaskar <- function (sol="la04",verbose=T)
{
i=0
if(sol=="la04" || sol=="la10a" || sol=="la10b" || sol=="la10c" || sol=="la10d" || sol=="la11" || sol=="insolation") i=1
if(i==0)
{
cat("\n**** ERROR: specified solution is not available.\n")
stop(" TERMINATING NOW!")
}
tempLaskar <- tempfile()
# download data
if(sol=="la04")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2004) astronomical solution: La2004\n\n")
cat(" Please cite: Laskar, J., Robutel, P., Joutel, F., Gastineau, M.,\n")
cat(" Correia, A.C.M., Levrard, B., 2004, A long term numerical solution \n")
cat(" for the insolation quantities of the Earth: Astron. Astrophys., Volume 428, 261-285.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la04.txt.bz2",tempLaskar)
}
if(sol=="la10a")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2011) astronomical solution: La2010a\n\n")
cat(" Please cite: Laskar, J., Fienga, A., Gastineau, M., Manche, H., 2011,\n")
cat(" La2010: A new orbital solution for the long-term motion of the Earth:\n")
cat(" Astron. Astrophys., Volume 532, A89.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la10a.txt.bz2",tempLaskar)
}
if(sol=="la10b")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2011) astronomical solution: La2010b\n\n")
cat(" Please cite: Laskar, J., Fienga, A., Gastineau, M., Manche, H., 2011,\n")
cat(" La2010: A new orbital solution for the long-term motion of the Earth:\n")
cat(" Astron. Astrophys., Volume 532, A89.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la10b.txt.bz2",tempLaskar)
}
if(sol=="la10c")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2011) astronomical solution: La2010c\n\n")
cat(" Please cite: Laskar, J., Fienga, A., Gastineau, M., Manche, H., 2011,\n")
cat(" La2010: A new orbital solution for the long-term motion of the Earth:\n")
cat(" Astron. Astrophys., Volume 532, A89.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la10c.txt.bz2",tempLaskar)
}
if(sol=="la10d")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2011) astronomical solution: La2010d\n\n")
cat(" Please cite: Laskar, J., Fienga, A., Gastineau, M., Manche, H., 2011,\n")
cat(" La2010: A new orbital solution for the long-term motion of the Earth:\n")
cat(" Astron. Astrophys., Volume 532, A89.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la10d.txt.bz2",tempLaskar)
}
if(sol=="la11")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2011) astronomical solution: La2011\n\n")
cat(" Please cite: Laskar, J., Fienga, A., Gastineau, M., Manche, H., 2011,\n")
cat(" La2010: A new orbital solution for the long-term motion of the Earth:\n")
cat(" Astron. Astrophys., Volume 532, A89.\n")
cat(" AND:\n")
cat(" Laskar, J., Gastineau, M., Delisle, J.-B., Farres, A., Fienga, A.: 2011,\n")
cat(" Strong chaos induced by close encounters with Ceres and Vesta:\n")
cat(" Astron. Astrophys., Volume 532, L4.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/la11.txt.bz2",tempLaskar)
}
if(sol=="insolation")
{
if(verbose)
{
cat(" * Downloading Laskar et al. (2004) 65 deg North summer insolation\n\n")
cat(" Please cite: Laskar, J., Robutel, P., Joutel, F., Gastineau, M.,\n")
cat(" Correia, A.C.M., Levrard, B., 2004, A long term numerical solution \n")
cat(" for the insolation quantities of the Earth: Astron. Astrophys., Volume 428, 261-285.\n")
}
download.file("http://www.geology.wisc.edu/~smeyers/astrochron/insolation.txt.bz2",tempLaskar)
}
if(verbose) cat(" * Decompressing solution\n")
la <- read.table(bzfile(tempLaskar),header=T)
unlink(tempLaskar)
if(sol != "la11" && sol != "insolation" ) la=data.frame(cbind((1:249001)-1,la))
if(sol == "la11") la=data.frame(cbind((1:100000)-1,la))
if(sol == "insolation") la=data.frame(cbind((1:20001)-1,la))
colnames(la)[1] = c("Time_ka")
return(la)
### END function getLaskar
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.