fig02c.tex

\documentclass[svgnames]{standalone}
\usepackage[T1]{fontenc}
\usepackage[frenchb]{babel}
\usepackage{bera}
\usepackage{pstricks-add}
\usepackage{xcolor}
\usepackage[moon,parallaxe]{xint-meeus}
\usepackage{pst-caelum}
 
\usepackage{readarray}
\readdef{lune02.txt}\lignes%
\readarray{\lignes}\lune[120,2]%
 
 
%% =============================================================================
\xintdeffloatvar lat  := 80.815;% 
\xintdeffloatvar lon  := 66.805;%
\xintdeffloatvar alt  := 0.0;%
\JulianDay(2021,6,10,10,41,51)%
\xintdeffloatvar ras,des := xmRAS(TS),xmDS(TS);%
\FCoordsMoon{para}{parac}
\xintdeffloatvar DTL  := pararv/149597870.7;%
\xintdeffloatvar dt := 10*15+41/4+51/240;%
\JulianDay(2021,6,10,0,0,0)%
\xintdeffloatvar ahs := STG0(TS,dt)-lon-ras;%
\xintdeffloatvar rasa := xmPNRA(ras,des,1,lat,ahs,alt);%
\xintdeffloatvar desa := xmPNDE(ras,des,1,lat,ahs,alt);%
\def\CA#1{%
  \xintdeffloatvar ram,dem := \lune[#1,1],\lune[#1,2];%
  \xintdeffloatvar ahm := STG0(TS,dt+15(#1-61))-lon-ram;% 
  \xintdeffloatvar rama := xmPNRA(ram,dem,DTL,lat,ahm,alt);%
  \xintdeffloatvar dema := xmPNDE(ram,dem,DTL,lat,ahm,alt);%
}%
%% ==============================================================================
\def\RASun{\FloatToFix{rasa}}% Ascension droite du Soleil
\def\DESun{\FloatToFix{desa}}% Déclinaison du Soleil 
 
\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{(#5)ax/20}}%
  \global\edef\ay{\FloatToFix{(#5)ay/20}}%
  \global\edef\bx{\FloatToFix{(#5)bx/20}}%
  \global\edef\by{\FloatToFix{(#5)by/20}}%
  \xdef\axm{\FloatToFix{(#5)ax/20-0.5}}%
  \xdef\aym{\FloatToFix{(#5)ay/20-0.5}}%
  \xdef\bxm{\FloatToFix{(#5)bx/20+0.5}}%
  \xdef\bym{\FloatToFix{(#5)by/20+0.5}}%
  \xdef\rotation{\FloatToFix{rot}}%
  \xdef\scale{#5}%
}
\def\GetVariables {
  \xintdeffloatvar ax,ay,bx,by,rot,scale := \ax,\ay,\bx,\by,\rotation,\scale;%
}
 
\AtBeginDocument{\sffamily\boldmath}
\begin{document}
 
\CoinDeCiel(\ftf8{\RASun-30},\ftf8{\DESun-20},\ftf8{\RASun+30},\ftf8{\DESun+20},4)
% ==============================================================================
\psset[caelum]{scale=\scale,linewidth=0.1pt,rotate=\rotation}%
\begin{pspicture}(\axm,\aym)(\bxm,\bym)%
\begin{psclip}{\psframe[linecolor=MidnightBlue,fillstyle=solid,fillcolor=black](\ax,\ay)(\bx,\by)}%
  \caelum{Sky}%
  \caelum[linestyle=dotted,linecolor=DimGray]{Boundaries}%
  \caelum[fillcolor=LightGoldenrodYellow]{Stars}%
  \caelum[linecolor=Crimson,linewidth=0.2pt]{Ecliptic}%
  \caelum[linecolor=LightSkyBlue,fillcolor=Goldenrod]{Asterisms}%
  %% ---------------------------------------------------------------------------
  \psdots[linecolor=Yellow,dotsize=12pt](\ftf8{-(\DESun-90)\scale/20};\ftf8{-(\RASun+90)+rot})%
  \multido{\i=1+6}{20}{\CA{\i}%
  \psdots[linecolor=DimGray,dotsize=11pt](\ftf8{-(dema-90)\scale/20};\ftf8{-(rama+90)+rot})}%
  %% ---------------------------------------------------------------------------
\end{psclip}%
\rput[l](\ax,\ftf8{\ay-0.25}){Éclipse du 10 juin 2021 à 10\textsuperscript{h}41'51'' TU}%
\rput[r](\ftf8{\bx-0.1},\ftf8{\ay+0.25}){\color{White}\textbf{PST-Caelum}}%
\rput[r](\bx,\ftf8{\ay-0.25}){Lieu d'observation: Nord de la Mer de Baffin (Arctique)}%
\end{pspicture}%
\end{document}