Failure to account for periodicity in set! on TripolarGrid

[glwagner]

Consider

using Oceananigans
using GLMakie

Nx = Ny = Nz = 128
Nz = 16
grid = TripolarGrid(size=(Nx, Ny, Nz), z=(0, 1))
c = CenterField(grid)
cᵢ(λ, φ, z) = exp(-(λ^2 + φ^2) / 200)
set!(c, cᵢ)
heatmap(interior(c, :, :, 1))

produces

I'm not entirely sure what the solution is. The problem is that one has to know that domain the function is defined on / what direction we should apply periodicity.