grads把两个水平分辨率不一样的数据插值为1度*1度，然后求空间差，遇到问题，求指导？

lgg00828

grads把两个水平分辨率不一样的数据插值为1度*1度，然后求空间差，
'reinit'
'open F:/678/6.ctl'
'sdfopen f:/678/t2m0.75.nc'
'set lat 18 55'
'set lon 73 135'
'set lev 1'
'set grads off'   
'set t 2'
**关掉图的时间下标**
'set mpdset cnworld'
'd re(t2.1,62,linear,73,1,37,linear,18,1,bs)-re(t2m.2,62,linear,73,1,37,linear,18,1,bs)'，


6.ctl控制文件为（垂直为19层）

xdef  822 linear   50.96148   0.12162163
ydef  387 linear   12.52175   0.12162163
zdef   19 levels  
1000.00000 950.00000 900.00000..........100.00000



t2m0.75.nc控制文件为
……
zdef 1  linear 0 1

垂直方向set lev 1时，可以单独画出re(t2m.2,62,linear,73,1,37,linear,18,1,bs）
set z 1时，可以单独画出re(t2,62,linear,73,1,37,linear,18,1,bs），
垂直方向怎么设置才能做两个之差呢？求指导！

lgg00828

mofangbao 发表于 2013-12-30 11:16

                               
登录/注册后可看大图

t2.1(z=1) 写变量的时候类似这样写就行了

谢谢，我试试

mofangbao

t2.1(z=1) 写变量的时候类似这样写就行了

lgg00828

mofangbao 发表于 2013-12-30 11:16

                               
登录/注册后可看大图

t2.1(z=1) 写变量的时候类似这样写就行了

'd re(t2.1(z=1),62,linear,73,1,37,linear,18,1,bs)-re(t2m.2(z=1),62,linear,73,1,37,linear,18,1,bs)'，,,就可以了，单独设置一下z就可以，把两个文件 的z一块设置就不行！非常感谢，困扰我一天的问题终于解决了！

晚晴山

re是个什么函数，GrADS命令大全里怎么没有？

mofangbao

晚晴山 发表于 2013-12-30 12:45

                               
登录/注册后可看大图

re是个什么函数，GrADS命令大全里怎么没有？

opengrads的扩展函数 去opengrads的页面找

甜狗狗

mofangbao 发表于 2013-12-30 16:06
opengrads的扩展函数 去opengrads的页面找

请问可以给一个opengrads页面的链接吗~

mofangbao

甜狗狗 发表于 2014-10-14 10:44
请问可以给一个opengrads页面的链接吗~

多学会搜索~