NOAA-GFDL / CEFI-regional-MOM6

A repository containing essential tools, XML files, and source codes for collaborators of the Climate, Ecosystems, and Fisheries Initiative (CEFI) to conduct simulations.
Other
19 stars 16 forks source link

Update/g interface #89

Closed yichengt900 closed 2 months ago

yichengt900 commented 2 months ago

This PR updates the MOM6 and ocean_BGC tags, incorporating the latest changes to the generic tracer interface. With the removal of eqn_of_STATE, the mixed layer depth (MLD) used for photoacclimation, as calculated by MOM6, will now be correctly applied in COBALTv3. The model configurations have also been updated to reflect these changes.

These changes also introduce additional flexibility in how the photoacclimation depth is calculated at runtime through a number of new parameters.

PHA_MLD_CALC (default false) if false, use a fixed value for the photoacclimation mixed layer depth. This option should be used to skip calculating a MLD. Likely this will be true for people who are not running COBLATv3 and do not want to calculate an unused depth. PHA_MLD_VAL (default 0.0) This is the fixed value for PHA_MLD_CALC. The purpose of using 0.0 as a default is to lead to a crash or obviously wrong solution when using COBALTv3. A warning message has been added to COBALT to give users some insight into the problem. PHA_MLD_USE_DELTA_RHO (default false) If true, use a density difference to calculate the photoacclimation MLD. PHA_MLD_HREF (default 0.0) Similar to the diagnostic MLD, the reference depth can be adjusted. A value of 10.0m would be most similar to de Boyer Montegut et al. 2004 PHA_MLD_DRHO (default 0.03) The target difference in the density at h_ref and the bottom of photoacclimation MLD. PHA_MLD_USE_DELTA_ENG (default false) If true, use the energy based MLD to find the photoacclimation depth (not tested). PHA_MLD_DENG (default 25.0) The threshold for the energy based MLD (not tested)

While the photoacclimation MLD calculated by MOM6 is quite similar to the previous MLD from COBALT, it is not identical, which will lead to changes in the baselines.