高斯格点插值均匀网格

wode_q_x

我记得在论坛里好像下载过，不过由于对程序不太懂，所以没再继续，惭愧啊。你的程序是论坛里的还是其他地方的？可否发给我看下

wxxxren

wode_q_x 发表于 2012-3-21 10:23

                               
登录/注册后可看大图

我记得在论坛里好像下载过，不过由于对程序不太懂，所以没再继续，惭愧啊。你的程序是论坛里的还是其他地方 ...

我也是在论坛上找的，我程序也不好，刚入门，可以给你发，你留邮箱吧~

wode_q_x

pangjing582@126.com  谢谢啦，一起学习！

wxxxren

wode_q_x 发表于 2012-3-21 10:55

                               
登录/注册后可看大图

谢谢啦，一起学习！

发过去了。

花楹

您好，能把你手头那份程序传给我下吗？QQ：573165079 传我邮箱里也行 多谢啊

d1268y

您好，能那程序传给我一份吗？d1268y@126.com. 多谢!

vibgyor

ztftom 发表于 2012-3-13 17:09

                               
登录/注册后可看大图

NCL有现成的程序g2fsh，可以参考下面网址：
http://www.ncl.ucar.edu/Document/Functions/Built-in/g2fsh. ...

我提取出来了grib2中的0.371*0.371°Gaussian数据，想转为经纬度中的1°*1°的数据
用g2fsh试了一试，结果出来还是0.371*0.371

do nmonth=1,12
if((nmonth.eq.1).or.(nmonth.eq.3).or.(nmonth.eq.5).or.(nmonth.eq.7).or.(nmonth.eq.8).or.(nmonth.eq.10).or.(nmonth.eq.12)) then
k=123
end if
if((nmonth.eq.4).or.(nmonth.eq.6).or.(nmonth.eq.9).or.(nmonth.eq.11)) then
k=119
end if
if(nmonth.eq.2) then
k=111
end if
if(nmonth.lt.10) then
filename="/mnt/hgfs/data/grib2/dswr/dswsfc.gdas.20030"+nmonth+".grb2"
else
filename="/mnt/hgfs/data/grib2/dswr/dswsfc.gdas.2003"+nmonth+".grb2"
end if
in=addfile(filename,"r")
do n=0,k
do m=0,5
a=in->DSWRF_P8_L1_GGA0_avg(n,m,{-89.761:89.761},:)
b=g2fsh(a,(/180,360/))
fbindirwrite("/mnt/hgfs/data/2003fixeddswr2.dat",b)
end do
end do
end do


程序如上，不知道错在哪里，请高人指教，感谢不已~

tiantianyixin

同求这个程序，版主发给我一份吧。545171101@qq.com
多谢多谢

wode_q_x

楼主解决了没？可不可以请教一下？

wxxxren

wode_q_x 发表于 2012-7-3 14:50

                               
登录/注册后可看大图

楼主解决了没？可不可以请教一下？

放弃这个程序了。