[gradsusr] netcdf control file error

Roberto Mera rjmera at coas.oregonstate.edu
Wed Oct 6 12:28:03 EDT 2010


  Thanks, Jennifer. Here goes.

ncdump -c Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
netcdf nhqhga.pdg6dec.pnw {
dimensions:
     X = 110 ;
     Y = 117 ;
     Z = 1 ;
variables:
     int time(Z) ;
         time:long_name = "time" ;
         time:units = "month" ;
     float latitude(Y) ;
         latitude:long_name = "latitude" ;
         latitude:units = "degrees_north" ;
     float longitude(X) ;
         longitude:long_name = "longitude" ;
         longitude:units = "degrees_east" ;
     float Tmin(Z, Y, X) ;
         Tmin:long_name = "Monthly average of daily Tmin" ;
         Tmin:units = "C" ;
         Tmin:Fill_Value = -9999.f ;
     float Frost(Z, Y, X) ;
         Frost:long_name = "Number of frost days per month" ;
         Frost:units = "days" ;
         Frost:Fill_Value = -9999.f ;
     float Tmin_1(Z, Y, X) ;
         Tmin_1:long_name = "First monthly minimum of Tmin" ;
         Tmin_1:units = "C" ;
         Tmin_1:Fill_Value = -9999.f ;
     float Tmin_2(Z, Y, X) ;
         Tmin_2:long_name = "Second montlhy minimum of Tmin" ;
         Tmin_2:units = "C" ;
         Tmin_2:Fill_Value = -9999.f ;
     float Tmin_3(Z, Y, X) ;
         Tmin_3:long_name = "Third monthly minimum of Tmin" ;
         Tmin_3:units = "C" ;
         Tmin_3:Fill_Value = -9999.f ;
     float Tmax(Z, Y, X) ;
         Tmax:long_name = "Monthly average of daily Tmax" ;
         Tmax:units = "C" ;
         Tmax:Fill_Value = -9999.f ;
     float TT30(Z, Y, X) ;
         TT30:long_name = "Number of days with Tmax > 30 degC" ;
         TT30:units = "days" ;
         TT30:Fill_Value = -9999.f ;
     float TT35(Z, Y, X) ;
         TT35:long_name = "Number of days with Tmax > 35 degC" ;
         TT35:units = "days" ;
         TT35:Fill_Value = -9999.f ;
     float TT40(Z, Y, X) ;
         TT40:long_name = "Number of days with Tmax > 40 degC" ;
         TT40:units = "days" ;
         TT40:Fill_Value = -9999.f ;
     float Tmax_1(Z, Y, X) ;
         Tmax_1:long_name = "First monthly maximum of Tmax" ;
         Tmax_1:units = "C" ;
         Tmax_1:Fill_Value = -9999.f ;
     float Tmax_2(Z, Y, X) ;
         Tmax_2:long_name = "Second monthly maximum of Tmax" ;
         Tmax_2:units = "C" ;
         Tmax_2:Fill_Value = -9999.f ;
     float Tmax_3(Z, Y, X) ;
         Tmax_3:long_name = "Third monthly maximum of Tmax" ;
         Tmax_3:units = "C" ;
         Tmax_3:Fill_Value = -9999.f ;
     float U-std(Z, Y, X) ;
         U-std:long_name = "Monthly standard deviation of U surface wind" ;
         U-std:units = "m s-1" ;
         U-std:Fill_Value = -9999.f ;
     float V-std(Z, Y, X) ;
         V-std:long_name = "Monthly standard deviation of V surface wind" ;
         V-std:units = "m s-1" ;
         V-std:Fill_Value = -9999.f ;
     float Max_wind_speed(Z, Y, X) ;
         Max_wind_speed:long_name = "Monthly maximum of daily wind speed" ;
         Max_wind_speed:units = "m s-1" ;
         Max_wind_speed:Fill_Value = -9999.f ;
     float Dry_max_1(Z, Y, X) ;
         Dry_max_1:long_name = "Max number of consecutive dry days at 
start of the month" ;
         Dry_max_1:units = "days" ;
         Dry_max_1:Fill_Value = -9999.f ;
     float Dry_max_2(Z, Y, X) ;
         Dry_max_2:long_name = "Max number of consecutive dry days at 
end of the month" ;
         Dry_max_2:units = "days" ;
         Dry_max_2:Fill_Value = -9999.f ;
     float Dry_max_3(Z, Y, X) ;
         Dry_max_3:long_name = "Max number of consecutive dry days 
during the month" ;
         Dry_max_3:units = "days" ;
         Dry_max_3:Fill_Value = -9999.f ;
     float Wet_days(Z, Y, X) ;
         Wet_days:long_name = "Number of wet days during the month" ;
         Wet_days:units = "days" ;
         Wet_days:Fill_Value = -9999.f ;
     float Precip_1(Z, Y, X) ;
         Precip_1:long_name = "Total precipitation within the 1st day of 
the month" ;
         Precip_1:units = "mm" ;
         Precip_1:Fill_Value = -9999.f ;
     float Precip_2(Z, Y, X) ;
         Precip_2:long_name = "Total precipitation on the 21st days of 
the month" ;
         Precip_2:units = "mm" ;
         Precip_2:Fill_Value = -9999.f ;
     float Precip_3(Z, Y, X) ;
         Precip_3:long_name = "Total precipitation within 3 days" ;
         Precip_3:units = "mm" ;
         Precip_3:Fill_Value = -9999.f ;
     float Precip_4(Z, Y, X) ;
         Precip_4:long_name = "Total precipitation within the 2 last 
days of the month" ;
         Precip_4:units = "mm" ;
         Precip_4:Fill_Value = -9999.f ;
     float Precip_5(Z, Y, X) ;
         Precip_5:long_name = "Total precipitation within the last day 
of the month" ;
         Precip_5:units = "mm" ;
         Precip_5:Fill_Value = -9999.f ;
     float Pge40(Z, Y, X) ;
         Pge40:long_name = "Number of days with daily P > 40 mm" ;
         Pge40:units = "days" ;
         Pge40:Fill_Value = -9999.f ;
     float Pge80(Z, Y, X) ;
         Pge80:long_name = "Number of days with daily P > 80 mm" ;
         Pge80:units = "days" ;
         Pge80:Fill_Value = -9999.f ;
     float Pge120(Z, Y, X) ;
         Pge120:long_name = "Number of days with daily P > 120 mm" ;
         Pge120:units = "days" ;
         Pge120:Fill_Value = -9999.f ;
     float Pmax_1(Z, Y, X) ;
         Pmax_1:long_name = "First max daily total precipitation rate 
over the month" ;
         Pmax_1:units = "mm" ;
         Pmax_1:Fill_Value = -9999.f ;
     float Pmax_2(Z, Y, X) ;
         Pmax_2:long_name = "Second max daily total precipitation rate 
over the month" ;
         Pmax_2:units = "mm" ;
         Pmax_2:Fill_Value = -9999.f ;
     float Pmax_3(Z, Y, X) ;
         Pmax_3:long_name = "Third max daily total precipitation rate 
over the month" ;
         Pmax_3:units = "mm" ;
         Pmax_3:Fill_Value = -9999.f ;
     float Z(Z, Y, X) ;
         Z:long_name = "Monthly mean geopotential height at 500mb" ;
         Z:units = "m" ;
         Z:Fill_Value = -9999.f ;
     float Z-std(Z, Y, X) ;
         Z-std:long_name = "Monthly standard deviation of geopotential 
height at 500mb" ;
         Z-std:units = "m" ;
         Z-std:Fill_Value = -9999.f ;
     float Snow_mass(Z, Y, X) ;
         Snow_mass:long_name = "Snow mass on the first day of the month" ;
         Snow_mass:units = "kg m-2" ;
         Snow_mass:Fill_Value = -9999.f ;
     float DlastSM(Z, Y, X) ;
         DlastSM:long_name = "Date of last snow mass" ;
         DlastSM:units = "day" ;
         DlastSM:Fill_Value = -9999.f ;
     float Max_runoff(Z, Y, X) ;
         Max_runoff:long_name = "Maximum daily runoff" ;
         Max_runoff:units = "kg m-2 s-1" ;
         Max_runoff:Fill_Value = -9999.f ;
     float Date_max_runoff(Z, Y, X) ;
         Date_max_runoff:long_name = "Date of maximum daily runoff" ;
         Date_max_runoff:units = "day" ;
         Date_max_runoff:Fill_Value = -9999.f ;
data:
}

--------------------
in grads:

GX Package Initialization: Size = 11 8.5
Command line history in /home/rjmera/.grads.log
ga-> sdfopen Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
Scanning self-describing file:  Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
gadsdf: SDF file has no discernable X coordinate.

ga-> open Desktop/CPDN/pdg6dec.ctl
Scanning description file:  Desktop/CPDN/pdg6dec.ctl
Data file Desktop/CPDN/nhqhga.pdg6dec.pnw.nc is open as file 1
LON set to 347.79 371.77
LAT set to -11.99 13.31
LEV set to 1 1
Time values set: 1967:12:6:0 1967:12:6:0
E set to 1 1
ga-> q file
File 1 : Pacific northwest
   Descriptor: Desktop/CPDN/pdg6dec.ctl
   Binary: Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
   Type = Gridded
   Xsize = 110  Ysize = 117  Zsize = 1  Tsize = 1  Esize = 1
   Number of Variables = 36
      tmin  1  t,z,y,x  Tmin (C)
      frost  1  t,z,y,x  Frost (days)
      tmin1  1  t,z,y,x  Min1 (C)
      tmin2  1  t,z,y,x  Min2 (C)
      tmin3  1  t,z,y,x  Min3 (C)
      tmax  1  t,z,y,x  Tmax (C)
      tt30  1  t,z,y,x  30Ctmax (days)
      tt35  1  t,z,y,x  35Ctmax (days)
      tt40  1  t,z,y,x  40Ctmax (days)
      tmax1  1  t,z,y,x  Max1 (C)
      tmax2  1  t,z,y,x  Max2 (C)
      tmax3  1  t,z,y,x  Max3 (C)
      ustd  1  t,z,y,x  monstdU (m s-1)
      vstd  1  t,z,y,x  monstdU (m s-1)
      mxwspeed    1  1  t,z,y,x  maxwindsped (m s-1)
      drymax1    1  1  t,z,y,x  daysmaxdrystart (days)
      drymax2    1  1  t,z,y,x  daysmaxdrymid (days)
      drymax3    1  1  t,z,y,x  daysmaxdryend (days)
      wetdays    1  1  t,z,y,x  numwetdays (days)
      p1    1  1  t,z,y,x  1stdayprec (mm)
      p2    1  1  t,z,y,x  21stdayprec (mm)
      p3    1  1  t,z,y,x  3daysprec (mm)
      p4    1  1  t,z,y,x  preclast2days (mm)
      p5    1  1  t,z,y,x  preclastday (mm)
      pge40    1  1  t,z,y,x  40mmprec (days)
      pge80    1  1  t,z,y,x  80mmprec (days)
      pge120    1  1  t,z,y,x  120mmprec (days)
      pmax1    1  1  t,z,y,x  1stmaxprec (C)
      pmax2    1  1  t,z,y,x  2ndmaxprec (C)
      pmax3    1  1  t,z,y,x  3rdmaxprec (C)
      z    1  1  t,z,y,x  hgtmean500 (m)
      zstd    1  1  t,z,y,x  hgtstdv (m)
      snowmass    1  1  t,z,y,x  snowmass (kg m-2)
      lastsnow    1  1  t,z,y,x  daylastsnw (day)
      mxrunoff    1  1  t,z,y,x  maxrunoff (kg m-2 s-1)
      datemxrunoff    1  1  t,z,y,x  maxrunoffdate (day)
ga-> d p1
gancgrid error: nc_get_vara_double failed; NetCDF: Start+count exceeds 
dimension bound
Data Request Error:  Error for variable 'p1'
   Error ocurred at column 1
DISPLAY error:  Invalid expression
   Expression = p1
ga->


On 10/6/10 5:17 AM, Jennifer Adams wrote:
> Roberto,
> You need to give the full ncdump -c output (just the time and latitude 
> coordinate metadata is not adequate), as well as the xdfopen-style 
> descriptor, and the full text of your GrADS session that delivers the 
> error message.
> --Jennifer
>
> On Oct 5, 2010, at 4:59 PM, Roberto Mera wrote:
>
>>
>>
>> Regarding the questions below: I'm using opengrads 2.0.a7.oga.3. and 
>> grads 1.9. on linux.
>> First, using sdfopen I get this message:
>>
>> sdfopen Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
>> Scanning self-describing file:  Desktop/CPDN/nhqhga.pdg6dec.pnw.nc
>> gadsdf: SDF file has no discernable X coordinate.
>>
>> Subsequently, I tried xdfopen with a control file. That did not work
>> either. Then I created another ctl file that would read its variables
>> and coordinate system. It opened fine but when I try to just draw the
>> map I get this message no matter which variable I'm working with:
>>
>> gancgrid error: nc_get_vara_double failed; NetCDF: Start+count exceeds
>> dimension bound
>> Data Request Error:  Error for variable 'p1'
>>    Error ocurred at column 1
>>
>> ncdump -h :
>>
>> netcdf nhqhga.pdg6dec.pnw {
>> dimensions:
>> X = 110 ;
>> Y = 117 ;
>> Z = 1 ;
>> variables:
>> int time(Z) ;
>> time:long_name = "time" ;
>> time:units = "month" ;
>> float latitude(Y) ;
>> latitude:long_name = "latitude" ;
>> latitude:units = "degrees_north" ;
>>
>>
>> .....plus other variables (37)
>>
>> Any thoughts?
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> Jennifer M. Adams
> IGES/COLA
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705
> jma at cola.iges.org <mailto:jma at cola.iges.org>
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101006/e6cbb278/attachment-0003.html 


More information about the gradsusr mailing list