# Nitrogen deposition ancillary file for SSP370 for ACCESS-ESM1.5 **URL:** \ **Category:** Earth System\ **Tags:** help, um-ancillary, access-esm\ **Created:** [10 November 2023 00:27 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583 "2023-11-10T00:27:50Z")\ **Posts on this page:** 20\ **Page:** 1
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [10 November 2023 00:27 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/1 "2023-11-10T00:27:50Z")
Hi, we require help with generating the nitrogen deposition ancillary file for our future Scenario SSP3-7.0 with ACCESS-ESM1.5. I turns out we used the SSP2-4.5 nitrogen deposition for those runs in the past and we don’t actually have the correct SSP3-7.0 data set. The nitrogen deposition fields are available from input4MIPS and the UKMO has a script to get this in the right format for the UM. This has also been briefly discussed with @MartinDix and Martin already provided some suggestions on how to get this done. Thank you. Cheers, Tilo ---
**Author:** ![MartinDix](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/martindix/32/32_2.png) [@MartinDix](https://forum.access-hive.org.au/u/MartinDix)\ **Post date:** [10 November 2023 00:47 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/2 "2023-11-10T00:47:35Z")
[https://code.metoffice.gov.uk/trac/ancil/wiki/CMIP6/ForcingData/NitrogenDeposition](https://code.metoffice.gov.uk/trac/ancil/wiki/CMIP6/ForcingData/NitrogenDeposition) describes creation of the ancillaries and points to a python script used to generate them from the input4MIPS source. It says only 1-2.6, 2-4.5 and 5-8.5 versions were created. However the 3-7.0 source files are available from input4MIPS (e.g. [Data Information View](https://handle-esgf.dkrz.de/lp/21.14100/b6aeb46a-f58b-3b6f-913d-ce09113d3ddc)) so we should be able to do it. ---
**Author:** ![MartinDix](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/martindix/32/32_2.png) [@MartinDix](https://forum.access-hive.org.au/u/MartinDix)\ **Post date:** [10 November 2023 00:49 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/3 "2023-11-10T00:49:27Z")
I think this would be useful way to get more familiar with ants as well as helping Tilo. ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [10 November 2023 01:16 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/4 "2023-11-10T01:16:29Z")
I will have a look at it! ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [19 November 2023 22:57 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/5 "2023-11-19T22:57:08Z")
Hi @atteggiani. Do you have an update on this? Anything I can help with? Cheers, Tilo ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [20 November 2023 04:38 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/6 "2023-11-20T04:38:13Z")
Hi @tiloz, Sorry I have starting looking at it more in details only today. However, I already have almost all the steps sorted out. I need to check that the script I wrote is effectively producing what you need, and I will then send you the path to the Nitrogen deposition ancillary file for SSP370 (together with the script used to create it). Cheers Davide ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [20 November 2023 04:46 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/7 "2023-11-20T04:46:59Z")
Great. Thank you @atteggiani! Cheers, Tilo ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [23 November 2023 05:27 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/8 "2023-11-23T05:27:21Z")
Hello @tiloz, I need a bit more information on some specifics of what you want, to produce the ancillary file with all correct headers. Can you please point me to the SSP2-4.5 file you used in the past for the experiment? Thank you, Davide ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [23 November 2023 07:24 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/9 "2023-11-23T07:24:14Z")
Hi @atteggiani. You can find the nitrogen deposition file that we used for our SSP245 simulations here: /g/data1a/p66/txz599/data/ancil/CMIP6/Ndep\_2014\_2101\_ssp245.anc Cheers, Tilo ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [27 November 2023 08:07 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/10 "2023-11-27T08:07:32Z")
Hi @tiloz, Thank you for providing your file. You can find the Nitrogen deposition for SSP370 here: `/scratch/public/dm5220/tilo/ndep_2015_2099_ssp370_esm1.5.anc` The script to create the file is `make_nitrogen_deposition_ancil.py`. I have also created a public gist → [Example of how to create the Nitrogen deposition ancillary file for ACCESS-ESM1.5 from a netCDF source. · GitHub](https://gist.github.com/atteggiani/797585a6ff215a590174dadd3526bcef) Please let me know if you have problems running ACCESS-ESM1.5 with the file. Cheers Davide ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [27 November 2023 22:42 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/11 "2023-11-27T22:42:41Z")
Thanks @atteggiani! I will have a look at the file. Cheers, Tilo ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [28 November 2023 03:49 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/12 "2023-11-28T03:49:47Z")
Hi @atteggiani, I had a look at the file and I think the numbers are too small. From your script it looks like you are only including the “drynhx” component. I think there are four components: dry and wet nitrate (drynoy and wetnoy) and dry and wet ammonium (drynhx and wetnhx) that go into the nitrogen deposition field. This should all be all available from input4mips. Also, the model requires an additional year of data at either end (for interpolation), which means we need the nitrogen deposition data set from 2014 to 2101 to run our simulation from 2015 to 2100. The 2014 data can be obtained from the last year of the historical period and the 2100/2101 I think are just a repeat of the last year available. Cheers, Tilo ![input4mips](https://us1.discourse-cdn.com/flex020/uploads/access1/original/2X/4/42b87883b32216be6a4712e3a52bea36d0748f45.png) ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [28 November 2023 04:43 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/13 "2023-11-28T04:43:16Z")
Hi @tiloz, > [@tiloz](#): > > I had a look at the file and I think the numbers are too small. From your script it looks like you are only including the “drynhx” component. I think there are four components: dry and wet nitrate (drynoy and wetnoy) and dry and wet ammonium (drynhx and wetnhx) that go into the nitrogen deposition field. This should all be all available from input4mips. In the ancillary file you provided (`/scratch/public/dm5220/tilo/ndep_2015_2099_ssp370_esm1.5.anc`) there is only one field (STASH 884) for the nitrogen deposition. Where are the other fields you mention? For the values being too small, they are sensible with the [source file](http://aims3.llnl.gov/thredds/fileServer/user_pub_work/input4MIPs/CMIP6/ScenarioMIP/NCAR/NCAR-CCMI-ssp370-1-0/atmos/mon/drynhx/gn/v20181018/drynhx_input4MIPs_surfaceFluxes_ScenarioMIP_NCAR-CCMI-ssp370-1-0_gn_201501-209912.nc) I used. I calculated the units for the Nitrogen deposition ancillary to be `gN m-2 d-1/`, while the source file has `kgN m-2 s-1`. Do you think that makes sense? Do you think the source file is the correct one? > [@tiloz](#): > > the model requires an additional year of data at either end (for interpolation), which means we need the nitrogen deposition data set from 2014 to 2101 to run our simulation from 2015 to 2100. The 2014 data can be obtained from the last year of the historical period and the 2100/2101 I think are just a repeat of the last year available. I did not know about it. In that case, where is the historical source file (or SSP370 ancil file) from where I can take year 2014? Also, which website/search did you do to get the outputs of the image you shared? Cheers Davide ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [28 November 2023 05:36 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/14 "2023-11-28T05:36:41Z")
Hi @atteggiani, We only provide one nitrogen deposition field (STASH 884) in our model which I think is the sum of the four components (drynoy, wetnoy, drynhx, wetnhx). The units of g/m2/d are correct. You can download the individual components here: [https://esgf-node.llnl.gov/search/input4mips/](https://esgf-node.llnl.gov/search/input4mips/) The nitrogen deposition file for the historical period (1850-2014) is here: /g/data/p66/txz599/data/ancil/CMIP6/Ndep\_1849\_2015.anc Cheers, Tilo ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [30 November 2023 01:10 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/15 "2023-11-30T01:10:33Z")
Hi @tiloz, > [@tiloz](#): > > We only provide one nitrogen deposition field (STASH 884) in our model which I think is the sum of the four components (drynoy, wetnoy, drynhx, wetnhx). I see now. I will download the source files from the website you pointed out and try producing the field as sum of the four components. I will also prepend the last year from the historical run and then duplicate the final year so the new ancil file will span from 2014 to 2101 (end of 2100 techically). I should be able to provide you with the updated files within the end of the week. Cheers Davide ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [30 November 2023 08:59 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/16 "2023-11-30T08:59:24Z")
Hi @tiloz, Please find the modified version of the file and script in the same directory: `/scratch/public/dm5220/tilo` The updated gist is also available at the same link: > [@atteggiani](#): > > public gist → [Example of how to create the Nitrogen deposition ancillary file for ACCESS-ESM1.5 from a netCDF source. · GitHub](https://gist.github.com/atteggiani/797585a6ff215a590174dadd3526bcef) In the final ancil data: - Year 2014 has been taken from the historical data. Historical data’s last year is actually 2015, not 2014. However, year 2015 is present in the source files. Therefore, in the final data, year 2014 has been taken from the historical data’s second-last year, whereas year 2015 comes from the source files; - Years from 2015 to 2099 have been computed from the source files as the sum of the 4 fields you indicated: `wetnoy`, `drynoy`, `wetnhx` and `drynhx`; - Years 2100 and 2101 have been computed by repeating year 2099 twice (since the source data ends at the end of 2099). Please let me know if the processing done for the different years make sense to you, and if the new `ndep` ancil file has sensible values. Cheers Davide ---
**Author:** ![tiloz](https://avatars.discourse-cdn.com/v4/letter/t/ba8739/32.png) [@tiloz](https://forum.access-hive.org.au/u/tiloz)\ **Post date:** [1 December 2023 04:41 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/17 "2023-12-01T04:41:44Z")
Hi @atteggiani, the data is looking good and in line with the other scenarios. I have just set-up a new run for testing and it seems to run fine with the new ndep ancillary file. Thanks again for your help with this! Cheers, Tilo ![Screenshot 2023-12-01 at 3.39.28 pm](https://us1.discourse-cdn.com/flex020/uploads/access1/original/2X/4/4df3f434a025d0840fbad96448e2f6193b932104.png) ---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [3 December 2023 04:13 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/18 "2023-12-03T04:13:19Z")
---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [3 December 2023 04:13 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/19 "2023-12-03T04:13:42Z")
---
**Author:** ![atteggiani](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/atteggiani/32/212_2.png) [@atteggiani](https://forum.access-hive.org.au/u/atteggiani)\ **Post date:** [3 December 2023 04:13 UTC](https://forum.access-hive.org.au/t/nitrogen-deposition-ancillary-file-for-ssp370-for-access-esm1-5/1583/20 "2023-12-03T04:13:51Z")