add save velocity

py/modules/simulations/ctree.py

@@ -200,6 +200,7 @@ class CTree:
              v=v[shuffle_idx_threshold, ]
              v=np.reshape(v,(-1,))
         cxv=CSnapshot(None).center(p,v,m*rho)
+
         if v is None: # no velocities
             cxv[3]=0.0
             cxv[4]=0.0
@@ -208,6 +209,8 @@ class CTree:
             unso=CunsOut(outfile,"nemo")
             unso.setValueF("time",self.__time)
             unso.setArrayF("all","pos",p)
+            if self.__vel is not None:
+                unso.setArrayF("all","vel",v)
             unso.setArrayF("all","mass",m)
             unso.setArrayF("all","rho",rho)
             unso.setArrayF("all","hsml",hsml)