yearday

 YEARDAY: converts calender month and day into yearday.
 yd = YEARDAY(mon,day,leapyr) converts calender month and day into yearday.
 If year is not a leap year, then omit the third input variable. If year
 is a leap year, then enter leapyr=1.

 INPUT:   mon    - month
          day    - day
          leapyr - set to 1 if it is a leap year, otherwise omit

 OUTPUT:  yd - year day

   Examples:  15 March 1995 = yearday(3,15) = 74
              15 March 1996 = yearday(3,15,1) = 75

0001 function yd=yearday(mon,day,leapyr)
0002 % YEARDAY: converts calender month and day into yearday.
0003 % yd = YEARDAY(mon,day,leapyr) converts calender month and day into yearday.
0004 % If year is not a leap year, then omit the third input variable. If year
0005 % is a leap year, then enter leapyr=1.
0006 %
0007 % INPUT:   mon    - month
0008 %          day    - day
0009 %          leapyr - set to 1 if it is a leap year, otherwise omit
0010 %
0011 % OUTPUT:  yd - year day
0012 %
0013 %   Examples:  15 March 1995 = yearday(3,15) = 74
0014 %              15 March 1996 = yearday(3,15,1) = 75
0015 
0016 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0017 % 8/19/98: version 1.1
0018 % 8/5/99: version 2.0
0019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0020 
0021 m=mon;
0022 d=day;
0023 
0024 % compute yd for non-leap year
0025 if m==1,y=d;end;
0026 if m==2,y=31+d;end;
0027 if m==3,y=59+d;end;
0028 if m==4,y=90+d;end;
0029 if m==5,y=120+d;end;
0030 if m==6,y=151+d;end;
0031 if m==7,y=181+d;end;
0032 if m==8,y=212+d;end;
0033 if m==9,y=243+d;end;
0034 if m==10,y=273+d;end;
0035 if m==11,y=304+d;end;
0036 if m==12,y=334+d;end;
0037 end
0038 yd=y;
0039 % adjust for leap year
0040 if nargin > 2,
0041   if (leapyr==1 & m>=3)
0042      yd=yd+1;
0043   end 
0044 end
0045 
0046 
0047