# ACCESS-ESM1.5 AMIP run with payu: model output volume

**URL:** <https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853>\
**Category:** Earth System Model\
**Tags:** help-needed, payu, esm\
**Created:** [14 June 2023 03:06 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853 "2023-06-14T03:06:40Z")\
**Posts on this page:** 15\
**Page:** 1

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [14 June 2023 03:06 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/1 "2023-06-14T03:06:40Z")

</div>

I have just done a 1-year AMIP test run of ACCESS-ESM1.5 using payu ( [GitHub - coecms/access-esm-amip](https://github.com/coecms/access-esm-amip)) and appreciated the clear instructions including the ‘Excessive Output Warning’. Years ago I used to be able to select the output using the UMUI but what is the best way to do it now? Editing the STASHC? Or does anyone else running ESM1.5 AMIP already have a set-up for a smaller set of output that they are able to share?  
Thanks, Rachel

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [14 June 2023 04:49 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/2 "2023-06-14T04:49:52Z")

</div>

@holger or @Scott might have some suggestions, or know of others who have run this configuration.

---

<div class="post-metadata">

**Author:** ![spencerwong](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/spencerwong/32/1894_2.png) [@spencerwong](https://forum.access-hive.org.au/u/spencerwong)\
**Post date:** [21 June 2023 23:58 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/3 "2023-06-21T23:58:30Z")

</div>

Hi Rachel,

I’ve attached a STASHC file that @ariaan has put together which I’ve been using for my AMIP runs. The file gets rid of the daily and sub daily output and some of the monthly fields as well, reducing the amount of output each year from 157G to ~14G. Hive won’t let me upload a txt file so I’ve had to paste the text into a pdf.  
[STASHC\_APmod2.pdf](https://forum.access-hive.org.au/uploads/short-url/jZcPxnbklVBnqUTbdwaGOttl8TX.pdf) (32.1 KB)

Cheers,  
Spencer

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [23 June 2023 00:05 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/4 "2023-06-23T00:05:10Z")

</div>

> [@spencerwong](#):
>
> Hive won’t let me upload a txt file so I’ve had to paste the text into a pdf.

I have added `.txt` as a valid upload format. Thanks for the heads-up for that as an issue. If you, or anyone else, has a similar problem uploading a file please feel free to get in touch.

---

<div class="post-metadata">

**Author:** ![spencerwong](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/spencerwong/32/1894_2.png) [@spencerwong](https://forum.access-hive.org.au/u/spencerwong)\
**Post date:** [23 June 2023 05:21 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/5 "2023-06-23T05:21:14Z")

</div>

Thanks for changing that @Aidan!

I’ve just attached the txt file here  
[STASHC\_APmod2.txt](https://forum.access-hive.org.au/uploads/short-url/29YwMDALwGmdpSoU3zyVCGlU8YY.txt) (20.3 KB)

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [23 June 2023 07:23 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/6 "2023-06-23T07:23:20Z")

</div>

@spencerwong - Thanks for the help. It may take me a while to test this but hopefully it will also be useful for others doing AMIP ACCESS\_ESM1.5 runs. Perhaps worth an update to the standard payu case to use the smaller output STASHC?

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [23 June 2023 08:53 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/7 "2023-06-23T08:53:53Z")

</div>

> [@RachelLaw](#):
>
> Perhaps worth an update to the standard payu case to use the smaller output STASHC?

If you want to do that it is best to start an issue in the repo and then make a pull request that fixes it.

> **[Issues · coecms/access-esm-amip](https://github.com/coecms/access-esm-amip/issues)**
>
> Contribute to coecms/access-esm-amip development by creating an account on GitHub.

---

<div class="post-metadata">

**Author:** ![spencerwong](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/spencerwong/32/1894_2.png) [@spencerwong](https://forum.access-hive.org.au/u/spencerwong)\
**Post date:** [29 June 2023 23:48 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/8 "2023-06-29T23:48:49Z")

</div>

I just realised that I really should have included a list of the variables that it saves! I’ve just included a list of the names from the output here.  
[monthly\_vars.txt](https://forum.access-hive.org.au/uploads/short-url/lrzviDr8Qy70qVO0iSb8ZdjgKuE.txt) (7.2 KB)

---

<div class="post-metadata">

**Author:** ![clairecarouge](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/clairecarouge/32/380_2.png) [@clairecarouge](https://forum.access-hive.org.au/u/clairecarouge)\
**Post date:** [4 July 2023 07:03 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/9 "2023-07-04T07:03:07Z")

</div>

It’s tricky to update the STASH file for the standard payu case. For the moment, it has been designed to reproduce the case used for CMIP submissions exactly, as far as I know. That would mean keeping the STASH file as required by CMIP.

So I guess the question is how to best handle configurations with different levels of output so that people can find and identify them easily.

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [5 July 2023 01:12 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/10 "2023-07-05T01:12:03Z")

</div>

> [@clairecarouge](#):
>
> So I guess the question is how to best handle configurations with different levels of output so that people can find and identify them easily.

As long as it is well documented then it should be possible to keep configuration tweaks in separate branches.

If the main use case of this model is _not_ CMIP6 style output, then the CMIP configuration could be in a branch. I’m in favour of this, as we know the model outputs a huge amount of data no-one ever uses, and causes issues down the track with data storage and not knowing what to retain and what to discard, so everything is kept by default.

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [6 July 2023 00:46 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/11 "2023-07-06T00:46:47Z")

</div>

@spencerwong , @clairecarouge , @Aidan - thanks for all your contributions. Spencer - it was a good idea to post the list of variables, as your list didn’t include the carbon variables that I needed so your stashc doesn’t suit my needs. Instead I’ve taken the stashc from the 4xco2 run. It produces output that is intermediate in size between the amip case and Spencer’s example. While it may be what we used for CMIP6, it also has output that wouldn’t be of general interest - the pgh files are timeseries of CO2 output at CO2 monitoring sites, leftover from when we were doing more CO2 transport focussed work. Perhaps this is worth a discussion in the ESM working group sometime - whether we want a ‘library’ of STASHC files for a range of output applications and/or some documentation around how to change the output. Is using the UMUI still the easiest way to do this?

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [6 July 2023 05:44 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/12 "2023-07-06T05:44:48Z")

</div>

> [@RachelLaw](#):
>
> whether we want a ‘library’ of STASHC files for a range of output applications

Given how much trouble people seem to have with them, that sounds like a really good idea.

> [@RachelLaw](#):
>
> Is using the UMUI still the easiest way to do this?

I’m not sure about that. @holger would know, but not sure if he is about right now.

---

<div class="post-metadata">

**Author:** ![holger](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/holger/32/320_2.png) [@holger](https://forum.access-hive.org.au/u/holger)\
**Post date:** [24 July 2023 04:05 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/13 "2023-07-24T04:05:23Z")

</div>

Hi

Sorry for answering so late. I don’t find the `STASHC` all that confusing, you have to mostly look at the `STREQ` lines. I think I have created a document on how to read it a while ago. I’ll have a look.

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [3 May 2024 12:27 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/14 "2024-05-03T12:27:35Z")

</div>

> [@holger](#):
>
> I think I have created a document on how to read it a while ago. I’ll have a look.

Any luck finding that @holger?

---

<div class="post-metadata">

**Author:** ![holger](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/holger/32/320_2.png) [@holger](https://forum.access-hive.org.au/u/holger)\
**Post date:** [8 May 2024 00:26 UTC](https://forum.access-hive.org.au/t/access-esm1-5-amip-run-with-payu-model-output-volume/853/15 "2024-05-08T00:26:12Z")

</div>

I did remember writing about how to read the `STASHC` but couldn’t find it.

I’ve now looked again, and I found what I wrote in an email that I sent, so no wonders I didn’t find it in any of the documents I wrote on my disk.

Here’s the gist of it:

The configuration on which fields to output is controlled by the `STASHC` file in the control directory. It’s a convoluted process, and it’s almost impossible not to mess this up if you’re not exactly sure what you’re doing.

The file begins with a namelist called `STASHNUM` which gives the number of entries for the requests, as well as the number of domain, time, and use profiles further down the list:

```auto
&STASHNUM NUM_REQ= 234, NUM_DOM=37, NUM_TIM=20, NUM_USE=12 /

```

Let’s go through this, for an example. “Fraction of Sea Ice after Timestep” (Timestep referring to a specific point in the calculation) has the stash code: Model 1, Section 0, Item 31. (You can look this up in the file `~access/umdir/vn7.3/cntldata/STASHmaster/STASHmaster_A` )

In the STASHC file, we see this namelist entry:

```auto
&STREQ IMOD= 1, ISEC=0, ITEM=31, IDOM=1, ITIM=3, IUSE=2 /

```

This says that we have a Stash Request for model 1, Section 0, Item 31, with Domain Profile 1, Time Profile 3, and Use profile 2.  
The numbers mean the profile in the order in the same file. Below the Stash Requests, there are the `TIME` domain namelists, and since we have Time Profile 3, we need to look at the 3rd entry:

```auto
&TIME NAME="TDMPMN", ITYP=3
INTV=1,UNT1="DU",IOFF=0,ISAM=1,UNT2="T ",
IOPT=1,
ISTR=1,IEND=-1,IFRE=1,UNT3="DU",
/

```

This gives us a human-readable name (which is good), and then a list of parameters that I don’t know from memory. But the name is TDMPMN , where the T stands for time profile, and DMPMN probably means ‘Dump Mean’, i.e. a mean over the dump period (which I think is a month).

After the time profiles we get the entries for the domain profiles, and we need the first one:

```auto
&DOMAIN NAME="DIAG", IOPL=5
PLT=0,
IOPA=1,
IMSK=1, IMN=0, IWT=0,
TS="N",TSNUM=0,
/

```

This describes which levels to output, whether theta or rho levels, etc.  
In this case, it’s a surface field, and `DIAG` is just that.  
For other fields, you might see names like `DALLTH` (all theta levels) or `DSOIL` (soil levels).  
Finally, we get to the Use Profiles, i.e. what to do with these fields. Usually it’s which files to output them to.

```auto
&USE NAME="Upamonth", LOCN=3
IUNT=60,
/

```

It goes into the PA file, which is meant for monthly outputs.

As far as I know, the `NAME` attributes of the profiles have zero impact on the model, and are just there so that we can read them easier.

If you do need to make changes, be careful. I would advise against changing the individual profiles, and only edit the STREQ entries. If you add or remove entries, make sure that at the beginning of the file in the `STASHNUM` section, you update the correct total number of entries.
