Opened 3 years ago

Closed 3 years ago

Last modified 2 years ago

#2000 closed help (fixed)

making ancillary file for ENDGame grid with Xancil

Reported by: marcus Owned by: jeff
Component: UM Model Keywords: Xancil, ENDGame
Cc: Platform: ARCHER
UM Version: 10.4

Description

Hi, I am trying to create ancillary files to replace fields in a start dump. After generating the file with Xancil it turned out that the vertical level structure is incorrect, i.e. it starts with level 1 at 20.0 whereas level 86 is 0.0.
Please how can I create a correct ancillary file for UM 10.4?

All the relevant data can be found in /work/n02/n02/marcus/um/startdumps/ag292a.19600901/:

The netcdf files *_1960.nc contain the five input data fields for the ancillary file.

The Xancil job file is ag9292a.job and the created ancillary file is ag929a.long-lived_1960_rescale.anc.

I have specified the atmospheric level parameters to be taken from /work/n02/n02/hum/vn10.4/ctldata/vert/vertlevs_L85_50t_35s_85km.

I loaded the *.anc file in Python-Iris and saved the first cube as iris_out.nc.

Please if you could help this would be great.
Many thanks,

Marcus

Change History (11)

comment:1 Changed 3 years ago by jeff

  • Owner changed from um_support to jeff
  • Status changed from new to accepted

Hi Marcus

Just to let you know I'm looking into your problem. It seems that endgame vertical coordinates have changed slightly from the new dynamics ones, so xancil needs to be updated to handle this. I'll let you know when there is a new version of xancil available.

Jeff.

comment:2 Changed 3 years ago by jeff

Hi Marcus

I've built a test version of xancil (/home/n02/n02/jwc/bin/xancil_test) which will hopefully solve your problem.

In the Grid Configuration panel you need to set the number of model levels to 85, to match the levels in file /work/n02/n02/hum/vn10.4/ctldata/vert/vertlevs_L85_50t_35s_85km. Then in the Generalised Ancillary Files / Ancillary file 1 panel, for each field, you need to set the "Number of levels" to 86 and set "Enter first model level" to 0.

Let me know if this works and I will update xancil to this version.

Jeff.

comment:3 Changed 3 years ago by marcus

Hi Jeff,

Many thanks for this. I have tried this but I am encountering an error which I don't understand how to interpret:

/home2/n02/n02/jwc/bin/mkancil0.55 failed
forrtl: severe (19): invalid reference to variable in NAMELIST input, unit -4, file stdin, line 115, position 18
    while executing
"exec $exec_arg < $namelist > $output"
    (procedure "run_exec" line 24)
    invoked from within
"run_exec $exec"
    (procedure "bb_createancil" line 5)
    invoked from within
"bb_createancil"
    invoked from within
".create invoke"
    ("uplevel" body line 1)
    invoked from within
"uplevel #0 [list $w invoke]"
    (procedure "tk::ButtonUp" line 22)
    invoked from within
"tk::ButtonUp .create"
    (command bound to event)

Perhaps I have made an error in the configuration. The job file is stored here:
/work/n02/n02/marcus/um/startdumps/ag929a.19600901/ag929a_new.job

Best regards,
Marcus

comment:4 Changed 3 years ago by jeff

Hi Marcus

That was my fault, it was using the old mkancil0.55 executable. Try it again.

Jeff.

comment:5 Changed 3 years ago by jeff

  • Status changed from accepted to pending

comment:6 Changed 3 years ago by marcus

Hi Jeff,

Apologies for the delay, I had to leave in a rush and was on sick leave for a week.

As far as I can see your new version of the xancil executable now creates an ancillary file with 85 levels. That's great.

I was under the impression the the ENDGame grid had an additional surface level, i.e. 86 levels in total? Should the ancillary file have 85 levels? Could you perhaps shed a bit of light on this for me? Thanks a lot.

Best regards,
Marcus

comment:7 Changed 3 years ago by marcus

Hi Jeff,

Sorry my mistake: I had not carried out your instructions above with "number of levels" and "enter first model level". After I did those changes it now creates and ancillary file with 86 levels.

Sorry for the confusion. I think it's all fine now as far as I can see.

Many thanks,
Marcus

comment:8 Changed 3 years ago by jeff

  • Resolution set to fixed
  • Status changed from pending to closed

comment:9 Changed 2 years ago by luke

Hi Jeff,

Can I check, is this version of Xancil on Monsoon2, and if so, where?

Thanks,
Luke

comment:10 Changed 2 years ago by jeff

Hi Luke

This fix was in xancil0.56, that version and version 0.57 are on monsoon in my directory /home/d00/jecole/bin

Jeff.

comment:11 Changed 2 years ago by luke

Great! Many thanks.

Luke

Note: See TracTickets for help on using tickets.