Opened 3 months ago

Last modified 6 weeks ago

#2536 new help

vn10.9 UKCA STASH issues

Reported by: admg26 Owned by: um_support
Priority: normal Component: UKCA
Keywords: Stash, Monsoon, nudging Cc:
Platform: Monsoon2 UM Version: 10.9

Description (last modified by willie)

Hi,

Monsoon

Rose suite: u-az303 running UKCA vn10.9

Output after 1 month run: moose:/crum/u-az303/

I am trying to output u,v and T on all model levels and every 6 hours so I can use the output file to nudge a different model run.

I am confused about how domains work.

In STASH, I changed the dom_name of these 3 variables in stream UPC to be 'DALLRH'. After the run, I seems to be getting all horizontal grid points (nx=192, ny=145) but only 17 vertical levels.

On the other hand,

V COMPNT OF WIND AFTER TIMESTEP, output to UPG 
dom_name=DALLRH 
nx=192, ny=145 and nz=85

but

U COMPNT OF WIND AFTER TIMESTEP, output to UPG 
dom_name=DALLRH 
nx=192, ny=144 (!) and nz=85

I am using xconv to look at the pp files.

Cheers,
Alison

Change History (10)

comment:1 Changed 2 months ago by willie

  • Description modified (diff)

Hi Alison,

Are you running this on Monsoon?

If you change the STASH then you need to run the STASH TidyStashTransform macro afterwards to re index the table.

I notice your suite has a large number of errors - these are the little red triangles in the left window. This might be because I'm looking on the wrong computer. But it is not a good idea to run the model with errors.

Regards
Willie

comment:2 Changed 2 months ago by admg26

Hi,

Yes this in on monsoon. I had forgotten to STASH TidyStashTransform?. Thanks!

The large number of errors is because I have added a number of stash items and they go away when I open the suite with -M path/to/correct/stashfile.

Thanks,
Alison

comment:3 Changed 2 months ago by admg26

Hi,

The STASH TidyStashTransform? has made the sizes of u and v consistent. I still have a problem of them being on 17 levels rather than 85.

On Monsoon, to edit the suite, I use

rose edit -M /home/d01/almin/vn10.9_cl2_photolysis/rose-meta

Cheers,
Alison

comment:4 Changed 2 months ago by grenville

Stream UPC doesn't output the same U,V,T as does UPG; it outputting

TEMPERATURE on P LEV/UV GRID
V COMPONENT of WIND on P LEV/UV GRID etc

are those what you want?

comment:5 Changed 2 months ago by admg26

Hello,

I am after T, U and V output every 6 hours on the right grid for use as input files to the nudging code. I am not sure what the right grid is.

Cheers,
Alison

comment:6 Changed 2 months ago by admg26

Grenville,

Your comment makes sense. I think what I want are the 1st 3 stash items

U COMPONENT of WIND after timestep
V COMPONENT of WIND after timestep
Theta after timestep

I have made a new stream and am outputting these every 6 hours.

Hopefully, this will work.

Cheers,
Alison

Last edited 2 months ago by admg26 (previous) (diff)

comment:7 Changed 2 months ago by luke

Hi Alison,

If you're making your own nudging files they will need to be identical (in terms of grid, metadata, filenames etc.) as the files contained here:

/projects/ukca-admin/analyses/era-in/

I've never tried nudging with UM model data, although it should be technically possible. You won't get identical dynamics however.

Thanks,
Luke

comment:8 Changed 2 months ago by admg26

Hi Luke,

I am looking at the nudging code in vn10.9 and there is an option to nudge to UM analysis data.

It looks like it reads in T and then converts it to theta. I am going to change that so it read in theta directly because that is what I get as UKCA output.

Cheers,
Alison

comment:9 Changed 2 months ago by willie

  • Platform set to Monsoon2

comment:10 Changed 6 weeks ago by admg26

Hi,

Just a quick update.

I have managed to output all the U, V and T fields I need from the model to do the nudging (in the right format and in the right netcdf files!).

I have also changed the nudging code so that UM analysis files can be read in.

The model is currently falling over but that is a different issue I think. Shall I open a different ticket and close this one?

Cheers,
Alison

Note: See TracTickets for help on using tickets.