利用grads描述文件中的pdef做插值的实用方法，转模式的同学看过来~~

dannybear

东风急流 发表于 2016-3-19 15:33
skHGT表示的是实况高度场，就是fnl资料的，本来就没有pdef哈

请问一下，为什么在wrf场之后提取的hgt的变量，把这个提取的变量的经纬度 这么设置，ydef 46 linear 10 1
xdef 81 linear 60 1，它在原来wrf的文件里面的分辨率不是0.0810吗，
xdef 1080 linear   56.71197   0.08108108
ydef  646 linear    7.71816   0.08108108
这样修改，会不会出错，为什么我也是按照你的方法错了，是我哪里理解错误了吗？

东风急流

dannybear 发表于 2016-3-19 20:12
请问一下，为什么在wrf场之后提取的hgt的变量，把这个提取的变量的经纬度 这么设置，ydef 46 linear 10 1 ...

你没有理解我的意思。。。我写这个的目的就是弄成分辨率一样的去画图，改成xdef 81 linear 60 相当于插值

dannybear

东风急流 发表于 2016-3-20 10:32
你没有理解我的意思。。。我写这个的目的就是弄成分辨率一样的去画图，改成xdef 81 linear 60 相当于插值

那不要在gs文件里面写差值的函数吗？只要把x,y方向的格点数修改一下吗？

liuchunchun

非常有用，十分感谢。

东风急流

dannybear 发表于 2016-3-20 20:38
那不要在gs文件里面写差值的函数吗？只要把x,y方向的格点数修改一下吗？

真的不用，自动差值的，你再仔细看下。

dannybear

东风急流 发表于 2016-3-21 08:51
真的不用，自动差值的，你再仔细看下。

请问一下，楼主，pdef的作用是什么？是不是仅仅是对画图的时候起作用，对于数据排列规则没有作用，数据的排列规则是不是按照正常的格点来排列的

当未来

楼主你好，前面你说的我明白了。如果想用cdo直接把ctl转成nc，去掉pdef的话，xdef和ydef该怎么改呢？
我原来的文件是
pdef  351 252 lcc  22.878  113.323  176.000  126.500  40.00000  20.00000  113.50000   3000.000   3000.000
xdef  809 linear  107.84959   0.01351351
ydef  532 linear   19.23486   0.01351351
如果用列举法的话，xdef是不是要列举出351×252个值。。

dannybear

东风急流 发表于 2016-3-21 08:51
真的不用，自动差值的，你再仔细看下。

请问一下，楼主，原来数据的经纬度的范围和插值之后的经纬度的范围不重合，也没有关系吗？C:\Users\Administrator\Desktop\QQ截图20160321105156.png

东风急流

dannybear 发表于 2016-3-21 09:48
请问一下，楼主，pdef的作用是什么？是不是仅仅是对画图的时候起作用，对于数据排列规则没有作用，数据的 ...

恰好反了，pdef描述的是原始数据的排列方式，而xdef,ydef是描述的画图所用的数据(在画图时候自动插值成)，，grads画图得是格点数据。

东风急流

当未来 发表于 2016-3-21 10:41
楼主你好，前面你说的我明白了。如果想用cdo直接把ctl转成nc，去掉pdef的话，xdef和ydef该怎么改呢？
我原 ...

这个转换我没用过，但我觉得不能去pdef吧，因为这个是描述原始数据。你可以用fwrite重新写一下数据（相当于做插值了），这样再去掉即可。