This repository contains scripts for the post-processing of the output from EC-Earth AMIP experiments.
The code aims to calculate the atmospheric meridional energy transport (AMET) based on the direct output from EC-Earth simulation. The full pipeline includes steps as follows:
1.mass budget correction
2.vertical integral of zonally integrated meridional energy transport
More information about the implementation of the mass budget correction through the penalty on baratropic winds can be found in Trenberth's paper 'Climate Diagnostics from Global Analyses: Conservation of Mass in ECMWF Analyses'.
For the calculation of meridional transport, please refer to the work by Trenberth and Caron 'Estimates of Meridional Atmosphere and Ocean Heat Transports'.
Following fields are saved as netcdf files with a frequency of 3 hour:
- 3 hourly U at 850, 500, 200 hPa
- 3 hourly V at 850, 500, 200 hPa
- 3 hourly T at 850, 500, 200 hPa
- 3 hourly Z at 850, 500, 200 hPa
- 3 hourly Q at 850, 500, 200 hPa
- 3 hourly PT,T2M,U10M,V10M,SLHF,SP,MSL,LSP,CP,TCC,SSHF,SSR,STR,TSR,TTR,SRO
Each netCDF file containing all the fields above is approximately 7.6GB per month.
In addition, the monthly mean of each component of meridional energy transport is saved as well:
- monthly mean total AMET
- monthly mean internal energy transport
- monthly mean latent energy transport
- monthly mean geopotential transport
- monthly mean kinetic energy transport
- monthly mean meridional baratropic wind correction
- monthly mean zonal baratropic wind correction
The fields above are saved as two seperate files containing spatial information and zonal integral. The netCDF files are approximately 22MB per month.
The fields above are saved through the command line tool "CDO".
Please run the script in the following folder to perform the post-processing:
ece-postprocess/ece_postprocess/scripts/
$ python ece-postprocess --rundir <path> --expname <exp>
The details about the arguments can be found inside the script.
The output from EC-Earth must be provided with the standard name as:
ICMGG<exp>+<time> # output on gaussian grid
ICMSH<exp>+<time> # output on spectral coordinate, must be changed to gaussian grid
where is the experiment name, and is the time of the exp. They shoud be given as the argument.
The script is recommended to be used together with the workflow manager suite.rc.