fig07e.tex

\documentclass[svgnames,a4paper]{standalone}
\usepackage[T1]{fontenc}
\usepackage[frenchb]{babel}
\usepackage{bera}
\usepackage{xcolor}
\usepackage{pst-plot,pst-caelum}
 
\usepackage{readarray}
\readdef{202012_jupiter.dta}\rowa
\readarray{\rowa}\jupiter[40,3]
\readdef{202012_saturn.dta}\rowb
\readarray{\rowb}\saturn[40,3]
 
 
\def\toCaelum#1#2{%
  \ftf8{-4(#2-90)/20};\ftf8{-(#1+90)}%
}
 
\def\Saturn#1{\psdots[linecolor=Khaki,dotsize=2pt](\toCaelum{\saturn[#1,1]}{\saturn[#1,2]})}
\def\Jupiter#1{\psdots[linecolor=Peru,dotsize=2pt](\toCaelum{\jupiter[#1,1]}{\jupiter[#1,2]})}
 
%% == Calcul de la séparation angulaire
\def\separation#1{
  \xintdeffloatvar ja,jd := \jupiter[#1,1],\jupiter[#1,2];%
  \xintdeffloatvar sa,sd := \saturn[#1,1],\saturn[#1,2];%
  \xintdeffloatvar cd := sind(jd)sind(sd)+cosd(jd)cosd(sd)cosd(sa-ja);%
  \xintdeffloatvar sep := acosd(cd);%
}
 
 
\renewcommand{\pshlabel}[1]{\footnotesize\color{ForestGreen}#1}
\renewcommand{\psvlabel}[1]{\footnotesize\color{ForestGreen}#1}
 
\begin{document}
 
\def\Graphe{%
  \begin{pspicture}(0,0)(4,3)
  \psset{xunit=0.1cm,yunit=1cm,linecolor=Black}
  \psaxes[Dy=0.5,Ox=1,Dx=7,ticks=none,labelsep=-2pt]{->}(1,0)(40,3)
    \multido{\i=1+1}{40}{
      % Séparation en ascension droite
      \psdots[linecolor=Blue,dotsize=1.4pt](\i,\ftf8{abs(\jupiter[\i,1]-\saturn[\i,1])})
      % Séparartion en déclinaison
      \psdots[linecolor=Green,dotsize=1.4pt](\i,\ftf8{abs(\jupiter[\i,2]-\saturn[\i,2])})
      % Séparation réelle
      \separation{\i}\psdots[linecolor=Tomato,dotsize=1.4pt](\i,\ftf8{sep})
    }
  \end{pspicture}
}
 
\begin{pspicture}(16,-16)(24,-8)
%% Ciel de Paris le 21 décembre 2021 à 18h21m TU
\psset[caelum]{
  scale=4,
  lon=-2.3,lat=48.805,
  year=2020,month=12,day=21,hour=18.35}
\caelum{Sky}
\caelum[linecolor=MidnightBlue!90!LightBlue,linewidth=0.5pt]{Graticule}%
\caelum[linewidth=0.2pt,linecolor=DimGray]{Boundaries}
\caelum[magmax=8,fillcolor=White]{Stars}
\caelum[linewidth=0.3pt,linecolor=Crimson]{Ecliptic}
\psdots[linecolor=Yellow,dotsize=6pt](\Sun)
\multido{\i=1+1}{40}{\Saturn{\i}\Jupiter{\i}}
\rput[c](21.4,-13.9){\psframebox[%
   framesep=12pt,fillstyle=solid,fillcolor=Gainsboro,linecolor=DimGray]{\Graphe}}
\end{pspicture}
\end{document}