Opened 3 months ago

Closed 11 days ago

#3117 closed help (fixed)

Cannot create land surface fraction ancillary suitable for MORUSES

Reported by: nx902220 Owned by: jeff
Component: Land Surface Modelling Keywords:
Cc: Platform:
UM Version:

Description

Hi,

I am taking the qrparm.veg.frac file which contains 10 pseudo levels (corresponding to tiles) with information on the land surface fractions.

I convert it to netcdf using xconv. I then use Python to edit the urban land fraction (I'm doing an experiment). I then use xancil to recreate the ancillary (hopefully with exact same metadata just updated fractions).

My suite u-bc220 fails with error:

LOOKUP TABLE
640 64-bit words long
c_io ( 13):Close
IO: Close: /projects/phdcase/lblunn/ANCILS/MORUSES/LONDON_SUITE_ANCILS/500m/qrparm.veg.frac on unit 13
Land surface configuration does not match /projects/phdcase/lblunn/ANCILS/MORUSES/LONDON_SUITE_ANCILS/500m/qrparm.veg.frac
Ancillary configuration 1 2 3 4 5 6 7 8 9 10
Namelist configuration 1 2 3 4 5 7 8 9 601 602

????????????????????????????????????????????????????????????????????????????????
???!!!???!!!???!!!???!!!???!!! ERROR ???!!!???!!!???!!!???!!!???!!!
? Error code: 99216
? Error from routine: inancila_rcf_inancila
? Error message: Land surface configuration does not match ancillary
? Error from processor: 0
? Error number: 0
????????????????????????????????????????????????????????????????????????????????

This makes sense. The original qrparm.veg.ancil ancillary has configuration

1 2 3 4 5 7 8 9 601 602

and this is what the namelist is expecting (the suite then works).

My new qrparm.veg.ancil has configuration

1 2 3 4 5 6 7 8 9 10

and this is not what the namelist is expecting and fails.

The stash code of the ancillary is 216 i.e FRACTIONS OF SURFACE TYPES. When I use xancil I tell it to make the new ancillary have stash code 216. It seems the default is to output the ancillary with the old 1 tile scheme configuration (1 2 3 4 5 6 7 8 9 10). Please can you advise me on how to output an ancillary with the new MORUSES configuration (1 2 3 4 5 7 8 9 601 602)?

My xancil job is contained in:
/home/d04/lblunn/ancils/qrparm_no_urban_500m.job
and the output ancillary from the job is:
/projects/phdcase/lblunn/ANCILS/MORUSES/LONDON_SUITE_ANCILS/500m/qrparm.veg.frac

Best wishes,

Lewis

Change History (13)

comment:1 Changed 3 months ago by pmcguire

  • Owner pmcguire deleted
  • Status changed from new to assigned

comment:2 Changed 3 months ago by pmcguire

  • Owner set to jcole

comment:3 Changed 3 months ago by pmcguire

  • Owner changed from jcole to jeff

comment:4 Changed 3 months ago by jeff

Hi Lewis

To fix your problem I will need to modify the xancil source code. I'm working on a new release now so hopefully this will be done soon.

Jeff.

comment:5 Changed 2 months ago by nx902220

Hi Jeff,

Thanks. Do you think that will be within the next month?

Best wishes,

Lewis

comment:6 Changed 2 months ago by jeff

Yes, hopefully it won't take too long.

Jeff.

comment:7 Changed 2 months ago by nx902220

Hi Jeff,

Sorry for the slow reply, I was at AMS.

Thank you. Please let me know when the new release is ready.

Lewis

comment:8 Changed 2 months ago by jeff

Hi Lewis

I've put a development version of xancil on monsoon/nexcs (/home/d00/jecole/bin/xancil0.58). I've added the ability to either specify the pseudo levels or read them from the input netCDF file, in the Grid Configuration panel.

There are a few things to note when creating the ancil file.

Make sure you select the correct UM version in the Configuration panel, you had this set incorrectly before. It doesn't have to be an exact version but you had 4.5 which won't work.

It is best to create the input netCDF file using 64 bit reals, otherwise the UM complains about the ancillary file grid not matching its grid for certain fields.

You can use the Atmosphere Ancillary Files → Vegetation Fractions panel to create your ancillary file instead of using Generalised Ancillary Files.

Let me know if this works.

Jeff.

comment:9 Changed 8 weeks ago by nx902220

Hi Jeff,

Thanks.

I have tried doing what you say. I specify the pseudo levels to be 1,2,3,4,5,7,8,9,601,602.
My xancil job is in /home/d04/lblunn/ancils/qrparm_no_urban_500m_try2.job

I get the same type of error message:


Land surface configuration does not match /projects/phdcase/lblunn/ANCILS/MORUSES/LONDON_SUITE_ANCILS/500m/qrparm.veg.frac
Ancillary configuration 1 2 3 4 5 7 8 9 601 0
Namelist configuration 1 2 3 4 5 7 8 9 601 602


The ancillary configuration has updated from 1,2,3,4,5,6,7,8,9,10 to 1,2,3,4,5,7,8,9,601,0. For some reason the last pseudo level in the ancillary configuration is wrong i.e. 0.

Do you know why that might be? I have checked that my xancil job has 602 rather than 0 for the last pseudo level.

Thanks,
Lewis

comment:10 Changed 8 weeks ago by jeff

Hi Lewis

In the Vegetation Fractions panel you have set Number of surface types to 9, change this to 10.

Jeff.

comment:11 Changed 8 weeks ago by nx902220

Hi Jeff,

Yes I missed that one. It now runs past the error!

Thank you for sorting this.

I might need to modify my ancillaries a few more times. Is it ok if you leave /home/d00/jecole/bin/xancil0.58 where it is for a while or let me know if it moves?

Lewis

comment:12 Changed 8 weeks ago by jeff

Hi Lewis

Glad it's working now.

Yes I will leave xancil0.58 where it is, I might update it at some point but that shouldn't affect you.

Jeff.

comment:13 Changed 11 days ago by jeff

  • Resolution set to fixed
  • Status changed from assigned to closed
Note: See TracTickets for help on using tickets.