function S=NouvLune(Annee_calc,format)
% NL=NOUVELLELUNE(ANNEE,FORMAT) calcule les dates de la nouvelle lune
% pour ANNEE. FORMAT est 'DATE' pour la date en clair ou 'LONG'
% pour la date julienne.

if nargin==1
	format='DATE';
end

% Constantes
d1=1236.853086;
d0=0.827361;

% initialisation
Lunaison_0=fix(d1*(Annee_calc-2000)/100);
S=[];
for Lunaison_i=Lunaison_0:Lunaison_0+13
	T_Nouv_Lune=(Lunaison_i-d0)/d1;
	T_Nouv_Lune=ameliore(T_Nouv_Lune);
	T_Nouv_Lune=ameliore(T_Nouv_Lune);
	mjd_nouv_lune=36525*T_Nouv_Lune+51544.5;
	[jour,mois,annee,heure]=caldat(mjd_nouv_lune);
	minutes=fix(frac(heure)*60);
	heure=fix(heure);
	if annee==Annee_calc
		if upper(format)=='DATE'
			S=[S; sprintf('%2i/%2i/%4i %2i:%2i',jour,mois,annee,heure,minutes)];
		else 
			S=[S;mjd_nouv_lune];
		end
	end
end
