I need to get pressure of vertex in the 'NEWXYZ' subroutine, and use this to caculate the vertex movement.But I don't know how to get the vertex's pressure value in the 'NEWXYZ' subroutine.Please help me!

thanks Henry

You can pass cell center pressure to newxyz through posdat.f: creat a common block of your own, say:

common/user1/puser(10000)

(you can change 10000 into any number of cells).

Then this common block can be moved to newxyz.

Similar you can calculate vertex pressure in posdat.

What is name of the varible which means vertex pressure in the posdat.f? thanks