<div dir="ltr"><div style="font-family:arial,sans-serif;font-size:13px"><div><div><div><br></div><div>Hi NCLers<br></div>I am trying to plot of map of GOES IR Brightness temperature, but am getting a flat field. When I print the data they&#39;re fine, the values are correct.<br><br></div>The file I am using I created with <br>ncrcat<br></div>to put all of the 15 minutes GOES netcdf images into a daily file.  Maybe I erred here, maybe something wrong with the time dimension.<br><br></div><span style="font-family:arial,sans-serif;font-size:13px">Here is ncdump -h  of the file I want to plot</span><br style="font-family:arial,sans-serif;font-size:13px"><div style="font-family:arial,sans-serif;font-size:13px">-------------------------------------------------------------------------------<br>ncdump -h <a href="http://test.nc/" target="_blank">test.nc</a><br>netcdf test {<br>dimensions:<br>    yc = 221 ;<br>    xc = 141 ;<br>    time = UNLIMITED ; // (130 currently)<br>variables:<br>    float new_temp(time, yc, xc) ;<br>        new_temp:time = 1372896900 ;<br>        new_temp:units = &quot;K&quot; ;<br>        new_temp:coordinates = &quot;lon lat&quot; ;<br>        new_temp:type = &quot;GVAR&quot; ;<br>        new_temp:long_name = &quot;Temperature&quot; ;<br>    float new_lat(time, yc, xc) ;<br>        new_lat:units = &quot;degrees_north&quot; ;<br>        new_lat:long_name = &quot;lat&quot; ;<br>    float new_lon(time, yc, xc) ;<br>        new_lon:units = &quot;degrees_east&quot; ;<br>        new_lon:long_name = &quot;lon&quot; ;<br>    int time(time) ;<br>        time:units = &quot;seconds since 2013-185 00:15&quot; ;<br><br>// global attributes:<br>        :history = &quot;Mon Oct  6 18:26:36 2014: ncrcat <a href="http://goes13_4_2013_185_0015_subset.nc/" target="_blank">goes13_4_2013_185_0015_subset.nc</a> ............<br><a href="http://goes13_4_2013_185_2345_subset.nc/" target="_blank">goes13_4_2013_185_2345_subset.nc</a> <a href="http://test.nc/" target="_blank">test.nc</a>\n&quot;,<br>            &quot;Mon Oct  6 18:24:58 2014: ncatted -O -a units,time,c,c,seconds since 2013-185 00:15 <a href="http://out.nc/" target="_blank">out.nc</a>\n&quot;,<br>            &quot;Mon Oct  6 18:24:58 2014: ncap2 -O -v -s time=array(0,0,$time) <a href="http://out.nc/" target="_blank">out.nc</a>\n&quot;,<br>            &quot;Mon Oct  6 18:24:58 2014: ncrename -O -d record,time <a href="http://out.nc/" target="_blank">out.nc</a>&quot; ;<br>        :nco_openmp_thread_number = 1 ;<br>}<br>----------------------------------------------------------------------------------------------------------------------<br></div><div style="font-family:arial,sans-serif;font-size:13px"><br>Here is the map plotting code<br></div><div style="font-family:arial,sans-serif;font-size:13px">-----------------------------------------------------------------------------------------------------------------------<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br><br>;*********************************<br>; read variable<br>;*********************************<br><br>f    = addfile(&quot;<a href="http://test.nc/" target="_blank">test.nc</a>&quot;,&quot;r&quot;)<br>d    = f-&gt;new_temp          ; float data(time, yc, xc) ;<br>                           ; data:type = &quot;IR&quot; ;<br>   printVarSummary(d)<br>;*********************************<br>; Fix the variable so it has recognizable missing data<br>;*********************************<br>   new_lat  = f-&gt;new_lat<br>   new_lon  = f-&gt;new_lon<br><br>;*********************************<br>; create plot<br>;*********************************<br>   pltType = &quot;pdf&quot;                             ; &quot;ps&quot;, &quot;eps&quot;, &quot;pdf&quot;, &quot;png&quot;<br>   pltDir  = &quot;./&quot;<br><br>   wks = gsn_open_wks(pltType, &quot;test&quot;)<br>   gsn_define_colormap(wks,&quot;BlAqGrYeOrReVi200&quot;); choose 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 lines<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@mpMinLatF           = 23.0  ;    ; region to zoom in on<br>   res@mpMaxLatF           = 33.00  ; <br>   res@mpMinLonF           = -115.200   ; <br>   res@mpMaxLonF           = -102.20 ; <br>   res@mpFillOn            = False<br>   res@mpOutlineBoundarySets = &quot;National&quot;    ; turn on country boundaries<br>  res@trGridType          = &quot;TriangularMesh&quot;   ; Necessary b/c lat, lon<br>  res@lbLabelStride        = 5.0             ; every other label bar label<br><br>   plot = gsn_csm_contour_map_ce(wks,d(111,:,:), res)<br>----------------------------------------------------------------------------------------------------<br><br></div><div style="font-family:arial,sans-serif;font-size:13px">Also, I am getting these warnings,<br></div><div style="font-family:arial,sans-serif;font-size:13px"><br>0)    check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn&#39;t contain one at all.<br>(0)    A valid latitude coordinate array should have a &#39;units&#39; attribute equal to one of the following values: <br>(0)        &#39;degrees_north&#39; &#39;degrees-north&#39; &#39;degree_north&#39; &#39;degrees north&#39; &#39;degrees_N&#39; &#39;Degrees_north&#39; &#39;degree_N&#39; &#39;degreeN&#39; &#39;degreesN&#39; &#39;deg north&#39;<br>(0)    check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn&#39;t contain one at all.<br>(0)    A valid longitude coordinate array should have a &#39;units&#39; attribute equal to one of the following values: <br>(0)        &#39;degrees_east&#39; &#39;degrees-east&#39; &#39;degree_east&#39; &#39;degrees east&#39; &#39;degrees_E&#39; &#39;Degrees_east&#39; &#39;degree_E&#39; &#39;degreeE&#39; &#39;degreesE&#39; &#39;deg east&#39;<br><br></div><div style="font-family:arial,sans-serif;font-size:13px">Any ideas?<br><br>thanks in advance,<br></div><div style="font-family:arial,sans-serif;font-size:13px">Dave</div></div>