Project

General

Profile

struggling with intlevel3d

Added by Catherine Drinkorn about 5 years ago

Hi all,

I am trying to interpolate from a hybrid model output depths (i.e. 0 > depth > -1) to actual depths. The original horizontal grid is curvilinear but using remapbil I can easily convert to something regular/rectangular.
I tried to make use of intlevel3d but I got error like unsupported file structure, too many variables or no missing values allowed. The way I understood it intlevel3d is exactly what I need to do what I described above. If I am wrong I'd be happy for any hint about another possibility. If I am right can someone please give me further instructions on how to use this operator?

Here is the sinfo of my icoordinate file

File format : NetCDF4
-1 : Institut Source T Steptype Levels Num Points Num Dtype : Parameter ID
1 : unknown unknown c instant 35 1 1 1 F64 : -1
2 : unknown unknown v instant 1 2 1040041 2 F32 : -2
3 : unknown unknown v instant 1 2 1040041 2 F32 : -3
Grid coordinates :
1 : generic : points=1
2 : lonlat : points=1040041 (2881x361)
lon : 0 to 360 by 0.125 degrees_east
lat : 45 to 90 by 0.125 degrees_north
Vertical coordinates :
1 : generic : levels=35
s_rho : -0.9857143 to -0.01428571 by 0.02857143
2 : surface : levels=1
Time coordinate : 11 steps
RefTime = 1970-01-01 00:00:00 Units = seconds Calendar = proleptic_gregorian
YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss
1970-01-01 00:00:00 1970-01-01 00:45:00 1970-01-01 01:30:00 1970-01-01 02:15:00
1970-01-01 03:00:00 1970-01-01 03:45:00 1970-01-01 04:30:00 1970-01-01 05:15:00
1970-01-01 06:00:00 1970-01-01 06:45:00 1970-01-01 07:30:00
cdo sinfo: Processed 3 variables over 11 timesteps [0.04s 31MB]

My variable is a 3D variable depending on s_rho in the vertical. But how is the target vertical variable file supposed to look like? No vertical coordinate I extract from any other dataset seems to be supported by cdo...

Thank you for any hint!

// Cate


Replies (1)

RE: struggling with intlevel3d - Added by Karin Meier-Fleischer about 5 years ago

Hi Cate,

without your data I can only guess why you're getting an error.
Your data file contains variables depending on two different grids.

Grid coordinates :
1 : generic : points=1
2 : lonlat : points=1040041 (2881x361)

CDO documentation:
If there are more then one grid in the data file the grid description of the first variable will be used. Add the extension :N to the name of the datafile to select grid number N.

You can either choose the grid to be used or you select the wanted variable before using intlevel3d.

-Karin

    (1-1/1)