Easy forcing perturbation experiments in ACCESS-OM2 + Tutorial

COSIMA
,
1

I’m not sure how many people are aware that there is a system in ACCESS-OM2 (actually in YATM) to easily setup forcing perturbation experiments where the atmospheric forcing fields are modified. If you like running experiments like this you really should be!

The system was written by Nic Hannah, Andrew Kiss and others and is described here: Tutorials · COSIMA/access-om2 Wiki · GitHub

It’s quite flexible - a range of additive or multiplicative perturbations to any of the forcing fields that are spatially or temporally dependent/independent (or combinations of these) can be done, without the need to create an entire new copy of the JRA55 (or ERA5, or whatever you’re using) forcing files.

I’ve written up a very brief tutorial on how to do a simple perturbation experiment for a student here: cosima-scripts/ACCESS-OM2_forcing_perturbation_tutorial.ipynb at master · rmholmes/cosima-scripts · GitHub . Note that it was written in a rush since I’m off on leave for a while tomorrow. If someone has the inclination, it would be great to improve on it, expand it to some more complicated cases, make it consistent with current NRI best practice and put it somewhere more accessible.

@cbull maybe this is another candidate for your training series?

2

Great, thanks @rmholmes, I’ve added a link to your tutorial on the wiki Tutorials · COSIMA/access-om2 Wiki · GitHub

3

Writing this down before I forget:

  • We discovered a “bug” where in ERA5 simulations with “separable” forcing perturbations the forcing perturbations are not logged in the YATM log file atmosphere/log/matmxx.pe00000.log. I don’t know why, but the forcing anomalies are still applied.
  • When running forcing perturbations with ERA5 forced simulations there are a few additional issues with netcdf packing, fill values and units that could cause you problems if you’re unaware. See this tutorial for a minimally-documented example of how to do this kind of thing in ERA5 (as opposed to with JRA55 in the more documented tutorial).
4

Also worth mentioning that there was a training presentation on this last August - see these slides and video.