General relativistic MHD large eddy simulations with gradient subgrid-scale model