sellonlatbox issue
Added by Justin Mankin almost 9 years ago
Hi,
I've searched to identify the issue, but I'm stuck.
I've run the sellonlatbox command just fine on CESM land component output, but am running into trouble with the atmospheric component output. I have a sense it's a grid issue, but I can't diagnose it. I'll provide two examples below and my CDO version:
My CDO version is:
Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo)
Compiled: by unknown on c672c4f7589c (x86_64-unknown-linux-gnu) Mar 6 2016 21:40:54
Compiler: gcc -std=gnu99 -g -O2 -fopenmp
version: gcc (GCC) 4.8.5
Features: DATA PTHREADS OpenMP HDF5 NC4/HDF5 OPeNDAP Z PROJ.4 SSE2
Libraries: HDF5/1.8.15 proj/4.91
Filetypes: srv ext ieg grb nc nc2 nc4 nc4c
CDI library version : 1.7.0 of Mar 6 2016 21:40:00
CGRIBEX library version : 1.7.3 of Sep 14 2015 10:58:44
netCDF library version : 4.3.3.1 of Feb 22 2016 16:27:42 $
HDF5 library version : 1.8.15
SERVICE library version : 1.4.0 of Mar 6 2016 21:39:46
EXTRA library version : 1.4.0 of Mar 6 2016 21:39:40
IEG library version : 1.4.0 of Mar 6 2016 21:39:44
FILE library version : 1.8.2 of Mar 6 2016 21:39:41*
On the file QBOT.nc I run the following to select North America, which works beautifully:
cdo sellonlatbox,-167,-52,5,83 QBOT.nc QBOT.NA_selection.nc
However, when I run the same command on the file named TREFHT.nc, or any CAM (atmospheric model) output, as this:
cdo sellonlatbox,-167,-52,5,83 TREFHT.nc TREFHT.NA_selection.nc
I receive the following error:
Warning (genlonlatbox_reg) : Allocation of 0 bytes! [ line 350 file Selbox.c ] Segmentation fault (core dumped)
The SINFO command on the TREFHT.nc file gives:
cdo sinfo TREFHT.nc
File format : netCDF4
-1 : Institut Source Ttype Levels Num Points Num Dtype : Parameter ID
1 : unknown CAM instant 1 1 55296 1 F32 : -1
2 : unknown CAM instant 1 1 1 2 F64 : -2
3 : unknown CAM instant 1 1 1 2 F64 : -3
4 : unknown CAM instant 1 1 1 2 I32 : -4
5 : unknown CAM instant 1 1 1 2 I32 : -5
6 : unknown CAM instant 1 1 1 2 F64 : -6
7 : unknown CAM instant 1 1 1 2 F64 : -7
8 : unknown CAM constant 1 1 192 3 F64 : -8
9 : unknown CAM instant 1 1 1 2 F64 : -9
10 : unknown CAM instant 1 1 1 2 I32 : -10
11 : unknown CAM constant 1 1 192 3 I32 : -11
12 : unknown CAM instant 1 1 1 2 I32 : -12
13 : unknown CAM instant 1 1 1 2 I32 : -13
14 : unknown CAM instant 1 1 1 2 F64 : -14
15 : unknown CAM constant 1 1 191 4 F64 : -15
16 : unknown CAM constant 1 1 192 3 I32 : -16
Grid coordinates :
1 : lonlat : points=55296 (288x192)
lon : 0 to 358.75 by 1.25 degrees_east circular
lat : -90 to 90 by 0.942408 degrees_north
2 : generic : points=1
3 : lonlat : points=192
lat : -90 to 90 by 0.942408 degrees_north
4 : lonlat : points=191
slat : -89.5288 to 89.5288 by 0.942408 degrees_north
Vertical coordinates :
1 : surface : levels=1
Any guidance would be greatly appreciated!
thanks,
justin
Replies (1)
RE: sellonlatbox issue - Added by Jaison-Thomas Ambadan almost 9 years ago
Hi,
have a look at "ncdump -h TREFHT.nc"
netcdf TREFHT {
dimensions:
lon = 288 ;
lat = 192 ;
slat = 191 ;
time = UNLIMITED ; // (1 currently)
bnds = 2 ;
variables:
...
...
At the moment, CDO don't support non-standard dimensions like "slat". In your input file TREFHT.nc, there is only one variable TREFHT with lat-lon dimensions. so just select the variable before calling "sellonlatbox" - it should work.
cdo sellonlatbox,-167,-52,5,83 -selvar,TREFHT TREFHT.nc TREFHT.NA_selection.nc
Cheers,
J