Error activating new UKCA tracers

Hello,

I am trying to move a scheme I implemented into UM vn13.4 to UM vn13.9 for running on Monsoon3. The scheme includes 9 new UKCA tracers. I am getting the following error when running the new branch at vn13.9 with the scheme included (https://code.metoffice.gov.uk/trac/um/browser/main/branches/dev/isabellesangha/vn13.9_PSC_scheme) on Monsoon3 using a copy of the vn13.9 test suite (https://code.metoffice.gov.uk/trac/roses-u/log/d/r/8/8/0/):

???
???!!!???!!!???!!!???!!!???!!! ERROR ???!!!???!!!???!!!???!!!???!!!
? Error code: 1
? Error from routine: UKCA_MAKE_TR_MAP
? Error message: Missing 9 of 218 required tracers
? Error from processor: 95
? Error number: 48
???

I’ve added some print statements and can see that the tracers aren’t being activated, but I can’t figure out why. (option codes are currently set to zero for debugging, but I get the same error even when I adjust the option codes).

Any help would be greatly appreciated!

Thank you :slight_smile:

Hi Isabelle,

It’s a long shot but when you committed the branch, did you wait a few minutes for it to sync with the mirror on MOSRS before running the suite?

Cheers,

James

Hi James,

Yes, I did. I committed it last week when I was first trying to debug the error and it’s still occurring today.

Thanks though (I also thought of this when it first happened!)

Cheers,

Izzy

In the pe_output file, can you see your additional section 34 tracers listed with the other section 34 tracers?

James

Yeah they seem to be.
In the /cylc-run/u-dr880/run2/work/1/atmos_cristrat_rigorous_16x8/pe_output/atmos.fort6.pe000 file (for example), they are displayed with the other tracers in sections 34 when listing the prognostic fields:

0 86 27648 1 34260 21313 RU10NO3 MASS MIXING RATIO AFTER TS 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34271 21399 NAT BIN 1 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34272 21485 NAT BIN 2 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34273 21571 NAT BIN 3 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34274 21657 NAT BIN 4 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34275 21743 NAT BIN 5 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34276 21829 NAT BIN 6 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34277 21915 NAT BIN 7 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34278 22001 NAT BIN 8 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34279 22087 NAT BIN 9 (VMR) 2008: 9: 1: 0: 0: 0 -32768
0 86 27648 1 34863 22173 DHPR12O2 MASS MIXING RATIO AFTER T 2008: 9: 1: 0: 0: 0

Another long shot while I think about this further, if you did this all via the terminal, can you activate the GUI (if one exists in Monsoon3) and see if it spots any errors?

I have had the “Missing X of Y required tracers” before and it was because stuff wasn’t being picked up in the new stashmaster but your changes look fine.

If this is a test suite, are all the atoms tasks failing?

Cheers,

James

I’ve been using the cylc TUI on Monsoon3 to visualise what is failing and it seems to be all the atmos tasks (from looking at the job.err file of each atmos task).

I’ve also had this error previously when working with version 13.4, but it was an issue of a missing line in the rose-suite.conf file that was preventing the stashmaster file from being found. I’ve gone through and checked all files related to the stashmaster and can’t see anything wrong which is why I am a bit stumped :frowning:

Isabelle

Please allow us read access to your roses directory and u-dr880.

Grenville

Hi Grenville,

They should both have read access now. Sorry about that!

Thanks,
Isabelle

Not yet:

[grenville.lister.ext@cazccylc1 users]$ ls -ld /home/users/isabelle.sangha.ext/
drwx------. 26 isabelle.sangha.ext isabelle.sangha.ext 4096 Aug 27 15:31 /home/users/isabelle.sangha.ext/

run this:

chmod -R a+rX /home/users/isabelle.sangha.ext

reset permissions on files where required (things you want to be hidden)

Grenville

Sorry – it should be all set now for permissions!

what was the UM 13.4 suite?

There shouldn’t be a UM vn13.4 suite in my Monsoon3 directory. I previously (on Monsoon2) implemented this scheme in vn13.4 and had it working using suites such as u-dd630 (d/d/6/3/0 (log) – roses-u).

The error I am getting now with the tracers is using a copy of the test suite at vn13.9 on Monsoon3 (u-dr880). I am trying to implement the same scheme that I put in UM vn13.4 to UM vn13.9.

Isabelle

The Stashmaster file for the 9 tracers has value “0” for all option codes - that appears to be the reason why the tracers are not on. Please look at the explanation in

/home/users/isabelle.sangha.ext/cylc-run/u-dr880/run2/share/fcm_make/extract/um/src/control/ukca_interface/tstmsk_ukca_mod.F90

and see line 137 for what the option codes need to be for specific chemistry schemes.

This may be of help too Controlling UKCA tracers from vn9.0 - UKCA

Grenville

Hi Grenville,

I had set the stashcodes to zero for debugging, but I still get the error when I have them set to what I used in vn13.4 (see cylc-run/u-dr880/run4/log/job/1/atmos_cristrat_rigorous_16x8/01/job.err).

I’m not sure what else to test because in a side-by-side comparison my STASHmaster file looks the same as the STASHmaster in version 13.4 except that I have had to move the tracers up one so that they now start at 34271 instead of 34270.

Thanks for your help with this!

Best wishes,

Izzy

Hi Isabelle

Sorry if I’m being dense, but where can I see what was used in the successful UM 13.4 job?

Grenville

from here I suuppose: fcm:um.xm_br/dev/isabellesangha/vn13.4_nat_psc_mods

Hi Grenville,

Yes - that should be the up-to-date UM vn13.4 branch that successfully ran with the new tracers. It was run in many suites for different simulations, but an example suite that worked with the branch is u-dd630.

Thanks again for taking a look at this! Are there any other tests I could be running that would help diagnose the problem?

Isabelle

I’m obviously no UKCA expert, but notice that the list of the names of tracers in tracer_varnames (appearing in here for example tr_lookup(i) = arrayloc(nm_spec_active, tracer_varnames(i)) does not include BIN_1,… BIN_9.

These are the names it has:

O3 NO NO3 N2O5 HO2NO2 HONO2 H2O2 CH4 CO HCHO MeOOH HONO C2H6 EtOOH MeCHO PAN C3H8 i-PrOOH EtCHO Me2CO PPAN MeONO2 C5H8 MPAN HCOOH MeCO3H MeCO2H Cl ClO Cl2O2 OClO Br BrCl BrONO2 N2O HOCl HBr HOBr ClONO2 CFCl3 CF2Cl2 MeBr N O(3P) C2H4 C4H10 oXYLENE TBUT2ENE APINENE BPINENE C2H2 BENZENE H2 DMS SO2 H2SO4 MSA DMSO NH3 CS2 COS H2S H OH HO2 MeOH Sec_Org C3H6 SO3 MEK TOLUENE NO2 BrO HCl Nuc_SOL_N Nuc_SOL_SUAit_SOL_N Ait_SOL_SUAit_SOL_BCAit_SOL_OMAcc_SOL_N Acc_SOL_SUAcc_SOL_BCAcc_SOL_OMAcc_SOL_SSCor_SOL_N Cor_SOL_SUCor_SOL_BCCor_SOL_OMCor_SOL_SSAit_INS_N Ait_INS_BCAit_INS_OMNuc_SOL_OMEtOH i-PrOH n-PrOH HOCH2CHO HOC2H4OOH EtCO3H HOCH2CO3H NOA EtONO2 AGE OF AIRi-PrONO2 MeO2NO2 HOC2H4NO3 PHAN MeSCH2OO MeS MeSO MeSO2 MeSO3 MSIA CARB14 CARB17 CARB11A CARB7 CARB10 CARB13 CARB16 UCARB10 CARB3 CARB6 CARB9 CARB12 CARB15 UCARB12 NUCARB12 UDCARB8 UDCARB11 UDCARB14 TNCARB26 TNCARB10 RN10NO3 RN13NO3 RN16NO3 RN19NO3 RA13NO3 RA16NO3 RA19NO3 RTX24NO3 RN10OOH RN13OOH RN16OOH RN19OOH RN8OOH RN11OOH RN14OOH RN17OOH RU14OOH RU12OOH RU10OOH NRU14OOH NRU12OOH RN9OOH RN12OOH RN15OOH RN18OOH NRN6OOH NRN9OOH NRN12OOH RA13OOH RA16OOH N2O HOCl HBr HOBr ClONO2 CFCl3 CF2Cl2 MeBr N O(3P) C2H4 C4H10 oXYLENE TBUT2ENE APINENE BPINENE C2H2 BENZENE H2 DMS SO2 H2SO4 MSA DMSO NH3 CS2 COS H2S H OH HO2 MeOH Sec_Org C3H6 SO3 MEK TOLUENE NO2 BrO HCl Nuc_SOL_N Nuc_SOL_SUAit_SOL_N Ait_SOL_SUAit_SOL_BCAit_SOL_OMAcc_SOL_N Acc_SOL_SUAcc_SOL_BCAcc_SOL_OMAcc_SOL_SSCor_SOL_N Cor_SOL_SUCor_SOL_BCCor_SOL_OMCor_SOL_SSAit_INS_N Ait_INS_BCAit_INS_OMNuc_SOL_OMEtOH i-PrOH n-PrOH HOCH2CHO HOC2H4OOH EtCO3H HOCH2CO3H NOA EtONO2 AGE OF AIRi-PrONO2 MeO2NO2 HOC2H4NO3 PHAN MeSCH2OO MeS MeSO MeSO2 MeSO3 MSIA CARB14 CARB17 CARB11A CARB7 CARB10 CARB13 CARB16 UCARB10 CARB3 CARB6 CARB9 CARB12 CARB15 UCARB12 NUCARB12 UDCARB8 UDCARB11 UDCARB14 TNCARB26 TNCARB10 RN10NO3 RN13NO3 RN16NO3 RN19NO3 RA13NO3 RA16NO3 RA19NO3 RTX24NO3 RN10OOH RN13OOH RN16OOH RN19OOH RN8OOH RN11OOH RN14OOH RN17OOH RU14OOH RU12OOH RU10OOH NRU14OOH NRU12OOH RN9OOH RN12OOH RN15OOH RN18OOH NRN6OOH NRN9OOH NRN12OOH RA13OOH RA16OOH RA19OOH RTN28OOH NRTN28OOH RTN26OOH RTN25OOH RTN24OOH RTN23OOH RTN14OOH RTN10OOH RTX28OOH RTX24OOH RTX22OOH NRTX28OOH RAROH14 RAROH17 RTN26PAN CCARB12 TNCARB15 RCOOH25 TXCARB24 TXCARB22 RN9NO3 RN12NO3 RN15NO3 RN18NO3 RU14NO3 RTN28NO3 RTN25NO3 RTX28NO3 RTX22NO3 AROH14 ARNOH14 AROH17 ARNOH17 ANHY DHPCARB9 HPUCARB12 HUCARB9 IEPOX HMML DHPR12OOH DHCARB9 RU12NO3 RU10NO3

Maybe you know where tracer_varnames gets set?

Grenville

Isabelle

I can get the atmos_strat_col_16x8 case to run, but only if I change BIN_1, BIN_2 etc to BIN1, BIN2

I note that in

…/ukca_interface/ukca_nmspec_mod.F90

there is no reference to BIN_??

Grenville