load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

a = addfile("WRF/wrfhrly_d03_2010-03-01_00:00:00.nc","r")

lat2d  = a->XLAT(0,:,:)
lon2d = a->XLONG(0,:,:)

DATADir = "/scratch/y98/mmughal/BOM_nc_various-2014/"
fils = systemfunc (" ls -1 " + DATADir + "*.nc ")
numfils = dimsizes(fils)
do j = 0,numfils-1,1

f = addfile(fils(j) + ".nc","r")

lat1=f@station_lat
lon1=f@station_lon

lat=dimsizes(lat1)
lon=dimsizes(lon1)
;do p=0,dimsizes(lat1)-1
  


nm = getind_latlon2d (lat2d,lon2d, lat, lon)

  print(nm)

  do k=0,dimsizes(lat)-1
     n = nm(k,0)
     m = nm(k,1)
     print(lat2d(n,m)+"   "+lon2d(n,m))
   end do
 end do
end