<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 = &quot;./&quot;<br>   fili = systemfunc(&quot;ls goes1*.nc&quot;)<br><br>   print(fili)<br>   nfili= dimsizes(fili)<br>   print(&quot;nfili=&quot;+nfili)<br><br>   f    = addfiles(diri+fili,&quot;r&quot;)<br>   d    = f[:]-&gt;data              ; float data(time, yc, xc)  ; data:type = &quot;IR&quot; ;<br>   d@_FillValue = 330<br>   lat = f[:]-&gt;lat<br>   lon = f[:]-&gt;lon<br>   time= f[:]-&gt;time<br><br>   nl     = 600<br>   dhov = d(:,nl,:)    ; (time,longitude)  ; &lt;== plot<br>; Fix the variable so it has recognizable missing data<br>;*********************************<br><br>   lat2d  = f[:]-&gt;lat<br>   lon2d  =  f[:]-&gt;lon<br>   d@lat2d = lat2d        ; (yc, xc)<br>   d@lon2d = lon2d<br><br> crTime  = f[:]-&gt;crTime     ; assorted info<br>   crDate  = f[:]-&gt;crDate<br>   lineRes = f[:]-&gt;lineRes<br>   elemRes = f[:]-&gt;elemRes<br><br>;*********************************<br>; create plot<br>;*********************************<br>   fNameBase = str_get_field(fili,1,&quot;.&quot;) + str_get_field(fili,2,&quot;.&quot;)<br><br>  ;pltName = fNameBase<br>   pltName = &quot;goes_4&quot;<br>   pltType = &quot;pdf&quot;                              ; &quot;ps&quot;, &quot;eps&quot;, &quot;pdf&quot;, &quot;png&quot;,&quot;x11&quot;<br>   pltDir  = &quot;./&quot;<br><br>   wks = gsn_open_wks(pltType, pltDir+pltName)<br>   gsn_define_colormap(wks,&quot;BlAqGrYeOrReVi200&quot;); specify a color map<br><br>   res                     = True<br>   res@cnFillOn            = True            ; turn on color<br>  res@cnFillMode          = &quot;RasterFill&quot;    ; 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 = &quot;Always&quot;      ; use NCL default<br>   res@lbLabelAutoStride   = True            ; let NCL decide spacing<br>   res@cnLevelSelectionMode = &quot;ManualLevels&quot;  ;&quot;ExplicitLevels&quot;<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 = &quot;McIDAS area file version&quot; ;<br>        int sensorID ;<br>        sensorID:long_name = &quot;McIDAS sensor number&quot; ;<br>    int imageDate ;<br>        imageDate:long_name = &quot;image year and day of year (in ccyyddd format)&quot; ;<br>    int imageTime ;<br>        imageTime:long_name = &quot;image time in UTC (hours/minutes/seconds, in HHMMSS format)&quot; ;<br>    int startLine ;<br>        startLine:long_name = &quot;starting image line (in satellite coordinates)&quot; ;<br>    int startElem ;<br>        startElem:long_name = &quot;starting image element (in satellite coordinates)&quot; ;<br>    int time(time) ;<br>        time:long_name = &quot;seconds since 1970-1-1 0:0:0&quot; ;<br>        time:units = &quot;seconds since 1970-1-1 0:0:0&quot; ;<br>    int dataWidth ;<br>        dataWidth:long_name = &quot;number of 8-bit bytes per source data point&quot; ;<br>    int lineRes ;<br>        lineRes:long_name = &quot;resolution of each pixel in line direction&quot; ;<br>        lineRes:units = &quot;km&quot; ;<br>    int elemRes ;<br>        elemRes:long_name = &quot;resolution of each pixel in element direction&quot; ;<br>        elemRes:units = &quot;km&quot; ;<br>    int prefixSize ;<br>        prefixSize:long_name = &quot;line prefix size in 8-bit bytes&quot; ;<br>    int crDate ;<br>        crDate:long_name = &quot;image creation year and day of year in ccyyddd format&quot; ;<br>    int crTime ;<br>        crTime:long_name = &quot;image creation time in UTC in hhmmss format&quot; ;<br>    int bands ;<br>        bands:long_name = &quot;satellite channel number&quot; ;<br>    char auditTrail(auditCount, auditSize) ;<br>        auditTrail:long_name = &quot;audit trail&quot; ;<br>    float data(time, yc, xc) ;<br>        data:long_name = &quot;Temperature&quot; ;<br>        data:type = &quot;VISR&quot; ;<br>        data:coordinates = &quot;lon lat&quot; ;<br>        data:units = &quot;K&quot; ;<br>    float lat(yc, xc) ;<br>        lat:long_name = &quot;lat&quot; ;<br>        lat:units = &quot;degrees_north&quot; ;<br>    float lon(yc, xc) ;<br>        lon:long_name = &quot;lon&quot; ;<br>        lon:units = &quot;degrees_east&quot; ;<br><br>// global attributes:<br>        :Conventions = &quot;CF-1.4&quot; ;<br>        :Source = &quot;McIDAS Area File&quot; ;<br>        :Satellite\ Sensor = &quot;G-13 IMG    &quot; ;<br><br><br><br><br><br><br><br><br><div><div><div><br><br><br></div></div></div></div></div>