Concatenate AMIP Ancillaries - Revision history

S.wales at 03:27, 12 December 2019

2019-12-12T03:27:44Z

S.wales at 05:00, 12 March 2019

2019-03-12T05:00:04Z

ScottWales: Imported from Wikispaces

2015-07-24T02:34:34Z

Imported from Wikispaces

New page

The ancillary values used for AMIP model runs come in two sections - historical values and RCP projections. The crossover between these sections happens between the years 2000 and 2002, with the RCP files beginning in January 2000 and the historical values ending in December 2001.

For AMIP model runs going over this date range either the job need to be stopped, the ancillary files changed over and the the job restarted or the two ancillary sections need to be combined.

The AMIP ancillary files can be found in ~access/data/ancil/CMIP5. The following ancillary categories are available for the RCP scenarios:
<syntaxhighlight>
BC_hi: Soot
BIO: Biomass emissions
OCFF: Organic carbon & fossil fuels
sulp: Sulphyr emissions
</syntaxhighlight>

To create a combined file following the full range:

1. Convert both segments to netcdf (e.g. with xconv)

2. Trim the ancillary files so that the date ranges don't overlap
> <syntaxhighlight>
cdo selyear,2002/2100 OCFF_RCP85_2000_2100.nc OCFF_2002-2100.nc
</syntaxhighlight>

3. Concatenate the two netcdf segments
> <syntaxhighlight>
ncrcat OCFF_IPCC_1850_2000.nc OCFF_2002-2100.nc OCFF_full.nc
</syntaxhighlight>

4. Use xancil to convert the concatenated netcdf file to UM format

* Make sure to set the version in Configuration/General to 7.3
* Use the 'generalised ancillary' section. Add the netcdf file and select the variables, the STASH code should get filled in automatically.
* You may need to override the dates in the file, specifying the start date and interval (1 month)
** To get values on the 16th of each month with a Gregorian calendar generate the ancillary using a 360 day calendar, then post-process by running
>> <syntaxhighlight>
python ~access/apps/pythonlib/umfile_utils/change_calendar365.py ANCIL
</syntaxhighlight>
Check the ancillary file using xancil, making sure the field & dates are correct

← Older revision		Revision as of 05:00, 12 March 2019
Line 35:		Line 35:
	</syntaxhighlight>		</syntaxhighlight>
	Check the ancillary file using <span style="font-family:monospace">xancil</span>, making sure the field & dates are correct		Check the ancillary file using <span style="font-family:monospace">xancil</span>, making sure the field & dates are correct
		+
		+	[[Category:Unified Model]]

@@ Line 4: / Line 4: @@
 The AMIP ancillary files can be found in <span style="font-family:monospace">~access/data/ancil/CMIP5</span>. The following ancillary categories are available for the RCP scenarios:
-<syntaxhighlight>
+<syntaxhighlight lang=text>
 BC_hi: Soot
 BIO:   Biomass emissions
@@ Line 16: / Line 16: @@
 . Trim the ancillary files so that the date ranges don't overlap
-> <syntaxhighlight>
+> <syntaxhighlight lang=bash>
 cdo selyear,2002/2100 OCFF_RCP85_2000_2100.nc OCFF_2002-2100.nc
 </syntaxhighlight>
 . Concatenate the two netcdf segments
-> <syntaxhighlight>
+> <syntaxhighlight lang=bash>
 ncrcat OCFF_IPCC_1850_2000.nc OCFF_2002-2100.nc OCFF_full.nc
 </syntaxhighlight>
@@ Line 31: / Line 31: @@
 * You may need to override the dates in the file, specifying the start date and interval (1 month)
 ** To get values on the 16th of each month with a Gregorian calendar generate the ancillary using a 360 day calendar, then post-process by running
->> <syntaxhighlight>
+>> <syntaxhighlight lang=bash>
 python ~access/apps/pythonlib/umfile_utils/change_calendar365.py ANCIL
 </syntaxhighlight>