void sod(float *xx, float *vx, float *yy, float *vy, float *zz,
	 float *vz, float *mass, float *den, float *p, float *u,
	 float *hh, float gamma, float xmax, float xmin, int
	 *partno, float hfac, int n, int ientro);