<div dir="ltr"><div><div><div><div>Hi Sebastián,<br></div>Try moving your copy_VarCoords line down beneath this line: <br>sst_avg(0,:,:)=sst_avg(0,:,:)/times<br><br></div>Also, note that there is a function that will average over a dimension:<br><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/dim_avg_n_Wrap.shtml">https://www.ncl.ucar.edu/Document/Functions/Contributed/dim_avg_n_Wrap.shtml</a><br><br></div>You could change all this:<br><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/),"float",0)</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 <br> copy_VarCoords(sst,sst_avg)<br> printVarSummary(sst_avg)<br><br></div><div>to this:<br></div><div>sst_avg = dim_avg_n_Wrap(sst,0)<br></div><br></div><div>Finally, note that setting _FIllValue = 0 is completely fine, but you will likely want to change that before attempting to plot the sst_avg array (using your original averaging coding) as some of the NCL plotting functions do not accept _FillValue = 0..<br><br></div><div>Hope that all helps. If it does not, or if you have further questions, please respond to the ncl-talk email list.<br></div><div>Adam<br></div><div><div><br><div><br><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 6, 2015 at 3:36 AM, Sebastian Otarola-Bustos <span dir="ltr"><<a href="mailto:Sebastian.F.Otarola-Bustos.1@nd.edu" target="_blank">Sebastian.F.Otarola-Bustos.1@nd.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi everybody, <br><br>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:<br><br><br><div>sst = a->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/),"float",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 = "png"</div><div> wks_type@wkWidth = 1700</div><div> wks_type@wkHeight = 1536</div><div><div> wks = gsn_open_wks(wks_type,"animate")</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>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>