Skip to content

Commit aa057e7

Browse files
Remove assumption of radius == EARTH_RADIUS in triangular icosahedral grid (#156)
* fix (only) exports * replace EARTH_RADIUS with specified radius in tri ico mesh
1 parent a36b7af commit aa057e7

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

src/meshes/dgmulti_icosahedron_tri.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -42,7 +42,7 @@ function DGMultiMeshTriIcosahedron2D(dg::DGMulti{NDIMS, <:Tri}, radius;
4242
vk = ntuple(n -> Vxyz_old[n][k], NDIMS_AMBIENT)
4343
vl = ntuple(n -> Vxyz_old[n][l], NDIMS_AMBIENT)
4444
midpoint = 0.5 .* (vk .+ vl)
45-
midpoint = midpoint .* (EARTH_RADIUS / norm(midpoint)) # Normalize to outer radius
45+
midpoint = midpoint .* (radius / norm(midpoint)) # Normalize to outer radius
4646
for n in 1:NDIMS_AMBIENT
4747
push!(Vxyz[n], midpoint[n])
4848
end

0 commit comments

Comments
 (0)