[gradsusr] xdfopen with multiple files
Kiswendsida Hyacinthe Guigma
K.Guigma at sussex.ac.uk
Fri Apr 13 03:47:10 EDT 2018
Thank you Jeff for these ideas. But none of them works.I am allowed to use up to 15 characters while I'm actually using only 14. And I removed several names to keep only 2 of them, but I still get the same error message. In the past I used similar length of names (24) with binary files (then using open) and it worked fine. So I don't think this could be the cause of the issue.
GUIGMA
________________________________________
From: gradsusr [gradsusr-bounces at gradsusr.org] on behalf of Jeff Duda [jeffduda319 at gmail.com]
Sent: Thursday, April 12, 2018 8:52 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] xdfopen with multiple files
Could it be that there is a length or character limit for entries in a control file and that your edef entry is exceeding that limit? What happens if you remove several of those names from the list in the EDEF line? If it works after removing several entries then that would suggest it is indeed a length limit issue.
Jeff Duda
On Thu, Apr 12, 2018 at 12:11 PM, Kiswendsida Hyacinthe Guigma <K.Guigma at sussex.ac.uk<mailto:K.Guigma at sussex.ac.uk>> wrote:
Hi grads-users!
I have many Netcdf files into one folder that I want to plot using grads. Each file contains 111 variables and the X and Y dimensions are the same for them all. They don't contain time nor level dimension. ncdump gives me the following result:
netcdf IND2_ERA_1_1_HImax_4d_95perc_stat {
dimensions:
lon = 71 ;
lat = 31 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
lon:long_name = "lon" ;
lon:axis = "X" ;
double lat(lat) ;
lat:units = "degrees_north" ;
lat:long_name = "lat" ;
lat:axis = "Y" ;
float int(lat, lon) ;
int:units = "deg_C" ;
int:_FillValue = 1.e+32f ;
int:long_name = "Mean intensity" ;
float int1(lat, lon) ;
int1:units = "deg_C" ;
int1:_FillValue = 1.e+32f ;
int1:long_name = "Mean intensity in January" ;
float int2(lat, lon) ;
int2:units = "deg_C" ;
int2:_FillValue = 1.e+32f ;
int2:long_name = "Mean intensity in February" ;
float int3(lat, lon) ;
int3:units = "deg_C" ;
int3:_FillValue = 1.e+32f ;
int3:long_name = "Mean intensity in March" ;
float int4(lat, lon) ;
.......}
It is possible for me to read each one of them using sdfopen. But when it comes to open them with xdfopen using the TEMPLATE option, it gives me this error message:
gadxdf error: Missing or invalid dimension size.
--> The invalid description file record is:
--> edef 24 names min_3d_90perc min_4d_90perc min_5d_90percperc min_6d_90perc min_3d_95percperc min_4d_95perc min_5d_95perc min_6d_95perc min_3d_99perc min_4d_99perc min_5d_99perc min_6d_99perc max_3d_90perc max_4d_90perc max_5d_90perc max_6d_90perc max_3d_95perc max_4d_95perc max_5d_95perc max_6d_95perc max_3d_99perc max_4d_99perc max_5d_99perc max_6d_99perc
The data file was not opened.
This is my ctl script:
DSET /research/geog/data2/DATA/kg312/Fortran/IND2_ERA_1_1_%e_stat.nc<http://e_stat.nc>
DTYPE NETCDF
OPTIONS template
EDEF 24 NAMES min_3d_90perc min_4d_90perc min_5d_90percperc min_6d_90perc min_3d_95percperc min_4d_95perc min_5d_95perc min_6d_95perc min_3d_99perc min_4d_99perc min_5d_99perc min_6d_99perc max_3d_90perc max_4d_90perc max_5d_90perc max_6d_90perc max_3d_95perc max_4d_95perc max_5d_95perc max_6d_95perc max_3d_99perc max_4d_99perc max_5d_99perc max_6d_99perc
VARS 111
int 1 99 int
int1 1 99 int1
int2 1 99 int2
int3 1 99 int3
....
ENDVARS
I really don't understand that. Is is because my netCDF files are 2D (if this was a problem I don't think I would be able to open them using sdfopen)?
Your help is more than welcome.
Cheers
GUIGMA
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr
--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Division
Boulder, CO
More information about the gradsusr
mailing list