#!/bin/sh

cat << *** | gp -q 
read("saisons.gp");
read("format.gp");
read("trigo.gp");
a = TempsPerigeeSoleilM($1);
b = JD(1,1,$1,0,0,0);
print(format_date(CD(a+b)));
print(format_date(printemps($1)));
print(format_date(ete($1)));
print(format_date(automne($1)));
print(format_date(hiver($1)));
***