@zender thanks for detailed feedback. I will run more tests next week and get back with you on the results. My guess in relation to fill and miss value is that you could have an ocean or terrestrial model that doesn't simulate anything on land or viceversa. Maybe it's for the ocean/land mask. Just a guess.
Hi again, I managed to install 5.3.6. I tried to skip removing the string variables but I still get a segmentation fault. Caught signal 11 (Segmentation fault: invalid permissions for mapped object at address 0x7ff80903b39f) ==== backtrace (tid: 8619) ==== 0 /home/a/miniconda3/envs/ncoenv/bin/../lib/./././libucs.so.0(ucs_handle_error+0x2fd) [0x7f91fe33e84d] 1 /home/a/miniconda3/envs/ncoenv/bin/../lib/./././libucs.so.0(+0x2fa3f) [0x7f91fe33ea3f] 2 /home/a/miniconda3/envs/ncoenv/bin/../lib/./././libucs.so.0(+0x2fc0a)...
Hi Charlie, Thanks for looking into this. I have version 5.2.1 sudo apt install nco Reading package lists... Done Building dependency tree... Done Reading state information... Done nco is already the newest version (5.2.1-1build2). 0 upgraded, 0 newly installed, 0 to remove and 207 not upgraded. $ ncks --version NCO netCDF Operators version 5.2.1 "Shabu Shabu" built by buildd on lcy02-amd64-080 at Apr 1 2024 07:01:37 ncks version 5.2.1 I am on WSL on Windows 10. Does this mean I have to build my...
Hi All, I am trying to perform several operations with .nc files associated to the latest IPCC report AR6 and that are used in the IPCC Atlas. The files are here https://digital.csic.es/handle/10261/332744 and the problem is the same on all the files. What I am trying to do is for example selecting only March over a period of consecutive year, and averaging over these years and over the ensemble models. Another operation is averaging over consecutive years and ensemble. If I use the command ncwa...
Hi All, I am trying to perform several operations with .nc files associated to the latest IPCC report AR6 and that are used in the IPCC Atlas. The files are here https://digital.csic.es/handle/10261/332744 and the problem is the same on all the files. What I am trying to do is for example selecting only March over a period of consecutive year, and averaging over these years and over the ensemble models. Another operation is averaging over consecutive years and ensemble. If I use the command ncwa...
Thanks that works!
I can't find in the -d documentation if there is a way to subset data at irregular interval. I would like to do something like this: ncrcat -d depth,0-10,13,14 tmp1.nc tmp2.nc Or do I have to a for loop and use a manually built index and than concatenate? Thanks
I was thinking that the z axis represents levels, is that not so? If you have x,y,z,time then you can also have levels? So a 5 dimension file? Because when I try to use "cdo -selzaxis,3/10 myfile.nc test.nc" it also doesn't work. So how can I select certain slices on the z axis? I am really confused. EDIT: I think I see the problem now, 'veget' is the time counter, so my question is how do I set 'veget' as vertical level instead of time counter?
I was thinking that the z axis represents levels, is that not so? If you have x,y,z,time then you can also have levels? So a 5 dimension file? Because when I try to use "cdo -selzaxis,3/10 myfile.nc test.nc" it also doesn't work. So how can I select certain slices on the z axis? I am really confused.
Sorry Charlie I didn't explain myself clearly. It seems CDO is not recognizing 'veget' as the z dimension. For example when I try to select one level I get: cdo -sellevel,1 myfile.nc test.nc Warning (cdf_set_dim): Inconsistent dimension definition for veget! dimid = 0; type = 4; newtype = 3 Warning (cdf_read_mapping_atts): Text attribute crs:grid_mapping_name missing! setBaseTime : Unsupported TIMEUNIT: -! Warning: Level 1 not found! Normally my dimension are like this: dimensions: lat = 720 ; lon...
I have the following netcdf file dimensions: lon = 1440 ; lat = 720 ; veget = UNLIMITED ; // (15 currently) variables: int crs ; crs:proj4 = "+proj=longlat +datum=WGS84 +no_defs" ; double lon(lon) ; lon:units = "degrees_east" ; lon:long_name = "Longitude" ; lon:valid_min = -180. ; lon:modulo = 360. ; lon:axis = "X" ; lon:valid_max = 180. ; lon:topology = "circular" ; double lat(lat) ; lat:units = "degrees_north" ; lat:long_name = "Latitude" ; lat:valid_min = -90. ; lat:valid_max = 90. ; lat:axis...
Thanks Henry, that this logic apply to all index aside from lat/long? If I want the equator I use lat=0?
I am confused about the indexing of levels and band, does it starts at zero or 1. Where can I find this explained in the help guide? It seems it starts at 0 because I get this error "message ncap2: ERROR ncap_lmt_evl(): Lower limit 15 for dim veget is outside range 0-14" Does this apply to all nco commands? Is there a reason for this, being used to all other netcdf tools I find it confusing. I find it confusing also because when I look at the file info the levels start at 1. Thanks
I am confused about the indexing of levels and band, does it starts at zero or 1. Where can I find this explained in the help guide? It seems it starts at 0 because I get this error "message ncap2: ERROR ncap_lmt_evl(): Lower limit 15 for dim veget is outside range 0-14" Does this apply to all nco commands? I find it confusing because when I look at the file info the levels start at 1. Thanks