- 积分
- 32
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-5-4
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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)'
|
|