Ahh. I was afraid this would come back to bite me. This is where it gets complicated (even more than it currently is). Very many apologies for the long message which follows.
So the orography that your version of the ancillary-creation suite (u-dn772) creates, at /work/n02/n02/grenvill/cylc-run/u-dn772/share/data/n96e_orca025_go6/orography/globe30/qrparm.orog_from_herold is, as you can easily see, not using a modern land sea mask, but rather one that is appropriate for about 50 million years ago (the early Eocene), which is what I was originally working on several years ago. The suite uses as its input /work/n02/n02/grenvill/wilro.2021-03-31/TGRES/herold_orog.new which, as you can see, is using the same land sea mask but is at very very high spatial resolution. This is to account for any sub-grid scale variability in the orography. So somehow the suite is designed to read this in and interpolate it to the standard model grid of 192,144.
The problem is this file was created several years before I started working on this, even before I moved to Bristol. It basically reads in very high resolution (satellite derived, ARC-30 second) data showing the modern orography, and modifies it to be appropriate for the Eocene. Although I have the original Python, I have no idea where the original data comes from - the person who did this has long since left Bristol, and their home directory (and all others) has been deleted. So I don’t know how to do this step. I have contacted the person, but so far waiting for a reply.
In order to run my version of the ancillary-creation suite (u-do273), instead of using this high resolution orography as input to the suite, I instead used the standard/original orography file (at /work/y07/shared/umshared/ancil/atmos/n96e/orca1/orography/globe30/v6/qrparm.orog) which, as you can see, is already on the standard model grid. Although this works, and the suite generates the equivalent orography file to yours (at /work/n02/n02/cjrw09/cylc-run/u-do273/share/data/n96e_orca025_go6/orography/globe30/qrparm.orog), all but 2 of the fields are completely blank i.e. full of zeros. The only fields in this file that look okay are “Silhouette orographic roughness” and “Half of (peak to trough)…”. And a couple of the other fields make xconv crash. So clearly something is going wrong with the ancillary-creation suite, and it doesn’t like using as it’s input a file that is already at low resolution.
Instead, when running my actual suite (u-do321) I pointed the orography to the standard/original orography file (at /work/y07/shared/umshared/ancil/atmos/n96e/orca1/orography/globe30/v6/qrparm.orog). The reason I thought this would be okay is because the orography file is a filled field i.e. it does not contain any land sea mask. So I was hoping that, for testing purposes, I could get away with using this original file, because it shouldn’t conflict with my new land sea mask. In the same way, a lot of my other ancillaries (e.g. the aerosols) are also currently set to the standard/original versions, because they are also filled fields and therefore don’t conflict with my new land sea mask.
Obviously, eventually, when it comes to running this properly, I will need to solve the above problem and create my own orography. Although my new land sea mask is very similar to the standard/original version, it is not the same. But for now, I was just wanting to see if the model would run with all of these new ancillaries.
But clearly I was wrong in using the standard/original orography. Is it possible that, even though it is a filled field, that the model is crashing because e.g. it is seeing a value for orography where the land sea mask says it is ocean, or vice versa? So it is seeing a high value in the orography (e.g. a mountain) but then finding no land in the new land sea mask? As I said, my new land sea mask is very similar to the standard (which would match the orography), but not exactly the same e.g. Australia is attached to Papua New Guinea.
If this is the case, I can’t even try running with your version of the orography, because it will have the same problem i.e. your version is early Eocene, so again it will find e.g. high values for orography where there is no land, or vice versa.
So clearly I have 2 problems here. Either the ancillary-creation suite has a problem, which is why it is creating a load of zeros when making the orography. Or it was never designed to read in orography that was already on the model grid, but only works with a much higher resolution version of the orography. In which case, I need to go back and try harder to get hold of the person who originally created this high resolution file for the Eocene, in order to replicate it for my new land sea mask.
Does that make sense? Are you able to see anywhere in my ancillary-creation suite (u-do273) where it specifies that the input orography needs to be at very high resolution?
Thank you,
Charlie