//input is date in form mmddyyyy; output is julian day number & fraction 
  function toJul(dd, mm, yyyy) { 
    if ( mm==1 || mm==2 ) { 
      yyyyp = yyyy - 1;  
      mmp = mm + 12; 
    } else { 
      yyyyp=yyyy; 
      mmp=mm; 
    } 
 
    yyyymmdd=(yyyy*10000)+(mm*100)+dd; 
    if(yyyymmdd >=15821015 ) { 
      a = Math.floor(yyyyp/100); 
      b = 2 - a + Math.floor(a/4); } 
    else b=0; 
 
    t2 = 365.25 * yyyyp; 
    if(yyyyp<0) 
      t2 = t2 - 0.75; 
 
    c = Math.floor(t2); 
    t1 = 30.6001 * (mmp + 1); 
    d = Math.floor(t1); 
    return b+c+d+dd+1720994.5; 
  }