<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>[1] <br>Please subscribe to <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br></div><div>No further 'applications' related questions will be answered on this list.<br><br>[2] <br></div><div>I have placed a version of the CRU data that spans all years at:<br><br></div><div> ftp <a href="http://ftp.cgd.ucar.edu">ftp.cgd.ucar.edu</a><br></div><div> anonymous<br></div><div> your_email<br></div><div> cd pub/shea<br></div><div> ls<br></div><div> prompt<br></div><div> get cru_ts3.23.1901.2014.pet.dat.nc.gz<br></div><div> quit<br><br>[3]<br></div><div>NCL will not unzip the data. You must do that!<br><br></div><div>%> gzip -d cru_ts3.23.1901.2014.pet.dat.nc.gz<br><br>[4]<br></div><div>BEFORE DOING ANYTHING: ***LOOK AT THE DATA FILE! ***<br><br></div><div>%> ncdump -h <a href="http://cru_ts3.23.1901.2014.pet.dat.nc">cru_ts3.23.1901.2014.pet.dat.nc</a> | less<br><br></div><div>or<br><br>%> ncl_filedump <a href="http://cru_ts3.23.1901.2014.pet.dat.nc">cru_ts3.23.1901.2014.pet.dat.nc</a> | less<br><br></div><div>Look at dimension sizes; look at variable names; look at units; look at variable types; look at global attributes.<br><br>After doing that .... look some more!!!!<br></div><div><br>[5] <br></div><div>Extract **and*** examine extracted the variable. Use NCL syntax for coordinate variables: {...}<br></div><div><br></div> latS = ... ; latitudes span -89.87 to 89.75<br></div> latN = ...<br></div> lonL = ... ; longitudes span -179.75 to 179.75<br></div> lonR = ...<br><br> f = addfile("<a href="http://cru_ts3.23.1901.2014.pet.dat.nc">cru_ts3.23.1901.2014.pet.dat.nc</a>" ,"r")<br></div> x = f->pet(:,{latS:latN},{lonL:lonR}) <br></div> printVarSummary(x)<br></div> printMinMax(x)<br></div> print("----")<br><br>; <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml">http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml</a> <br><pre> opt = True
opt@PrintStat = True
stats = <strong>stat_dispersion</strong>(x, opt )</pre> </div>[6] <br>Please read NCL documentation: <a href="http://www.ncl.ucar.edu/Document/">http://www.ncl.ucar.edu/Document/</a><br><br></div>In particular: the following tutorial may be of some help to you.<br> <a href="http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/">http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/</a><br></div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 1, 2017 at 10:24 AM, Dawit Abebe <span dir="ltr"><<a href="mailto:dawit.kan@gmail.com" target="_blank">dawit.kan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi Rick,<div><br></div><div>Thanks for trying to help but it is frustrating! </div><div><br></div><div>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">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">)</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">http://data.ceda.ac.uk/badc/<wbr>cru/data/cru_ts/cru_ts_3.23/<wbr>data/pre/</a> 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>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><br></div><div>Do you have any advice how I should do this?</div><div><br></div><div>I would appreciate your help so much.</div><div><br></div><div>Dawit</div><div> </div><div><br></div><div>On Apr 1, 2017, at 7:03 PM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<div><blockquote type="cite"><br class="m_5095800277789372094Apple-interchange-newline"><div><div dir="ltr"><div>Hi Dawit,<br><br></div><div>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><br>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. I was also able to open just that one file with the multi "addfiles()" command as well. Whats a bit puzzling is that the error message does not contain the leading "./" on the filename as it should. 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><br></div><div>print(diri)<br></div><div>print(files)<br></div><div>print(diri + files)<br><br></div><div>I'm not sure what else to tell you. BTW: I'm using NCL 6.4.0; what version are you using?<br><br></div><div>Wish I could be of more help...<br></div><div>Rick<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 1, 2017 at 9:34 AM, Dawit Abebe <span dir="ltr"><<a href="mailto:dawit.kan@gmail.com" target="_blank">dawit.kan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Rick,<div><br></div><div>I think uploaded the file</div><div><br></div><div>Thanks for your help</div><span class="m_5095800277789372094HOEnZb"><font color="#888888"><div><br></div><div>Dawit</div></font></span><div><div class="m_5095800277789372094h5"><div><br><div><blockquote type="cite"><div>On Apr 1, 2017, at 6:09 PM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:</div><br class="m_5095800277789372094m_-3709410322213034287Apple-interchange-newline"><div><div dir="ltr"><div><div>Yes, put refers to the filename you want to send; the command:<br><br></div> put <a href="http://pre_199501_mondat.nc/" target="_blank">pre_199501_mondat.nc</a><br><br></div>will upload the file to the ftp server.<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 1, 2017 at 8:47 AM, Dawit Abebe <span dir="ltr"><<a href="mailto:dawit.kan@gmail.com" target="_blank">dawit.kan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear Rick,<div><br></div><div>Sorry I am bothering you, what is put…put…refers to, is it the files?</div><span class="m_5095800277789372094m_-3709410322213034287HOEnZb"><font color="#888888"><div><br></div><div>Dawit</div><div><br></div></font></span><div><br><div><blockquote type="cite"><span class=""><span><div>On Apr 1, 2017, at 5:13 PM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:</div><br class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887Apple-interchange-newline"></span></span><div><div class="h5"><div><div class="m_5095800277789372094m_-3709410322213034287h5"><div><div dir="ltr"><div><div>Hi,<br><br></div>It sounds like that one file in particular, <span style="font-variant-ligatures:no-common-ligatures">"<a href="http://pre_199501_mondat.nc/" target="_blank">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? If you want to send it to our ftp site, I'll take a look at it (instructions below).<br><br></span></div><span style="font-variant-ligatures:no-common-ligatures">Rick<br></span><div><span style="font-variant-ligatures:no-common-ligatures"><br></span><br><b>You can ftp larger datasets and files to:</b><div>
<br class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887webkit-block-placeholder"></div><pre>ftp <a href="http://ftp.cgd.ucar.edu/" target="_blank">ftp.cgd.ucar.edu</a>
anonymous
<use your email address for the password>
cd incoming
put ...
put ...
.
.
.
quit
</pre><p>
<font color="red">Note that for security reasons, you cannot list the
contents of this directory, and neither can we. We need to know
the <i>exact</i> names of the files in order to retrieve them.</font>
</p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Apr 1, 2017 at 7:40 AM, Dawit Abebe <span dir="ltr"><<a href="mailto:dawit.kan@gmail.com" target="_blank">dawit.kan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Dear NCL users,<div><br></div><div>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. </div><div><br></div><div><font color="#0433ff">Command used</font></div><div><br></div><div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"><span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap"> </span>diri = "./" </span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 29> files = systemfunc("ls *.nc") </span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 30> f = addfiles(diri+files, "r")</span></div></div><div><br></div><div><font color="#0433ff">Return message for all files</font></div><div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><br></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">warning:_NclOpenFile: Can not open file <<a href="http://pre_199501_mondat.nc/" target="_blank">pre_199501_mondat.nc</a>>; file format not supported or file is corrupted</span></div><div><br></div><div><br></div><div><font color="#0433ff">The command and return for single file seems working as showed bellow</font>: </div><div><br></div><div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 10> diri = "./"</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 11> fili = "<a href="http://pre_201412_mondat.nc/" target="_blank">pre_201412_mondat.nc</a>" </span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 12> f = addfile(diri+fili, "r")</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 13> prc = f->pre </span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">ncl 14> printVarSummary(prc)</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Variable: prc</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Type: double</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Total Size: 288 bytes</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> 36 values</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Number of Dimensions: 3</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">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"><span style="font-variant-ligatures:no-common-ligatures">Coordinates: </span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> time: [41987..41987]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> lat: [3.75..5.25]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> lon: [36.75..40.75]</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures">Number Of Attributes: 5</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> correlation_decay_distance :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap"> </span> 450</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="font-variant-ligatures:no-common-ligatures"> 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"><span style="font-variant-ligatures:no-common-ligatures"> 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"><span style="font-variant-ligatures:no-common-ligatures"> 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"><span style="font-variant-ligatures:no-common-ligatures"> _FillValue :<span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887m_-6484431340853010663Apple-tab-span" style="white-space:pre-wrap"> </span>9.96920996839e+36</span></div><div><br></div></div><div><font color="#0433ff">Where is my problem?</font></div><div><font color="#0433ff"><br></font></div><div><font color="#0433ff">Thanks</font></div><span class="m_5095800277789372094m_-3709410322213034287m_-1755242649157655887HOEnZb"><font color="#888888"><div><font color="#0433ff">Dawit</font></div></font></span></div></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><br></div>
</div></div></div></div></div></blockquote></div><br></div></div></blockquote></div><br></div>
</div></blockquote></div><br></div></div></div></div></blockquote></div><br></div>
</div></blockquote></div><br></div></div><br>______________________________<wbr>_________________<br>
ncl-install mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-install" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-install</a><br>
<br></blockquote></div><br></div>