Description Usage Arguments Details Value Author(s) References See Also Examples

Compute the Right Ascension and Declination of the Moon at specified Julian date(s)

1 | ```
moonpos(jd, radian=F)
``` |

`jd` |
Julian ephemeris date, scalar or vector |

`radian` |
if =TRUE, then all output variables are given in radians rather than degrees (default=FALSE) |

The method is derived from the Chapront ELP2000/82 Lunar Theory (Chapront-Touze and Chapront, 1983), as described by Jean Meeus (1998, Chpt. 47). Meeus quotes an approximate accuracy of 10" in longitude and 4" in latitude, but he does not give the time range for this accuracy. Comparison of this procedure with the example in “Astronomical Algorithms” reveals a very small discrepancy (~1 km) in the distance computation, but no difference in the position calculation.

`ra` |
apparent right ascension of the moon, referred to the true equator of the specified date(s), in degrees |

`dec` |
declination of the moon, in degrees |

`dis` |
Earth-moon distance between the center of the Earth and the center of the Moon, in km |

`geolong` |
apparent longitude of the moon referred to the ecliptic of the specified date(s), in degrees |

`geolat` |
apparent longitude of the moon referred to the ecliptic of the specified date(s), in degrees |

Written W. Landsman

R adaptation by Arnab Chakraborty June 2013

Chaprint-Touze, M, and Chapront, J. 1983, The lunar ephemeris ELP 2000, Astron. Astrophys. 124, 50-62.

Meeus, J., 1998, “Astronomical Algorithms”, 2nd ed., Willmann-Bell

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
# Find the position of the moon on April 12, 1992
# Result: 08 58 45.23 +13 46 6.1
# This is within 1" from the position given in the Astronomical Almanac
jd = jdcnv(1992,4,12,0) # get Julian date
pos = moonpos(jd) # get RA and Dec of moon
adstring(pos$ra,pos$dec,1)
# Plot the Earth-moon distance for every day at 0 TD in July, 1996
jd = jdcnv(1996,7,1,0) # get Julian date of July 1
pos = moonpos(jd+rep(0,31)) # Moon position at all 31 days
plot(seq(1,31), pos$dis, xlab="July 1996 (day)", ylab="Lunar distance (km)")
``` |

```
[1] " 8 58 45.34 +13 46 5.3"
Warning messages:
1: In d * d_lng :
longer object length is not a multiple of shorter object length
2: In m * m_lng :
longer object length is not a multiple of shorter object length
3: In mprime * mp_lng :
longer object length is not a multiple of shorter object length
4: In f * f_lng :
longer object length is not a multiple of shorter object length
5: In omega * om_lng :
longer object length is not a multiple of shorter object length
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.