Nothing
knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
To compute the position of celestial body or star with SE (Swiss Ephemeris), you do the following steps:
swephR
:library(swephR)
swe_set_ephe_path("C:\\sweph\\ephe")
year <- 2000 month <- 1 day <- 1 hour <- 12 jdut <- swe_julday(year, month, day, hour, SE$GREG_CAL) jdut
ipl <- SE$SUN iflag <- SE$FLG_MOSEPH + SE$FLG_SPEED result <- swe_calc_ut(jdut, ipl, iflag) result
or a fixed star (in below example: Sirius):
starname = "sirius" result <- swe_fixstar2_ut(starname, jdut, iflag) result
options(digits=15) result <- swe_heliacal_ut(jdut,c(0,50,10),c(1013.25,15,50,0.25),c(25,1,1,1,5,0.8),starname, SE$HELIACAL_RISING,SE$HELFLAG_HIGH_PRECISION+SE$FLG_MOSEPH) result
options(digits=6) swe_set_ephe_path(NULL) iflag = SE$FLG_SPEED + SE$FLG_MOSEPH { #get year jyear <- 2000 #get month jmon <- 1 #get day jday <- 1 #get time jhour <- 12 #determine julian day number (at 12:00 GMT) tjd_ut <- swe_julday(jyear, jmon, jday, jhour, SE$GREG_CAL) cat("Julian day number (UT) :", tjd_ut, "(",jyear,",",jmon,",",jday,"; proleptic Gregorian calendar)\n") cat("planet :", c("longitude", "latitude", "distance", "long. speed", "lat. speed"), "\n") cat("===========================================================\n") # loop over all planets for (p in SE$SUN:SE$OSCU_APOG) { # get the name of the planet p objectname = swe_get_planet_name(p) # do the coordinate calculation for this planet p i = swe_calc_ut(tjd_ut, p, iflag) if (i$return < 0) { cat("Error :", i$err, "(", objectname, ")\n") } else { # print data cat (objectname, ":", i$xx[0:5], "\n") } } }
swe_close()
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.