obsids_spirephot_txt = simpleAsciiTableReader(file="/Users/bmerin/0_demo/obsids_spirephot.txt")

outdir="/Users/bmerin/0_demo/Results/HIPE6/"

for i in range(len(obsids_spirephot_txt["c0"].data)):
     obsId = obsids_spirephot_txt["c0"].data[i]
     poolname = obsids_spirephot_txt["c1"].data[i]
     print "Reprocessing ",obsId," from pool ",poolname
     print "======================================================="
     # Retrieve the observation from a pool
     obs = getObservation(obsId,poolName=poolname,instrument="SPIRE")
     # reprocessing with current HIPE
     #spirePhotLargeMapPipeline(obs=obsId,mapping="madmap",level="level1")
     # Extract the 250 micron map
     level2_map_psw = obs.refs["level2"].product.refs["PSW"].product
     #save the map as a FITS file :
     outfile = outdir+ "map_"+ str(obsId) + "_PSW.fits"
     print "Saving file: " + outfile
     simpleFitsWriter(level2_map_psw,outfile)  
     # Extract the 350 micron map
     level2_map_pmw = obs.refs["level2"].product.refs["PMW"].product
     #save the map as a FITS file :
     outfile = outdir+ "map_"+ str(obsId) + "_PMW.fits"
     print "Saving file: " + outfile
     simpleFitsWriter(level2_map_pmw,outfile)  
     # Extract the 500 micron map
     level2_map_plw = obs.refs["level2"].product.refs["PLW"].product
     #save the map as a FITS file :
     outfile = outdir+ "map_"+ str(obsId) + "_PLW.fits"
     print "Saving file: " + outfile
     simpleFitsWriter(level2_map_plw,outfile)  
pass