Skip to content

Conversation

@trontrytel
Copy link
Member

What we discussed today - apply "moisture fixer" inside dss to fix small tracer at the expense of water vapor

@trontrytel
Copy link
Member Author

Will test in aquaplanet 1M and let you know how it looks.

if tracer_name in propertynames(Y.c)
ρq_tot = Y.c.ρq_tot
tracer = getproperty(Y.c, tracer_name)
@. tracer += min(max(0, -tracer), max(0, ρq_tot))
Copy link
Member

@szy21 szy21 Jan 13, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Note that the sum of the change in tracers (e.g. if all four are negative) can still be larger than ρq_tot because it is constrained individually. I guess we don't want to have an artificial order for the clipping. But then, do we need the constraint of ρq_tot at all?

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

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants