站点数据处理ctl报错

river

hello-lyn 发表于 2014-5-24 09:38
好像还没到那一步诶，不是要先生成站点映射文件才能继续处理的吗？

好吧，是没到那一步。问题是grads已经提示你错误了啊，ctl 里面缺少维数的设置啊。你仔细去看清风的帖子，对比你缺的是哪一项。你要是还找不出来那就回复我，我再告诉你

lqouc

hello-lyn 发表于 2014-5-24 09:41
就是这样的，不知是哪里出了问题

我说的不是这个，是fortran的报错。
你不是说改成real就报错么？
那个地方根据我的经验必须用real，这个你不修改之后肯定都不对。

hello-lyn

lqouc 发表于 2014-5-24 09:55
我说的不是这个，是fortran的报错。
你不是说改成real就报错么？
那个地方根据我的经验必须用real，这 ...

换了CVF6.5运行以后又可以了，也许是之前fortran版本问题，现在就是grads那边报错。这两天重新改了下程序，先处理成各站点月平均，再转为grd。现在是这样的：Program sta2grd
character*8 stid(9)
real rh(9,12),lon(9),lat(9),tim
integer nlev,flag
open(3,file='d:\0512\RHave.txt')
open(4,file='d:\0512\stid.txt')
open(5,file='d:\0512\jw.txt')
open(7,file='d:\0512\RHnew.grd',form='binary')
  read(3,*) ((rh(i,j),j=1,12),i=1,9)
  do i=1,9
read(4,*) stid(i)
  print*,stid(i)
enddo
do i=1,9
read(5,*) lat(i),lon(i)
enddo
do i=1,12
tim=0.0
nlev=1
flag=1
do j=1,9
write(7)stid(j),lat(j),lon(j),tim,nlev,flag,rh(j,i)
!print*,stid(j),lat(j),lon(j),tim,nlev,flag,rh(j,i)
write(8,*)stid(j),lat(j),lon(j),tim,nlev,flag,rh(j,i)
enddo
nlev=0
write(7)stid(j-1),lat(j-1),lon(j-1),tim,nlev,flag
enddo
close(7)
end
程序可以运行，但stid输出的时候只有一个问号，不知道为什么；然后就会出现grads的报错。

lqouc

hello-lyn 发表于 2014-5-24 10:34
换了CVF6.5运行以后又可以了，也许是之前fortran版本问题，现在就是grads那边报错。这两天重新改了下程序 ...

不明白你说的输出问号是什么意思，我觉得还是先用整形读取stid，然后在输出之前转换char就行了。
我之前做其他程序的时候出过类似的问题，你说stid（1）是指stid数组的第一个数，还是stid字符串的第一个字符呢？
还有之前那个报错的截图是stnmap的时候出现的？

hello-lyn

river 发表于 2014-5-24 09:48
好吧，是没到那一步。问题是grads已经提示你错误了啊，ctl 里面缺少维数的设置啊。你仔细去看清风的帖子 ...

一直是按照教程来的，菜鸟年幼无知呀，请大大指点一二

river

hello-lyn 发表于 2014-5-24 10:52
一直是按照教程来的，菜鸟年幼无知呀，请大大指点一二

你在输入ctl 的路径的时候不要用  \  改成  /  试一下

sun_shine_Xia

两位大神都在指导，那我就不来掺合了~~

hello-lyn

lqouc 发表于 2014-5-24 10:45
不明白你说的输出问号是什么意思，我觉得还是先用整形读取stid，然后在输出之前转换char就行了。
我之前 ...

之前的问题都解决了，不知道是不是电脑的问题，新建了个文件把原来的内容贴进去就能读了o(╯□╰)o，有时候用do循环也提示错误，必须用隐循环来做⊙﹏⊙b汗

hello-lyn

sun_shine_Xia 发表于 2014-5-24 15:51
两位大神都在指导，那我就不来掺合了~~

嘻嘻，大大们多多赐教呀

hello-lyn

river 发表于 2014-5-24 11:00
你在输入ctl 的路径的时候不要用  \  改成  /  试一下

这样也不行，但最后新建了文件就可以了，不知道这算什么问题。。囧 多谢指导！