Hi Edward, The Jacobians are calculated in subroutine glmapm1 in coef.f and stored in array JACM1. Jacobian for each GLL point is essential calculated as the determinant of the Jacabi matrix which accounts for coordinates transformation. Ketan