<div dir="ltr"><div><div>The given script is basically doing the same thing, It extracting the grid element of the selected grid box and then averaging of it. <br><br></div>The next step I need to link the avg value with each time. That's what I am looking for. <br><br></div>But Still, I need to confirm that am I doing the right process for getting grid box mean value. <br><br><br><br>;----------------------------------------------------------------------<br>; Main code<br>;----------------------------------------------------------------------<br> dir = "./"<br> fnames = systemfunc("ls " + dir + "maiactaot*.nc")<br> print(fnames)<br> a = addfile(fnames,"r")<br> data = short2flt(a->Optical_Depth_055_grid1km)<br>;printVarSummary (data)<br>printMinMax (data, 0)<br>print("-----")<br> nFill = num(ismissing(data))<br>print("nFill="+nFill)<br>print("-----")<br><br> lat2d = a->GridLat_grid1km<br> lon2d = a->GridLon_grid1km<br>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/Document/Functions/Contributed/getind_latlon2d.shtml</a><br> ; latv = (/22, 24.0, 24.0, 22.0, 22.0/)<br> ; lonv = (/76, 76.0, 78.0, 78.0, 76.0/)<br><br> latv = (/ 13, 25/) ; for extracting grid box data <br> lonv = (/ 72, 85 /)<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> aot_subset = aot_subset@_FillValue<br> aot_subset(:,ilt1:ilt2,iln1:iln2) = data(:,ilt1:ilt2,iln1:iln2)<br> <br> <br> <br> printVarSummary(aot_subset)<br> print(aot_subset)<br> print("ilt1 = " + ilt1)<br> print("ilt2 = " + ilt2)<br> print("iln1 = " + iln1)<br> print("iln2 = " + iln2)<br> <br> <br> aot_subset_avg = avg(aot_subset)<br> print(aot_subset_avg)<br> print("-----")<br><br> do k=0,dimsizes(latv)-1<br> n = nm(k,0)<br> m = nm(k,1)<br> print(latv+" "+lonv) <br> print(" "+aot_subset_avg) ; nearest grid point at all time steps<br> print("-----")<br> end do<br><br><br> <br><div class="gmail_extra">regards<br clear="all"></div><div class="gmail_extra"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left: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 Wed, Oct 25, 2017 at 3:46 PM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@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">Kunal:<div><br></div><div>Why don't you extract the grid elements you want and then use the NCL average function to find the mean of those grid points? You can extract a grid of the size you want or just the points.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Barry</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 25, 2017 at 10:07 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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><div><div>Thanks for the information. I did not include the time variable. So<br><br></div>I changed aot_subset(ilt1:ilt2,iln1:iln2<wbr>) = data(ilt1:ilt2,iln1:iln2) to aot_subset(<span style="color:rgb(255,0,0)"><b>:,</b></span>ilt1:ilt2,iln1:il<wbr>n2) = data(<span style="color:rgb(255,0,0)"><b>:,</b></span>ilt1:ilt2,iln1:iln2)<br><br></div>AND it worked. <br><br><br></div>But from the given script, I am not getting the mean value of the grid box. It producing the values of every grid separately. As given below<br><br>(0) 22.17971308103524 69.17124237598317 <br>(0) -28672<br>(1) 0.208<br>(0) -----<br>(0) 24.55800304046942 79.47472902537179<br>(0) -28672<br>(1) -28672<br>(0) -----<br><br><br></div>I need to extract the data of square box and mean it. so that I can get the values with every time steps. <br><br> latv = (/ 13, 25/) ; for extracting grid box data <br> lonv = (/ 72, 85 /)<br><div><br></div><div>The script is</div><div><br></div><div><br></div><div><span>;-----------------------------<wbr>------------------------------<wbr>-----------<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>printVarSummary (data)<br>printMinMax (data, 0)<br>print("-----")<br> nFill = num(ismissing(data))<br>print("nFill="+nFill)<br>print("-----")<br><br> lat2d = a->GridLat_grid1km<br> lon2d = a->GridLon_grid1km<br>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" target="_blank">http://www.ncl.ucar.edu/Do<wbr>cument/Functions/Contributed/<wbr>getind_latlon2d.shtml</a><br> ; latv = (/22, 24.0, 24.0, 22.0, 22.0/)<br> ; lonv = (/76, 76.0, 78.0, 78.0, 76.0/)<br><br></span> latv = (/ 13, 25/) ; for extracting grid box data <br> lonv = (/ 72, 85 /)<span><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></span> aot_subset = aot_subset@_FillValue<br> aot_subset(:,ilt1:ilt2,iln1:il<wbr>n2) = data(:,ilt1:ilt2,iln1:iln2)<span><br> <br> <br> <br> printVarSummary(aot_subset)<br> print("ilt1 = " + ilt1)<br> print("ilt2 = " + ilt2)<br> print("iln1 = " + iln1)<br> print("iln2 = " + iln2)<br> <br></span> ; aot_avg = avg(data)<br> aot_subset_avg = avg(data)<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></span></div><span class="m_872981164938646875HOEnZb"><font color="#888888"><div> <br><div><br><br><br></div></div></font></span><div class="gmail_extra"><span class="m_872981164938646875HOEnZb"><font color="#888888"><br clear="all"><div><div class="m_872981164938646875m_3760947518708914871gmail_signature" data-smartmail="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-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></font></span><div><div class="m_872981164938646875h5">
<br><div class="gmail_quote">On Wed, Oct 25, 2017 at 12:25 PM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@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">Kunal:<div><br></div><div>The error message says that you have used three subscripts on a variables with only two dimensions.</div><div><br></div><div>You need to fix this error for your script to proceed to the next step.</div><div><br></div><div>Barry</div></div><div class="gmail_extra"><div><div class="m_872981164938646875m_3760947518708914871h5"><br><div class="gmail_quote">On Wed, Oct 25, 2017 at 8:59 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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">It also shows an error if I replaced := with = in 38 line of the given script (attached earlier), even though variable is in 2D array. <br><span><br>fatal:["Execute.c":7950]:Numbe<wbr>r of subscripts on right-hand-side do not match<br> number of dimensions of variable: (2), Subscripts used: (3)<br><br></span><span>fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 38 in file extract1.ncl<br><br><br></span><div class="gmail_extra"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254HOEnZb"><font color="#888888"><br clear="all"><div><div class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><div><br></div><div><br><br></div><div><br></div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></font></span><div><div class="m_872981164938646875m_3760947518708914871m_-1490809836611793254h5">
<br><div class="gmail_quote">On Tue, Oct 24, 2017 at 4:46 PM, Helen Macintyre <span dir="ltr"><<a href="mailto:Helen.Macintyre@phe.gov.uk" target="_blank">Helen.Macintyre@phe.gov.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="blue" vlink="purple" lang="EN-GB">
<div class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">You are trying to delete a variable that has not been assigned.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Using “:=” tries to delete aot_subset when it doesn’t exist yet.
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">More info here:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><a href="http://www.ncl.ucar.edu/Document/Language/varb_assign.shtml" target="_blank">http://www.ncl.ucar.edu/Docume<wbr>nt/Language/varb_assign.shtml</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"" lang="EN-US">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"" lang="EN-US"> ncl-talk [mailto:<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.<wbr>edu</a>]
<b>On Behalf Of </b>Kunal Bali<br>
<b>Sent:</b> 24 October 2017 06:28<br>
<b>To:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] grid_area_avg<u></u><u></u></span></p><div><div class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I am not able to sort out the query regarding
<br>
<br>
fatal:Variable (aot_subset) is undefined<br>
fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 38 in file extract1.ncl<u></u><u></u></p>
</div>
<p class="MsoNormal">I have attached my script. please let me know how to solve it (the file regarding this script is attached in the previous mail).
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-st">sorry for the inconvenience</span><i>.
</i><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">regards<br clear="all">
<u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Kunal Bali<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif""><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Sat, Oct 21, 2017 at 1:53 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">thanks for providing this information. I am positive that the given variable in the data is a 2D array and the error is coming on 38 line as I mentioned earlier.
<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I am trying to debug this problem accordingly as you mentioned.
<u></u><u></u></p>
</div>
<p class="MsoNormal">However, I have attached the data (size =30mb, tried via ftp but facing some error while uploading) and the full script for extracting the grid box mean value.
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thank You for all the support. <u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u><u></u></p>
<div style="border:solid #dddddd 1.0pt;padding:4.0pt 4.0pt 4.0pt 4.0pt;max-height:18px">
<p class="MsoNormal" style="background:whitesmoke"><b><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:#222222"><a href="https://drive.google.com/file/d/0B6MS-JLMHECwaVh3QVVHUU1Wazg/view?usp=drive_web" target="_blank"><span style="border:none windowtext 1.0pt;padding:0cm;text-decoration:none"><img id="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417_x0000_i1025" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png" border="0"> </span><span style="color:#1155cc;border:none windowtext 1.0pt;padding:0cm;text-decoration:none">maiactaot.h00v02.20010040635-<wbr>20011550640.nc.bz2</span></a><u></u><u></u></span></b></p>
</div>
<p class="MsoNormal"><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><br clear="all">
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-hoenzb"><u></u><u></u></span></span></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Fri, Oct 20, 2017 at 9:14 PM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Kunal,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Multiple people have tried to help you, but you continue to post error messages without trying to debug them yourself. If you don't make an effort, people are likely to stop trying to help. You are including code snippets, but this is not
always enough for us to help you debug the problem. You need to include full scripts, output from print/printVarSummary commands, AND all error messages.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">For example, you originally included this code:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> latv = (/ 18, 24/) ; for extracting grid box data </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> lonv = (/ 74, 78 /)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> nm = getind_latlon2d (lat2d, lon2d, latv, lonv)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> ilt1 = nm(0,0) ; start lat index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> ilt2 = nm(1,0) ; start lon index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> iln1 = nm(0,1) ; end lat index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> iln2 = nm(1,1) ; end lon index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> lat_subset := lat2d(ilt1:ilt2,iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> lon_subset := lon2d(ilt1:ilt2,iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset := data</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset(ilt1:ilt2,iln1:iln2<wbr>) := data(ilt1:ilt2, iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> </span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> ;printVarSummary(aot_subset)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_avg = avg(data)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset_avg = avg(aot_subset)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">and then said that you were getting an error about aot_subset being undefined. From the perspective of somebody who is just looking at the code snippet above and doesn't have access to your full script, we are going to wonder about this
line:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset := data</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Since this appears to be the first time aot_subset is set, our next question is going to be "is data defined"? In the code snippet you included, we don't see where "data" is defined, so we have no way to help. This is where a "printVarSummary(data)"
line would have helped.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">In your last email, in which you did include some information about "data", things are still not making any sense. You said you were subsetting data as a 2D array with:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">aot_subset(ilt1:ilt2,iln1:iln2<wbr>) = data(ilt1:ilt2,iln1:iln2)<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">which you said was causing this error:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">fatal:["Execute.c":7950]:Numbe<wbr>r of subscripts on right-hand-side do not match<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> number of dimensions of variable: (2), Subscripts used: (3)<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">This error doesn't seem possible, because it's saying that you are using 3 subscripts, but the line only has 2 subscripts being used in both arrays. Are you positive that this error is coming from this line?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I have to ask you again to please include your <b>full script</b> next time and
<b>all error messages</b> and <b>all printed output</b>, so people can see exactly what line numbers in your script the errors correspond with and what your data looks like. But don't do this until you have shown that you tried to debug the problem yourself.
Use "print" and "printVarSummary" *before* the locations where the error messages are occurring, so you can see what might be causing the error.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">For example, if the error says <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> number of dimensions of variable: (2), Subscripts used: (3)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">then this means you are subscripting a variable as if it were a 2D array, when it is actually a 3D array. Add
<b>printVarSummary</b> statements on ALL subscripted variables that occur in that line, so you can be sure you are doing things correctly. For example, if the error occurs on this line:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset(ilt1:ilt2,iln1:iln2<wbr>) := data(ilt1:ilt2, iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal">then do this:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><br>
<span style="font-family:"Courier New""> printVarSummary(data)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> printVarSummary(aot_subset)</span><u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset(ilt1:ilt2,iln1:iln2<wbr>) := data(ilt1:ilt2, iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">It may also help to print the il* variables:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> printVarSummary(data)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> printVarSummary(aot_subset)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> print("ilt1 = " + ilt1)</span><u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> print("ilt2 = " + ilt2)</span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> print("iln1 = " + iln1)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> print("iln2 = " + iln2)</span><u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> aot_subset(ilt1:ilt2,iln1:iln2<wbr>) := data(ilt1:ilt2, iln1:iln2)</span><u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Moving forward, you *must* show that you are making some effort to debug your problem.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thank you for your consideration.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888">--Mary<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Wed, Oct 18, 2017 at 12:23 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">printVarSummary(aot_subset), shows that aot_subset is undefined, which I am unable to fix this problem yet.<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal">If I used aot_subset(ilt1:ilt2,iln1:iln2<wbr>) <b><span style="color:red">=</span></b> data(ilt1:ilt2,iln1:iln2)
<b><span style="color:#cc0000">instead of </span></b>aot_subset(ilt1:ilt2,iln1:iln2<wbr>)<span style="color:red">
<b>:=</b></span><b><span style="color:blue"> </span></b>data(ilt1:ilt2,iln1:iln2)<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
then the error appeared as <br>
<br>
<b><span style="color:red">fatal:["Execute.c":7950]:Numbe<wbr>r of subscripts on right-hand-side do not match<br>
number of dimensions of variable: (2), Subscripts used: (3)</span></b><u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:black">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" target="_blank">maiactaotE.h00v03.2000055052<wbr>5-20003660700.nc</a><br>
<br>
Variable: data<br>
Type: float<br>
Total Size: <a href="tel:(254)%20592-0000" target="_blank">2545920000</a> 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_s<wbr>ubset is undefined, can not subscript an undefined variable<br>
fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 38 in file test_extract.ncl</b></span><u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:red">38 line is
</span></b><br>
<b><span style="color:red">aot_subset(ilt1:ilt2,iln1:iln2<wbr>) = data(ilt1:ilt2,iln1:iln2)</span></b><u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:black">fatal:Variable (aot_subset) is undefined</span></b><u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal">regards<span style="color:#888888"><br clear="all">
</span><u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali<u></u><u></u></span></p>
</div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Tue, Oct 17, 2017 at 3:05 PM, Guido Cioni <<a href="mailto:guidocioni@gmail.com" target="_blank">guidocioni@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<p class="MsoNormal">Why did you comment this line ?<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal"> ;printVarSummary(aot_subset)<u></u><u></u></p>
</div>
</div>
</div>
</blockquote>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">That would have helped you in debugging the error on your own.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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 <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal"> <b> <span style="color:red">aot_subset</span>(ilt1:ilt2,iln1:il<wbr>n2) := data(ilt1:ilt2, iln1:iln2)</b><u></u><u></u></p>
</div>
</div>
</div>
</blockquote>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">as NCL has no way to know the shape of <span style="font-family:"Menlo","serif"">
aot_subset</span> beforehand. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">The following should be used instead (notice the assignment with "=" in the second line).<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Menlo","serif""> aot_subset := data</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Menlo","serif""> aot_subset = aot_subset@_FillValue<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-"> aot_subset(ilt1:ilt2,iln1:il<wbr>n2)
<b>=</b> data(ilt1:ilt2, iln1:iln2)</span></span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Menlo","serif"">aot_subset_avg=avg(data(ilt1:i<wbr>lt2, iln1:iln2))</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">If you have the right coordinate attributes you can even jump the first 7 lines and do <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Menlo","serif"">aot_subset_avg=avg(data({18:24<wbr>},{74:78})</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">but I don't know your data. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">As it was already suggested by Mary and others, PLEASE take time to debug the errors on your own. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Regards,<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal">On 17. Oct 2017, at 11:25, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">I did the following changes but still getting an error of <b>
<span style="color:red">aot_subset</span></b> which is undefined. <u></u><u></u></p>
</div>
<p class="MsoNormal">please let me know that how to fix or defined the <b><span style="color:red">aot_subset
</span></b>this problem.<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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:red">aot_subset</span>(ilt1:ilt2,iln1:iln2<wbr>) := data(ilt1:ilt2, iln1:iln2)</b><br>
<br>
;printVarSummary(aot_subset)<br>
aot_avg = avg(data)<br>
aot_subset_avg = avg(<b><span style="color:red">aot_subset</span></b>)<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>
<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal">regards<br clear="all">
<u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Kunal Bali<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif""><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Tue, Oct 17, 2017 at 12:41 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<div>
<div>
<p class="MsoNormal">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<u></u><u></u></p>
</div>
<p class="MsoNormal">ncl_script = extract_iiserb.ncl<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I am not able to define the var_subset in the given script.
<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">please let me know if you have received the data and script?<span style="color:#888888"><br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-hoenzb"><u></u><u></u></span></span></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><br clear="all">
</span><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-hoenzb"><u></u><u></u></span></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p>
</div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Oct 16, 2017 at 11:07 PM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Kunal,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/ftp_files.shtml" target="_blank">http://www.ncl.ucar.edu/ftp_fi<wbr>les.shtml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Document/Language/error_messages.shtml" target="_blank">http://www.ncl.ucar.edu/Docume<wbr>nt/Language/error_messages.sht<wbr>ml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">The particular error you're seeing can be found here:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Document/Language/error_messages.shtml#NumSubRHS" target="_blank">http://www.ncl.ucar.edu/Docume<wbr>nt/Language/error_messages.sht<wbr>ml#NumSubRHS</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">As a side note, you are running a 3-year old version of NCL. I suggest that you upgrade to NCL V6.4.0.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888">--Mary<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Oct 16, 2017 at 12:52 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">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>-----------<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>
printVarSummary (data)<br>
printMinMax (data, 0)<br>
print("-----")<br>
nFill = num(ismissing(data))<br>
print("nFill="+nFill)<br>
print("-----")<br>
<br>
lat2d = a->GridLat_grid1km<br>
lon2d = a->GridLon_grid1km<br>
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" target="_blank">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/)<br>
; lonv = (/76, 76.0, 78.0, 78.0, 76.0/)<br>
<br>
<b><span style="color:red">latv = (/ 22, 24/)<br>
lonv = (/ 76, 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>
var_subset = var(ilt1:ilt2,iln1:iln2)<br>
var_subset_avg = avg(var_subset)</span></b><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>
<br>
<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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/" target="_blank">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/" target="_blank">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:red">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) -----<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.26559325846936 77.93996773594735<br>
(0) -28672<br>
(1) 0.393<br>
(2) -28672<br>
(3) 0.156<br>
(4) 0.313<br>
(5) -28672<br>
(6) -28672<br>
(0) -----<u></u><u></u></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><br>
<br>
<br>
<br>
<br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961hoenzb">
</span></span><u></u><u></u></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><br clear="all">
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961hoenzb"><u></u><u></u></span></span></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p>
</div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Sat, Oct 14, 2017 at 12:24 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Kunal,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Please include the code where you are actually taking the average of your data using the information returned from getind_latlon2d.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> lats = (/ min_lat, max_lat /)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> lons = (/ min_lon, max_lon /)</span><u></u><u></u></p>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> nm = getind_latlon2d (lat2d, lon2d, lats, lons)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> ilt1 = nm(0,0) ; start lat index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> ilt2 = nm(1,0) ; start lon index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> iln1 = nm(0,1) ; end lat index</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> iln2 = nm(1,1) ; end lon index<br>
var_subset = var(ilt1:ilt2,iln1:iln2)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New""> var_subset_avg = avg(var_subset)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New";color:#888888">--Mary</span><span style="color:#888888"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Thu, Oct 12, 2017 at 1:19 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
<div>
<p class="MsoNormal">Dear Haley Ma'am,<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks for providing the information and <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
example. </span><u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">But if I follow the getind_latlon2d function as I did, <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">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" target="_blank">
http://www.ncl.ucar.edu/Docume<wbr>nt/Functions/Contributed/getin<wbr>d_latlon2d.shtml</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">But still I am not getting the mean of the grid box on every time step or am I reading incorrectly the result ?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">regards<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Fri, Oct 13, 2017 at 12:11 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<u></u><u></u></p>
<div>
<div>
<p class="MsoNormal">Kunal,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">You have to be careful with <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
subsetting</span> curvilinear data, because specifying a <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
lat</span>/<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">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.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">See example wrf_gsn_10<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">.ncl</span>
at:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/wrfgsn.shtml#ex10" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/wrfgsn.shtml#ex10</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">--Mary<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<p class="MsoNormal">On Wed, Oct 11, 2017 at 12:03 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p>
</div>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Dear NCL users,<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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 class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
dir</span> = "<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">.</span>/"<br>
fnames = systemfunc("ls " + dir + "maiactaot*.nc")<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">print(fnames</span>)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
a</span> = <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
addfile(fnames</span>,"r")<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
data</span> = short2flt<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">(</span>a->Optical_Depth_055<wbr>_grid1km)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
nFill</span> = <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
num(ismissing(</span>data))<br>
lat2d = a->GridLat_grid1km<br>
lon2d = a->GridLon_grid1km<br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">;</span>---grid
points of interest<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
latv</span> = (/23.0/)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
lonv</span> = (/77.0 /)<br>
<br>
nm = getind_latlon2d (lat2d,lon2d, latv, lonv)<br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
do</span> k=0<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,dimsizes(latv)</span>-1<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
n</span> = <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
nm(</span>k<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>0)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
m</span> = <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
nm(</span>k<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>1)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
print(</span>lat2d<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">(</span>n<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>m<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">)</span>+"
"+lon2d<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">(</span>n<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>m))<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
;</span> grid point location<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
print(</span>" "+data<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">(</span>:<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>n<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">,</span>m))<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566"> <wbr>
;</span> nearest grid point at all time steps<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
print(</span>"-----")<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
end</span> do<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"> <br>
NOW, I want to extract the square grid box of 100 km around the coordinates (23N,77E) such as
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
latv</span> = (/22, 24.0, 24.0, 22.0, 22.0/)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
lonv</span> = (/76, 76.0, 78.0, 78.0, 76.0 /)<br>
<br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">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.<u></u><u></u></p>
</div>
<p class="MsoNormal">So when I <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
replacing</span> <br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
latv</span> = (/23.0/)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
lonv</span> = (/77.0/)<u></u><u></u></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">with</span><br>
<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
latv</span> = (/22, 24.0, 24.0, 22.0, 22.0/)<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
lonv</span> = (/76, 76.0, 78.0, 78.0, 76.0 /)<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">then</span>
It's not giving <span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
mean</span> of the square grid box. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">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 class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">
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<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">etc</span>.<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">.etc</span>..<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<br>
<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">regards</span><span style="color:#888888"><br clear="all">
</span><u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#888888">Kunal Bali<u></u><u></u></span></p>
</div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#888888"><u></u> <u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">______________________________<wbr>_________________<br>
<span class="m_872981164938646875m_3760947518708914871m_-1490809836611793254m_8536465301339776101m_-2618359072133921417gmail-m1857073434174617002m-3452838036563629326m-452582749034241339m3908073124649922015m7028704873234973401gmail-m475446440518907222gmail-m-5786444228229381634m3806601695242320961m7618211217714691404gmail-m-7893384156443960248m168316882208566">ncl</span>-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" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal">______________________________<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" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><u></u><u></u></p>
</div>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black"><br>
Guido Cioni<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black"><a href="http://guidocioni.altervista" target="_blank">http://guidocioni.altervista</a>.o<wbr>rg<u></u><u></u></span></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div></div></div>
<br clear="all">
******************************<wbr>******************************<wbr>**************<br>
The information contained in the EMail and any attachments is confidential and intended solely and for the attention and use of the named addressee(s). It may not be disclosed to any other person without the express authority of Public Health England, or the intended recipient, or both. If you are not the intended recipient, you must not disclose, copy, distribute or retain this message or any part of it. This footnote also confirms that this EMail has been swept for computer viruses by Symantec.Cloud, but please re-sweep any attachments before opening or saving. <a href="http://www.gov.uk/PHE" target="_blank">http://www.gov.uk/PHE</a><br>
******************************<wbr>******************************<wbr>**************<br>
</div>
</blockquote></div><br></div></div></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><br clear="all"><div><br></div></div></div>-- <br><div class="m_872981164938646875m_3760947518708914871m_-1490809836611793254gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:091443%2023108" value="+919144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_872981164938646875gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:091443%2023108" value="+919144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</div></div></blockquote></div><br></div></div>