求双线性插值的fortran程序

fjk

急需要一个双线性插值的fortran程序，可以实现不同分辨率下格点资料的转换。自己正在尝试编，但感觉转换后有些不对。哪位大神可以把代码贴一下让我学习一下吗！不胜感激！

lqouc

论坛有，可以搜搜
话说如果是grads的话貌似可以直接用re函数，其他的软件也都有类似功能

fjk

那个双线性插值的程序真的需要原来各格点的坐标吗？好麻烦呀

fjk

我想用grads中自带的re函数将2.5*2.5的二进制格点资料重新插值成1*1的格点资料，资料的经纬度范围为东经30～160，北纬12.5～80
ctl文件为：dset   f:\height.dat
undef  1.e+36
xdef   53 linear 30   2.5
ydef   28 linear 12.5 2.5
zdef   17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef   480 linear 00:00Z1jan2009 6hr
vars   1
h    17   99     height
endvars
gs文件为：'reinit'
'open f:\1.ctl'
'set grads off'
'set lon -180 180'
'set lat -90 90'
'set z 1'
'set t 1'
'set gxout fwrite'
'set fwrite f:\xin.grd'
'd re(h,131,linear,30.0,1.0,69,linear,12.5,1.0,b1)'
'disable fwrite'
;
运行时报re函数有错。
请问哪位大神可以帮忙解答一下呀？