grads多时次格点插站点结果提取

HAHAHAH

格点数据ctl文件如下：
dset    d:\gr2stn\1981.grd
title   1981 grid data
undef   -9999.0
xdef    141 linear 70 0.5
ydef     81 linear 15 0.5
zdef      1 levels 1000
tdef    365 linear 1jan1981 1dy
vars    9
temp 0 99 average temperature data
tmax 0 99 maxinum temperature data
tmin 0 99 minimum temperature data
rain 0 99 rain data
fsw  0 99 net absorbed solar energy flux(W/m2)
flw  0 99 net infrared energy flux(W/m2)
q    0 99 specific humidity
u10m 0 99 meridianal wind
v10m 0 99 zonal wind
endvars

站点数据ctl文件如下：
dset    d:\gr2stn\sta.grd
dtype   station
stnmap  d:\gr2stn\sta.map
undef   -32766
title   station information
tdef    1 linear 1jan1981 1dy
vars    1
g     0  99  rain data
endvars

gs文件如下
'reinit'
'i=1'
while(i<=365)
'open d:\gr2stn\1981.ctl'
'set lon 70 140'
'set lat 15 55'
'set t 'i''
'define a'i'=temp'
'close1'
i=i+1
endwhile
i=1
while(i<=365)
'open d:\gr2stn\sta.ctl'
'grads2ascii gr2stn(a'i',g) d:\gr2stn\1981_a'i'.txt'
i=i+1
endwhile

HAHAHAH

求大神指点我真的用尽洪荒之力了啊

男紫汗

提取哪个变量就用哪个变量，直接用temp

HAHAHAH

男紫汗 发表于 2017-6-26 14:23
提取哪个变量就用哪个变量，直接用temp

谢谢回复，但我还是不太懂，能否说的细致些，谢谢您

男紫汗

HAHAHAH 发表于 2017-7-2 16:23
谢谢回复，但我还是不太懂，能否说的细致些，谢谢您

'reinit'

'open d:\gr2stn\1981.ctl''open d:\gr2stn\sta.ctl'
'set lon 70 140'
'set lat 15 55'

i=1
while(i<=365)

'grads2ascii gr2stn(temp,g.2) d:\gr2stn\1981_'i'.txt'
i=i+1
endwhile

;
试试看

HAHAHAH

男紫汗 发表于 2017-7-3 08:38
'reinit'

'open d:\gr2stn\1981.ctl''open d:\gr2stn\sta.ctl'

太感谢您了，成功了！可是为什么输出的数据全为负呢?

HAHAHAH

HAHAHAH 发表于 2017-7-3 10:16
太感谢您了，成功了！可是为什么输出的数据全为负呢?

刚看了下，对应站点的temp数据确实为付，太感谢您啦！

HAHAHAH

男紫汗 发表于 2017-7-3 08:38
'reinit'

'open d:\gr2stn\1981.ctl''open d:\gr2stn\sta.ctl'

刚看了下，对应站点的temp数据确实为付，太感谢您啦！

HAHAHAH

男紫汗 发表于 2017-6-26 14:23
提取哪个变量就用哪个变量，直接用temp

您好，我想问一下为什么我这样插出来每天的数据是一样的呢？

男紫汗

HAHAHAH 发表于 2017-7-8 20:40
您好，我想问一下为什么我这样插出来每天的数据是一样的呢？

你检查下你的每天的站点数据一样不一样