更新:
NCL试试。它它有WRF-related包功能,包括双向的坐标转换例程。
转换坐标对网格指标,遵循这里的文档:wrf_user_ll_to_ij。向后(即准备协调)例程可以在这里找到:wrf_user_ij_to_ll
他们。他们很方便的。
对于你的场景,您可以加载父域文件(“geo_em”geo_em
或“wrfout”数控wrfout
NetCDF文件),然后指定子域坐标开始,你会得到结果。
这里有一些例子从官方文档:
从坐标转换到网格索引:
负载”美元NCARG_ROOT / lib / ncarg / nclscripts / csm / gsn_code。ncl美元“负载”NCARG_ROOT / lib / ncarg nclscripts / wrf / WRFUserARW。ncl”= addfile (“wrfout_d01_2000-01-24_12:00:00.nc”、“r”) loc = wrf_user_ij_to_ll (10。,30岁。,真正的)打印(“经度位置是:“+ loc(0))打印(“纬度位置是:”+ loc (1))
从坐标转换到网格索引:
负载”美元NCARG_ROOT / lib / ncarg / nclscripts / csm / gsn_code。ncl美元“负载”NCARG_ROOT / lib / ncarg nclscripts / wrf / WRFUserARW。ncl”= addfile (“wrfout_d01_2000-01-24_12:00:00.nc”、“r”) loc = wrf_user_ij_to_ll (10。,30岁。,真正的)打印(“经度位置是:“+ loc(0))打印(“纬度位置是:”+ loc (1))
2.转换转换从电网协调指数:
负载”美元NCARG_ROOT / lib / ncarg / nclscripts / csm / gsn_code。ncl美元“负载”NCARG_ROOT / lib / ncarg nclscripts / wrf / WRFUserARW。ncl“= addfile (“wrfout_d01_2000-01-24_12:00:00.nc”、“r”) res = True res@returnInt = False;返回实际值loc = wrf_user_ll_to_ij (-100.0, 40.0, res)打印(“X / Y位置是:”+ loc)