diff --git a/scripts/plot_one_degree.jl b/scripts/plot_one_degree.jl index ae9c091..4ad6224 100644 --- a/scripts/plot_one_degree.jl +++ b/scripts/plot_one_degree.jl @@ -1,6 +1,6 @@ using Oceananigans, CairoMakie, GeoMakie -P = FieldTimeSeries("bgc/near_global_360_150_48_ConvectiveAdjustmentVerticalDiffusivity_fields.jld2", "P"); +P = FieldTimeSeries("bgc/near_global_360_150_48_ConvectiveAdjustmentVerticalDiffusivity_fields.jld2", "P") n = Observable(1) @@ -10,13 +10,13 @@ fig = Figure() ax = GeoAxis(fig[1, 1];title, dest = "+proj=natearth2") -xc, yc, zf = nodes(P.grid, Center(), Center(), Face()) +x, y, z = nodes(P) -hm = heatmap!(ax, xc, yc, P_plt, colorrange = (0.0001, 5), colorscale = log10) +hm = heatmap!(ax, x, y, P_plt, colorrange = (0.0001, 5), colorscale = log10) Colorbar(fig[2, 1], hm, vertical= false, label = "Phytoplankton (mmolN/m³)") -CairoMakie.record(fig, "P.mp4", 1:10:length(P.times)) do i; - n[] = i; - @info "$(n.val)"; -end \ No newline at end of file +CairoMakie.record(fig, "P.mp4", 1:10:length(P.times)) do i + n[] = i + @info "$(n.val)" +end