站点资料插值成为格点资料用grads作图，还可否固定lat，lon画时间变化曲线图，或纬向

311wen

最近把站点资料作图弄清楚，
想问问大家，除了画某范围内的图，还可否固定lat，lon等画时间变化曲线图，或者径向/纬向环流指数

我试了一下，貌似是不行的。也可能是我的oacres(g,k.2)用法还没弄懂，问问大家

river

上脚本看看

311wen

river 发表于 2015-11-23 10:45
上脚本看看

以下是想话3天的k指数平均场，可以运行
'reinit'
'open c:\msdev\projects\k_9\grid.ctl'
'open c:\msdev\projects\k_9\k_4.ctl'
'set lon 70 140'
'set lat 15 55'
'set lev 1000'
'define k1=oacres(g,k.2(t=1))'
'define k2=oacres(g,k.2(t=2))'
'define k3=oacres(g,k.2(t=3))'
'd (k1+k2+k3)/3'
'd (k1+k2+k3)/3'
'printim c:\msdev\projects\k_9\avek_4.png white'
;

但是如果求时间序列的话，假如固定某个经纬度
'reinit'
'open c:\msdev\projects\k_9\grid.ctl'
'open c:\msdev\projects\k_9\k_4.ctl'
'set lon 70'
'set lat 30'
'set lev 1000'
'set t 1 4 '
'define k=oacres(g,k.2)'感觉这里不对，但是又不知道怎么处理才能画出时序图
'd k'

river

311wen 发表于 2015-11-23 21:03
以下是想话3天的k指数平均场，可以运行
'reinit'
'open c:\msdev\projects\k_9\grid.ctl'

不感觉不对，肯定是错的啊。经纬度都固定在一个点上了，还怎么插值啊。倒过来啊，先插值，在固定经纬度画图

311wen

river 发表于 2015-11-23 21:52
不感觉不对，肯定是错的啊。经纬度都固定在一个点上了，还怎么插值啊。倒过来啊，先插值，在固定经纬度画 ...

非常感谢您~
您好换了下顺序，grads可以运行了。但是出来的图是空白
格点插值的是1°的，三个时次的。时间不用管么
program  griddata
integer,parameter::nx=71,ny=41,nt=3
real grid
integer i,j,k
open(11,file='grid.grd',form='binary')
grid=1.0
do k=1,nt
    do i=1,nx
        do j=1,ny
            write(11)grid
        enddo
    enddo
enddo
close(11)
stop
end

----------------------------------------
'reinit'
'open c:\msdev\projects\k_9\grid.ctl'
'open c:\msdev\projects\k_9\k_4.ctl'
'define kk=oacres(g,k.2)''set lon 100'
'set lat 30'
'set lev 1000'
'set t 1 3'
'd kk'
'printim c:\msdev\projects\k_9\tk_4.png white'
;

river

311wen 发表于 2015-11-23 22:06
非常感谢您~
您好换了下顺序，grads可以运行了。但是出来的图是空白
格点插值的是1°的，三个时次的。 ...

出不来图，那你要自己先排除一下啊。别的都不要，你先看看你插值出来的kk能不能出图，然后再一句一句的往下顺······

311wen

因为插值是三个时次生成的.map
所以画1个时次还好（t=1就可以），画连续多次时序的就不知道怎么画了。

现在不知道站点多时次插值出来，能不能画时间序列图。
因为家园里没看见有人做过同类的。

星雨

可以的，不过set lat 和lon的范围不要超过你插值区域的范围了。我之前画过雷达资料的，没问题。

311wen

星雨 发表于 2015-11-24 19:57
可以的，不过set lat 和lon的范围不要超过你插值区域的范围了。我之前画过雷达资料的，没问题。

嗯嗯，谢谢。现在一个站点的也ok了，不知道多个站点的时间平均也可以么。
比如格点可以用d ave(ave(k,lon=0,lon=360),t=1,t=3)，站点就不太清处了

小白菜99

{:5_275:}好高大上的感觉，看看能不能学习一下