<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>  dir    = "./"<br>  fnames = systemfunc("ls " + dir + "maiactaot*.nc")<br>print(fnames)<br>  a      = addfile(fnames,"r")<br>  data   = short2flt(a->Optical_Depth_<wbr>055_grid1km)<br>  nFill  = num(ismissing(data))<br>  lat2d  = a->GridLat_grid1km<br>  lon2d  = a->GridLon_grid1km<br><br>;---grid points of interest<br>  latv = (/23.0/)<br>  lonv = (/77.0 /)<br> <br>  nm = getind_latlon2d (lat2d,lon2d, latv, lonv)<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))   ; grid point location<br>     print("     "+data(:,n,m))                <wbr>   ; nearest grid point at all time steps<br>     print("-----")<br>  end 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>   latv = (/22, 24.0, 24.0, 22.0, 22.0/)<br>   lonv = (/76, 76.0, 78.0, 78.0, 76.0 /)<br><br><br>after 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 replacing <br>  latv = (/23.0/)<br>  lonv = (/77.0/)<div><div><div><div><br></div><div>with<br><br>   latv = (/22, 24.0, 24.0, 22.0, 22.0/)<br>   lonv = (/76, 76.0, 78.0, 78.0, 76.0 /)<br><br></div><div>then It's not giving mean 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 asciiwrite 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>etc..etc..<br></div><div><br><br><br></div><div>regards<br clear="all"></div><div><div><div><div class="gmail-m_-4699656085058546873gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><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>
</div></div></div></div></div></div>