Alberto F. Martin amartinhuertas

## mwe_grad_perp.jl
using Gridap

# Distorted mesh
function rndm(p::Point)
  r, s = p
  x = r + 0.1*sin(2π*r)*sin(2π*s)
  y = s + 0.1*sin(2π*r)*sin(2π*s)
  Point(x,y)
end

## assemble_L2MM_invL2MM.jl
using Gridap
using GridapGeosciences

order=1
degree=4

model = CubedSphereDiscreteModel(1; radius=rₑ)


# Forward integration of the shallow water equations
	using Gridap

	# Distorted mesh
	function rndm(p::Point)
	r, s = p
	x = r + 0.1sin(2πr)sin(2πs)
	y = s + 0.1sin(2πr)sin(2πs)
	Point(x,y)
	end
	using Gridap
	using GridapGeosciences

	order=1
	degree=4

	model = CubedSphereDiscreteModel(1; radius=rₑ)


	# Forward integration of the shallow water equations