[gradsusr] netcdf control file error
Jennifer Adams
jma at cola.iges.org
Wed Oct 6 12:52:31 EDT 2010
You forgot to include your descriptor pgd6dec.ctl, and you omitted the
critical part from the beginning of the GrADS session that includes
the version you are running. I don't ask for these things for trivial
reasons!!! Based on what you did send, it looks like you might have
something like this in your descriptor:
Precip_1=>p1 1 t,z,y,x precip in 1st day of month
But the variable Precip_1 has only Z,Y,X dimensions. Replace
"t,z,y,x" with "z,y,x" and try again.
--Jennifer
On Oct 6, 2010, at 12:28 PM, Roberto Mera wrote:
> 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
>>> 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
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101006/f044f1f2/attachment-0003.html
More information about the gradsusr
mailing list