[ncl-talk] check_for_y_lat_coord

[Ehsan Taghizadeh]

 I've just want to inform about the solution. The problem was that I've written:plot = gsn_csm_contour_map(wks,tp(nt,:,:)*1000, res)
and the problem is "*1000". So I added following line:tp = (/ tp*1000 /)
plot = gsn_csm_contour_map(wks,tp(nt,:,:), res)
and it's just fixed.
SincerelyEhsan
    On Tuesday, February 26, 2019, 10:45:23 AM GMT+3:30, Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com> wrote:  
 
  Dear Rick,Thank you so much for your helpful reply. I've uncommented the following lines, however it didn't change the result:
   tp&latitude at units="degrees_north"   tp&longitude at units="degrees_east"
I've also attached .tar.gz file of data, if it could help.SincerelyEhsan    On Tuesday, February 26, 2019, 1:40:51 AM GMT+3:30, Rick Brownrigg <brownrig at ucar.edu> wrote:  
 
 Hi,
I don't have a means to open the .rar file you sent, but it sounds like its the latitude and longitude variables in the file are what do not have the units attribute.  From your printFileSummary(tp):
Variable: tpType: floatTotal Size: 3755712 bytes            938928 valuesNumber of Dimensions: 3Dimensions and sizes:   [time | 1262] x [latitude | 24] x [longitude | 31]Coordinates:            time: [1008072..1038336]            latitude: [41.25..24]            longitude: [42.75..65.25]
we can see there's already coordinates attached to tp. Perhaps try verifying that longitude and latitude in your file are lacking units, and if so, add them to those variables -- I don't think you want to do what you found on the link you cited.
Hope that helps...Rick

On Mon, Feb 25, 2019 at 10:47 AM Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com> wrote:

Dears,I've tried to plot precipitation data from ERA-Interim with the attached script.However after trying this I encountered the following errors:
Variable: tpType: floatTotal Size: 3755712 bytes            938928 valuesNumber of Dimensions: 3Dimensions and sizes:   [time | 1262] x [latitude | 24] x [longitude | 31]Coordinates:            time: [1008072..1038336]            latitude: [41.25..24]            longitude: [42.75..65.25]Number Of Attributes: 8  _FillValue :  -32767  missing_value :       -32767  units :       m  long_name :   Total precipitation  _FillValue_original : -32767  missing_value_original :      -32767  time :        1008060  NCL_tag :     calculate_daily_values: arith=sum(0)     ==================================(0)     check_for_y_lat_coord: Warning: Data either does not contain(0)     a valid latitude coordinate array or doesn't contain one at all.(0)     A valid latitude coordinate array should have a 'units'(0)     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(0)     a valid longitude coordinate array or doesn't contain one at all.(0)     A valid longitude coordinate array should have a 'units'(0)     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'
I've found this link (https://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml), and based on it I tried something like these:   tp!0="time"   tp!1="lat"   tp!2="lon"   tp&lat at units="degrees_north"   tp&lon at units="degrees_east"
But unfortunately they doesn't help!May I ask your help to plot precipitation data from ERA-Interim. I've just to say I did a processing to original ERA-Interim data to obtain daily precipitation (It has been done before with ncl-talk help). Attached are the script and data.Any help would be appreciated.

SincerelyEhsan_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

    
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190226/900e1263/attachment.html>