Higher-order finite volume differential operators with selective upwinding on the icosahedral spherical grid

These operators compute the cell-average divergence, curl, and gradient with accuracy between second and fourth order for test functions. When applied to shallow water test cases, the resulting method is between third and fourth order accurate for the transport equation and between second and third order accurate for the full shallow water system.