Discrete boundary

[phyo-wai-thaw]

Hi all,

I’m wondering if this is the correct way to define the discrete boundary conditions at the top.

The full code is below.
https://github.com/phyo-wai-thaw/melting/blob/3952cd908d546a3be76fdfbf45c09de25eca049f/ISOBL_coarse_for_github.ipynb

I also have a question about the data in the halo cells.

When I check the average temperature (Tg = interior(slice_time_series.Tg, :, :, :, 151)) using the slicer at grid.Nz+1, I get -2.2763659. However, when I check directly (Tt1 = avg_time_series.Tavg.data[:, :, grid.Nz+1, 151]) without using the slicer, the temperature at the first ghost cell is zero.

How can I extract the data of the halo cells? I believe the default JLD2 output writer stores halo data.

Many thanks,
Phyo Wai