dirI='/obs/bcarry/data/mining/omegacam/' dirO= dirI+'skybot/' fileIN = 'listFOV-2013-10-19.dat' fileEX = 'progress.dat' fovBOX= '1.1x1.1' if file_exist( dirO+fileEX ) then begin readfmt, dirO+fileEX, $ '(F10.6,2x,F10.6,2x,A-20,2x,D16.8,2x,I2)', $ raArr, deArr, isoArr, jdArr, foundArr, $ /SILENT endif else begin readfmt, dirI+fileIN, 'A18,1x,I2,1x,I2,1x,F6.3,1x,'+$ 'A1,I2,1x,I2,1x,F5.2,1x,A22', $ name, rah,ram,ras, $ decpm, dech,decm,decs, isoArr, /SILENT nbDate = n_elements( isoArr ) foundArr = intarr( nbDate ) jdArr = dblarr( nbDate ) raArr = fltarr(nbDate) deArr = fltarr(nbDate) for kD=0, nbDate-1 do begin jdArr(kD) = date_conv( isoArr(kD), 'JULIAN' ) raArr(kD) = ten( rah(kd),ram(kD),ras(kd) )*15. deArr(kD) = ten( dech(kd),decm(kD),decs(kd) ) if decpm(kD) eq '-' then deArr(kD)*=-1. ; print, kD, rah(kd),ram(kD),ras(kd), raArr(kD) ; print, kD, dech(kd),decm(kD),decs(kd), deArr(kD) endfor endelse toSearch = where( foundArr eq 0, nbAst ) ; nbAst = 1000L jdStart = SYSTIME( /JULIAN, /UTC ) for kAst=0L, nbAst-1 do begin kCur = toSearch(kAst) jdArr(kCur) = date_conv( isoArr(kCur), 'JULIAN' ) print, '-- '+strtrim(string(kAst+1, format='(I6)'),2)+'/'+$ strtrim(string(nbAst , format='(I6)'),2) fileOut = dirO+isoArr(kCur)+'.xml' epoch = strtrim(isoArr(kCur),2) fovRA = strtrim(string(raArr(kCur),format='(F10.6)'),2) fovDEC= strtrim(string(deArr(kCur),format='(F10.6)'),2) spawn, 'wget "http://vo.imcce.fr/webservices/skybot/skybotconesearch_query.php?'+$ '-ep='+epoch+$ ;-Observing epoch '&-ra='+fovRA+$ ;-Center of the FoV RA '&-dec='+fovDEC+$ ;-Center of the FoV DEC '&-bd='+fovBOX+$ ;-Box of the FoV in degree '&-loc=568'+$ '&-mime=votable'+$ '&-objFilter=110'+$ '&-from=BCarry'+$ ;-User '" -O '+fileOut+' -q' testString = 'QUERY_STATUS' spawn, 'grep '+testString+' '+fileOut, output split1 = strsplit( output(0), '<=>', /EXTRACT ) status= strsplit(split1(3),'"/', /EXTRACT ) ;-- The request worked -- Proceed if status eq 'OK' then begin split2 = strsplit( output(1), '<=>', /EXTRACT ) result= strsplit(split2(2),'"/', /EXTRACT ) ;-- Object found or not if result(0) eq 'OK' then foundArr(kCur)= 1 $ else foundArr(kCur)=-1 endif else begin wait, 60. endelse endfor jdEnd = SYSTIME( /JULIAN, /UTC ) forprint, raArr, deArr, isoArr, jdArr, foundArr, $ format='(F10.6,2x,F10.6,2x,A-20,2x,D16.8,2x,I2)', $ textout=dirO+fileEX, /SILENT, /NOCOMMENT end