function [eldia,elmes]=quefecha_fun(dia,any)




%disp('Calcula la data respecte al dia-del-any')
%any=input('Any?=');
%dia=input('Dia de l.any?=');





ndias=[31,28,31,30,31,30,31,31,30,31,30,31,0];

%if any==1948|any==1952|any==1956|any==1960|any==1964|any==1968|any==1972|any==1976|any==1980|any==1984|any==1988|any==1992|any==1996|any==2000|any==2004|any==2008|any==2012|any==2016|any==2020
anyobisiesto=[1936:4:2100];
clear si,si=find(any==anyobisiesto);
if length(si)==1  %Anyo bisiesto

ndias=[31,29,31,30,31,30,31,31,30,31,30,31,0];

end



if dia>sum(ndias),disp('Wrong Julian Day'),return,end


k=1;

if dia<ndias(1:1),
	elmes=1;
	eldia=dia;
else

while dia>=sum(ndias(1:k)) & k<=12;
clear elmes,clear eldia
elmes=k+1;
eldia=dia-sum(ndias(1:k));
k=k+1;
end

end %if dia<ndias(1:1)








if eldia==0, elmes=elmes-1;eldia=ndias(elmes);end

if isempty(elmes)==1 &dia<=31; eldia=dia;elmes=1;end









%if isempty(elmes)==1,elmes=1;eldia=dia;end







%disp(['dia=',int2str(eldia)])

%disp(['mes=',int2str(elmes)])



