<div dir="ltr"><div class="gmail_default" style="font-size:small">I believe the error is telling you that you have some duplicate lat/lon points.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Try removing any duplicate points before you call cssgrid. You can also average the points if they are different values, but this is something that you need to determine what works best for your data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"> </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 9, 2015 at 7:40 AM, Adv <span dir="ltr"><<a href="mailto:advita6@gmail.com" target="_blank">advita6@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 dir="ltr"><div><div>Hi,<br></div>Could someone help me in this? Its quite confusing. I would like to compute objective analysis of all the station data sets to plot them spatially..<br><br></div>Thank you<br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 8, 2015 at 2:15 PM, Adv <span dir="ltr"><<a href="mailto:advita6@gmail.com" target="_blank">advita6@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 dir="ltr"><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Adv</b> <span dir="ltr"><<a href="mailto:advita6@gmail.com" target="_blank">advita6@gmail.com</a>></span><br>Date: Tue, Dec 8, 2015 at 2:10 PM<br>Subject: sub: css grid and objective analysis...<br>To: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br><br><br><div dir="ltr"><div><div>Hi,<br></div>I have been trying to compute objective analysis of station data sets. <br><br>begin<br> strll = asciiread("stationmac.txt", -1, "string") ;<br> print(strll) ; STATION","STATION_NAME", "LAT", "LON"<br> print("=====")<br><br>sloc = str_get_field(strll , 2, ",")<br> lat = tofloat( str_get_field(strll , 3, ",") )<br> lon = tofloat( str_get_field(strll , 4, ",") )<br> print(sloc +" : "+lat+" "+lon)<br><br> NOBS = 65<br>diri="/home/dailysimms/prec-obs/01/"<br>z1 = asciiread(diri+"PrecJan1965-2014",-1,"float")<br>print(z1)<br>printVarSummary(z1)<br>;a = (/1,2,3,4,5,6,7,8/)<br>;391 rows by 50 columns<br>a0 = onedtond(z1,(/391,50/))<br>;print(a0)<br>x=dim_avg(a0)<br><br>;------------------------------------------------------------------<br>; Read data<br>;------------------------------------------------------------------<br><br> diri = "./"<br> fili = "Ecearth320x160x1961-2005.nc"<br> f = addfile(diri+fili, "r")<br> Z = f->t2(0,{40:49},{242:264})<br> Z1 = f->t2(0,{40:49},{242:264})*86400<br>printVarSummary(Z)<br><br>;------------------------------------------------------------------<br>; Generate N random values from the above<br>;------------------------------------------------------------------<br> dimz = dimsizes( Z )<br>nlat = dimz(0)<br> mlon = dimz(1)<br> nx = 8<br> ny = 19<br> plat = 90.*fspan(-1.,1.,nx)<br> plon = 180.*fspan(-1.,1.,ny)<br> N = nlat*mlon<br> iZ = generate_unique_indices( N )<br> trip = cssgrid(lat, lon, x,nlat,mlon)<br>dlon=0.1<br> reps = random_uniform (-dlon, dlon, N) ; minor location perturbations<br><br> IOBS = iZ(0:NOBS) ; convenience<br> ; clarity only ... create explicit variables<br> rlon = trip(0,IOBS) + reps(IOBS)<br> rlat = trip(1,IOBS) + reps(IOBS(::-1))<br> rZ = trip(2,IOBS)<br>;******************************************************************<br>; call Wrap version of obj_anal_ic<br>;******************************************************************<br><br> opt = True<br> opt@timing = True<br><br> dcrit = (/45,35,25, 15,10 /)<br> grid = obj_anal_ic_Wrap(lon, lat, x, rlon, rlat, dcrit, opt)<br>printVarSummary(grid)<br>print(grid)<br>return<br><br>ERROR:<br><br>CSSGRID - error -319, coordinates 319 and M coincide for some M > 319 >= 1 (coordinate numbering starting at 1).<br>Error number -319 returned from c_cssgridd<br>fatal:cssgrid: Error number -319.<br>fatal:["Execute.c":8575]:Execute: Error occurred at or near line 45 in file oban.ncl<br><br></div><div>Could someone help me to fix this bug? <br></div><div>Thank you<span><font color="#888888"><span><font color="#888888"><br></font></span></font></span></div><span><font color="#888888"><span><font color="#888888"><div>Adv<br></div></font></span></font></span></div>
</div><br></div>
</blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>