请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13122|回复: 8

[图形美化] grads的tcorr函数可以用于做站点相关吗?

[复制链接]

新浪微博达人勋

发表于 2020-6-2 18:22:48 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x

                               
登录/注册后可看大图

这是我从气象家园找的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做站点相关,那么应该怎么做?



相关结果

相关结果
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-6-2 23:57:41 | 显示全部楼层
本帖最后由 river 于 2020-6-3 00:00 编辑

就是一个点和一个场的相关,我认为整体的做法是对的。细节上有些小问题。
'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
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-3 03:31:33 | 显示全部楼层
river 发表于 2020-6-2 23:57
就是一个点和一个场的相关,我认为整体的做法是对的。细节上有些小问题。
①'set gxout contour'不需要这 ...

非常感谢,终于给我试出来了,加set t 1还是会出现warning,站点出不来,能出图能出站点的是在tcorr前加一句reset,然后把所有经纬度,地图全部重新set一遍,这次把时间设成set t 1,然后就能出图了,站点也能出来了,但是很奇怪的是还是会warning,不过无所谓,结果是对的就行,感谢大神的耐心解答。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-3 03:38:04 | 显示全部楼层
本帖最后由 云一林 于 2020-6-3 03:59 编辑

现在把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不知道为啥,但是能出相关图,站点也能显示了。

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-14 09:13:16 | 显示全部楼层
楼主好棒!求q1summer.ctl
3.ctl
grid.ctl
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-14 10:09:48 | 显示全部楼层
求3个ctl文件!!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-14 20:43:44 | 显示全部楼层
qxtjls 发表于 2021-7-14 10:09
求3个ctl文件!!!

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

已经不用grads一年了,不太清楚grid.ctl找的对不对
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-7-14 20:46:23 | 显示全部楼层
njzqxt 发表于 2021-7-14 09:13
楼主好棒!求q1summer.ctl
3.ctl
grid.ctl

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

已经不用grads一年了,不太清楚grid.ctl找的对不对
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-14 23:11:30 | 显示全部楼层
本帖最后由 njzqxt 于 2021-7-16 22:47 编辑

图出来了,感谢楼主!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表