fortran处理micaps第四类数据————模式雨量运算

oh嘟嘟小千

本人fortran+Grads菜鸟一枚，最近突然被分配了一个任务，其中有画模式雨量图，但是给的资料是3小时的累计雨量，必去处理成24小时的。初学fortran编了一个程序，自己也知道有错，但是不知道怎么改，由于时间紧迫，所以求助各位大神！
现在有2个文件，EC雨量数据，14050120.48和14050120.24，需要的是14050120.48减去14050120.24的数据。
求各位指点修改
program main
integer,parameter::x=361,y=201
integer i,j
real one(y,x),two(y,x),th(y,x)
open(3,file='f:\de\EC\14050120.048',status='old')
read(1,*)
read(1,*)
read(1,*)((one(j,i),i=1,x),j=1,y)
close(3)
open(4,file='f:\de\EC\14050120.024',status='old')
read(1,*)
read(1,*)
read(1,*)((two(j,i),i=1,x),j=1,y)
close(4)
do j=1,y
  do i=1,x
     th(j,i)=one(j,i)-two(j,i)
  end do
end do
open(7,file='F:\bi\EC\2014050120.txt',status='replace',form='formatted')
do j=1,y
  do i=1,x
   write(3)th(j,i)
  end do
end do
close(3)
end

MeteoInfo

参考此贴：MeteoInfoLab脚本示例：micaps4类数据
http://bbs.06climate.com/forum.p ... 105&fromuid=106
(出处: 气象家园)

547168768

虽然我看不懂，但我知道楼主说的很高端

547168768

open(7,file='F:\bi\EC\2014050120.txt',status='replace',form='formatted')
请问这个文件在哪里？

wuxing3838

{:eb502:}{:eb502:}

船头尺

很好，很强大，多谢。

tulalang

楼主你这个程序不就是已经做好了么。。剩下就是grads出图了。。应该比较好完成的了

qxtlyf

建议直接写成二进制文件在Grads中，相减不更好吗?

whw334455

论坛上有第四类数据的读写方式，望参考