; Specifically for date='2009:03:25_09h:17m:21s'
pro fdoy_date,date,fdoy,year

    year  = fix(strmid(date,0,4))
    month = fix(strmid(date,5,2))
    day   = fix(strmid(date,8,2))
    hh=FIX(STRMID(date,11,2))
    mm=FIX(STRMID(date,15,2))
    ss=FIX(STRMID(date,19,2))
    dmy2doy,year,month,day,doy
    fdoy=doy+hh/24.+mm/(24.*60)+ss/86400.

end