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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3951|回复: 3

[脚本编辑] 关于海温数据逐月变季节性以及各种格式数据转换处理的一些问题

[复制链接]

新浪微博达人勋

发表于 2018-10-17 17:46:28 | 显示全部楼层 |阅读模式

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

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

x
初学grads,然后老师今天布置了一道题目,整个宿舍做了半天都报错,想问问大家怎么处理

贴一下题目:
把海温数据HadISST_sst.nc中海冰的数据设置为缺测值,并把逐月的数据转换成季节平均的数据,并输出至HadISST.seasonal.1870-2016.dat文件中,写出相应的ctl(注意,变量为:sstapr,sstsum,sstaut,sstwin)文件,并转换成HadISST.seasonal.1870-2016.nc文件。

贴一下gs:
'reinit'
'sdfopen D:\grads\data\HadISST_sst.nc'
'set gxout fwrite'
'set fwrite D:\grads\data\HadISST.seasonal.1870-2016.dat'
t0=4
while(t0<=1772)
'set t 't0
'set z 1'
'set x 1 360'
'set y 1 180'
'define sstapr=ave(sst,t-1,t+1)'
'd maskout(sstapr,sstapr)'
t0=t0+12
endwhile
t1=7
while(t1<=1772)
'set t 't1
'set z 1'
'set x 1 360'
'set y 1 180'
'define sstsum=ave(sst,t-1,t+1)'
'd maskout(sstsum,sstsum)'
t1=t1+12
endwhile
t2=10
while(t2<=1772)
'set t 't2
'set z 1'
'set x 1 360'
'set y 1 180'
'define sstaut=ave(sst,t-1,t+1)'
'd maskout(sstaut,sstaut)'
t2=t2+12
endwhile
t3=13
while(t3<=1772)
'set t 't3
'set z 1'
'set x 1 360'
'set y 1 180'
'define sstwin=ave(sst,t-1,t+1)'
'd maskout(sstwin,sstwin)'
t3=t3+12
endwhile
'disable fwrite'

*

贴一下ctl:
dset D:\grads\data\HadISST.seasonal.1870-2016.dat
title Seasonal version of HadISST sea surface temperature component
undef -9.99e+08
dtype netcdf
xdef 360 linear -179.5 1
ydef 180 linear -89.5 1
zdef 1 linear 0 1
tdef 590 linear Mar1870 1yr
vars 4
sstapr=>sstapr 0 t,y,x  mean sst of spring
sstsum=>sstsum 0 t,y,x  mean sst of summer
sstaut=>sstaut 0 t,y,x  mean sst of autumn
sstwin=>sstwin 0 t,y,x  mean sst of winter

endvars

现在dat已经输出来了,但是再转nc文件一直报错,我觉得是ctl有问题吧,但是实在找不到怎么解决,求大佬帮忙看看,谢谢!

微信图片_20181017002727.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-10-19 16:14:18 | 显示全部楼层
大佬们帮帮忙吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-10-20 10:44:39 | 显示全部楼层
我猜 应该是ctl的dtype不能写成netcdf,因为是二进制数据,要写成grid吧
但是我的gs文件跟你一样啊为什么我死循环了啊啊啊啊
(能认出我是谁吗hhh)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-10-20 19:50:07 | 显示全部楼层
哈哈哈认出来了,grads这磨人的小妖精,我改一下试一下
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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