Forums » Support »
enlarge command not giving correct result
Added by S BR about 5 years ago
Dear Users,
I am restricting (sellonlatbox) the north Atlantic region from a global data. Later, I am trying to expand (enlarge) that domain to the global grid, but it is not giving correct result. Following is my syntax.
cdo sellonlatbox,-80,0,0,65 SST_grid_r128x64.nc SST_grid_r128x64_80W-0E_0-65N.nc
cdo enlarge,SST_grid_r128x64.nc SST_grid_r128x64_80W-0E_0-65N.nc SST_grid_r128x64_80W-0E_0-65N_enlarge.nc
The input data is attached here.
The CDO version is:
Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo)
Compiled: by root on e0002 (x86_64-unknown-linux-gnu) Nov 11 2015 14:54:22
Compiler: icc -std=gnu99 -fPIC -fopenmp
version: icc (ICC) 15.0.3 20150407
Features: DATA PTHREADS OpenMP4 HDF5 NC4/HDF5 Z PROJ.4 SSE2
Libraries: HDF5/1.8.15 proj/4.8
Filetypes: srv ext ieg grb grb2 nc nc2 nc4 nc4c
CDI library version : 1.7.0 of Nov 11 2015 14:52:45
CGRIBEX library version : 1.7.3 of Sep 14 2015 10:58:44
GRIB_API library version : 1.14.2
netCDF library version : 4.3.3.1 of Oct 23 2015 17:15:28 $
HDF5 library version : 1.8.15
SERVICE library version : 1.4.0 of Nov 11 2015 14:52:34
EXTRA library version : 1.4.0 of Nov 11 2015 14:52:29
IEG library version : 1.4.0 of Nov 11 2015 14:52:32
FILE library version : 1.8.2 of Nov 11 2015 14:52:30
Hi S,
when I understood you right then you want to generate a data set containing values of a sub-region and outside the sub-region the values are missing values.
Why don't you use one of the remap operators?
cdo remapycon,SST_grid_r128x64.nc SST_grid_r128x64_80W-0E_0-65N.nc SST_grid_r128x64_80W-0E_0-65N_global.nc
-Karin
Hi Karin,
Yes, you understood my point right but the given syntax gives some strange value north of 65N.
Please check the attached image.
Thank you.
Best,
S
With the given data it doesn't give strange values north of 65N. Which CDO version are you using?
cdo sellonlatbox,-80,0,0,65 SST_grid_r128x64.nc SST_grid_r128x64_80W-0E_0-65N.nc
cdo remapycon,$infile SST_grid_r128x64_80W-0E_0-65N.nc tmp_SST.nc
cdo grfill,device="png" -selvar,tos -seltimestep,1 tmp_SST.nc plot_SST
-Karin
Dear Karin,
I am using Version 1.7.0 (I have already mentioned in my first message). However, when I tried the given syntax with version 1.9.3 in another computer, it worked fine and it didn't show any strange value north of 65N. It looks like there was some issue with this syntax in Version 1.7.0.
However, the enlarge syntax fails in both the version 1.7.0 and 1.9.3. I was trying to use this syntax to solve the problem in different way.
Thank you for your support.
S
