<div dir="ltr">Hi everybody, <br><br>I&#39;m trying to plot a contour with week average of sst anomalies, but I&#39;m having that error. I read that could be because of the coordinates, but I&#39;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&#39;y need all of them( I don&#39;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 &amp; sst_avg), and I&#39;m showing the essential part of the script:<br><br><br><div>sst = a-&gt;anom(:,:,:)    ; times, lat ,lon</div><div> sst@missing_value = 0.</div><div> sst@_FillValue = 0.</div><div><br></div><div> printVarSummary(sst)   </div><div> times=dimsizes(sst(:,0,0))</div><div> row  = dimsizes(sst(0,:,0))</div><div> col  = dimsizes(sst(0,0,:))</div><div> sst_avg = new((/times,row,col/),&quot;float&quot;,0)</div><div> copy_VarCoords(sst,sst_avg)</div><div> printVarSummary(sst_avg)</div><div>         do it =0,times-1</div><div>             sst_avg(0,:,:)=sst_avg(0,:,:)+sst(it,:,:)</div><div>         end do</div><div> sst_avg(0,:,:)=sst_avg(0,:,:)/times    ; Calculo anomalia semanal promedio</div><div><br></div><div> wks_type            = &quot;png&quot;</div><div> wks_type@wkWidth    = 1700</div><div> wks_type@wkHeight   = 1536</div><div><div> wks = gsn_open_wks(wks_type,&quot;animate&quot;)</div><div> mpid    = gsn_csm_map(wks,mpres)</div><div> ST  = gsn_csm_contour(wks,sst_avg(0,:,:),Tres)</div><div>;Overlaps maps in order</div><div> overlay(mpid,ST)</div><div> maximize_output(wks,True)</div></div><div><br></div><div><br></div> Any suggestion to try new direction would be really appreciated.<br><br>Best regards,<br>Sebastián.</div>