grads和ncl的纬度对应问题

奚子惠

请教各位大神，grads和ncl的纬度是怎么对应的？
我下载的ncep资料发现纬度写的是90~-90
但是grads里面在编写ctl的时候不是-90~90吗？
现在我有个grads现成的程序，但输入的数据我是先用ncl处理出来的gr ，那么ncl在生成grd格式的文件时：
          system("/bin/rm -f " + filenc)                 
            fout                      = addfile(filenc,"c")
          fout->p_ano               = p_ano(:,{0:90},:)
p_ano是写(:,{0:90},:)  还是(:,{90:0},:)呢？



另外，如果ncl需要读入的数据是grd格式的，那么读进去之后需要为数组行列命名，
是写lat=fspan(-90,90,73)  还是lat=fspan90,-90,73)呢~~


grads和ncl的纬度对应问题困扰我很久了~~大神们给菜鸟讲讲~~

天玖Y

我也遇到过这种问题，如果ncep资料用ncl处理过的话，纬度跟原来是反的，所以在写出来时应该p_ano（：，：：-1，：）。  这样纬度就对了~

奚子惠

天玖Y 发表于 2015-4-10 09:07
我也遇到过这种问题，如果ncep资料用ncl处理过的话，纬度跟原来是反的，所以在写出来时应该p_ano（：，：： ...

我师兄在做论文时 他先用grads处理ncep数据的 生成grd 然后用ncl画图 自己在配属性的时候就是写的-90 90，但是原本ncep是90 -90啊，所以我也搞不清楚了~~~

niuda

应该是一样的吧，看经纬度坐标的单位是degrees_north的吧，应该是从南到北的吧

天玖Y

奚子惠 发表于 2015-4-10 11:15
我师兄在做论文时 他先用grads处理ncep数据的 生成grd 然后用ncl画图 自己在配属性的时候就是写的-90 90 ...

具体原因我也不清楚，反正ncep的nc资料用ncl处理成grd资料的时候要这样设置的~

yang69can

Necp再分析资料一般都是从90N到-90N，但是从人的习惯来看喜欢从小到大，而且ncl中有些函数要求单调递增，所以读取的时候可以从小到大，比如：a=necp->lat(::-1),这样就倒过来了！你可以print(lat）来查看你的数据是否处理完成。另外如果要写成二进制格式，建议用fbindirwrite函数，是fortran喜欢用的函数，当然用哪个函数获得二进制文件都无所谓，只要你自己知道存储方式就可以了！~

奚子惠

niuda 发表于 2015-4-10 20:23
应该是一样的吧，看经纬度坐标的单位是degrees_north的吧，应该是从南到北的吧

好的谢谢哈 我研究研究

奚子惠

yang69can 发表于 2015-4-13 07:54
Necp再分析资料一般都是从90N到-90N，但是从人的习惯来看喜欢从小到大，而且ncl中有些函数要求单调递增，所 ...

谢谢