关于站点数据多时次处理和插值保存问题

麦兜兜

mark...........

麦兜兜

楼主，为什么我用Fortran输出的变量如果是二维的话，在grads里生成。map的时候grads就崩溃了呢。。

孤峤蟠烟

帮大忙了。

shows

在最近数据处理并作图中，遇到这样的问题，比如降水数据是60年12个月，160个点，用fortran读出站点资料再画图。看到楼主提供的程序只是多年，160个站点，楼主的程序如下，
cccccccccccccc写站点数据
      do i=1,62cccccccccccccccccccccccccccc多时次读取，这里是62年，自己改值
do j=1,160cccccccccccccccccccccccccccc160站，使用其他站的资料，也可以照样改值即可
id(j)=char(j)
tim=0.0
nlev=1
nflag=1
write(9)id(j),lat(j),lon(j),tim,nlev,nflag,x(j,i)ccccccccccccccccccccccccccc单时次为x(i)即可
enddoccccccccccccccccccccccccccccccccccccccccc
tim=0.0
nlev=0
nflag=1
       write(9)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflag
       end doccccccccccccccccccccccccccccccccccccccccccc若不是多时次，则此do循环,enddo则不用了
end
这样只能处理多年资料，但如何处理上述多年且12个月资料呢，还请楼主指点为盼，谢谢了。

mollycar

Sylvia 发表于 2012-9-17 08:37
我之前试过，是不可以的，但是这也许是每个人电脑和grads的不同造成的。之前也遇到过这样的问题，可能在 ...

你好，看到你写的这个多时次站点资料处理，我按照这种方法写了，在grads中运行的时候，还是出现reading error错误，想请教一下这个该如何改呢？

Jackluo

这么好的东西，先收藏了。

细水长流1

楼主你好，请问你没有做过冬季平均呀？冬季平均, 用12、1、2月进行平均， ave(i)=ave(i)+((a12,j)+a(i,j+1）+ a2（i,j+1)/3)/62
这样改对吗？Fortran一直报错，没法生成grd文件

丫丫君

是不是应该写成tim=0.0
nlev=1
nflag=1
write(9)id(j),lat(j),lon(j),tim,nlev,nflag,x(i,j)ccccccccccccccccccccccccccc单时次为x(j)即可
enddoccccccccccccccccccccccccccccccccccccccccc

轻拢烟云

非常感谢，正好需要

lincf

楼主  顶一个！ghhhhhhhh