在使用NCL處理nc文件時冲茸，
如果是wrfout結(jié)果，可以用“wrf_ij_to_ll”、“wrf_ll_to_ij ”或者“wrf_user_xy_to_ll”函數(shù)找到距離特定位置（如觀測點）對應(yīng)的網(wǎng)格文件中最近的網(wǎng)格點毅舆；
如果是普通的netcdf文件，可能用distance和nearest關(guān)鍵字找不到想要的函數(shù)愈腾，ncl在介紹時用了“closest”關(guān)鍵字憋活。我們簡單看一下兩個函數(shù)的差異：
wrf_ll_to_ij ： Finds the nearest model grid indices (i,j) to the specified location(s) in longitude and latitude (deprecated). 翻譯為：查找經(jīng)度和緯度到指定位置最近的模型網(wǎng)格索引（i，j）虱黄。（不建議使用悦即，6.6版本ncl建議使用更新的wrf_user_xy_to_ll等函數(shù)）。
getind_latlon2d ：Returns the indices (subscripts) of two-dimensional latitude/longitude arrays closest to a user-specified latitude/longitude coordinate pair.? 翻譯為：返回最接近用戶指定的緯度/經(jīng)度坐標對的二維緯度/經(jīng)度數(shù)組的索引（下標）。
還有個特殊的重要問題：有的矩形網(wǎng)格的nc文件中辜梳，lat和lon是線性變化的粱甫，所以是一緯數(shù)組。而上述函數(shù)使用的是二維數(shù)組作為被查閱的數(shù)據(jù)作瞄，這時候就要將一維的數(shù)據(jù)轉(zhuǎn)換為二維數(shù)據(jù)茶宵，方法如下：
? ? ? ? nlat_1d = dimsizes(lat_1d)
nlon_1d = dimsizes(lon_1d)
lat2d = conform_dims((/nlat_1d(0),nlon_1d(0)/),lat_1d,0) ;構(gòu)建netcdf文件二維經(jīng)緯度
lon2d = conform_dims((/nlat_1d(0),nlon_1d(0)/),lon_1d,1)
構(gòu)建好了就可以使用上述找最近距離的函數(shù)了。
