函数g2fsh的用法

vibgyor

RT~~
不知道哪位有用过g2fsh这个函数？
我提取出来了grib2中的0.371*0.371°Gaussian数据，想转为经纬度中的1°*1°的数据，ncl中的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


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

hzwjy

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)

问题可能是在fbindirwrite上面了。
因为fbindirwrite这个函数写数据会接着已有的文件写。如果你在运行这个脚本前“2003fixeddswr2.dat”这个文件已存在，画图的又没考虑到这个问题，那么得到的可能就不是你期望的那个数据。

另外前面的时间循环写得比较麻烦，如果电脑性能够，可以用g2fsh函数一次就把一月的数据都做格点转换了。