In [11], the authors developed a class of high-order numerical schemes for the Hamilton-Jacobi (H-J) equations, which are unconditionally stable, yet take the form of an explicit scheme. This paper extends such schemes, so that they are more effective at capturing sharp gradients, especially on nonuniform meshes. In particular, we modify the weighted essentially non-oscillatory (WENO) methodology in the previously developed schemes by incorporating an exponential basis and adapting the previously developed nonlinear filters used to control oscillations. The main advantages of the proposed schemes are their effectiveness and simplicity, since they can be easily implemented on higher-dimensional nonuniform meshes. We perform numerical experiments on a collection of examples, including H-J equations with linear, nonlinear, convex and non-convex Hamiltonians. To demonstrate the flexibility of the proposed schemes, we also include test problems defined on non-trivial geometry.Recent work on the MOL T has involved extending the method to solve more general nonlinear PDEs, for which an integral solution is generally not applicable. This work includes the nonlinear degenerate convection-diffusion equations [10], as well as the H-J equations [11]. The key idea of these papers involved exploiting the linearity of a given differential operator, rather than requiring linearity in the underlying equations. This allowed derivative operators in the problems to be expressed through kernel representations developed for linear problems. Formulating applicable derivative operators in this way ultimately facilitated the stability of the schemes, since a global coupling was introduced through the integral operator. As part of this embedding process, a kernel parameter β was introduced, and through a careful selection, was shown to yield schemes which are A-stable. Remarkably, it was shown that one could couple these representations for the derivative operators with an explicit time-stepping method, such as the strong-stability-preserving Runge-Kutta (SSP-RK) methods [17] and still obtain schemes which maintain unconditional stability [10,11]. To address shock-capturing and control non-physical oscillations, the latter two papers introduced quadrature formulas based on WENO reconstruction, along with a nonlinear filter.This paper seeks to extend the work in [10,11] to the H-J equations (1.1) defined on non-uniformly distributed spatial domains. In particular, several improvements are given. First, we develop the MOL T for mapped grids using a general coordinate transformation function, which allows for a non-uniform distribution of grid points. We show that, with this mapping, our numerical scheme is able to preserve the conservation property for the derivative of the solution to the H-J equation. We also describe a novel WENO-based quadrature for the spatial discretization, which uses a basis consisting of exponential polynomials, to improve the shock capturing capabilities of the method. Another difference in this paper, compared to...