<div dir="ltr"><div><div><div><div><div>printVarSummary(aot_subset), shows that aot_subset is undefined, which I am unable to fix this problem yet.<br><br></div></div>If I used aot_subset(ilt1:ilt2,iln1:iln2) <span style="color:rgb(255,0,0)"><b>=</b></span> data(ilt1:ilt2,iln1:iln2) <span style="color:rgb(204,0,0)"><b>instead of </b></span>aot_subset(ilt1:ilt2,iln1:iln2)<span style="color:rgb(255,0,0)"> </span><span style="color:rgb(0,0,255)"><b><span style="color:rgb(255,0,0)">:=</span> </b></span>data(ilt1:ilt2,iln1:iln2)<br></div><br>then the error appeared as <br><br><span style="color:rgb(255,0,0)"><b>fatal:["Execute.c":7950]:Number of subscripts on right-hand-side do not match<br> number of dimensions of variable: (2), Subscripts used: (3)<br><br></b></span></div><span style="color:rgb(255,0,0)"><font color="#000000">basically, the data summary is <br><br>Variable: fnames<br>Type: string<br>Total Size: 8 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates: <br>(0) ./<a href="http://maiactaotE.h00v03.20000550525-20003660700.nc">maiactaotE.h00v03.20000550525-20003660700.nc</a><br><br>Variable: data<br>Type: float<br>Total Size: 2545920000 bytes<br> 636480000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: <b>[time | 442] x [YDim_grid1km | 1200] x [XDim_grid1km | 1200]</b><br>Coordinates: <br> time: [1301.416666666686..8767]<br>Number Of Attributes: 8<br> coordinates : GridLat_grid1km, GridLon_grid1km<br> hdfeos_name : Optical_Depth_055<br> projection : Albers Conical Equal_Area<br> unit : None<br> long_name : AOT at 0.55 micron<br> _FillValue_original : -28672<br> _FillValue : -28672<br> valid_range : ( -0.1, 5 )<br>(0) latitude (degrees_north) : min=11.2940312768897 max=24.54921980985786<br>(0) longitude (degrees_east) : min=68.08003670367275 max=81.63171349927691<br><b>fatal:["Execute.c":7776]:aot_subset is undefined, can not subscript an undefined variable<br>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 38 in file test_extract.ncl<br><br></b></font></span></div><span style="color:rgb(255,0,0)"><font color="#000000"><b><span style="color:rgb(255,0,0)">38 line is </span></b></font></span><br><span style="color:rgb(255,0,0)"><font color="#000000"><b><span style="color:rgb(255,0,0)">aot_subset(ilt1:ilt2,iln1:iln2) <b>=</b> data(ilt1:ilt2,iln1:iln2)</span><br><br></b></font></span><div><span style="color:rgb(255,0,0)"><font color="#000000"><b>fatal:Variable (aot_subset) is undefined<br></b><br></font></span><div><div><div><div><span style="color:rgb(255,0,0)"><b></b></span><div><div><div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div></div><div class="gmail_extra">regards<br clear="all"></div><div class="gmail_extra"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><p style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div>
<br><div class="gmail_quote">On Tue, Oct 17, 2017 at 3:05 PM, Guido Cioni <span dir="ltr"><<a href="mailto:guidocioni@gmail.com">guidocioni@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Why did you comment this line ?<div><br></div><div><blockquote type="cite"><div dir="ltr"><div><div> ;printVarSummary(aot_subset)<br></div></div></div></blockquote><div><br></div><div>That would have helped you in debugging the error on your own.</div><div><br></div><div>You are using the re-assignment operator ":=". This means that every time you make an assignment the variable is created again. Thus, you cannot make this assignment </div><span class="gmail-"><div><br></div><div><blockquote type="cite"><div dir="ltr"><div><div> <b> <span style="color:rgb(255,0,0)">aot_subset</span>(ilt1:ilt2,iln1:<wbr>iln2) := data(ilt1:ilt2, iln1:iln2)</b><br></div></div></div></blockquote><br></div></span><div>as NCL has no way to know the shape of <font face="Menlo">aot_subset</font> beforehand. </div><div><br></div><div>The following should be used instead (notice the assignment with "=" in the second line).</div><div><br></div><div><font face="Menlo"> aot_subset := data</font></div><div><font face="Menlo"> aot_subset = aot_subset@_FillValue<span class="gmail-"><br> aot_subset(ilt1:ilt2,iln1:<wbr>iln2) <b>=</b> data(ilt1:ilt2, iln1:iln2)</span></font></div><div><br></div><div>Notice that this will (1) create an array which is shaped like data (2) fill it with fill values as you only want to compute the average over the subset (3) populate the array only in the lat-lon box provided. However I don't see why you would want to use such a complicated syntax when the following code would produce the same result</div><div><br></div><div><font face="Menlo">aot_subset_avg=avg(data(ilt1:<wbr>ilt2, iln1:iln2))</font></div><div><br></div><div>If you have the right coordinate attributes you can even jump the first 7 lines and do </div><div><br></div><div><font face="Menlo">aot_subset_avg=avg(data({18:<wbr>24},{74:78})</font></div><div><br></div><div>but I don't know your data. </div><div><br></div><div>As it was already suggested by Mary and others, PLEASE take time to debug the errors on your own. </div><div><br></div><div>Regards,</div><div><div class="gmail-h5"><br><div><blockquote type="cite"><div>On 17. Oct 2017, at 11:25, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com">kunal.bali9@gmail.com</a>> wrote:</div><br class="gmail-m_475446440518907222Apple-interchange-newline"><div><div dir="ltr"><div><div><br></div>I did the following changes but still getting an error of <span style="color:rgb(255,0,0)"><b>aot_subset</b></span> which is undefined. <br></div>please let me know that how to fix or defined the <span style="color:rgb(255,0,0)"><b>aot_subset </b></span>this problem.<br><div><div><br><br><br> latv = (/ 18, 24/) ; for extracting grid box data <br> lonv = (/ 74, 78 /)<br> nm = getind_latlon2d (lat2d, lon2d, latv, lonv)<br> ilt1 = nm(0,0) ; start lat index<br> ilt2 = nm(1,0) ; start lon index<br> iln1 = nm(0,1) ; end lat index<br> iln2 = nm(1,1) ; end lon index<br> lat_subset := lat2d(ilt1:ilt2,iln1:iln2)<br> lon_subset := lon2d(ilt1:ilt2,iln1:iln2)<br> aot_subset := data<br> <b> <span style="color:rgb(255,0,0)">aot_subset</span>(ilt1:ilt2,iln1:<wbr>iln2) := data(ilt1:ilt2, iln1:iln2)</b><br> <br> ;printVarSummary(aot_subset)<br> aot_avg = avg(data)<br> aot_subset_avg = avg(<span style="color:rgb(255,0,0)"><b>aot_subset</b></span>)<br> <br> ;print(nm)<br> print("-----")<br><br> do k=0,dimsizes(latv)-1<br> n = nm(k,0)<br> m = nm(k,1)<br> print(lat2d(n,m)+" "+lon2d(n,m)) <br> print(" "+data(:,n,m)) <wbr> ; nearest grid point at all time steps<br> print("-----")<br> end do<br><br><br><br></div><div><br></div></div><div class="gmail_extra">regards<br clear="all"></div><div class="gmail_extra"><div><div class="gmail-m_475446440518907222gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><br></div><div><div style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Tue, Oct 17, 2017 at 12:41 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com">kunal.bali9@gmail.com</a>></span> wrote:<br><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><div>I have sent the ncl script and data to the incoming dir with the given names as <br><br>data = maiactaotE.h00v03.20003660525-<wbr>20003660700.nc.bz2<br></div>ncl_script = extract_iiserb.ncl<br></div><div><br></div><div>I am not able to define the var_subset in the given script. <br></div>please let me know if you have received the data and script?<span class="gmail-m_475446440518907222gmail-HOEnZb"><font color="#888888"><br><br><br><div><br><br></div></font></span><div class="gmail_extra"><span class="gmail-m_475446440518907222gmail-HOEnZb"><font color="#888888"><br clear="all"><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><br></div><div><div style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></div></div></div></div></div></div></div></font></span><div><div class="gmail-m_475446440518907222gmail-h5">
<br><div class="gmail_quote">On Mon, Oct 16, 2017 at 11:07 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>></span> wrote:<br><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 style="font-size:small">Kunal,</div><div style="font-size:small"><br></div><div style="font-size:small">You must make an effort to debug some of these problems on your own. The error is telling you that something is wrong around line 34 and it is telling you exactly what the problem is. Please try to fix this error.</div><div style="font-size:small"><br></div><div style="font-size:small">If you continue to have problems after fixing this error, then I need to ask you to provide the full script and data via ftp:</div><div style="font-size:small"><br></div><div><a href="http://www.ncl.ucar.edu/ftp_files.shtml">http://www.ncl.ucar.edu/ftp_fi<wbr>les.shtml</a><br></div><div><br></div><div>However, make sure you debug any other new problems that come up, instead of just emailing them to ncl-talk. We have a page to help you understand some of NCL's errors, and what causes them:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Document/Language/error_messages.shtml">http://www.ncl.ucar.edu/Docume<wbr>nt/Language/error_messages.sht<wbr>ml</a><br></div><div><br></div><div>The particular error you're seeing can be found here:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Document/Language/error_messages.shtml#NumSubRHS">http://www.ncl.ucar.edu/Docume<wbr>nt/Language/error_messages.sht<wbr>ml#NumSubRHS</a><br></div><div><br></div><div style="font-size:small">As a side note, you are running a 3-year old version of NCL. I suggest that you upgrade to NCL V6.4.0.</div><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634HOEnZb"><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div><div style="font-size:small"><br></div></font></span></div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634HOEnZb"><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 16, 2017 at 12:52 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com">kunal.bali9@gmail.com</a>></span> wrote:<br><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 simply include the code (which you mentioned) were I was taking the avg of the data as given below (please correct me if I did wrong)<br><br><br>;-----------------------------<wbr>------------------------------<wbr>-----------<span><br>; Main code<br>;-----------------------------<wbr>------------------------------<wbr>-----------<br> dir = "./"<br> fnames = systemfunc("ls " + dir + "maiactaot*.nc")<br>print(fnames)<br> a = addfile(fnames,"r")<br> data = short2flt(a->Optical_Depth_055<wbr>_grid1km)<br></span>printVarSummary (data)<br>printMinMax (data, 0)<br>print("-----")<br> nFill = num(ismissing(data))<br>print("nFill="+nFill)<br>print("-----")<span><br><br> lat2d = a->GridLat_grid1km<br> lon2d = a->GridLon_grid1km<br></span>printMinMax (lat2d, 0)<br>printMinMax (lon2d, 0)<br>print("-----")<br><br>;---grid points of interest<br>;---<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/getind_latlon2d.shtml">http://www.ncl.ucar.edu/Do<wbr>cument/Functions/Contributed/g<wbr>etind_latlon2d.shtml</a><br> ; latv = (/22, 24.0, 24.0, 22.0, 22.0/)<span><br> ; lonv = (/76, 76.0, 78.0, 78.0, 76.0/)<br><br></span> <span style="color:rgb(255,0,0)"><b>latv = (/ 22, 24/)<br> lonv = (/ 76, 78 /)<br> nm = getind_latlon2d (lat2d, lon2d, latv, lonv)<span><br> ilt1 = nm(0,0) ; start lat index<br> ilt2 = nm(1,0) ; start lon index<br> iln1 = nm(0,1) ; end lat index<br> iln2 = nm(1,1) ; end lon index<br> var_subset = var(ilt1:ilt2,iln1:iln2)<br> var_subset_avg = avg(var_subset)</span></b></span><br><br> print(nm)<br> print("-----")<span><br><br> do k=0,dimsizes(latv)-1<br> n = nm(k,0)<br> m = nm(k,1)<br> print(lat2d(n,m)+" "+lon2d(n,m)) <br></span><span> print(" "+data(:,n,m)) <wbr> ; nearest grid point at all time steps<br> print("-----")<br> end do<br><br><br><br><br><br></span></div><b>But still not getting the desired output as a single square box mean value . It is showing the errors as</b><br><br>kunal@kunal-Vostro-1015:/media<wbr>/Local Disk_/IIT_DELHI/2001/MAIACTAOT<wbr>_TIME$ ncl extract1.ncl <br> Copyright (C) 1995-2014 - All Rights Reserved<br> University Corporation for Atmospheric Research<br> NCAR Command Language Version 6.2.1<br> The use of this software is governed by a License Agreement.<br> See <a href="http://www.ncl.ucar.edu/">http://www.ncl.ucar.edu/</a> for more details.<br><br><br>Variable: fnames<br>Type: string<br>Total Size: 4 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates: <br>(0) ./<a href="http://maiactaot.h00v02.20010040635-20011550640.nc/">maiactaot.h00v02.20010040635<wbr>-20011550640.nc</a><br><br>Variable: data<br>Type: float<br>Total Size: 40320000 bytes<br> 10080000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 7] x [YDim_grid1km | 1200] x [XDim_grid1km | 1200]<br>Coordinates: <br> time: [8862.583333333314..12486.6666<wbr>6666669]<br>Number Of Attributes: 8<br> coordinates : GridLat_grid1km, GridLon_grid1km<br> hdfeos_name : Optical_Depth_055<br> projection : Albers Conical Equal_Area<br> unit : None<br> long_name : AOT at 0.55 micron<br> _FillValue_original : -28672<br> _FillValue : -28672<br> valid_range : ( -0.1, 5 )<br>(0) AOT at 0.55 micron: min=0 max=3.924<br>(0) -----<br>(0) nFill=6285271<br>(0) -----<br>(0) latitude: min=21.86737357453027 max=34.99855492029426<br>(0) longitude: min=64.23634233372005 max=79.4747290253718<br>(0) -----<br><span style="color:rgb(255,0,0)">fatal:["Execute.c":7863]:Numbe<wbr>r of subscripts on right-hand-side do not match<br> number of dimensions of variable: (2), Subscripts used: (3)<br><br>fatal:["Execute.c":8578]:Execu<wbr>te: Error occurred at or near line 34 in file extract1.ncl<br><br>fatal:Variable (var_subset) is undefined<br>fatal:["Execute.c":8578]:Execu<wbr>te: Error occurred at or near line 35 in file extract1.ncl</span><br><br><br><br><br>Variable: nm<br>Type: integer<br>Total Size: 16 bytes<br> 4 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [2] x [2]<br>Coordinates: <br>Number Of Attributes: 1<br> long_name : indices closest to specified LAT/LON coordinate pairs<br>(0,0) 1199<br>(0,1) 792<br>(1,0) 1199<br>(1,1) 1041<br>(0) -----<span><br>(0) 23.76603307545667 75.53995171540016<br>(0) 0.179<br>(1) 0.142<br>(2) -28672<br>(3) 0.169<br>(4) -28672<br>(5) -28672<br>(6) -28672<br>(0) -----<br></span>(0) 24.26559325846936 77.93996773594735<br>(0) -28672<br>(1) 0.393<br>(2) -28672<br>(3) 0.156<br>(4) 0.313<span><br>(5) -28672<br>(6) -28672<br>(0) -----<br><br></span><div><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961HOEnZb"><font color="#888888"><br><br><br><br><br><br> <br></font></span><div class="gmail_extra"><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961HOEnZb"><font color="#888888"><br clear="all"><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><div style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></div></div></div></div></div></div></div></font></span><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961h5">
<br><div class="gmail_quote">On Sat, Oct 14, 2017 at 12:24 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>></span> wrote:<br><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 style="font-size:small">Kunal,</div><div style="font-size:small"><br></div><div style="font-size:small">Please include the code where you are actually taking the average of your data using the information returned from getind_latlon2d.</div><div style="font-size:small"><br></div><div style="font-size:small">It's important to note that getind_latlon2d simply returns the individual index values that are closest to each of the lat/lon pairs that you gave it. It doesn't work in such a way that you give it a lat/lon box and it returns all the values in that box. </div><div style="font-size:small"><br></div><div style="font-size:small">That's what I was trying to illustrate with the wrf_gsn_10.ncl I created. I gave the function two corners of the box, and then using the two pairs of index values I got back, I was able to get all the values inside the box area of interest:</div><div style="font-size:small"><br></div><div style="font-size:small"><div><font face="monospace, monospace"> lats = (/ min_lat, max_lat /)</font></div><div><font face="monospace, monospace"> lons = (/ min_lon, max_lon /)<br></font><div><font face="monospace, monospace"> nm = getind_latlon2d (lat2d, lon2d, lats, lons)</font></div><div><font face="monospace, monospace"> ilt1 = nm(0,0) ; start lat index</font></div><div><font face="monospace, monospace"> ilt2 = nm(1,0) ; start lon index</font></div><div><font face="monospace, monospace"> iln1 = nm(0,1) ; end lat index</font></div><div><font face="monospace, monospace"> iln2 = nm(1,1) ; end lon index<br> var_subset = var(ilt1:ilt2,iln1:iln2)</font></div><div><font face="monospace, monospace"> var_subset_avg = avg(var_subset)</font></div><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-HOEnZb"><font color="#888888"><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">--Mary</font></div><div><br></div></font></span></div></div></div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-HOEnZb"><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 12, 2017 at 1:19 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com">kunal.bali9@gmail.com</a>></span> wrote:<br><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">Dear Haley Ma'am,<div><br></div><div>Thanks for providing the information and <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0bbd9b95-cb2b-4f0a-b8b2-4edcfc107596" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">example. </span></div><div class="gmail_extra"><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail_signature"><div dir="ltr"><div dir="ltr"><div>But if I follow the getind_latlon2d function as I did, <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-1c477331-4a11-4c84-ad39-3a8ba33ef414" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">then</span> it seems like I followed the same procedure as given in the example <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/getind_latlon2d.shtml">http://www.ncl.ucar.edu/Docume<wbr>nt/Functions/Contributed/getin<wbr>d_latlon2d.shtml</a></div><div><br></div><div>But still I am not getting the mean of the grid box on every time step or am I reading incorrectly the result ?</div><div><br></div><div><br></div><div>regards</div><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248HOEnZb"><font color="#888888"><div>Kunal Bali<br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></div></div></font></span></div></div></div></div><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248h5">
<br><div class="gmail_quote">On Fri, Oct 13, 2017 at 12:11 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>></span> wrote:<br><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 style="font-size:small">Kunal,</div><div style="font-size:small"><br></div><div style="font-size:small">I created an example for you, showing three ways you can subset a WRF grid (since I already had a similar example). It shows how to use the getind_latlon2d function.</div><div style="font-size:small"><br></div><div style="font-size:small">You have to be careful with <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-1289eca5-8b9d-4751-8b41-203824a2d70e" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">subsetting</span> curvilinear data, because specifying a <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-f9850e60-a155-4ea7-8bbb-7b44d4bb976b" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lat</span>/<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-a9f8cbe1-6e90-406d-9d33-36e9eb5603ab" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lon</span> box to take an average over can mean two different things. You will be able to see the difference when you look at the graphic.</div><div style="font-size:small"><br></div><div style="font-size:small">See example wrf_gsn_10<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-400a1526-0159-44ad-9b85-2943a7748b6b" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">.</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-58d9264d-b4fb-44f7-a8d1-5273e9437d65" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">ncl</span> at:</div><div style="font-size:small"><br></div><div><a href="http://www.ncl.ucar.edu/Applications/wrfgsn.shtml#ex10">http://www.ncl.ucar.edu/Applic<wbr>ations/wrfgsn.shtml#ex10</a><br></div><div><br></div><div>--Mary</div><div><br></div><div style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-h5">On Wed, Oct 11, 2017 at 12:03 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com">kunal.bali9@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-h5"><div dir="ltr"><div><div><div><div>Dear NCL users,<br></div><br>Single point coordinates values of 23.0 N and 77.0 E is extracted or print by using the script below<br><br>------------------------------<wbr>------------------------------<wbr>----------<br>; Main code<br>;-----------------------------<wbr>------------------------------<wbr>-----------<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0eab2372-bf07-4566-8be4-194c036a7be5" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">dir</span> = "<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289eb66f378-df2f-4d2a-9b36-192d3bd52243" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289GINGER_SOFTWARE_mark gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">.</span>/"<br> fnames = systemfunc("ls " + dir + "maiactaot*.nc")<br><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-01b044f6-38bf-485e-acdf-ba1d22d492a4" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">print</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-adb5574c-48f3-4308-b709-4f2feea381d7" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0f74ff8f-0ce4-420a-81e6-bff096448b14" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">fnames</span>)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-29a33cf6-7d8b-42d2-a7a2-cc759938745b" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">a</span> = <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-41dd69d4-630b-4b7a-99f2-beefd157e052" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">addfile</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-ad398bb2-cf92-44ae-ab19-1f5afe3a7dd7" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-87206339-7931-4ae3-b2a1-079285d695c3" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">fnames</span>,"r")<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-bcae822f-609a-47c2-8824-3665cfd41e71" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">data</span> = short2flt<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-ac3b2eee-102a-4ed1-b9df-b3257ab6a9a6" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>a->Optical_Depth_055<wbr>_grid1km)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-219c3815-4e7b-4e5f-a3e6-d6a4a9b6a08f" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">nFill</span> = <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-bfaa5817-9b42-44ca-88ab-eff43b22b436" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">num</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-3659f5e6-ebcc-4fac-8c31-0bc9c70c1960" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-cf814eb0-0ab8-457e-86c8-ea3605715cb8" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">ismissing</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289efa1c965-9013-41b9-9c42-5be1289f675a" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289GINGER_SOFTWARE_mark gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>data))<br> lat2d = a->GridLat_grid1km<br> lon2d = a->GridLon_grid1km<br><br><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-04f30fc0-40ae-4d18-8ba4-5e77b3bf2f04" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">;</span>---grid points of interest<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-f3b93122-e518-4ae2-9fb9-0fd4b214dc3e" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">latv</span> = (/23.0/)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-115643a3-d152-4114-ba20-f1a60b04afd3" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lonv</span> = (/77.0 /)<br> <br> nm = getind_latlon2d (lat2d,lon2d, latv, lonv)<br> <br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-099d53de-f551-4764-9543-8221a12abf0c" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">do</span> k=0<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-a69fa0cf-7759-4145-bdcf-d978798a71e9" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-3ba935a2-1b90-46bd-b638-0a071abbfc6b" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">dimsizes</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-76de012b-7e94-4223-96c8-d5547fdb1cc1" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-9a02ff64-99e8-4db4-aa9b-7a77268ca367" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">latv</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-a4c57e6b-18dc-4a07-a369-56bf54f86b40" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">)</span>-1<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0f8c0fe5-1980-431c-820a-15f372f351c0" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">n</span> = <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-7bb9a579-6d60-4292-8ff9-9f4e310f9d9c" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">nm</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-4d516442-f36a-4c47-a61c-94b85425d47e" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>k<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-1e0c19a9-d4ca-4a78-8688-9aa3a0e68224" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>0)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-446a2187-9e79-4bde-850c-b105ae793c20" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">m</span> = <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289ea6e60ca-f646-4f45-a9b9-9f9264aa5877" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289GINGER_SOFTWARE_mark gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">nm</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-074d96f4-e799-43f9-b539-3c1a4fdfac89" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>k<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-8363988c-9f57-4aaa-bf35-7d0be783383e" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>1)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-35df3b6d-3bb3-446b-a73f-f7189cccf843" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">print</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-9a0da512-5789-45e2-93b3-ed87c6139deb" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>lat2d<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-5dfb6b38-c317-46b6-a8d2-78e14ee6a2c1" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>n<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-1137e035-7e87-4678-888a-8d0a89811411" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>m<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-5058d776-df51-4303-bfc8-b4191bad89eb" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">)</span>+" "+lon2d<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-d2d5e97e-676f-4384-abe0-898083cc7a94" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>n<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-276f3f05-2f4d-44b3-9e93-68e031c9fd57" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>m))<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-4519e924-0bd5-4ac0-977e-76c89a450d87" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark"> ;</span> grid point location<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-f7f19858-5fb9-4ca5-8e98-4431cf8a937a" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">print</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-dfb984a4-43ea-4e52-bcaf-e1123af92a28" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>" "+data<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-61e5e32a-6ac6-482b-a3fe-3aaeb0328a9e" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>:<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-acebda5f-cde4-4701-84c7-9c54bdc02f2f" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>n<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-332c9976-ed1b-45b3-a10c-447af1bff740" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">,</span>m))<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-7fedb156-b839-4bc1-b32e-1ca36913b12d" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark"> <wbr> ;</span> nearest grid point at all time steps<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0ceadbb8-787e-4a29-a45d-ad7c4143403f" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">print</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-88590e2e-6ab2-4f48-a268-30abf27307c0" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">(</span>"-----")<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0cca89c0-44cc-4b16-85cf-8350baa58bfe" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">end</span> do<br><br></div></div> <br>NOW, I want to extract the square grid box of 100 km around the coordinates (23N,77E) such as <br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-56a37e98-4e37-4030-99df-0773ca663e8d" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">latv</span> = (/22, 24.0, 24.0, 22.0, 22.0/)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-40988c35-c7dd-41f2-8a8e-342f5982734a" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lonv</span> = (/76, 76.0, 78.0, 78.0, 76.0 /)<br><br><br><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0ff724fa-e29b-4eab-bbe1-b2de3975a31f" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">after</span> that I want to do the area mean of the grid box. So that I can get one single point value of that grid box.<br></div>So when I <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-a1b9144a-5584-426f-b36b-f421fd1a6205" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">replacing</span> <br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-a8e90c83-2c3b-4487-9623-659592d903b2" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">latv</span> = (/23.0/)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-81e03fdc-a361-4f48-98ce-b7b19919e539" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lonv</span> = (/77.0/)<div><div><div><div><br></div><div><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-523ff320-986b-48b4-9de2-d2b31c2fc288" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">with</span><br><br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-b72de5f5-11c2-4b2a-9b07-c6b5f0a43729" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">latv</span> = (/22, 24.0, 24.0, 22.0, 22.0/)<br> <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-23cf30b5-cd6b-4845-a2e5-656528aeaf35" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">lonv</span> = (/76, 76.0, 78.0, 78.0, 76.0 /)<br><br></div><div><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-533657e8-e50c-43eb-94bd-79f4eeacfce0" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">then</span> It's not giving <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-4560b2ea-4a62-4752-a123-7fad0b6e299f" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">mean</span> of the square grid box. <br><br></div><div>It shows the something like that (given below), which is incorrect. So please let me know how can I get the square grid box mean and then extract as <span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-caeb93ef-ff68-4e03-b87c-edd6fe4dd3d6" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">asciiwrite</span> or print?<br><br>(0) 23.76603307545667 75.53995171540016<br>(0) 0.179<br>(1) 0.142<br>(2) -28672<br>(3) 0.169<br>(4) -28672<br>(5) -28672<br>(6) -28672<br>(0) -----<br>(0) 24.00322458681621 75.99551564353467<br>(0) 0.171<br>(1) 0.16<br>(2) -28672<br>(3) 0.17<br>(4) -28672<br>(5) -28672<br>(6) -28672<br><br></div><div><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-606478b4-8f99-4ef2-b02b-40b38fce9195" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">etc</span>.<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-4e9d351b-658a-4db2-8526-0d74705d8994" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">.</span><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-1eb96f87-eea9-4539-b3f0-b5bd67743d13" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">etc</span>..<br></div><div><br><br><br></div><div><span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-0241a938-1e80-41a1-a7fc-22ceaffd8f72" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">regards</span><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-m_2067925330722622683HOEnZb"><font color="#888888"><br clear="all"></font></span></div><span class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-m_2067925330722622683HOEnZb"><font color="#888888"><div><div><div><div class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-m_2067925330722622683m_-3126016792024081524gmail-m_-4699656085058546873gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><div style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></div></div></div></div></div></div></div>
</div></div></font></span></div></div></div></div>
<br></div></div>______________________________<wbr>_________________<br>
<span id="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-c5e5f52d-823e-4d15-8982-b06388b22bd1" class="gmail-m_475446440518907222gmail-m_-5786444228229381634m_3806601695242320961m_7618211217714691404gmail-m_-7893384156443960248m_1683168822085668289gmail-GINGER_SOFTWARE_mark">ncl</span>-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">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div></div></div></div>
<br>______________________________<wbr>_________________<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">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div></div>
<br>______________________________<wbr>_________________<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">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div>
</blockquote></div><br></div></div>
______________________________<wbr>_________________<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">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br></div></blockquote></div><br></div></div><span class="gmail-HOEnZb"><font color="#888888"><div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br class="gmail-m_475446440518907222Apple-interchange-newline">Guido Cioni</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="http://guidocioni.altervista">http://guidocioni.altervista</a>.<wbr>org</div>
</div>
<br></font></span></div></div></blockquote></div><br></div></div></div></div></div></div></div></div>