diff --git a/Tools/Py_util/amrex_particles_to_vtp/write_pview_file.py b/Tools/Py_util/amrex_particles_to_vtp/write_pview_file.py
index 253b862db1a..57264745432 100644
--- a/Tools/Py_util/amrex_particles_to_vtp/write_pview_file.py
+++ b/Tools/Py_util/amrex_particles_to_vtp/write_pview_file.py
@@ -1,4 +1,6 @@
 import numpy as np
+import struct
+import os
 
 def write_paraview_file_structmesh(fname,xp,yp,ccdata,ncdata):
 
@@ -184,7 +186,7 @@ def write_paraview_file_particles(fname,pts,ncdata):
     Npts=pts.shape[0]
 
     outfile.write("<?xml version=\"1.0\"?>\n")
-    outfile.write("<VTKFile type=\"PolyData\" version=\"0.1\" byte_order=\"LittleEndian\">\n")
+    outfile.write("<VTKFile type=\"PolyData\" version=\"0.1\" byte_order=\"LittleEndian\" header_type=\"UInt64\">\n")
     outfile.write("<PolyData>\n")
     outfile.write("<Piece NumberOfPoints=\"%d\" NumberOfVerts=\"0\" NumberOfLines=\"0\" NumberOfStrips=\"0\" NumberOfPolys=\"0\">\n"%(Npts))
 
@@ -199,14 +201,25 @@ def write_paraview_file_particles(fname,pts,ncdata):
     outfile.write("</PointData>\n")
 
     outfile.write("<Points>\n")
-    outfile.write("<DataArray type=\"Float32\" Name=\"Points\" NumberOfComponents=\"3\" format=\"ascii\">\n")
-    for i in range(Npts):
-        outfile.write("%e\t%e\t%e\t"%(pts[i][0],pts[i][1],pts[i][2]))
-    outfile.write("\n</DataArray>\n")
-    outfile.write("</Points>\n")
-
+    outfile.write("<DataArray type=\"Float32\" Name=\"Points\" NumberOfComponents=\"3\" format=\"appended\" offset=\"0\">\n")
+    outfile.write("</DataArray>\n")
+    outfile.write("</Points>\n")    
     outfile.write("</Piece>\n")
     outfile.write("</PolyData>\n")
-    outfile.write("</VTKFile>\n")
 
+    # write AppendedData in special format
+    # https://vtk.org/Wiki/VTK_XML_Formats#Uncompressed_Data
+    outfile.write("<AppendedData encoding=\"raw\">\n")
+    outfile.write("\t_")
+    outfile.close()
+    with open(fname, 'ab') as binfile:
+        arrays_to_write = [pts]
+        for arr in arrays_to_write:
+            float_arr = arr.astype(np.float32)
+            binfile.write(struct.pack('<L', arr.nbytes))
+            binfile.write(float_arr.tobytes('C'))
+
+    outfile = open(fname, 'a')
+    outfile.write("\n</AppendedData>")
+    outfile.write("</VTKFile>\n")
     outfile.close()