Hey,


I would like to run a POD on a section of my 3D field on a coarser grid. For this, I would like the section to be equidistributed cartesian box and this is exactly what the '-m interppoints:box' utility does. However, I'm having trouble running the utility in parallel. When I run the command:

FieldConvert -m interppoints:fromxml=base/jet.xml:fromxml=base/mesh.xml:fromfld=base/mesh_1.chk:\
box=10,10,10,-0.001,0.001,-0.001,0.001,-0.003,0.003 box1.dat

the utility does exactly what I want.


When I run it in parallel:

mpirun -np 3 FieldConvert \
-m interppoints:fromxml=base/jet.xml:fromxml=base/mesh.xml:fromfld=base/mesh_1.chk:\
box=10,10,10,-0.001,0.001,-0.001,0.001,-0.003,0.003 box2.dat


box2.dat only contains one third of the points in the box. I also tried adding box2.dat:dat:writemultiplefiles, in which case I get the correct number of points split in three different files, but the file is still somehow corrupted when I open it with paraview. What am I doing wrong?

Petteri Peltonen, (M.Sc.)
Aalto University, School of Engineering
Department of Mechanical Engineering
PO Box 14300
FI-00076 Aalto, Finland
+358 44 0554984
petteri.peltonen@aalto.fi

_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users