This subroutine calculates the volume fraction of phase i. The call to dassig is:
call dassig(nx,ny,nz,nphase,prob,pix)