grads站点实况插值

guarenyouji

fortran.f
Program sta2grd
implicit none
character*8 stid
real lat,lon,rain,tim
integer nlev,flag
tim=0.0
nlev=1
flag=1
open(1,file='jiangshui.txt',status='old')
open(2,file='sta.grd',status='replace',form='binary')
10 read(1,*,end=100)stid,lat,lon,rain
write(2)stid,lat,lon,tim,nlev,flag,rain
goto 10
100 continue
close(1)
nlev=0
write(2)stid,lat,lon,tim,nlev,flag
close(2)

end



sta.ctl
dset  G:\sta.grd
dtype station
stnmap G:\sta.map
TITLE grid data
UNDEF -9999.0
tdef  1 linear 06jul2011 1mo
VARS  1
p     0 99 rainfall data
ENDVARS


fortran
parameter(nx=5001,ny=5001)
real lat(ny),lon(nx)
real s(nx,ny)
open(1,file='G:\zhengli\grid.grd',form='binary')
lat(1)=18.0
lon(1)=100.0
do j=1,ny-1
  lat(j+1)=lat(j)+0.01
enddo
do i=1,nx-1
  lon(i+1)=lon(i)+0.01
enddo

do i=1,nx
do j=1,ny
    s(i,j)=1
enddo
enddo
write(1) s

end

grid.ctl
dset G:\grid.grd
UNDEF -9999.0
title sample GRID Data
xdef  5001 linear 100 1
ydef  5001 linear 18  1
zdef  1    linear 500 1
tdef  1 linear 06jul2011 1mo
VARS 1
g    0 99 grid data
ENDVARS

gs
'reinit'
'open G:\grid.ctl'
'open G:\sta.ctl'
'set lon 100 130'
'set lat 25 38'
'set xlint 1'
'set ylint 1'
*'set mpdset jiangs'
'set grads off'
'set grid off'
'draw map'
'set gxout shaded'
'define a=oacres(g,p.2,0.18)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set cmin 30'
'd aa'
'cbarn 1 0 4.5 3.0'
'enable print G:\rain.gmf'
'print'
'disable print'
;

斥鷃

这个问题问得太抽象了点吧~要么贴张图给个gs？

guarenyouji

斥鷃 发表于 2014-4-17 18:08
这个问题问得太抽象了点吧~要么贴张图给个gs？

贴了，大哥

river

差不多啊，差值以后肯定有误差的，不可能一模一样

guarenyouji

river 发表于 2014-4-18 09:03
差不多啊，差值以后肯定有误差的，不可能一模一样

gs出图，最大降水才40mm，可在站点东西方向是有一个100mm落区的

river

guarenyouji 发表于 2014-4-18 10:10
gs出图，最大降水才40mm，可在站点东西方向是有一个100mm落区的

那你改改插值半径试试，先不设置看看grads默认插值的效果怎么样，好的话就不用自己设置了

guarenyouji

river 发表于 2014-4-18 14:50
那你改改插值半径试试，先不设置看看grads默认插值的效果怎么样，好的话就不用自己设置了

设置了呀，大哥，1.5 ，1， 2，3，4 ,0.5都设置过了

river

guarenyouji 发表于 2014-4-18 20:33
设置了呀，大哥，1.5 ，1， 2，3，4 ,0.5都设置过了

这本来就是一个需要耐心的东西，没有耐心就到不到最好的效果。你看看你的数据写进去的对不对，数据没问题，那就是插值的问题了，觉得不好你就换一个插值方法

风浪

不错不错，谢谢！！

daohua201314

谢谢楼主分享