Physical processes in the ocean span multiple spatial and temporal scales. Simultaneously resolving all these in a simulation is computationally challenging. Here we develop a more efficient technique to better study the interactions across scales, particularly focusing on the ocean surface turbulent mixing, by coupling a global ocean circulation model MPAS-Ocean and a large eddy simulation model PALM. The latter is customized and ported on a GPU to further accelerate the computation.