grads的tcorr函数可以用于做站点相关吗？

云一林

                               
登录/注册后可看大图

这是我从气象家园找的tcorr的用法，我做的是青藏高原夏季东部大气热源的57年区域平均，即对应图中的elnino，为57年的时间序列，slp对应的是西南地区夏季降水，是站点数据，并不是格点，所以我将57年降水数据插值完再做的相关，以下是我的相关gs，修改了一些地方：
'reinit'
'open d:\lw\q1summer.ctl'
'open d:\lw\3.ctl'
'open d:\lw\grid.ctl'

'enable print d:\lw\tcorr.png'
'set grads off'
'set grid off'
'set x 1 1'
'set y 1 1'
'set z 1'
'set t 1 57'
'define mdiv2=aave(q1,lon=90,lon=100,lat=30,lat=35)'

'set mpdset xinan2'
'set lon 97 110.5'
'set lat 21 34.5'
'set z 1'
'set t 1 57'为了插值57年将t修改
'define qx=oacres(g.3,sd.2)'
'define mm=tcorr(mdiv2,qx,t=1,t=57)'与插值结果做相关
*'set gxout contour'不需要这句话也能画出来相关图为什么？
'd mm'

'q w2xy  102.97  33.57'
x=subwrd(result,3)
y=subwrd(result,6)
'set strsiz 0.07'  
'draw string  'x' 'y' Ruoergai   '
'q w2xy  102.97  33.57'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.03'

'print'
'disable print'

grads运行结果如下：

ga-> D:\lw\tcorr.gs
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits

出来的图片不能双击打开，只能用gv32来看，请问这样用tcorr做出来的相关结果是正确的吗？如果不正确那么应该怎么做高原格点数据的热源和西南地区站点格式的降水的相关？需要相关的显示区域为西南地区。还有一个问题就是为什么grads直接某一年降水数据插值，即set t 1这样，grads运行就正常，后面绿色字体的站名标注就正常显示，而相关时warning能出图但是后面的站点就不能显示了？有没有大神可以给我指导下。如果不能用tcorr做站点相关，那么应该怎么做？

river

就是一个点和一个场的相关，我认为整体的做法是对的。细节上有些小问题。
①'set gxout contour'不需要这句话也能画出来相关图为什么？  因为你不设置出图类型，那它默认就给你出等值线图，这个不存在任何问题。
②出来的图片不能双击打开，只能用gv32来看。之所以这样是因为你用的GrADS是旧版，你使用了enable print 这个命令出来的是矢量图，只能用gv32来查看，无论你用什么后缀都一样。 要想直接双击打开的你可以直接用 printim 命令。然而新版的GrADS已经舍弃这两种命令了，全部统一为gxprint 命令了。
③为什么grads直接某一年降水数据插值，即set t 1这样，grads运行就正常，后面绿色字体的站名标注就正常显示，而相关时warning能出图但是后面的站点就不能显示了。 这个问题因为你对GrADS的运行机制不够清楚，你为了把57年降水场全部插值而设置了set t 1 57 ，这是一个隐循环，你后面求相关的语句是放在它后面的，相当于求相关的计算是在循环内的，也就是它会循环求57次。而你的相关场求出来就是一个场，没有多时次，所以从第二次循环开始，全部提示为数据超出文件限制。最终也就导致GrADS不会执行最后标注站点的命令。   解决这个问题，应该在求相关的语句'define mm=tcorr(mdiv2,qx,t=1,t=57)' 前面加一句set t 1

云一林

river 发表于 2020-6-2 23:57
就是一个点和一个场的相关，我认为整体的做法是对的。细节上有些小问题。
①'set gxout contour'不需要这 ...

非常感谢，终于给我试出来了，加set t 1还是会出现warning，站点出不来，能出图能出站点的是在tcorr前加一句reset，然后把所有经纬度，地图全部重新set一遍，这次把时间设成set t 1，然后就能出图了，站点也能出来了，但是很奇怪的是还是会warning，不过无所谓，结果是对的就行，感谢大神的耐心解答。

云一林

现在把gs附上，大家可以直接做格点数据时间序列和降水站点场的相关，生成的图片也能双击打开，不用gv32了，以下是gs：

'reinit'
'reset'
'open d:\lw\q1summer.ctl'夏季大气热源，格点数据
'open d:\lw\3.ctl' 站点降水数据
'open d:\lw\grid.ctl'站点插值所需的背景场
三个文件的ctl描述时间那行的是一模一样的，就复制粘贴使他们一样就行
'set x 1 1'
'set y 1 1'
'set z 1'
'set t 1 57'
'define mdiv2=aave(q1,lon=90,lon=100,lat=30,lat=35)'夏季大气热源的区域平均，57年，数据是先把夏季提取为dat文件再做区域平均，得到57年的时间序列

'set lon 97 110.5'
'set lat 21 34.5'
'set z 1'
'set t 1 57'
'define qx=oacres(g.3,sd.2)'57年插值
'reset'重点！！重新设置环境变量
'set grads off'
'set grid off'
'set mpdset xinan2'
'set lon 97 110.5'
'set lat 21 34.5'
'set z 1'
'set t 1'
'define mm=tcorr(mdiv2,qx,t=1,t=57)'grads自带的tcorr函数
'd mm'

'q w2xy  102.97  33.57'
x=subwrd(result,3)
y=subwrd(result,6)
'set strsiz 0.07'  
'draw string  'x' 'y' Ruoergai   '
'q w2xy  102.97  33.57'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3  'x' 'y' 0.03'站点标注

'printim d:\lw\tcorr.png x1600 y1200 white'
'disable print'
这个gs在grads运行时仍旧是warning不知道为啥，但是能出相关图，站点也能显示了。

njzqxt

楼主好棒！求q1summer.ctl
3.ctl
grid.ctl

qxtjls

求3个ctl文件！！！

云一林

qxtjls 发表于 2021-7-14 10:09
求3个ctl文件！！！

链接：https://pan.baidu.com/s/1RAjs7IALe8d9YvL8Qb-QHQ
提取码：erxf

已经不用grads一年了，不太清楚grid.ctl找的对不对

云一林

njzqxt 发表于 2021-7-14 09:13
楼主好棒！求q1summer.ctl
3.ctl
grid.ctl

链接：https://pan.baidu.com/s/1RAjs7IALe8d9YvL8Qb-QHQ
提取码：erxf

已经不用grads一年了，不太清楚grid.ctl找的对不对

njzqxt

图出来了，感谢楼主！