Merge pull request #180 from jbisits/jib/isotropic-kedissipation

Move to isotropic KE dissipation rate

Author: jbisits

src/StaircaseShenanigans.jl

@@ -14,7 +14,7 @@ using SeawaterPolynomials: thermal_expansion, haline_contraction, total_density
 using GibbsSeaWater: gsw_alpha, gsw_beta
 using NCDatasets, JLD2
 using Statistics
-using Oceanostics: KineticEnergyDissipationRate, KineticEnergy, PotentialEnergy, BuoyancyProductionTerm
+using Oceanostics: IsotropicKineticEnergyDissipationRate, KineticEnergy, PotentialEnergy, BuoyancyProductionTerm
 
 import SeawaterPolynomials.SecondOrderSeawaterPolynomials: RoquetSeawaterPolynomial
 

src/staircase_simulation.jl

@@ -318,7 +318,7 @@ function save_computed_output!(simulation, sdns, schedule, save_file, output_dir
     σ_ha = Average(σ, dims = (1, 2))
     N² = buoyancy_frequency(model)
     N²_ha = Average(N², dims = (1, 2))
-    ε = KineticEnergyDissipationRate(model)
+    ε = IsotropicKineticEnergyDissipationRate(model)
     ∫ε = Integral(ε)
     ε_maximum = Reduction(maximum!, ε, dims = (1, 2, 3))
     Eₖ = KineticEnergy(model)