<div dir="ltr"><div dir="ltr">Hi all,</div><div>I wanted to extract rainfall within 500 km from certain lat and lon. I am not getting the right one. any suggestions?</div><div dir="ltr"><div><br></div><div><div>f  = addfile("<a href="http://jan500.nc">jan500.nc</a>", "r")</div><div> prc =f->pre(:,{75:100},{0:30})</div><div> lon = f->longitude</div><div> lat =f->latitude</div><div> prc_reorder = prc(time|:,latitude|:,longitude|:)</div><div><br></div><div> prc_reorder&latitude = lat</div><div> prc_reorder&longitude = lon</div><div><br></div><div>lat1 =(/5.5/)</div><div>lon1 = (/87.0/)</div><div><br></div><div>do i = 0,dimsizes(lat)-1</div><div>      do j = 0,dimsizes(lon)-1</div><div><span style="white-space:pre">        </span>      if gc_latlon(lat1,lon1,i,j,2,4) .lt.500 then</div><div><span style="white-space:pre"> </span>          prc_reorder(:,i,j) =  prc_reorder(:,i,j)<span style="white-space:pre">     </span></div><div>              end if<span style="white-space:pre">           </span> </div><div>      end do</div><div> end do</div><div>;;;;</div><div>Jalal</div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font color="#000000" size="2"><b>Md. Jalal Uddin</b></font></div><div dir="ltr">MSc in Applied Meteorology (English Language)<br>Nanjing University of Information, Science and Technology, China</div><div dir="ltr"><span style="letter-spacing:0.2px">B.Sc. in Disaster Management (Hons.) </span><br></div><div dir="ltr"><span style="letter-spacing:0.2px">Patuakhali Science and Technology University, Bangladesh.</span></div><div dir="ltr">Cell: +8613260859092, <span style="letter-spacing:0.2px">+8801792052662</span></div><div dir="ltr"><font size="2">Web: <a href="http://www.dmjalal90.weebly.com" target="_blank">www.dmjalal90.weebly.com</a>  <br>Facebook: jalal.hossen.39  <br>LinkedIn: <a href="https://bd.linkedin.com/in/md-jalal-uddin-80a026b0" target="_blank">https://bd.linkedin.com/in/md-jalal-uddin-80a026b0</a>   <br>Twitter: dmjalal90  <br>Skype: dmjalal90</font><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>