On modification of pressure gradient operator in integrated ISPH for multifluid and porous media flow with free-surface