请问ncl利用PopLatLonV语句转换pop数据格点总报错是怎么回事？

Cyric

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"
;------------------------------------------------------------
begin
f = addfile("picon_nowa.pop.h.0860-01.nc","r")
depth = f->z_t(0:6)
fils = systemfunc("ls *.nc")
printVarSummary(depth)
fs = addfiles(fils,"r")
u1 = fs[:]->UVEL(:,0:6,:,:)
v1 = fs[:]->VVEL(:,0:6,:,:)
time = fs[:]->time
lon = fspan(0.5,359.5,192)
lat = fspan(-80,80,94)
u = new( (/dimsizes(time),dimsizes(depth),dimsizes(lat),dimsizes(lon)/),double)      
;v = new( (/dimsizes(time),dimsizes(depth),dimsizes(lat),dimsizes(lon)/),double)      
uv = PopLatLonV(u1, v1, "gx1v5", "T62", "aave", "da","061230")
u = uv(0,:,:,:,:)

大约就到这里会出错，显示：
ls： 无法访问rot_*:没有那个目录或文件
（0）  create_rotpopFileName:file not recognized/found: rot_gx1v5.nc

求大神指点，感激不尽

Cyric

已解决，需要角度文件，找个文件把名改了就行……

嘟嘟珍珍

我也遇到了和你一样的问题，请问你是在哪儿找的角度文件呢？

teleporter

嘟嘟珍珍 发表于 2017-11-25 23:29
我也遇到了和你一样的问题，请问你是在哪儿找的角度文件呢？

角度文件可以自己創，挑一個有海洋網格角度資訊的nc文件 (例如某個CESM輸出的pop文件)

執行以下NCO指令: (xxx為input角度的nc文件，自行替換成相應名稱)
ncks -v ANGLE,ANGLET xxx.pop.h.0201-01.nc rot_gx1v6.nc

再把角度文件放在同位置就可以了

xigua

嘟嘟珍珍 发表于 2017-11-25 23:29
我也遇到了和你一样的问题，请问你是在哪儿找的角度文件呢？

http://www.ncl.ucar.edu/Applications/popwgtfiles.shtml
https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/ncl/

xigua

嘟嘟珍珍 发表于 2017-11-25 23:29
我也遇到了和你一样的问题，请问你是在哪儿找的角度文件呢？

http://www.ncl.ucar.edu/Applications/popwgtfiles.shtml
https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/ncl/