爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4862|回复: 4

[脚本编辑] 求区域平均后如何画line型图

[复制链接]

新浪微博达人勋

发表于 2015-5-4 16:07:26 | 显示全部楼层 |阅读模式

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

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

x
画厄尔尼诺典型年份的1-12月的变化图,做区域平均,之后不会画line型的图,求大神帮忙,横坐标是1-12月,纵坐标是sst
下面是gs文件,还没有写完的,cp1表示CP型厄尔尼诺所有年份1月的平均,jun表示1月的北极海冰的区域平均,以此类推,想做line型的图,
'reinit'
'sdfopen F:/biye/HadISST_ice-187001_201312.nc'
'set lev 0'
'set lon 0 360'
'set lat 65 89.5'
'define cp1=(sst(t=1045)+sst(t=1117)+sst(t=1141)+sst(t=1357)+sst(t=1405)+sst(t=1453)+sst(t=1669))/7'
'define cp2=(sst(t=1046)+sst(t=1118)+sst(t=1142)+sst(t=1358)+sst(t=1406)+sst(t=1454)+sst(t=1670))/7'
'define cp3=(sst(t=1047)+sst(t=1119)+sst(t=1143)+sst(t=1359)+sst(t=1407)+sst(t=1455)+sst(t=1670))/7'
'define cp4=(sst(t=1048)+sst(t=1120)+sst(t=1144)+sst(t=1360)+sst(t=1408)+sst(t=1456)+sst(t=1670))/7'
'define cp5=(sst(t=1049)+sst(t=1121)+sst(t=1145)+sst(t=1361)+sst(t=1409)+sst(t=1457)+sst(t=1670))/7'
'define cp6=(sst(t=1050)+sst(t=1122)+sst(t=1146)+sst(t=1362)+sst(t=1410)+sst(t=1458)+sst(t=1670))/7'
'define cp7=(sst(t=1051)+sst(t=1123)+sst(t=1147)+sst(t=1363)+sst(t=1411)+sst(t=1459)+sst(t=1670))/7'
'define cp8=(sst(t=1052)+sst(t=1124)+sst(t=1148)+sst(t=1364)+sst(t=1412)+sst(t=1460)+sst(t=1670))/7'
'define cp9=(sst(t=1053)+sst(t=1125)+sst(t=1149)+sst(t=1365)+sst(t=1413)+sst(t=1461)+sst(t=1670))/7'
'define cp10=(sst(t=1054)+sst(t=1126)+sst(t=1150)+sst(t=1366)+sst(t=1414)+sst(t=1462)+sst(t=1670))/7'
'define cp11=(sst(t=1055)+sst(t=1127)+sst(t=1151)+sst(t=1367)+sst(t=1415)+sst(t=1463)+sst(t=1670))/7'
'define cp12=(sst(t=1056)+sst(t=1128)+sst(t=1152)+sst(t=1368)+sst(t=1416)+sst(t=1464)+sst(t=1670))/7'
'define jun=aave(cp1,lon=0,lon=360,lat=65,lat=89.5)'
'define feb=aave(cp2,lon=0,lon=360,lat=65,lat=89.5)'
'define mat=aave(cp3,lon=0,lon=360,lat=65,lat=89.5)'
'define api=aave(cp4,lon=0,lon=360,lat=65,lat=89.5)'
'define may=aave(cp5,lon=0,lon=360,lat=65,lat=89.5)'
'define june=aave(cp6,lon=0,lon=360,lat=65,lat=89.5)'
'define jul=aave(cp7,lon=0,lon=360,lat=65,lat=89.5)'
'define aug=aave(cp8,lon=0,lon=360,lat=65,lat=89.5)'
'define sep=aave(cp9,lon=0,lon=360,lat=65,lat=89.5)'
'define oce=aave(cp10,lon=0,lon=360,lat=65,lat=89.5)'
'define nov=aave(cp11,lon=0,lon=360,lat=65,lat=89.5)'
'define dec=aave(cp12,lon=0,lon=360,lat=65,lat=89.5)'


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-4 16:39:33 | 显示全部楼层
把经纬度坐标固定在某一点,任意点都行
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-5-4 17:07:45 | 显示全部楼层
传说中的谁 发表于 2015-5-4 16:39
把经纬度坐标固定在某一点,任意点都行

那d后面写什么呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-4 17:10:33 | 显示全部楼层
ai果汁小鱼 发表于 2015-5-4 17:07
那d后面写什么呢?

你问我我问谁呢。。。。那是你的数据你的变量啊亲
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-5-4 17:14:20 | 显示全部楼层
ai果汁小鱼 发表于 2015-5-4 17:07
那d后面写什么呢?

但是我新定义了呀,我想画的是我新定义的值,有12个
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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