coordinates units

Dear NCL users,

With the following script

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
   setfileoption("bin","ReadByteOrder","BigEndian")
   prec=fbindirread("gpcp_1dd_p1d.200001",1,(/360,180/),"float")

   lat=latGlobeFo(180,"lat","latitude","degrees_north")
   lat=lat(::-1)
   lon=lonGlobeFo(360,"lon","longitude","degrees_east")
  
   prec!0="lon"
   prec!1="lat"
   prec&lon=lon
   prec&lat=lat

   prec@_FillValue=-99999.
   prec_at_long_name="GPCP 1DD Data"
   prec_at_units="mm/day"

   xwks = gsn_open_wks("pdf","mapa")
   plot = gsn_csm_contour_map_ce(xwks,prec,False)
end

NCL gives me the following error messages:

(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'

I tried to put the lat and lon units by hand and with the latGlobeFo()
and lonGlobeFo(), but the messages were the same.
Can anyone help me?

Thanks
Mateus

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 10 2007 - 09:38:56 MST