【已解决】参考清风大侠的《fortran+grads站点作图详解》,生成站点映射文件（*map）时

beiji66

这是为什么啊，真郁闷

lysx

mofangbao 发表于 2013-4-21 19:09
那本书里面本身就有很多实例，看不懂是因为你没自己敲代码，把书上的代码一句句的敲进去，在这个过程中你 ...

fortran关于文件打开关闭中的关键字很多。我在编写程序时，即使使再简单的文件，我也严格按照书上的要求来做，一说明自己对这些关键字的理解程度，二是它可以帮助你养成良好的变成习惯，三是可以帮助你检查到一些因为粗心大意造成的错误。

beiji66

还是不明白啊，看来学习得不够

ni三色堇

留个爪印~~

jojingwj

楼主  你这个问题解决了么 我遇到了同样的问题

shirly

jojingwj 发表于 2014-10-7 20:20
楼主  你这个问题解决了么 我遇到了同样的问题

参考前面别人的回复，尤其清风大侠的回复

jojingwj

shirly 发表于 2014-10-8 09:22
参考前面别人的回复，尤其清风大侠的回复

如果只是一年的话，是可以画出来的，但是好几年就画不出来了。

绯宜物语

mofangbao 发表于 2013-4-21 20:18
我是说你最内的一个循环每次结束之后的那个write要改，in=1,n里面的是不用改的
因为内层的do循环结束之 ...

你好。我按照你说的这样写出来了，生成了map文件。但是sta.ctl打开后画图是这样的
这是我的程序。因为r1601.txt中只有降水数据，所以又读了station.txt读入站点号、经纬度。
program precititation
integer*4,parameter:: n=160,m=64
Character*8 stid(n)
Real lon(n),lat(n),tim
integer rain(n,m)
integer ::nlev,flag

open(1,file='f:\r1601.txt',status='old')
     read(1,*) ((rain(i,j),i=1,n),j=1,m)

open(2,file='f:\station.txt',status='old')
  do i=1,n
     read(2,*) stid(i),lon(i),lat(i)
  enddo

open(3,file='f:\sta.grd',status='replace',form='binary')

  do j=1,m
     tim=0.0
     nlev=1
     flag=1
      do i=1,n
        write(3) stid(i-1),lon(i-1),lat(i-1),tim,nlev,flag,rain(i,j)
      enddo
            nlev=0
        write(3) stid(n),lon(n),lat(n),tim,nlev,flag
  enddo
  close(1)
  close(2)
  close(3)
end

可以帮我看看是哪出错了吗？谢谢。

mofangbao

绯宜物语 发表于 2014-10-31 15:49
你好。我按照你说的这样写出来了，生成了map文件。但是sta.ctl打开后画图是这样的
这是我的程序。因为r ...

看看是不是经纬度写反了

绯宜物语

mofangbao 发表于 2014-10-31 20:51
看看是不是经纬度写反了

不是。之前写单时次的时候反了，后来改过来了。还专门注意了一下。