Stash code outputs in UKESM

Hi there,

I added a new species, Br2 (34261), along with some diagnoses to the UM-UKCA model (vn11.0). The model ran well (suite u-dv236). However, when I checked the outputs, I found that the added species were shown as “stash code =34261” rather than their descriptions or names given in the description, though their values are OK. Below are the relevant changes made in files STASHmaster_A:

1| 1 | 34 | 261 |Br2 MASS MIXING RATIO AFTER TS |

2| 2 | 0 | 1 | 1 | 2 | 40 | 11 | 0 | 0 | 0 | 0 |

3| 000000000000000000110000011000 | 00000000000000000001 | 1 |

4| 1 | 0 | -99 -99 -99 -99 -99 -99 -99 -99 -99 -99 |

5| 0 | 2240 | 0 | 65 | 0 | 0 | 0 | 0 | 0 |

and in STASHmaster-meta.conf:

[stashmaster:code(34261)]

description=Br2 MASS MIXING RATIO AFTER TSTEP

help=Br2 Mass Mixing Ratio in kg/kg(Air)

In addition, some corresponding changes were also made in section 51 for P-level output.

Please tell me where I need to change? Thanks!

Bests,

Xin

Hi Xin,

This is happening because Xconv is using a different STASHmaster_A file than you are in your suite. You can get Xconv to use your one by telling it to via the settings menu, or you can ignore it. The names in Xconv only match if the STASHmaster is correct - the thing that is important is the STASH code themselves, which are unique to your diagnostics.

I always advise people to be careful of the listings in Xconv, as while the names might be correct, the units might not etc. The best thing is to be sure that the STASH code has the correct data, and then treat it accodingly.

Best wishes,
Luke

Hi Luke,

Thanks for the prompt reply. It is a relief, as I thought I might have missed something. OK, I will ignore it, as any changes in Xconv will not affect the output files. I may have to manually update their descriptions for future use (mainly for others), but this is a pain if many tracers are added.
Best,

Xin

Hi Luke,

Sorry to bother you, but this is related to your previous reply about potential problems with Xconv. I have now encountered a strange issue when using Xconv to plot model outputs. For example, in u-dv236 outputs, all 3-hourly outputs look fine; however, the units of the same species in the monthly output file are incorrect (i.e. 1000 times larger!). The X1000 factor seems to apply to all variables. I am wondering whether this is only in visualisation or in the corresponding raw PP output files. If it is in the raw data, something (in STASHmaster?) must have gone wrong, and should be corrected to avoid misleading. Could you please provide a bit more detailed instructions on how to fix the code? Thanks!

Best,

Xin

Hi Xin,

By units being x1000 larger, do you mean that the values in the field are fine, but Xconv lists them as being e.g. kg rather than g, or that the values themselves are x1000 larger than expected?

In the STASH panel, what output stream (e.g. UPA etc.) are these going to and what is the time profile (e.g. TMONMN, T1HR etc.)?

Many thanks,
Luke

Hi Luke,

It seems the problem is on my side; in another run (with a “tidyup” of Metadata), the outputs do not exhibit the same problem in Xconv. So it is likely this is just a false alarm.

Best,

Xin