[ncl-talk] no valid values in scalar field; ContourPlot not possible:[errno=1101]

[Sebastian Otarola-Bustos]

Hi everybody,

I'm trying to plot a contour with week average of sst anomalies, but I'm
having that error. I read that could be because of the coordinates, but I'd
copied all the cordinates from sst(daily values (time, lat,lon)) to my
auxiliary sst_avg(week average(time,lat,lon)), despite I don'y need all of
them( I don't need time), just to be sure. I also tried setting to 0 all
fill_Values and missing values. I attached you an image where you can see
the summary of my two variables(sst & sst_avg), and I'm showing the
essential part of the script:


sst = a->anom(:,:,:)    ; times, lat ,lon
 sst at missing_value = 0.
 sst at _FillValue = 0.

 printVarSummary(sst)
 times=dimsizes(sst(:,0,0))
 row  = dimsizes(sst(0,:,0))
 col  = dimsizes(sst(0,0,:))
 sst_avg = new((/times,row,col/),"float",0)
 copy_VarCoords(sst,sst_avg)
 printVarSummary(sst_avg)
         do it =0,times-1
             sst_avg(0,:,:)=sst_avg(0,:,:)+sst(it,:,:)
         end do
 sst_avg(0,:,:)=sst_avg(0,:,:)/times    ; Calculo anomalia semanal promedio

 wks_type            = "png"
 wks_type at wkWidth    = 1700
 wks_type at wkHeight   = 1536
 wks = gsn_open_wks(wks_type,"animate")
 mpid    = gsn_csm_map(wks,mpres)
 ST  = gsn_csm_contour(wks,sst_avg(0,:,:),Tres)
;Overlaps maps in order
 overlay(mpid,ST)
 maximize_output(wks,True)


 Any suggestion to try new direction would be really appreciated.

Best regards,
Sebastián.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150806/2827446d/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: error_sst.PNG
Type: image/png
Size: 26655 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150806/2827446d/attachment.png