xancil wants formula_terms @ vertical coordinate

[lboljka]

Hi

I was trying to create an ancillary file from my ozone netCDF file but received error:

doesn't contain the formula_terms attribute needed 
to calculate vertical level values.

At first it gave me an error about standard name but I fixed that, however I am not sure how to specify this "formula_terms". I have checked in cf-convention/xancil documentation and it only says I should specify formula_terms and in documentation it says:

formula_terms = "a: var1 b: var2 orog: var3".

But this does not help me much as I am not sure how to add it in the file. I have used ancillary ozone file from $UMDIR:

/work/n02/n02/hum/ancil/atmos/n96e/ozone/sparc/1994-2005/v2/qrclim.ozone_L85_O85

but apparently it does not have specified standard_name nor formula_terms.
All I did was converted that ancillary file to netCDF via xconv, created annual average and symmetrised data about Equator. Now I am trying to convert it back to ancillary file using xancil.

The netCDF file I created is:

archer$ /work/n02/n02/lboljka/startdump_ancillary/ozone_symmetric.nc

Then I added appropriate standard_name and created ozone_symmetric1.nc from that file. In that same folder there should also be the job I was using to create ancillary that failed with the above error.

Do you have any suggestions as to what is the best way to go about this problem?

[ros]

Hi Lina,

When you create a ticket please leave the owner as <default> otherwise we don't get alerted to the ticket being created.

I have now assigned this to the team so one of us will now pick it up.

Cheers,
Ros.

[lboljka]

Hi Ros

Oh, I didn't know that! Thank you very much!!!!

Best wishes

Lina

[jeff]

Hi Lina

It would be easier to use one of xancil's alternative methods for specifying vertical levels. In panel
Xancil→Configuration→Grid Configuration, select "Specify atmosphere vertical levels" and fill in the details of your model run. You will be able to find the location of the vertical namelist file in your umui output directory on puma. Look in file INITHIS under the variable name VERT_LEV. Further details on how to do this can be found in the xancil documentation

http://cms.ncas.ac.uk/documents/xancil/index.html#document-config

Jeff.

[lboljka]

Hi Jeff

Thank you. This worked like a charm! Though I could not get the VERT_LEV variable in INITHIS, I took the file that the model is using from the "/work/n02/n02/hum/vn8.6/ctldata/vert/" (vertlevs_L85_50t_35s_85km).

Thank you again!

Best wishes

Lina

[jeff]

Hi Lina

Good to hear it worked. It looks like they moved the VERT_LEV variable to file INITFILEENV for vn8.3 and later versions of the UM.

Jeff.