Skip to content

WOMBATlite: Treat key biological loss terms semi-implicitly#97

Merged
dougiesquire merged 3 commits intomainfrom
semi-implicit-loss
Mar 25, 2026
Merged

WOMBATlite: Treat key biological loss terms semi-implicitly#97
dougiesquire merged 3 commits intomainfrom
semi-implicit-loss

Conversation

@dougiesquire
Copy link
Copy Markdown
Collaborator

@dougiesquire dougiesquire commented Mar 25, 2026

This PR implements the changes described in #96

To summarise, there are a few loss terms in the biological tendency equations that are particularly susceptible to instability with Forward Euler time-stepping:

  • phy equation: phymorq = f(phy^2)
  • zoo equation: zoomorq = f(zoo^2)
  • det equation: detremi = f(det^2)
  • caco3 equation: caldiss = f(caco3), aradiss = f(caco3), pocdiss = f(caco3), zoodiss = f(caco3)

This PR rewrites these update equations using a semi-implicit scheme for these terms.

Note, I haven't looked at WOMBATmid, but this presumably has similar issues.

@dougiesquire dougiesquire self-assigned this Mar 25, 2026
@dougiesquire dougiesquire requested a review from pearseb March 25, 2026 03:49
@dougiesquire
Copy link
Copy Markdown
Collaborator Author

@pearseb are you okay with these changes after our discussion today?

I've added your suggestions on the comments. I've also added another commit that guards against some potential overflows - see a443314

@dougiesquire dougiesquire merged commit 9341996 into main Mar 25, 2026
9 checks passed
@dougiesquire dougiesquire deleted the semi-implicit-loss branch March 25, 2026 22:11
@pearseb
Copy link
Copy Markdown
Collaborator

pearseb commented Mar 25, 2026

Nice :) Great contribution Dougie. ACCESS-OM2-0.1º-WOMBATlite here we come.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants