import os

convertvtk2vti = "/home/dfuentes/DDDAS/trunk/utilities/WriteVTIFile"
#joblist = os.listdir( "Processed" )
#joblist.sort()

jobid = "dog%d_980"
for canine,dt,nstep in [ (1,6.0,119),(2,6.0,122),(3,6.0,74),(4,5.0,61)]:
  print canine,dt,nstep
  for filebase,varname in [("temperature","mrti"),("tmapstd","dmrti")]:
    vti=open("%s/%s00/mrivis/%s.pvd" % (jobid%canine,jobid%canine,filebase),"w")
    vti.write('<?xml version="1.0"?>\n')
    vti.write('<VTKFile type="Collection" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">\n')
    vti.write('  <Collection>\n')
    for i in range(nstep+1):
         vti.write('   <DataSet timestep="%f" part="0" file="./%s.%04d.vti"/>\n' % (i*dt,filebase,i) ) 
    vti.write('  </Collection>\n')
    vti.write('</VTKFile>\n')
    for i in range(nstep+1):
      template= "%s %s/%s00/mrivis/%s.%04d.vtk %s/%s00/mrivis/%s.%04d.vti %s" 
      cmd =  template % (convertvtk2vti, 
                         jobid%canine,jobid%canine,filebase,i,
                         jobid%canine,jobid%canine,filebase,i,
                         varname)
      print cmd
      os.system( cmd )