General Structure

L0 - Concat different timesteps and generate a NETcdf file

This generates an AerData class which converts all matched iris files within a simulation into a NetCDF class of the same name. Combination of 4 monthly UM files with 270ish variable totalling ~2GB takes between 5 and 15 minutes to complete on a desktop.

If a NetCDF of already exists this is opened without needing to run the L0 script again.

import aervis
aervis.AerData('bk2783jk')

L1 - calculate useful metrics

d = AerData('bk2783jk') # load dataset
d.getL1(override=True) # run the L1 code

Wiki: https://github.com/wolfiex/AerVis/wiki/L1-Processing

L2 - Plotting prep and presentation

colourmaps

opearation functions

The basic plot function allows for a custom operation function to be employed on the data.

benchmark = AerData('normal')
perturbed = AerData('perturbed')


def plotfn(a,b):
   ''' a custom plotting function '''
   return a.data - b.data / a.data

aervis.plot(a=benchmark, b=perturbed, customfn = plotfn)