fig02a.tex

\documentclass[svgnames,a4paper]{standalone}
\usepackage[T1]{fontenc}
\usepackage[frenchb]{babel}
\usepackage[garamond]{mathdesign}
\usepackage{pstricks-add}
\usepackage{pst-caelum}
\usepackage{xint-meeus}
\usepackage{readarray}
 
\input{xint-sun.tex}
\readdef{lune01.txt}\lignes
\readarray{\lignes}\lune[120,2]
 
\JulianDay(2021,6,10,11,0,58)
\newcommand\calcTfor[1]{\xintdeffloatvar T := (JD+#1-2451545.0)/36525;}
\calcTfor{0}%
\def\RPSun{\FloatToFix{(90-xmDS(T))/20}}
\def\APSun{\FloatToFix{-90-xmADS(T)}}
 
\begin{document}
\begin{pspicture}(-7,-7)(7,7)
\caelum{Graticule}%
\caelum[linewidth=0.2pt,linecolor=DimGray,linestyle=dotted]{Boundaries}
\caelum[linecolor=Crimson]{Ecliptic}%
\psdots[linecolor=Crimson,dotsize=4pt](\RPSun;\APSun)%
\multido{\i=1+6}{20}{%
  \def\dec{\lune[\i,2]}
  \def\asc{\lune[\i,1]}
  \def\RPMoon{\FloatToFix{(90-\dec)/20}}
  \def\APMoon{\FloatToFix{-90-\asc}}
  \psdots[linecolor=DimGray,dotsize=4pt](\RPMoon;\APMoon)
}
\end{pspicture}
\end{document}