@@ -14,6 +14,7 @@ using ClimaOcean.ECCO
1414using ClimaOcean. JRA55
1515using ClimaOcean. OceanSimulations
1616using Oceananigans
17+ using Dates
1718using CairoMakie
1819
1920# # Computing fluxes on the ECCO2 grid
@@ -24,10 +25,10 @@ grid = ECCO_immersed_grid()
2425
2526# We visualize the bottom height of the ECCO grid using CairoMakie.
2627
27- fig = Figure ( )
28- ax = Axis (fig[1 , 1 ] )
29- heatmap! (ax, interior (grid . immersed_boundary . bottom_height, :, :, 1 ))
30- save (" ECCO_continents.png" , fig) # hide
28+ fig, ax, hm = heatmap (grid . immersed_boundary . bottom_height )
29+ Colorbar (fig[1 , 2 ], hm, height = Relative ( 3 / 4 ), label = " Depth (m) " )
30+
31+ save (" ECCO_continents.png" , fig)
3132
3233# 
3334
@@ -67,10 +68,10 @@ set!(ocean.model; T=T_metadata, S=S_metadata)
6768
6869coupled_model = OceanSeaIceModel (ocean; atmosphere, radiation= Radiation (eltype))
6970
70- # # Now that the surface fluxes are computed, we can extract and visualize them.
71- # # The turbulent fluxes are stored in `coupled_modelinterfaces .atmosphere_ocean_interface.fluxes`.
71+ # Now that the surface fluxes are computed, we can extract and visualize them.
72+ # The turbulent fluxes are stored in `coupled_model.interfaces .atmosphere_ocean_interface.fluxes`.
7273
73- fluxes = coupled_model. interfaces. atmosphere_ocean_interface. fluxes
74+ fluxes = coupled_model. interfaces. atmosphere_ocean_interface. fluxes
7475λ, φ, z = nodes (fluxes. sensible_heat)
7576
7677fig = Figure (size = (800 , 800 ), fontsize = 15 )
@@ -90,6 +91,6 @@ heatmap!(ax, λ, φ, interior(fluxes.y_momentum, :, :, 1); colormap = :bwr)
9091ax = Axis (fig[3 , 1 ], title = " Water vapor flux (kg m⁻² s⁻¹)" , xlabel = " Longitude" , ylabel = " Latitude" )
9192heatmap! (ax, λ, φ, interior (fluxes. water_vapor, :, :, 1 ); colormap = :bwr )
9293
93- save (" fluxes .png" , fig)
94+ save (" surface_fluxes .png" , fig)
9495
95- # 
96+ # 
0 commit comments