fig02b.tex

\documentclass[svgnames]{standalone}
\usepackage[T1]{fontenc}
\usepackage[frenchb]{babel}
\usepackage{bera}
\usepackage{pstricks-add}
\usepackage{xcolor}
\usepackage{pst-caelum}
\usepackage{readarray}
\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\RASun{\FloatToFix{xmRAS(T)}} % Ascension droite du Soleil le 10/6/21 à 11h00m58s
\def\DESun{\FloatToFix{xmDS(T)}} % Déclinaison du soleil le 10/6/21 à 11h00m58s 
 
\def\CoinDeCiel(#1,#2,#3,#4,#5){%
  \xintdeffloatvar d1,d2 := -(#2-90),-(#4-90);%
  \xintdeffloatvar aa := (#3-#1)/2;%
  \xintdeffloatvar rot := (#3+#1)/2;%
  \xintdeffloatvar ax,ay := -sind(aa)d1,-d1;%
  \xintdeffloatvar bx,by := sind(aa)d1,-cosd(aa)d2;%
  \global\edef\ax{\FloatToFix{#5ax/20}}%
  \global\edef\ay{\FloatToFix{#5ay/20}}%
  \global\edef\bx{\FloatToFix{#5bx/20}}%
  \global\edef\by{\FloatToFix{#5by/20}}%
  \global\edef\axm{\FloatToFix{(#5ax/20)-0.5}}%
  \global\edef\aym{\FloatToFix{(#5ay/20)-0.5}}%
  \global\edef\bxm{\FloatToFix{(#5bx/20)+0.5}}%
  \global\edef\bym{\FloatToFix{(#5by/20)+0.5}}%
  \global\edef\axr{\FloatToFix{#5(cosd(rot)ax+sind(rot)ay)/20}}%
  \global\edef\ayr{\FloatToFix{#5(-sind(rot)ax+cosd(rot)ay)/20}}%
  \global\edef\cxr{\FloatToFix{#5(cosd(rot)ax+sind(rot)by)/20}}%
  \global\edef\cyr{\FloatToFix{#5(-sind(rot)ax+cosd(rot)by)/20}}%
  \global\edef\bxr{\FloatToFix{#5(cosd(rot)bx+sind(rot)by)/20}}%
  \global\edef\byr{\FloatToFix{#5(-sind(rot)bx+cosd(rot)by)/20}}%
  \global\edef\dxr{\FloatToFix{#5(cosd(rot)bx+sind(rot)ay)/20}}%
  \global\edef\dyr{\FloatToFix{#5(-sind(rot)bx+cosd(rot)ay)/20}}%
  \global\edef\rotation{\FloatToFix{rot}}
  \global\edef\scale{#5}
}
\def\GetVariables {
  \xintdeffloatvar ax,ay,bx,by,rot,scale := \ax,\ay,\bx,\by,\rotation,\scale;%
}
\def\ObjetDuCiel(#1,#2){
 \xintdeffloatvar pa := -(#1+90)+rot;%
 \xintdeffloatvar pr := -(#2-90)scale/20.0;%
 \xintdeffloatvar px,py := cosd(pa)*pr,sind(pa)*pr;%
}
 
\AtBeginDocument{\sffamily\boldmath}
\begin{document}
 
\CoinDeCiel(\ftf5{\RASun-30},\ftf5{\DESun-20},\ftf5{\RASun+30},\ftf5{\DESun+20},4)
% ==============================================================================
\begin{pspicture}(\axm,\aym)(\bxm,\bym)
\psrotate(0,0){\rotation}{%
  \begin{psclip}{\psline[linecolor=MidnightBlue](\axr,\ayr)(\cxr,\cyr)(\bxr,\byr)(\dxr,\dyr)(\axr,\ayr)}%
    \psset[caelum]{scale=\scale,linewidth=0.1pt}%
    %% -------------------------------------------------------------------------
    %% Construction du Ciel
    %% Envisager la commande \CaelumTheme{<nom du thème>}
    %% -------------------------------------------------------------------------
    \caelum{Sky}%
    \caelum[linestyle=dotted,linecolor=DimGray]{Boundaries}%
    % \caelum[abrev=lyr,fillcolor=DodgerBlue]{Constellation}
    \caelum[fillcolor=LightGoldenrodYellow]{Stars}%
    \caelum[linecolor=Crimson,linewidth=0.2pt]{Ecliptic}%
    \caelum[linecolor=LightSkyBlue,fillcolor=Goldenrod]{Asterisms}%
    %% -------------------------------------------------------------------------
    \psdots[linecolor=Yellow,dotsize=12pt]%
       (\ftf5{-(\DESun-90)\scale/20};\ftf5{-(\RASun+90)})%
    \multido{\i=1+6}{20}{%
      \def\dec{\lune[\i,2]}
      \def\asc{\lune[\i,1]}
      \def\RPMoon{\ftf5{-(\dec-90)\scale/20}}
      \def\APMoon{\ftf5{-(\asc+90)}}
      \psdots[linecolor=DimGray,dotsize=12pt](\RPMoon;\APMoon)
    }%
    %% -------------------------------------------------------------------------
  \end{psclip}%
}%
%% -----------------------------------------------------------------------------
%% Objets à placer
%% -----------------------------------------------------------------------------
% \GetVariables%                       % Retrouver les variables du Coin de Ciel
%% -----------------------------------------------------------------------------
\rput[l](\ax,\ftf5{\ay-0.25}){Éclipse du 10 juin 2021 (conjonction à 11\textsuperscript{h}00'58'' TU)}
\rput[r](\bx,\ftf5{\ay-0.25}){\textbf{PST-Caelum}}
\end{pspicture}
\end{document}