<div dir="ltr"><div>Hi Dave,</div>yes, I've been plotting the 1st time step only for test purposes. Maybe I phrased it unluckily, I call it tice because in the end I am plotting on the t-grid, but the variable I am reading in is aice. How do you think could I avoid reading in missing values then?<div>Thanks for getting back to me,</div><div>Anne</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 15 February 2018 at 17:50, David Bailey <span dir="ltr"><<a href="mailto:dbailey@ucar.edu" target="_blank">dbailey@ucar.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"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi Anne,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I can try to help with this. What if you just plot "aice" as it is? Your variables tice and var_dtrend are three dimensional. Are you just plotting a single time slice? I think you are introducing some missing values with your var_dtrend computation.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Dave</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Thu, Feb 15, 2018 at 9:30 AM, Anne <span dir="ltr"><<a href="mailto:anne.seidenglanz@unive.it" target="_blank">anne.seidenglanz@unive.it</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hello,<div>I have the problem that I can perfectly plot <i>tice</i> from CICE output (based on <i>aice,</i> I have used one of the special ncl scripts for this purpose), but not anymore after I detrended it (dtrend; along with it using both <span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><i>copy_VarMeta</i> and <i>copy_VarCoord</i> of the original in order to keep the metadata). The plot (attached) looks distorted.</span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><br></div><div>I think the error is in the <b>lat2d</b> and <b>lon2d</b>, when doing printVarSummary:</div><div><br></div><div><div>Variable: <b>tice</b></div><div>Type: float</div><div>Total Size: 846028800 bytes</div><div> 211507200 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[time | 144] x [y | 1020] x [x | 1440]</div><div>Coordinates:</div><div> time: [727779..736325]</div><div>Number Of Attributes: 8</div><div> long_name :<span style="white-space:pre-wrap"> </span>ice area (aggregate)</div><div> units :<span style="white-space:pre-wrap"> </span>%</div><div> coordinates :<span style="white-space:pre-wrap"> </span>TLON TLAT</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>1e+30</div><div> missing_value :<span style="white-space:pre-wrap"> </span>1e+30</div><div> comment :<span style="white-space:pre-wrap"> </span>none</div><div> cell_methods :<span style="white-space:pre-wrap"> </span>time: mean</div><div> time_rep :<span style="white-space:pre-wrap"> </span>averaged</div></div><div><br></div><div>and printVarSummary od <i>var_dtrend:</i></div><div><i><br></i></div><div><div>Variable: <b>var_dtrend</b></div><div>Type: float</div><div>Total Size: 846028800 bytes</div><div> 211507200 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[144] x [1020] x [1440]</div><div>Coordinates:</div><div>Number Of Attributes: 10</div><div> missing_value :<span style="white-space:pre-wrap"> </span>1e+30</div><div><b> lat2d :<span style="white-space:pre-wrap"> </span><ARRAY of 1468800 elements></b></div><div><b> lon2d :<span style="white-space:pre-wrap"> </span><ARRAY of 1468800 elements></b></div><div> long_name :<span style="white-space:pre-wrap"> </span>ice area (aggregate)</div><div> units :<span style="white-space:pre-wrap"> </span>%</div><div> coordinates :<span style="white-space:pre-wrap"> </span>TLON TLAT</div><div> comment :<span style="white-space:pre-wrap"> </span>none</div><div> cell_methods :<span style="white-space:pre-wrap"> </span>time: mean</div><div> time_rep :<span style="white-space:pre-wrap"> </span>averaged</div><div> _FillValue :<span style="white-space:pre-wrap"> </span>1e+30</div></div><div><br></div><div>lat2d and lon2d are extra listed here again, and I don't know why. I have also tried to add </div><div><div>var_dtrend@lon2d = tlon2d</div><div>var_dtrend@lat2d = tlat2d</div></div><div><br></div><div>but this didn't solve the issue. I also attached the script. </div><div><br></div><div>Is the simplest thing just to interpolate it anyway?</div><div><br></div><div>Thanks, any help appreciated!</div><div>Anne</div></div>
<br>
</div></div><div><font style="font-family:Arial,Helvetica,sans-serif" size="2">Nota automatica aggiunta dal sistema di posta.</font></div><div><img src="https://www.unive.it/media/banner_150.jpg"></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_5216819926298395372gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">
<p>David A Bailey email: <a href="http://dbailey_at_ucar.edu" target="_blank">dbailey_at_ucar.edu</a><br>National Center for Atmospheric Research phone: <a href="tel:(303)%20497-1737" value="+13034971737" target="_blank">303-497-1737</a><br>PO Box 3000 fax : <a href="tel:(303)%20497-1700" value="+13034971700" target="_blank">303-497-1700</a><br>Boulder, CO 80307-3000<br> <a href="http://www.cgd.ucar.edu/staff/dbailey" target="_blank">http://www.cgd.ucar.edu/<wbr>staff/dbailey</a></p></div></div></div></div>
</font></span></div>
</blockquote></div><br></div>
<br>
<div><font style="font-family:Arial,Helvetica,sans-serif" size="2">Nota automatica aggiunta dal sistema di posta.</font></div><div><img src="https://www.unive.it/media/banner_150.jpg"></div>