[ncl-talk] fatal:Dimension (lon) of (tprec) does not have an associated coordinate variable

[Jayant]

my bad! Thanks Dennis

On Sun, Apr 11, 2021 at 10:22 PM Dennis Shea <shea at ucar.edu> wrote:

> @ is for attributes
> & is for coordinates
>
> tprec at lat=lat
> tprec at lon=lon
>
> should be
>
> tprec&lat=lat
> tprec&lon=lon
>
>
> On Sun, Apr 11, 2021 at 2:28 PM Jayant via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Hi Adam,
>> I am reading time-lat-lon data from a binary file. Next, I create
>> coordinate variables based on the accompanying CTL file. The CTL file
>> indicates the binary to be sequential. I assign the coordinates to a
>> variable. But, I am getting an error. Below is a working example:
>>
>>  nlat=192
>>  nlon=384
>>  ntim=24
>>  ;
>>  lon = lonGlobeF(nlon, "lon", "longitude", "degrees_east")
>>  lat = latGau(nlat, "lat", "latitude", "degrees_north")
>>  ;
>>  tim=new(ntim,integer)
>>  tim!0         = "time"
>>  tim at long_name = "time"
>>  tim at units     = "hours"
>>  tim&time      = tim
>>  tim=0
>>  ;
>>  tprec=new((/ntim,nlat,nlon/),float)    ; surface variable
>>  tprec!0="time"
>>  tprec!1="lat"
>>  tprec!2="lon"
>>  tprec at time=tim
>>  tprec at lat=lat
>>  tprec at lon=lon
>>  tprec at long_name="total precipitation"
>>  tprec at units="kg m**-2 day**-1"
>>  tprec=0.0
>>  tprec=lonFlip(tprec)
>> ;--- executing the above
>> ncl 29>  print(tprec(16,{-5.0:-3.0},{-90:-85}))
>> fatal:Dimension (lon) of (tprec) does not have an associated coordinate
>> variable
>> fatal:["Execute.c":7741]:Execute: Error occurred at or near line 29
>> ;
>> ncl 30>  print(tprec(16,{-5.0:-3.0},:))
>> fatal:Dimension (lat) of (tprec) does not have an associated coordinate
>> variable
>> fatal:["Execute.c":7741]:Execute: Error occurred at or near line 30
>> ;
>> ncl 31> printVarSummary(tprec)
>> Variable: tprec
>> Type: float
>> Total Size: 7077888 bytes
>>             1769472 values
>> Number of Dimensions: 3
>> Dimensions and sizes: [time | 24] x [lat | 192] x [lon | 384]
>> Coordinates:
>> Number Of Attributes: 7
>>   lonFlip : longitude coordinate variable has been reordered via lonFlip
>>   time : <ARRAY of 24 elements>
>>   lat : <ARRAY of 192 elements>
>>   lon : <ARRAY of 384 elements>
>>   long_name : total precipitation
>>   units : kg m**-2 day**-1
>>   _FillValue : 9.96921e+36
>>
>> ncl 32> print(isdimnamed(prec,ispan(0,2,1)))
>>     True
>>     True
>>     True
>> ;
>> ncl 33>  if(all(iscoord(tprec,(/"time","lat","lon"/)))) then
>>     print("min/max tprec&lat = " + min(tprec&lat) + "/" + max(tprec&lat))
>>     print("min/max tprec&lon = " + min(tprec&lon) + "/" + max(tprec&lon))
>>          end if
>> ncl 34>
>>
>> What am I missing?
>>
>> Also, the warnings (below) that appears when plotting (ignoring the print
>> statements) related to the above issue with coordinate variables?
>> (0) check_for_y_lat_coord: Warning: Data either does not contain a valid
>> latitude coordinate array or doesn't contain one at all.
>> (0) A valid latitude coordinate array should have a 'units' attribute
>> equal to one of the following values:
>> (0)    'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
>> 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
>> (0) check_for_lon_coord: Warning: Data either does not contain a valid
>> longitude coordinate array or doesn't contain one at all.
>> (0) A valid longitude coordinate array should have a 'units' attribute
>> equal to one of the following values:
>> (0)    'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>
>> Best,
>> Jayant
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210411/aacc0202/attachment-0001.html>