pyqg Subgrid Forcing Datasets

pyqg Subgrid Forcing Datasets

This repository contains a set of zarr datasets which (when uncompressed) contain snapshots from different series of pyqg.QGModel simulations run with different sets of parameters (one which produces isotropic eddies, stored in eddy/, and one which produces strong jets, stored in jet/. These parameters are stored as JSON in each dataset's metadata as the pyqg_params attribute.

Descriptions of each dataset:

low_res.zarr contains snapshots saved every 1000 hours from 5 simulations run for 10 years at nx=64 (which is eddy-permitting but not resolving)

high_res.zarr contains similar snapshots but for simulations at nx=256 (which is fully eddy-resolving)

forcing{1,2,3}.zarr contain subgrid forcing variables (described below) and coarsened state variables which can be used to predict them for 275 simulations. Each forcing dataset is computed using a different filtering and coarse-graining operator:

forcing1 is computed by truncating high-frequency potential vorticity (PV) modes in spectral space, then applying the sharp, scale-selective spectral filter described here

forcing2 is computed by truncating high-frequency PV modes in spectral space, then applying a Gaussian spectral filter used here

forcing3 is computed by applying a diffusive filter from GCM-Filters to the PV in real space, then averaging in real space

Descriptions and units for the variables in each dataset are stored in metadata which can be viewed when the files are loaded. More details (and links to datasets in different configurations) can be found at https://github.com/m2lines/pyqg_parameterization_benchmarks.