secondsInYear = ndaysnoleap * secondsperday
days_in_year=ndaysnoleap
ndays = (ndaysnoleap - 1) + cal_IsLeap( yi, myThid ) -
ndays = ndays + (ndaysnoleap - 1) +