<div dir="ltr"><div>Dear Dennis,</div><div>So many thanks for your advice. I run the script on a remote server with at least <u><span style="background-color:rgb(217,234,211)">1 terabyte</span></u> of memory. I think this is not the case.</div><div>I used the below commands to check the files:</div><div>===========================================<br></div><div><span style="background-color:rgb(255,217,102)">fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")<br>   print(fils)<br><br>   f   = addfiles (fils, "r")<br><br>   x   = short2flt(f[:]->NDVI)<br><br>   printVarSummary(x)</span><br>===========================================</div><div>, and attached the result. Please look at the image, and kindly help me sort the issue out. I do not know what is wrong with the files I downloaded.<br></div><div>best regards,<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Based on the provided information, I have no idea why you would get "returned nothing to me"</div><div><br></div><div>I know that NDVI files are 0.1 degree resolution (3600,7200). You have 36 years (1986-2021) of daily data.</div><div>So, 36*365=13140   [+ a few Leap Year days)<br></div><div><br></div><div>   prc(13140,3600,7200)    => daily data<br></div><div><br></div><div>A lot of memory is needed. Maybe your computer can not handle this?</div><div><br></div><div>After the <br></div><div>      xM = calculate_monthly_values (prc, "avg", 0,opt)</div><div>      printVarSummary(xM)</div><div><br></div><div>   xM(432,3600,7200)<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I need to calculate the <span style="background-color:rgb(241,194,50)">monthly anomaly of <font size="2"><span style="font-family:arial,sans-serif"><span style="font-weight:normal">the Normalized Difference Vegetation Index</span></span></font> (NDVI), from Advanced Very High-Resolution Radiometer (AVHRR) sensor</span>. I have downloaded AVHRR data (1986-2021) from the link below:</div><div><font size="2"><span style="font-weight:normal"><a href="https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html" target="_blank">https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html</a></span></font></div><div><font size="2"><span style="font-weight:normal">I first tried to prepare a script to combine all the data, and then calculate the monthly average from daily data ( the link provides daily data). After having monthly data I tried to figure out anomalies (the script has been attached). But NCL stopped and returned nothing to me. I just wonder what is wrong. Please kindly advise me on this issue.</span></font></div><div><font size="2"><span style="font-weight:normal">Best wishes,<br></span></font></div><div><font size="2"><span style="font-weight:normal"><br></span></font></div><div><h1><font size="2"><span style="font-family:arial,sans-serif"><span style="font-weight:normal"></span></span></font></h1><br>-- <br><div dir="ltr">S.Rahimi<br><br></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature">S.Rahimi<br><br></div>