<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div class="" style="">Denis,<o:p class="" style=""></o:p></div><div class="" style="">Thank you very much for all the work you have done for me.
This will provide me with learning material.<o:p class="" style=""></o:p></div><div class="" style="">!. I have noted the very important points you make about the
dangers of <span style="font-size: 10pt; font-family: 'Segoe UI', sans-serif;" class="">extrapolation. This is
clearly illustrated by results from the different extrapolation methods you
used. Only the fifth (prc_objanal) of the six methods resembling what one would
expect from a look at the original data distribution, and to a lesser extent,
the sixth (prc_triple).<o:p class="" style=""></o:p></span></div><div class="" style=""><span style="font-size: 10pt; font-family: 'Segoe UI', sans-serif;" class="">2. yes Grads does have
interpolation functions, but I did not try them (not really sure they would do what I want as I had the impression the functions deal with data on an already
regular grid).<o:p class="" style=""></o:p></span></div><div class="" style=""><span style="font-size: 10pt; font-family: 'Segoe UI', sans-serif;" class="">I feel deeply indebted to
your effort.<o:p class="" style=""></o:p></span></div><div class="" style=""><span class="" style="">
</span></div><div class="" style=""><span style="font-size: 10pt; font-family: 'Segoe UI', sans-serif;" class="">Zilore Mumba</span><o:p class="" style=""></o:p></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;" class=""> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;" class=""> <div dir="ltr" class="" style=""> <font size="2" face="Arial" class="" style=""> On Monday, July 28, 2014 9:33 PM, Dennis Shea <shea@ucar.edu> wrote:<br class="" style=""> </font> </div> <br class="" style=""><br class="" style=""> <div class="" style=""><div id="yiv4554107861" class="" style=""><div class="" style=""><div dir="ltr" class="" style=""><div class="" style=""><div class="" style=""><div class="" style="">I ttook a look at
this.<br clear="none" class="" style=""><br clear="none" class="" style="">[1]<br clear="none" class="" style=""></div>Obviously, there are not many data points. Further, they are not well distributed. Unless there is a physical basis to guide extrapolating values beyond the data boundaries, you should go by the rule: All data extrapolation is dangerous and should be avoided. Use of extrapolated values on subsequent analysis tasks contaminates the results. In particular by people who do not know the original random data distribution.<br clear="none" class="" style="">
<br clear="none" class="" style="">[2]<br clear="none" class="" style=""></div><div class="" style="">Different interpolation algorithms, will yield different results depending of what method (eg, splines) are being used.<br clear="none" class="" style=""><br clear="none" class="" style="">[3]<br clear="none" class="" style=""></div><div class="" style="">You mention "<span class="" style="">writing to binary and then contour in grads with which am more familiar.". I am sure GrADS has interpolation alogrithms for randomly based data. Not sure why these were not used.<br clear="none" class="" style="">
<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">Anyway, 'writing in binary' ... GrADS can read flat binary (access=direct) or fortran sequential (fortran default).<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">Did you go to NCL's IO functions? [Under Functions: Click Category; then, File Input/Output] You would see NCL's suite of IO functions. The attached script uses fbinrecwrite (fortran sequential)<br clear="none" class="" style="">
<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">The attached script shows results from several NCL interpolating (regridding) functions. The ESMF is very strict in the sense it will not extrapolate beyond data boundaries. All othe points are set to _Fillvalue.<br clear="none" class="" style="">
<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">Respond only to ncl-talk.<br clear="none" class="" style=""><br clear="none" class="" style="">===<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style=""><a rel="nofollow" shape="rect" target="_blank" href="http://www.ncl.ucar.edu/Document/Manuals/" class="" style="">http://www.ncl.ucar.edu/Document/Manuals/</a><br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">Download/Read: Mini-Language Manual<br clear="none" class="" style="">
<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style="">Good luck<br clear="none" class="" style=""></span></div><div class="" style=""><span class="" style=""><br clear="none" class="" style=""><br clear="none" class="" style="">
</span></div><div class="" style=""><br clear="none" class="" style=""></div><br clear="none" class="" style=""></div></div><div class="" style=""><br clear="none" class="" style=""><br clear="none" class="" style=""><div class="" style="">On Sun, Jul 27, 2014 at 2:32 AM, zilore mumba <span dir="ltr" class="" style=""><<a rel="nofollow" shape="rect" ymailto="mailto:zmumba@yahoo.com" target="_blank" href="mailto:zmumba@yahoo.com" class="" style="">zmumba@yahoo.com</a>></span> wrote:<br clear="none" class="" style="">
<blockquote class="" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="" id="yiv4554107861yqt11141" style=""><div class="" style=""><div style="color:#000;background-color:#fff;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt;" class="">
<div class="" style=""><span class="" style="">Dave,</span></div><div style="color:rgb(0,0,0);font-size:16px;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;background-color:transparent;font-style:normal;" class=""><span class="" style="">Once more thanks very much for your effort and time. You have moved me closer to the solution. Please see, in the attachment, the output I get from the script you sent me.</span></div>
<div style="color:rgb(0,0,0);font-size:16px;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;background-color:transparent;font-style:normal;" class=""><span class="" style="">I will work with this and see if I can move. I have seen the ESMF gridding examples, they can take me ages to decipher what to do.</span></div>
<div class="" style=""><div class="" style=""> <div class="" style=""><br clear="none" class="" style=""><br clear="none" class="" style=""></div><div style="display:block;" class=""> <div style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt;" class=""> <div style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt;" class="">
<div dir="ltr" class="" style=""> <font face="Arial" class="" style=""> On Sunday, July 27, 2014 3:13 AM, Dave Allured - NOAA Affiliate <<a rel="nofollow" shape="rect" ymailto="mailto:dave.allured@noaa.gov" target="_blank" href="mailto:dave.allured@noaa.gov" class="" style="">dave.allured@noaa.gov</a>> wrote:<br clear="none" class="" style=""> </font> </div> <br clear="none" class="" style="">
<br clear="none" class="" style=""> <div class="" style=""><div class="" style=""><div class="" style=""><div dir="ltr" class="" style="">Zilore,<br clear="none" class="" style=""><br clear="none" class="" style="">Sorry, that example #1 is not set up correctly for your kind of data. Try the attached script. This is another triangular mesh example that does use need extra coordinate arrays.<br clear="none" class="" style="">
<br clear="none" class="" style="">If you still want to regrid, then I found that ESMF regridding example works easily with your sample data.<br clear="none" class="" style="">
<br clear="none" class="" style=""><a rel="nofollow" shape="rect" target="_blank" href="http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex21" class="" style="">http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex21</a><br clear="none" class="" style=""><br clear="none" class="" style="">--Dave<br clear="none" class="" style="">
<br clear="none" class="" style=""><div class="" style=""><br clear="none" class="" style="">On Sat, Jul 26, 2014 at 5:24 PM, Dave Allured - NOAA Affiliate <<a rel="nofollow" shape="rect" ymailto="mailto:dave.allured@noaa.gov" target="_blank" href="mailto:dave.allured@noaa.gov" class="" style="">dave.allured@noaa.gov</a>> wrote:<br clear="none" class="" style="">
><br clear="none" class="" style="">> Zilore,<br clear="none" class="" style="">><br clear="none" class="" style="">> Please look at the first example on the Station Data example page:<br clear="none" class="" style="">> <a rel="nofollow" shape="rect" target="_blank" href="http://www.ncl.ucar.edu/Applications/station.shtml" class="" style="">http://www.ncl.ucar.edu/Applications/station.shtml</a><br clear="none" class="" style="">
><br clear="none" class="" style="">> You said that you wanted to plot contours. This simple example plots contours WITHOUT first putting data on a regular grid. It uses a triangular mesh plotting technique. For irregular coordinates, this technique will actually make a more accurate contour map than one made from a regular grid.<br clear="none" class="" style="">
><br clear="none" class="" style="">> Will this example meet your needs? My apologies if this was already discussed previously.<br clear="none" class="" style="">><br clear="none" class="" style="">> --Dave<br clear="none" class="" style="">><br clear="none" class="" style="">><br clear="none" class="" style="">
> On Sat, Jul 26, 2014 at 4:40 PM, zilore mumba <<a rel="nofollow" shape="rect" ymailto="mailto:zmumba@yahoo.com" target="_blank" href="mailto:zmumba@yahoo.com" class="" style="">zmumba@yahoo.com</a>> wrote:<br clear="none" class="" style="">
>><br clear="none" class="" style="">>> I have x,y,z data on an irregularly spaced geographical grid, which I want to put on a regular grid to be able to plot contours.<br clear="none" class="" style="">>> My data is below (id, long., lat., rainfall).<br clear="none" class="" style="">
>> 01 28.92 -2.47 286.5<br clear="none" class="" style="">
>> 02 29.02 -2.68 168.8<br clear="none" class="" style="">>> 03 29.25 -1.67 243.4<br clear="none" class="" style="">>> 04 29.42 -2.07 204.9<br clear="none" class="" style="">>> 05 29.55 -1.58 286.9<br clear="none" class="" style="">>> 06 29.57 -2.48 230.<br clear="none" class="" style="">
>> 07 29.6 -1.5 177.5<br clear="none" class="" style="">>> 08 29.72 -2.18 290.<br clear="none" class="" style="">
>> 09 30.05 -1.6 206.5<br clear="none" class="" style="">>> 10 30.13 -1.97 178.1<br clear="none" class="" style="">>> <a rel="nofollow" shape="rect" class="" style="" href="">11 30.5 -2.15 156.7</a><br clear="none" class="" style="">
>> My sincere apologies because this is a question I already asked more than 6 months ago and I still have not found a solution.<br clear="none" class="" style="">
>> The question is: Is it possible to do what I want to do in ncl?<br clear="none" class="" style="">>> 1. I have tried 2dgrid.ncl which seems to but am not able to figure out how to write output to file, and the z values all seem to be above 370 (well above what I have in my file).<br clear="none" class="" style="">
>> 2. I have try the Fortran single precision linear interpolation at Natgrid Examples and run the fortran program with ncargf77 filename I get the errors below:<br clear="none" class="" style="">>> Advice will be appreciated.<br clear="none" class="" style="">
>><br clear="none" class="" style="">
>> ncargf77 natgrid.f<br clear="none" class="" style="">>> gfortran -fPIC -fno-second-underscore -fno-range-check -fopenmp -O natgrid.f -L/usr/local/lib -L/usr/X11R6/lib64 -lncarg -lncarg_gks -lncarg_c -lX11 -lXext -lcairo -lfontconfig -lpixman-1 -lfreetype -lexpat -lpng -lz -liconv -lXrender -lbz2<br clear="none" class="" style="">
>> /usr/bin/ld: cannot find -lX11<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -lXext<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -lcairo<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -lfontconfig<br clear="none" class="" style="">
>> /usr/bin/ld: cannot find -lpixman-1<br clear="none" class="" style="">
>> /usr/bin/ld: cannot find -lfreetype<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -lexpat<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -liconv<br clear="none" class="" style="">>> /usr/bin/ld: cannot find -lXrender<br clear="none" class="" style="">
>> /usr/bin/ld: cannot find -lbz2<br clear="none" class="" style="">
>> collect2: error: ld returned 1 exit status<br clear="none" class="" style="">>><br clear="none" class="" style="">>> <br clear="none" class="" style="">>> <br clear="none" class="" style="">>> <br clear="none" class="" style="">>> <br clear="none" class="" style="">>> <br clear="none" class="" style="">
>> <br clear="none" class="" style="">>> Natgrid Examples<br clear="none" class="" style="">>> Natgrid examples Overview This module contains source codes for all of the examples referenced by the other modules in the documentation. Links to all of the example plots are provided.<br clear="none" class="" style="">
>> View on <a rel="nofollow" shape="rect" target="_blank" href="http://w3.pppl.gov/" class="" style="">w3.pppl.gov</a><br clear="none" class="" style="">>> Preview by Yahoo<br clear="none" class="" style=""><br clear="none" class="" style=""></div></div></div></div><br clear="none" class="" style=""><br clear="none" class="" style=""></div>
</div> </div> </div> </div></div></div></div></div><br clear="none" class="" style="">_______________________________________________<br clear="none" class="" style="">
ncl-talk mailing list<br clear="none" class="" style="">
List instructions, subscriber options, unsubscribe:<br clear="none" class="" style="">
<a rel="nofollow" shape="rect" target="_blank" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" class="" style="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br clear="none" class="" style="">
<br clear="none" class="" style=""></blockquote></div><br clear="none" class="" style=""></div></div></div><br class="" style=""><br class="" style=""></div> </div> </div> </div> </div></body></html>