Volcanic Optical Depth (SAOD)

The volcanic pipeline calculates stratospheric aerosol optical depth (SAOD) by integrating extinction over latitude and stratospheric vertical layers, producing monthly timeseries and tapering extensions.

Detailed Module Reference

cmip7_HI_volcanic_generate

esm1p6_ancil.volcanic.cmip7_HI_volcanic_generate

save_hi_stratospheric_aerosol_optical_depth(args, dataset_path)

Calculate the average stratospheric aerosol optical depth (SAOD) for each historical month by averaging extinction over latitude, and summing over stratospheric layers. Save to the save file.

cmip7_PI_volcanic_generate

esm1p6_ancil.volcanic.cmip7_PI_volcanic_generate

average_stratospheric_aerosol_optical_depth(dataset_path)

Calculate the average stratospheric optical depth (SAOD) by averaging extinction over both time and latitude, and summing over stratospheric layers.

cmip7_pi_volcanic_patch(average_saod)

Patch the VOLCTS_val variable in the coupling namelist

mean_over_pi_months(cube)

Find the time average SAOD by averaging over months in the pre-industrial year, weighted by month length.

cmip7_SM_volcanic_generate

esm1p6_ancil.volcanic.cmip7_SM_volcanic_generate

save_sm_stratospheric_aerosol_optical_depth(args, dataset_path, pi_mean_saod)

Calculate the average stratospheric aerosol optical depth (SAOD) for each historical month by averaging extinction over latitude, and summing over stratospheric layers. Save to the save file.

cmip7_volcanic

esm1p6_ancil.volcanic.cmip7_volcanic

constrain_to_latitude_band(cube, band)

Constrain to one of four equal latitude bands.

constrain_to_wavelength(cube, wavelength)

Constrain to just the prescribed wavelength.

constrain_to_year_month(cube, year, month)

Constrain to a given year and month. See #iris.coords.Cell.point in scitools-iris.readthedocs.io/en/stable/generated/api/iris.coords.html

mean_over_latitudes(cube)

Find the mean over all latitude bands, weighted by area.

save_early_saod(volcanic_beg_year, save_file, pi_mean_saod)

Save the PI average SAOD for all years before volcanic_beg_year

save_stratospheric_aerosol_optical_depth(args, volcanic_beg_year, volcanic_end_year, dataset_path, save_dirpath, pi_mean_saod=None, save_beg_year=SAOD_BEG_YEAR, save_end_year=SAOD_END_YEAR, hold_constant=False)

Calculate the average stratospheric aerosol optical depth (SAOD) for each month by averaging extinction over latitude, and summing over stratospheric layers. Save to the save file.

sum_over_height_layers(cube)

Calculate the stratospheric aerosol optical depth by summing over stratospheric layers, weighted by layer height.

taper_saod(volcanic_end_year, save_end_year, saod_for_beg_year, saod_for_end_year)

Interpolate between the saod values in saod_for_beg_year and saod_for_end_year. The SAOD values taper from saod_for_end_year towards saod_for_beg_year for NBR_TAPER_YEARS, and remain at saod_for_beg_year afterwards.