drop max melting rate to 10 m/yr

grahamedwards · Feb 17, 2024 · 5330ef6 · 5330ef6
1 parent 3af7d3f
commit 5330ef6
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/metropolis.jl b/src/metropolis.jl
@@ -20,7 +20,7 @@ Evalute strict constraints on priors that will automatically reject a proposal w
 - Onset date beyond the climate record timespan (in ka) (`record_max_age`)
 - Nonphysical subglacial thresholds -- melting at lower benthic δ¹⁸O than freezing or values exceeding the record extrema (`climatelimits`).
 - Freezing rate is ≤ 0.4dt/dz (to prevent an error in a log-calculation in [`PorewaterDiffusion.boundaryconditions`](@ref)).
-- Annual melting rate must be less than the thickness of the Antarctic Ice Sheet (~5 km)
+- Annual melting rate must be less than that observed at the Thwaites grounding line (<10 m/yr, [Davis+ 2023](https://www.nature.com/articles/s41586-022-05586-0)).
 
 """
 function strictpriors(p::Proposal, record_max_age::Number, climatelimits::Tuple{Number,Number}, k::Constants)
@@ -33,7 +33,7 @@ function strictpriors(p::Proposal, record_max_age::Number, climatelimits::Tuple{
     x &= p.frz2mlt < climatelimits[2]
 
     x &= p.dfrz <= 0.4k.dz/k.dt
-    x &= p.dmlt < 5000.
+    x &= p.dmlt < 10.
 
     x
 end

diff --git a/test/metro.jl b/test/metro.jl
@@ -8,7 +8,7 @@ spargs = (1000., (1.,4.), Constants())
 @test !PorewaterDiffusion.strictpriors(Proposal(500,.1,1,2,5), spargs...) # frz2mlt too high
 @test !PorewaterDiffusion.strictpriors(Proposal(500,.1,1,3,2), spargs...) # sez2frz < frz2mlt 
 @test !PorewaterDiffusion.strictpriors(Proposal(500,.3,1,2,3), spargs...) # p.dfrz too high.
-@test !PorewaterDiffusion.strictpriors(Proposal(500,.1,1e5,2,3), spargs...) # p.dfrz too high.
+@test !PorewaterDiffusion.strictpriors(Proposal(500,.1,20.,2,3), spargs...) # p.dfrz too high.
 
 
 # proposaljump