3 proj_var2d(xDim, yDim, zDim, "data", "fieldz") 4 proj_2d(xDim, yDim, zDim,"data","wfc",50000) 5 #proj_2d(xDim, yDim, zDim,"data","wfc",1000000) 6 #proj_k2d(xDim, yDim, zDim,"data","wfc",500000) 8 item_wfc = wfc_density(xDim, yDim, zDim,"data","wfc",50000) 9 #item_wfc = wfc_density(xDim, yDim, zDim,"data","wfc",1000000) 10 item_phase = wfc_phase(xDim, yDim, zDim,"data","wfc",50000) 11 #item_phase = wfc_phase(xDim, yDim, zDim,"data","wfc",1000000) 12 #item_var = var(xDim, yDim, zDim,"data","By_0") 13 #item2_var = var(xDim, yDim, zDim,"data","Edges_0") 15 #to_bvox(item_wfc, xDim, yDim, zDim, 1, "test_wfc.bvox") 16 #to_bvox(item_phase, xDim, yDim, zDim, 1, "test_phase.bvox") 17 #to_bvox(item_var, xDim, yDim, zDim, 1, "test_var.bvox") 18 #to_bvox(item2_var, xDim, yDim, zDim, 1, "test_edges.bvox") 22 print("iterating through data...") 23 # Writing out the potential 24 #item_pot = var(xDim, yDim, zDim, "data", "V_0") 25 #to_bvox(item_pot, xDim, yDim, zDim, 1, "test_pot.bvox") 26 #item_gauge = var(xDim, yDim, zDim, "data", "Az_0") 27 #to_bvox(item_gauge, xDim, yDim, zDim, 1, "test_gauge.bvox") 28 #item_edges = var(xDim, yDim, zDim,"data","Edges_0") 29 #to_bvox(item_edges, xDim, yDim, zDim, 1, "test_edges.bvox") 30 #to_vtk(item_edges, xDim, yDim, zDim, 1, "test_edges.vtk") 31 for i in range(0,101): 32 #for i in [86, 143, 195, 95, 152, 204]: 35 proj_2d(xDim, yDim, zDim, "data", "wfc_ev", num) 36 #item = wfc_density(xDim, yDim, zDim, "data", "wfc_ev", num) 37 #item_ph = wfc_phase(xDim, yDim, zDim, "data", "wfc", num) 38 #to_bvox(item, xDim, yDim, zDim, 1, "wfc_%s.bvox" %num) 39 #to_vtk(item, xDim, yDim, zDim, 1, "wfc_%s.vtk" %num) 40 #to_bvox(item_ph, xDim, yDim, zDim, 1, "wfc_ph_%s.bvox" %num) 48 item = var(xDim, yDim, zDim, "data", "Edges_%s"%num) 49 to_bvox(item, xDim, yDim, zDim, 1, "EDGES_%s.bvox" %num) 53 proj_var2d(xDim, yDim, zDim, "data", "fieldz") 54 item = var(xDim, yDim, zDim, "data", "fieldz") 55 to_bvox(item, xDim, yDim, zDim, 1, "HE11check.bvox") 59 proj_2d(xDim, yDim, zDim, "data", "wfc_ramp", 100000) 60 #item_edges = var(xDim, yDim, zDim,"data","Edges_0") 61 #to_bvox(item_edges, xDim, yDim, zDim, 1, "test_edges.bvox") 66 item_edges = var(xDim, yDim, zDim,"data","Edges_std") 67 to_bvox(item_edges, xDim, yDim, zDim, 1, "edges_std.bvox") 68 to_vtk(item_edges, xDim, yDim, zDim, 1, "edges_std.vtk") 69 print("starting linpol") 70 item_edges = var(xDim, yDim, zDim,"data","Edges_linpol") 71 to_bvox(item_edges, xDim, yDim, zDim, 1, "edges_linpol.bvox") 72 to_vtk(item_edges, xDim, yDim, zDim, 1, "edges_linpol.vtk") 76 item = wfc_density(xDim, yDim, zDim, "data", "wfc", 100000) 77 to_vtk(item, xDim, yDim, zDim, 1, "wfc.vtk") 78 to_bvox(item, xDim, yDim, zDim, 1, "wfc_check.bvox") 92 item_gx = var(xDim, yDim, zDim,"data","Ax_0") 93 to_bvox(item_gx, xDim, yDim, zDim, 1, "Ax_0.bvox") 94 item_gy = var(xDim, yDim, zDim,"data","Ay_0") 95 to_bvox(item_gy, xDim, yDim, zDim, 1, "Ay_0.bvox") 96 item_gz = var(xDim, yDim, zDim,"data","Az_0") 97 to_bvox(item_gz, xDim, yDim, zDim, 1, "Az_0.bvox") 98 item_V = var(xDim, yDim, zDim,"data","V_0") 99 to_bvox(item_V, xDim, yDim, zDim, 1, "V_0.bvox") 100 item_K = var(xDim, yDim, zDim,"data","K_0") 101 to_bvox(item_K, xDim, yDim, zDim, 1, "K_0.bvox") 104 for i in range(0,11): 107 item = var(xDim, yDim, zDim, "data", "Edges_%s" %num) 108 to_bvox(item, xDim, yDim, zDim, 1, "EDGES_%s.bvox" %num) 113 proj_2d(xDim, yDim, zDim, "data", "wfc", num) 114 item = wfc_density(xDim, yDim, zDim, "data", "wfc", num) 115 item_edges = var(xDim, yDim, zDim, "data", "Edges_0") 116 to_bvox(item, xDim, yDim, zDim, 1, "wfc_%s.bvox" %num) 117 to_bvox(item_edges, xDim, yDim, zDim, 1, "Edges_0.bvox") 118 to_vtk(item, xDim, yDim, zDim, 1, "wfc_%s.vtk" %num) 122 proj_var2d(xDim, yDim, zDim, "data", "V_0") 123 #proj_2d(xDim, yDim, zDim, "data", "wfc_ev", 0) 126 item = var(xDim, yDim, zDim,"data","Bx_0") 127 to_vtk(item, xDim, yDim, zDim, 1, "test_Bx.vtk") 128 item = var(xDim, yDim, zDim,"data","By_0") 129 to_vtk(item, xDim, yDim, zDim, 1, "test_By.vtk") 130 item = var(xDim, yDim, zDim,"data","Bz_0") 131 to_vtk(item, xDim, yDim, zDim, 1, "test_Bz.vtk") 132 item = var(xDim, yDim, zDim,"data","Br_0") 133 to_vtk(item, xDim, yDim, zDim, 1, "test_Br.vtk") 134 item = var(xDim, yDim, zDim,"data","Bphi_0") 135 to_vtk(item, xDim, yDim, zDim, 1, "test_Bphi.vtk") 136 #item = var(xDim, yDim, zDim,"data","Edges_0") 137 #to_bvox(item, xDim, yDim, zDim, 1, "test_wfc.bvox") 139 #comx, comy = wfc_com(xDim, yDim, zDim, "data", "wfc", 0) 144 item =
var(xDim, yDim, zDim,
"data_lin",
"Edges_0")
145 to_vtk(item, xDim, yDim, zDim, 1,
"Edges.vtk")
def to_vtk(item, xDim, yDim, zDim, nframes, filename)