;********************************************
;  Subroutine to find the number of days in the month
;********************************************

pro find_number_of_days,year,month,numdays

;  Check if the year is a leap year
leap_flag=0  ;not a leap year
check4=(year) mod 4  ;years divisible by 4 are leap years, unless...........
check100=(year) mod 100  ;years also divisible by  100 are not leap years, except.......
check400=(year)mod 400  ;years divisible by 400 are leap years.
if (  (check4 eq 0 and check100 ne 0) or (check400 eq 0) ) then leap_flag=1  ;leap year

;  Make an array of the number of days
if leap_flag eq 0 then begin
           ; J  F  M  A  M  J  J  A  S  O  N  D
  num_days=[31,28,31,30,31,30,31,31,30,31,30,31]
endif else if leap_flag eq 1 then begin
  num_days=[31,29,31,30,31,30,31,31,30,31,30,31]
endif

;  This is the number of days in the month
numdays=num_days[month-1]

return
end