Some bugs in the near global example

[simone-silvestri]

I think there is something wrong in the fluxes with the near-global example. I have ran it locally after it was not running in PR #137 and the results were a little odd.

The simulation is running but the extrema of SST drift extremely quickly, reaching -35 and +47 degree celsius in a month.
This is the video that is produced by the example that shows that temperature increases quite steeply in the tropics.

near_global_ocean_surface.mp4

The output is the following:

[ Info: In-painting ECCO temperature
[ Info: In-painting ECCO salinity
[ Info: Initializing simulation...
[ Info: Iter: 0, time: 0 seconds, Δt: 1.500 minutes, max|u|: (0.00e+00, 0.00e+00, 0.00e+00) m s⁻¹, extrema(T): (31.25, -1.96) ᵒC, wall time: 2.681 minutes
[ Info:     ... simulation initialization complete (6.395 seconds)
[ Info: Executing initial time step...
┌ Warning: Simulation stopped during initialization.
└ @ Oceananigans.Simulations ~/.julia/packages/Oceananigans/s1DfC/src/Simulations/run.jl:129
[ Info:     ... initial time step complete (3.408 minutes).
[ Info: Iter: 4800, time: 5 days, Δt: 1.500 minutes, max|u|: (1.73e+00, 2.26e+00, 1.87e-02) m s⁻¹, extrema(T): (34.14, -10.23) ᵒC, wall time: 40.442 minutes
[ Info: Simulation is stopping after running for 1.289 hours.
[ Info: Simulation time 10 days equals or exceeds stop time 10 days.
[ Info: Iter: 9600, time: 10 days, Δt: 1.500 minutes, max|u|: (2.07e+00, 2.05e+00, 9.71e-03) m s⁻¹, extrema(T): (37.57, -15.24) ᵒC, wall time: 36.879 minutes
[ Info: Initializing simulation...
[ Info:     ... simulation initialization complete (1.970 ms)
[ Info: Executing initial time step...
[ Info:     ... initial time step complete (3.960 seconds).
[ Info: Iter: 10320, time: 15 days, Δt: 10 minutes, max|u|: (1.94e+00, 1.97e+00, 1.61e-02) m s⁻¹, extrema(T): (41.27, -17.67) ᵒC, wall time: 5.780 minutes
[ Info: Iter: 11040, time: 20 days, Δt: 10 minutes, max|u|: (1.79e+00, 1.87e+00, 1.83e-02) m s⁻¹, extrema(T): (42.93, -25.24) ᵒC, wall time: 5.736 minutes
[ Info: Iter: 11760, time: 25 days, Δt: 10 minutes, max|u|: (1.73e+00, 2.56e+00, 1.75e-02) m s⁻¹, extrema(T): (45.31, -31.84) ᵒC, wall time: 5.722 minutes
[ Info: Iter: 12480, time: 30 days, Δt: 10 minutes, max|u|: (2.90e+00, 2.25e+00, 2.19e-02) m s⁻¹, extrema(T): (46.06, -35.93) ᵒC, wall time: 5.718 minutes

[navidcy]

Seems that after 5 days we have temperatures of -10. I thought there was a freezing temperature being set at -1.8, no?

[navidcy]

Are we sure that the reported values are not from "land" regions?

[navidcy]

I confirm that the subfreezing temperatures are in the Arctic. It's not possible to see with the plot above.

[navidcy]

What's going on with the speed plot??

[navidcy]

OK, the speed plot was just a visualization issue; is fixed in #239.

[glwagner]

Don't we have some validation of the fluxes? Can we re-run the validation?