<div dir="ltr"><div><div>Hi NCLers, <br><br></div>I have managed to get GOES IR data from curvilinear to rectangular grids. Lat and Lon are still 2d, however.  But the code I have for the Hovmoller diagrams gives this error. I am sure I am just not reading the data correctly.<br><br></div>The error, ncl code and ncdump of the file follow:<br><div><br></div><div>thanks,<br></div><div>Dave<br></div><div><br></div><div>---------------------- error ------------------------------------------------------<br></div><div>warning:IrTransInitialize: error creating spline approximation for trYCoordPoints; defaulting to linear<br>fatal:ContourPlotDraw: CPCICA - ONE OF THE CORNER POINTS OF THE CELL ARRAY IS INCORRECT<br>fatal:ContourPlotDraw: draw error<br>fatal:ContourPlotDraw: draw error<br>fatal:PlotManagerDraw: error in plot draw<br>fatal:_NhlPlotManagerDraw: Draw error<br><br></div><div>------------------------ NCL Code Snippet -----------------------------------------------<br>diri = "./"<br>  fili = systemfunc("ls goes1*.nc")<br><br>  print(fili)<br>  nfili= dimsizes(fili)<br>  print("nfili="+nfili)<br><br>  f   = addfiles(diri+fili,"r")<br>  d   = f[:]->data             ; float data(time, yc, xc) ; data:type = "IR" ;<br>  d@_FillValue = 330<br>  lat = f[:]->lat<br>  lon = f[:]->lon<br>  time= f[:]->time<br><br>  nl    = 600<br>  dhov = d(:,nl,:)   ; (time,longitude) ; <== plot<br>; Fix the variable so it has recognizable missing data<br>;*********************************<br><br>  lat2d = f[:]->lat<br>  lon2d = f[:]->lon<br>  d@lat2d = lat2d       ; (yc, xc)<br>  d@lon2d = lon2d<br><br> crTime = f[:]->crTime    ; assorted info<br>  crDate = f[:]->crDate<br>  lineRes = f[:]->lineRes<br>  elemRes = f[:]->elemRes<br><br>;*********************************<br>; create plot<br>;*********************************<br>  fNameBase = str_get_field(fili,1,".") + str_get_field(fili,2,".")<br><br> ;pltName = fNameBase<br>  pltName = "goes_4"<br>  pltType = "pdf"                             ; "ps", "eps", "pdf", "png","x11"<br>  pltDir = "./"<br><br>  wks = gsn_open_wks(pltType, pltDir+pltName)<br>  gsn_define_colormap(wks,"BlAqGrYeOrReVi200"); specify a color map<br><br>  res                    = True<br>  res@cnFillOn           = True           ; turn on color<br> res@cnFillMode         = "RasterFill"   ; cell mode<br>  res@cnLinesOn          = False          ; Turn off contour line<br>;  Create longitude labels on the x-axis<br>  res@gsnSpreadColors    = True           ; use full colormap<br>  res@gsnAddCyclic       = False          ; data not cyclic<br>  res@gsnMaximize        = True           ; ps, pdf, pdf<br>  res@pmTickMarkDisplayMode = "Always"     ; use NCL default<br>  res@lbLabelAutoStride  = True           ; let NCL decide spacing<br>  res@cnLevelSelectionMode = "ManualLevels" ;"ExplicitLevels"<br>  res@cnMinLevelValF      = 180.0                           ; set the minimum contour level<br>  res@cnMaxLevelValF      = 310.                         ; set the maximum contour level<br>  res@cnLevelSpacingF     = 5.0                           ; set the contour interval<br>  res@cnRasterSmoothingOn = True<br>  res@lbLabelStride       = 5.0            ; every other label bar label<br><br> plot = gsn_csm_contour(wks,dhov(:,:), res)<br>---------------------------------------------------------------------------------------------------------<br><br>netcdf goes13_4_2013_210_2002_1byte_rect {<br>dimensions:<br>   xc = 2720 ;<br>   yc = 1301 ;<br>   time = 1 ;<br>   auditCount = 2 ;<br>   auditSize = 80 ;<br>variables:<br>   int version ;<br>      version:long_name = "McIDAS area file version" ;<br>      int sensorID ;<br>      sensorID:long_name = "McIDAS sensor number" ;<br>   int imageDate ;<br>      imageDate:long_name = "image year and day of year (in ccyyddd format)" ;<br>   int imageTime ;<br>      imageTime:long_name = "image time in UTC (hours/minutes/seconds, in HHMMSS format)" ;<br>   int startLine ;<br>      startLine:long_name = "starting image line (in satellite coordinates)" ;<br>   int startElem ;<br>      startElem:long_name = "starting image element (in satellite coordinates)" ;<br>   int time(time) ;<br>      time:long_name = "seconds since 1970-1-1 0:0:0" ;<br>      time:units = "seconds since 1970-1-1 0:0:0" ;<br>   int dataWidth ;<br>      dataWidth:long_name = "number of 8-bit bytes per source data point" ;<br>   int lineRes ;<br>      lineRes:long_name = "resolution of each pixel in line direction" ;<br>      lineRes:units = "km" ;<br>   int elemRes ;<br>      elemRes:long_name = "resolution of each pixel in element direction" ;<br>      elemRes:units = "km" ;<br>   int prefixSize ;<br>      prefixSize:long_name = "line prefix size in 8-bit bytes" ;<br>   int crDate ;<br>      crDate:long_name = "image creation year and day of year in ccyyddd format" ;<br>   int crTime ;<br>      crTime:long_name = "image creation time in UTC in hhmmss format" ;<br>   int bands ;<br>      bands:long_name = "satellite channel number" ;<br>   char auditTrail(auditCount, auditSize) ;<br>      auditTrail:long_name = "audit trail" ;<br>   float data(time, yc, xc) ;<br>      data:long_name = "Temperature" ;<br>      data:type = "VISR" ;<br>      data:coordinates = "lon lat" ;<br>      data:units = "K" ;<br>   float lat(yc, xc) ;<br>      lat:long_name = "lat" ;<br>      lat:units = "degrees_north" ;<br>   float lon(yc, xc) ;<br>      lon:long_name = "lon" ;<br>      lon:units = "degrees_east" ;<br><br>// global attributes:<br>      :Conventions = "CF-1.4" ;<br>      :Source = "McIDAS Area File" ;<br>      :Satellite\ Sensor = "G-13 IMG   " ;<br><br><br><br><br><br><br><br><br><div><div><div><br><br><br></div></div></div></div></div>