search

firedrakeproject / gusto_case_studies

A repository for Gusto case studies
1 stars 0 forks source link

Terminator #9

Closed ta440 closed 3 months ago

ta440 commented 1 year ago

This runs the Terminator Toy Test (Lauritzen et al. 2014) using the terminator_toy physics scheme and the SplitPrescribedTransport timestepper. A dry density field of two Gaussians is advected using a Nair Lauritzen non-divergent flow profile. There is also the interaction of two chemical species (X and X2) through very stiff ODEs which are solved using a Backwards Euler method.

The dry density profile is the Gaussian setup from Nair-Lauritzen. Here are the initial (t=0days) and final (t=12days) states at ref_level=4:

terminator_toy_rf4rho_d_t_day0 0 terminator_toy_rf4rho_d_t_day12 0

We can examine the distributions of the species during the simulation (t=0,6,12 days, ref_level=3)

Species X: terminator_toy_testX_t_day0 0 terminator_toy_testX_t_day6 0 terminator_toy_testX_t_day12 0

Species X2: terminator_toy_testX2_t_day0 0 terminator_toy_testX2_t_day6 0 terminator_toy_testX2_t_day12 0

There should be conservation of the total mixing ratio, Xt = X +2X2, between the two interacting species. Additionally, there should be conservation of the tracer density, which is the product of the dry density and the total species mass, rho_X = rho_d*X_t. This is examined using the 'TracerDensity' diagnostic.

Here is the normalised change in the total tracer density over time (ref_level=3):

tracer_density_change_reflev3

tommbendall commented 3 months ago

Merging in despite test failures, and will fix lint here in a following PR