<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Shea,<div class=""><br class=""></div><div class="">Thank you so much. I managed to&nbsp;<span style="color: rgb(51, 51, 51); font-family: 'Lucida Grande', 'Lucida Sans Unicode', arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255);" class="">concatenate my files into one big file using NCO.&nbsp;</span></div><div class=""><span style="color: rgb(51, 51, 51); font-family: 'Lucida Grande', 'Lucida Sans Unicode', arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255);" class=""><br class=""></span></div><div class=""><span style="color: rgb(51, 51, 51); font-family: 'Lucida Grande', 'Lucida Sans Unicode', arial, sans-serif; font-size: 13px; background-color: rgb(255, 255, 255);" class="">Dawit</span></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Apr 1, 2017, at 9:04 PM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu" class="">shea@ucar.edu</a>&gt; wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class="">[1] <br class="">Please subscribe to <a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class=""></div><div class="">No further 'applications' related questions will be answered on this list.<br class=""><br class="">[2] <br class=""></div><div class="">I have placed a version of the CRU data that spans all years at:<br class=""><br class=""></div><div class="">&nbsp; ftp <a href="http://ftp.cgd.ucar.edu/" class="">ftp.cgd.ucar.edu</a><br class=""></div><div class="">&nbsp; anonymous<br class=""></div><div class="">&nbsp; your_email<br class=""></div><div class="">&nbsp; cd pub/shea<br class=""></div><div class="">&nbsp; ls<br class=""></div><div class="">&nbsp; prompt<br class=""></div><div class="">&nbsp; get cru_ts3.23.1901.2014.pet.dat.nc.gz<br class=""></div><div class="">&nbsp; quit<br class=""><br class="">[3]<br class=""></div><div class="">NCL will not unzip the data. You must do that!<br class=""><br class=""></div><div class="">%&gt; gzip -d cru_ts3.23.1901.2014.pet.dat.nc.gz<br class=""><br class="">[4]<br class=""></div><div class="">BEFORE DOING ANYTHING: ***LOOK AT THE DATA FILE! ***<br class=""><br class=""></div><div class="">%&gt; ncdump -h <a href="http://cru_ts3.23.1901.2014.pet.dat.nc/" class="">cru_ts3.23.1901.2014.pet.dat.nc</a> | less<br class=""><br class=""></div><div class="">or<br class=""><br class="">%&gt; ncl_filedump <a href="http://cru_ts3.23.1901.2014.pet.dat.nc/" class="">cru_ts3.23.1901.2014.pet.dat.nc</a> | less<br class=""><br class=""></div><div class="">Look at dimension sizes; look at variable names; look at units; look at variable types; look at global attributes.<br class=""><br class="">After doing that .... look some more!!!!<br class=""></div><div class=""><br class="">[5] <br class=""></div><div class="">Extract **and*** examine extracted the variable.&nbsp; Use NCL syntax for coordinate variables:&nbsp; {...}<br class=""></div><div class=""><br class=""></div>&nbsp;&nbsp; latS = ...&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; latitudes span -89.87 to 89.75<br class=""></div>&nbsp;&nbsp; latN = ...<br class=""></div>&nbsp;&nbsp; lonL = ...&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; longitudes span -179.75 to 179.75<br class=""></div>&nbsp;&nbsp; lonR = ...<br class=""><br class="">&nbsp;&nbsp; f = addfile("<a href="http://cru_ts3.23.1901.2014.pet.dat.nc/" class="">cru_ts3.23.1901.2014.pet.dat.nc</a>" ,"r")<br class=""></div>&nbsp;&nbsp; x = f-&gt;pet(:,{latS:latN},{lonL:lonR})&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br class=""></div>&nbsp;&nbsp; printVarSummary(x)<br class=""></div>&nbsp;&nbsp; printMinMax(x)<br class=""></div>&nbsp;&nbsp; print("----")<br class=""><br class="">; <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml" class="">http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml</a> <br class=""><pre class="">  opt = True
  opt@PrintStat = True
  stats = <strong class="">stat_dispersion</strong>(x, opt )</pre>&nbsp;&nbsp; </div>[6] <br class="">Please read NCL documentation:&nbsp;&nbsp; <a href="http://www.ncl.ucar.edu/Document/" class="">http://www.ncl.ucar.edu/Document/</a><br class=""><br class=""></div>In particular: the following tutorial may be of some help to you.<br class="">&nbsp;&nbsp;&nbsp;&nbsp; <a href="http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/" class="">http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/</a><br class=""></div><br class=""></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Apr 1, 2017 at 10:24 AM, Dawit Abebe <span dir="ltr" class="">&lt;<a href="mailto:dawit.kan@gmail.com" target="_blank" class="">dawit.kan@gmail.com</a>&gt;</span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Hi Rick,<div class=""><br class=""></div><div class="">Thanks for trying to help but it is frustrating!&nbsp;</div><div class=""><br class=""></div><div class="">The data is from CRU version 3.23 precipitation downloaded by subsetting for the region I am interested with from the Ceda Oge Web Services/Web Processing Service (<a href="http://wps-web1.ceda.ac.uk/" target="_blank" class="">http://wps-web1.ceda.ac.uk/</a><a href="http://wps-web1.ceda.ac.uk/jobviewer?status_url=http://wps-web1.ceda.ac.uk/status/00fec680f736dbd74b73af2425e3966e" target="_blank" class="">)</a>. As a second option I download the entire dataset (Global) from <a href="http://data.ceda.ac.uk/badc/cru/data/cru_ts/cru_ts_3.23/data/pre/" target="_blank" class="">http://data.ceda.ac.uk/badc/<wbr class="">cru/data/cru_ts/cru_ts_3.23/<wbr class="">data/pre/</a>&nbsp; and add to NCL and subset the section I am interested using coordinate subsetting command. My problem with this is the file is in zip format (cru_ts3.23.1901.2014.pre.dat.<wbr class="">nc.gz) and I am wondering if it can be called to NCL as it is or need to be unzipped. If it should be unzipped, how to do it?</div><div class=""><br class=""></div><div class="">Do you have any advice how I should do this?</div><div class=""><br class=""></div><div class="">I would appreciate your help so much.</div><div class=""><br class=""></div><div class="">Dawit</div><div class="">&nbsp;</div><div class=""><br class=""></div><div class="">On Apr 1, 2017, at 7:03 PM, Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" target="_blank" class="">brownrig@ucar.edu</a>&gt; wrote:<div class=""><blockquote type="cite" class=""><br class="m_5095800277789372094Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Hi Dawit,<br class=""><br class=""></div><div class="">First, I apologize -- I overlooked that you state in your original message that you are getting the error for *all* the files, not just the 199501 file?<br class=""><br class="">In any case, I was able to open that file just fine, which would echo what you were saying about opening a file with the singular "addfile()" command.&nbsp; I was also able to open just that one file with the multi "addfiles()" command as well.&nbsp; Whats a bit puzzling is that the error message does not contain the leading "./" on the filename as it should.&nbsp; If the files are in your current working directory, that should not matter. Nonetheless it might be useful to print several things before the addfiles() statement, just to make sure the values are what we think they should be:<br class=""><br class=""></div><div class="">print(diri)<br class=""></div><div class="">print(files)<br class=""></div><div class="">print(diri + files)<br class=""><br class=""></div><div class="">I'm not sure what else to tell you. BTW: I'm using NCL 6.4.0; what version are you using?<br class=""><br class=""></div><div class="">Wish I could be of more help...<br class=""></div><div class="">Rick<br class=""></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Apr 1, 2017 at 9:34 AM, Dawit Abebe <span dir="ltr" class="">&lt;<a href="mailto:dawit.kan@gmail.com" target="_blank" class="">dawit.kan@gmail.com</a>&gt;</span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Dear Rick,<div class=""><br class=""></div><div class="">I think uploaded the file</div><div class=""><br class=""></div><div class="">Thanks for your help</div><span class="m_5095800277789372094HOEnZb"><font color="#888888" class=""><div class=""><br class=""></div><div class="">Dawit</div></font></span><div class=""><div class="m_5095800277789372094h5"><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Apr 1, 2017, at 6:09 PM, Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" target="_blank" class="">brownrig@ucar.edu</a>&gt; wrote:</div><br class="m_5095800277789372094m_-3709410322213034287Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class="">Yes, put refers to the filename you want to send; the command:<br class=""><br class=""></div>&nbsp;&nbsp;&nbsp; put <a href="http://pre_199501_mondat.nc/" target="_blank" class="">pre_199501_mondat.nc</a><br class=""><br class=""></div>will upload the file to the ftp server.<br class=""><br class=""><br class=""></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Apr 1, 2017 at 8:47 AM, Dawit Abebe <span dir="ltr" class="">&lt;<a href="mailto:dawit.kan@gmail.com" target="_blank" class="">dawit.kan@gmail.com</a>&gt;</span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Dear Rick,<div class=""><br class=""></div><div class="">Sorry I am bothering you, what is put…put…refers to, is it the files?</div><span class="m_5095800277789372094m_-3709410322213034287HOEnZb"><font color="#888888" class=""><div class=""><br class=""></div><div class="">Dawit</div><div class=""><br class=""></div></font></span><div class=""><br class=""><div class=""><blockquote type="cite" class=""><span class=""><span class=""><div class="">On Apr 1, 2017, at 5:13 PM, Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" target="_blank" class="">brownrig@ucar.edu</a>&gt; wrote:</div><br class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887Apple-interchange-newline"></span></span><div class=""><div class="h5"><div class=""><div class="m_5095800277789372094m_-3709410322213034287h5"><div class=""><div dir="ltr" class=""><div class=""><div class="">Hi,<br class=""><br class=""></div>It sounds like that one file in particular, <span style="font-variant-ligatures:no-common-ligatures" class="">"<a href="http://pre_199501_mondat.nc/" target="_blank" class="">pre_199501_mondat.nc</a>", is corrupt. I note that it is not the same file as the single file instance that you cited, and I *think* NCL would have ingested all the files from 1901 up to the 199501 before throwing the error. Are you able to open the file with the ncl_filedump command, or any of the NCO/CDO operators?&nbsp; If you want to send it to our ftp site, I'll take a look at it (instructions below).<br class=""><br class=""></span></div><span style="font-variant-ligatures:no-common-ligatures" class="">Rick<br class=""></span><div class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><br class=""></span><br class=""><b class="">You can ftp larger datasets and files to:</b><div class="">

<br class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887webkit-block-placeholder"></div><pre class="">ftp <a href="http://ftp.cgd.ucar.edu/" target="_blank" class="">ftp.cgd.ucar.edu</a>
anonymous
&lt;use your email address for the password&gt;
cd incoming
put ...
put ...



quit
</pre><p class="">

<font color="red" class="">Note that for security reasons, you cannot list the
contents of this directory, and neither can we.  We need to know
the <i class="">exact</i> names of the files in order to retrieve them.</font>
</p></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Apr 1, 2017 at 7:40 AM, Dawit Abebe <span dir="ltr" class="">&lt;<a href="mailto:dawit.kan@gmail.com" target="_blank" class="">dawit.kan@gmail.com</a>&gt;</span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Dear NCL users,<div class=""><br class=""></div><div class="">I’m having a challenge to add multiple files to NCL. I’m intending to work on monthly files ranging from 1901 - 2014. I tried the command below but failed. &nbsp;</div><div class=""><br class=""></div><div class=""><font color="#0433ff" class="">Command used</font></div><div class=""><br class=""></div><div class=""><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""><span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>diri = "./"&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 29&gt; files = systemfunc("ls *.nc")&nbsp; &nbsp;</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 30&gt; f &nbsp; &nbsp; = addfiles(diri+files, "r")</span></div></div><div class=""><br class=""></div><div class=""><font color="#0433ff" class="">Return message for all files</font></div><div class=""><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><br class=""></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">warning:_NclOpenFile: Can not open file &lt;<a href="http://pre_199501_mondat.nc/" target="_blank" class="">pre_199501_mondat.nc</a>&gt;; file format not supported or file is corrupted</span></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><font color="#0433ff" class="">The command and return for single file seems working as showed bellow</font>:&nbsp;</div><div class=""><br class=""></div><div class=""><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 10&gt; diri = "./"</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 11&gt; fili = "<a href="http://pre_201412_mondat.nc/" target="_blank" class="">pre_201412_mondat.nc</a>" &nbsp;</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 12&gt; f&nbsp; &nbsp; = addfile(diri+fili, "r")</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 13&gt; prc = f-&gt;pre &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">ncl 14&gt; printVarSummary(prc)</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;min-height:13px" class=""><span style="font-variant-ligatures:no-common-ligatures" class=""></span><br class=""></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Variable: prc</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Type: double</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Total Size: 288 bytes</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 36 values</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Number of Dimensions: 3</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Dimensions and sizes:<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>[time | 1] x [lat | 4] x [lon | 9]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Coordinates:&nbsp;</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; time: [41987..41987]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; lat: [3.75..5.25]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; lon: [36.75..40.75]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">Number Of Attributes: 5</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; correlation_decay_distance :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>&nbsp;450</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; long_name :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>precipitation</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; units :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>mm/month</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; missing_value :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>9.96920996839e+36</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo" class=""><span style="font-variant-ligatures:no-common-ligatures" class="">&nbsp; _FillValue :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap">        </span>9.96920996839e+36</span></div><div class=""><br class=""></div></div><div class=""><font color="#0433ff" class="">Where is my problem?</font></div><div class=""><font color="#0433ff" class=""><br class=""></font></div><div class=""><font color="#0433ff" class="">Thanks</font></div><span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887HOEnZb"><font color="#888888" class=""><div class=""><font color="#0433ff" class="">Dawit</font></div></font></span></div></div><br class="">______________________________<wbr class="">_________________<br class="">
ncl-talk mailing list<br class="">
<a href="mailto:ncl-talk@ucar.edu" target="_blank" class="">ncl-talk@ucar.edu</a><br class="">
List instructions, subscriber options, unsubscribe:<br class="">
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/mailma<wbr class="">n/listinfo/ncl-talk</a><br class="">
<br class=""></blockquote></div><br class=""></div>
</div></div></div></div></div></blockquote></div><br class=""></div></div></blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></div></div></div></div></blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></div></div><br class="">______________________________<wbr class="">_________________<br class="">
ncl-install mailing list<br class="">
List instructions, subscriber options, unsubscribe:<br class="">
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-install" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/<wbr class="">mailman/listinfo/ncl-install</a><br class="">
<br class=""></blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></div></body></html>