ITKeyword,专注技术干货聚合推荐

注册 | 登录

使用grads中的gr2stn做格点插站点

dawran 分享于 2011-08-19

2020腾讯云8月秒杀活动,优惠非常大!(领取2860元代金券),
地址https://cloud.tencent.com/act/cps/redirect?redirect=1040

2020阿里云最低价产品入口,含代金券(新老用户有优惠),
地址https://www.aliyun.com/minisite/goods

这几天尝试用grads中的gr2stn对格点资料进行插值。个人觉得完全是黑匣子操作。根据动力论坛高人的方法尝试了一下。

1.先用Fortran编了一个grads空的可识别的站点资料。grads站点资料格式如下:

struct reportheader { char id[8]; float lat; float lon; float t; int nlev; int flag; };
2.站点和格点的描述文件
dset     E:\data\T639\2006\gtost\st.grd
dtype    station
stnmap   E:\data\T639\2006\gtost\st.map
undef    -9999.90
title rains to station
tdef 1 linear 29Dec2006 1dy
vars 1
gtostr      0   99  station rain2
endvars
Dset E:\data\T639\2006\cutdata\20061229
undef    -9999.90
title 20061229 
xdef   225  linear    27.0     0.5625
ydef   129  linear    0.0    0.5625
zdef     1  levels 999
tdef 1  linear 29Dec2006 1dy
vars 1
LSP   0 99  rain
endvars
 
3.使用函数gr2stn插值
'reinit'
'open E:\data\T639\2006\gtost\g.ctl'
'open E:\data\T639\2006\gtost\st.ctl'
'set gxout print'
'set t 1'
 
'd gr2stn(LSP.1,gtostr.2)'
rc=write('E:\data\T639\2006\gtost\out.txt',result)
4.得到的结果
Printing Stations -- 160 Reports -- Undef = -9999.9
54836    118.2     36.18     -9.99e+33 
9.8324e-01
54904    115.5     35.58     -9.99e+33 
2.5427e+00
54905    115.9     35.6      -9.99e+33 
1.8573e+00
54906    115.4     35.25     -9.99e+33 
2.5794e+00
54907    116.7     34.98     -9.99e+33 
2.0643e+00
54908    115.1     35.28     -9.99e+33 
3.2358e+00
54909    115.6     35.1      -9.99e+33 
2.3791e+00
54910    116.1     35.8      -9.99e+33 
1.7107e+00
54912    116.5     35.72     -9.99e+33 
1.4634e+00
54913    116.8     35.75     -9.99e+33 
1.5159e+00
54914    116.1     35.42     -9.99e+33 
1.8098e+00

这几天尝试用grads中的gr2stn对格点资料进行插值。个人觉得完全是黑匣子操作。根据动力论坛高人的方法尝试了一下。 1.先用Fortran编了一个grads空的可识别的站点资料。grads站点资料格式如下:

相关阅读排行


相关内容推荐

最新文章

×

×

请激活账号

为了能正常使用评论、编辑功能及以后陆续为用户提供的其他产品,请激活账号。

您的注册邮箱: 修改

重新发送激活邮件 进入我的邮箱

如果您没有收到激活邮件,请注意检查垃圾箱。