<div dir="ltr"><div>Dear Dennis,</div><div>I put all NDVI files in one folder(like the attached image), and tried to run the attached script to plot the NDVI anomaly from 1986 to 2021.</div><div>However, I faced another issue(attached). May I ask you, please to advise me why this message appeared? How can I remove this issue, please?</div><div>Many thanks for considering my request.</div><div>Best wishes,<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 21, 2022 at 10:09 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">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>I downloaded 5 files for a test.</div><div><br></div><div>%> ncl NDVI_anomaly.ncl</div><div><br></div><div>Copyright (C) 1995-2017 - All Rights Reserved<br> University Corporation for Atmospheric Research<br> NCAR Command Language Version 6.5.0 <br> The use of this software is governed by a License Agreement.<br> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.<br><br><br>Variable: fils<br>Type: string<br>Total Size: 40 bytes<br> 5 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [5]<br>Coordinates: <br>(0) AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc<br>(1) AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc<br>(2) AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc<br>(3) AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc<br>(4) AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc<br><br>Variable: prc<br>Type: float<br>Total Size: 518400000 bytes<br> 129600000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 5] x [latitude | 3600] x [longitude | 7200]<br>Coordinates: <br> time: [1826..2190]<br> latitude: [89.975..-89.97499]<br> longitude: [-179.975..179.975]<br>Number Of Attributes: 7<br> long_name : NOAA Climate Data Record of Normalized Difference Vegetation Index<br> units : 1<br> grid_mapping : crs<br> standard_name : normalized_difference_vegetation_index<br> _FillValue_original : -9999<br> _FillValue : -9999</div><div>========================================================================</div><div>========================================================================</div><div>========================================================================</div><div>It seems to me that you have path or file problems.</div><div>I am not sure what else I can do.</div><div><br></div><div>Others?</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 21, 2022 at 5:49 AM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@gmail.com</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 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" target="_blank">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">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr">S.Rahimi<br><br></div>