c_grad_p replaced by velocity scaling term c_v and advection scaling term c_a
use default device