The bash script ocean_pp.bash performs diagnostics on monthly or daily ocean model outputs. It relies on a set of functions included in common_ocean_post.bash that can also be used within the Ec-Earth and NEMO post templates when using autosubmit. Many of the functions included in common_ocean_post.bash have been built based on Hui Du's post-processing scripts.
1. A README file describes the various diagnostics which are available as options through the listpost argument.
The outputs of ocean_pp.bash are of 3 types :
2. The list of arguments of ocean_pp.bash are the followings :
Fill up either ltime0/ltimef or year0/yearf below
If you fill up the year argument, complete years will be processed, year by year from moni
Valid options :
Ec2.3_O1L42 for Ec-Earth 2.3 ORCA1 L42 Ec3.0_O1L46 for Ec-Earth 3.0 ORCA1 L46 Ec3.0_O25L46 for Ec-Earth 3.0 ORCA0.25 L46 N3.2_O1L42 for Nemo 3.2 ORCA1 L42 N3.3_O1L46 for Nemo 3.3 ORCA1 L42 nemovar_O1L42 for Nemo COMBINE and ORAS4 ORCA1L42
Beware that with a given horizontal and vertical resolution, let's say ORCA1L46 (O1L46), the 3d land/sea mask is dependent even on the subversion of NEMO. For example, although NEMO3.3 (N3.3_O1L46) is included in Ec-Earth3.0 (Ec3.0_O1L46), N3.3_O1L46 and Ec3.0_O1L46 will not give the exact same result because the masks are different
If you post-process daily or monthly means the output directory is daily_mean or monthly_mean. daily means is only tested for nemo currently
3. Some more specifications:
BUT : You have already postprocessed for exp t001: 1) members 0-2 between leadtimes 1 and 36, 2) members 0 between leadtimes 37 and 60 If you process member 2 between leadtimes 25 and 60, you'll find the two above files plus a new one with members 2 between leadtimes 25 and 60 on the nas
1. ocean_pp.bash handles now any range of leadtimes instead of complete years only. However, the max_moc diagnostic can not be computed if you do not have complete years since it is a diagnostic computed from annual means.
2. ocean_pp.bash now post-processes daily means also, instead of only monthly means before. You only need to specify as output directory daily_mean instead of monthly_mean because the number of time steps and dates are read from the output. It has been tested with NEMO only simulations but this required some changes in the NEMO templates so that some changes in the EcEarth templates might be needed also before using this option for EcEarth simulations.
An update was produced in February 2013 by Pierre-Antoine. The main point was to extend the list of variables that could be post-processed. Ocean_pp_new.bash includes the following changes:
If “default” option is specified, it will treat the same variables as before (sosstsst, sosaline and somxl010 for ocean, isnowthi, iicethic, ileadfra, iicetemp, ice_pres for sea ice). If nothing is specified, it will extract all the 2D variables present in the input files. Else, it will treat exclusively the variables listed.
If you want to include a new diagnostic. Look for the @ signs in the common_ocean_post.bash then ocean_pp.bash and you'll find the locations where you have to include lines of code and some additional instructions so that everything works fine. Feel free to ask if something is not clear.
Some plans for volunteers:
If you want to make available the post-processing tools for new grids, you need to choose a flag for this new grid (such as Ec2.3_O1L42), and then prepare all the configuration files such as done here : /cfu/autosubmit/con_files/ for example for any other flags. Namely you need to prepare the files :
virginie.guemas@ic3.cat, pierre-antoine.bretonniere@ic3.cat, isabel.andreu-burillo@ic3.cat
The purpose is to standardize CFU EC-Earth outputs, to make it follow the common netcdf file conventions of the ENSEMBLE project. Most importantly, it makes CFU people who works on different date sets in a common framework. That means common diagnostic scripts can be used on different dataset. With this in mind, we save all members of one variable of one starting date in one file.
Surface variables| Short Name| Long name| Unit |
tas | near surface temperature | K |
tos | sea surface temperature | K |
prlr | total precipitation | m s-1 |
psl | mean sea level pressure | Pa |
rss | surface net download shortwave flux | W s-2 |
rls | surface net download longwave flux | W s-2 |
hfssd | surface downward sensible heat flux | W s-2 |
hflsd | surface downward latent heat flux | W s-2 |
Pressure level variables| Short Name| Long name| Unit |
ta | air temperature | K |
ua | zonal wind positive to the east | m s-1 |
va | meridional wind positive to the north | m s-1 |
g | geopotential | m2 s-2 |
/cfu/pub/scripts/postp_ecearth/postp_ecearth.sh b014 19601101 # “b014” is the experiment ID and “19601101” is the starting date. In case you only want to post-process subset of the variables, you only need to specify it on the header of the script
In February 2013, an update was added to this script to automate the the treatment of the options and to extend the variable list that could be treated. are the main changes:
pierre-antoine.bretonniere@ic3.cat, muhammad.asif@ic3.cat, huidu@ic3.cat
The purpose is to select daily outputs (every 6 hours for variables with one level and every 12 hours for pressure level variables) from the GRIB files IFS produces, to make it follow the common NetCDF file conventions of the ENSEMBLES project in the same way as it is done for the EC-Earth monthly outputs. 24 basic variables are saved for the purpose of performing diagnostics if intraseasonal variability.
Surface variables| Short Name| Long name| Unit |
tas | near surface temperature | K |
tos | sea surface temperature | K |
prlr | total precipitation | m s-1 |
psl | mean sea level pressure | Pa |
rss | surface net downward shortwave flux | W s-2 |
rls | surface net downward longwave flux | W s-2 |
rsds | solar radiation incident at the surface | W s-2 |
rlds | longwave radiation incident at the surface | W s-2 |
hfssd | surface downward sensible heat flux | W s-2 |
hflsd | surface downward latent heat flux | W s-2 |
tdps | near-surface dewpoint temperature | K |
uas | near-surface zonal wind positive to the east | m s-1 |
vas | near-surface meridional wind positive to the north | m s-1 |
rst | difference of solar radiation from above and from below at the top of the atmosphere | W s-2 |
rlut | difference of longwave radiation from above and from below at the top of the atmosphere | W s-2 |
snld | snow depth in metres of water equivalent | m |
evlwr | moisture flux into the surface from the atmosphere | m s-1 |
ts | surface temperature | K |
clt | cloud area fraction for the whole atmosphere column | none |
Pressure-level variables| Short Name| Long name| Unit |
ta | air temperature | K |
ua | zonal wind positive to the east | m s-1 |
va | meridional wind positive to the north | m s-1 |
g | geopotential | m2 s-2 |
hus | mass fraction of water vapor in moist air | kg kg-1 |
/cfu/pub/scripts/postp_ecearth/postp_daily.sh $expid 19601101 # where $expid is the experiment id and “19601101” is the starting date; the script will postprocess all the GRIB files for this start date
muhammad.asif@ic3.cat, huidu@ic3.cat