【已解决】关于fortran提取特定月份资料的问题

爱知天气

问了同学没问出来，于是求助。先谢谢大家了。
现在我有1958-1977年的V风月平均资料，一共20年，共240个月，这时的资料经过GrADS画图检查，都是正确的。现在我想提取每年6--10月的资料出来，也就是一共20*5=100个月的资料。于是我写了一个fortran程序，但是提取出来的资料有问题，20年中每一年的6月分画出的图都是一样的（7、8、9、10月同样如此），也就是最后一年的6月应该把之前年份的6月全部覆盖掉了。我传一张1958年6月的图和1959年6月的图，它们是一样的。下面是fortran程序：
program main                 !!!! !!! This program is choosing the  data in June to October and rewrite them in a new file for better using
implicit none
integer i,j,t,yr
integer,parameter::n=20     !20年时间，1958—1977年
real v(53,21,240)

!!!把grd文件读入数组!!!
open(10,file='E:\BY\figure2\data\v1958-1977.grd',form="binary")
read(10)(((v(i,j,t),i=1,53),j=1,21),t=1,240)
close(10)

!!!提取所需数据重新写入新二进制文件，便于处理
open(11,file='E:\BY\figure2\data6to10\v1958-1977JtoO.grd',form="binary")
do yr=1,n
   do t=(n-1)*12+6,(n-1)*12+10       ! 提取6-10月的资料
           do j=1,21
              do i=1,53
          if(v(i,j,t)<-999.) then
           v(i,j,t)=-9999.
        endif
                 write(11)v(i,j,t)
              end do
        end do
      end do
end do
close(11)
end



这是V的ct文件：
dset E:\BY\figure2\data6to10\omega1958-1977JtoO.grd
title vwind 1958-1977 JUN to OCT
undef -9.99e08
xdef 53 linear 70 2.5
ydef 21 linear -10 2.5
zdef 1 levels 500
tdef 100 linear 00Z01JUN1958 1mo
vars 1
v 1 0 Monthly v
endvars

爱知天气

把原来的数据和运行后的数据也挂上把

蓝莲花

你把gs文件也贴出来吧

爱知天气

蓝莲花 发表于 2014-4-26 15:55
你把gs文件也贴出来吧

测试我是直接输入的：现在我打一下
open ....\v1958-1977JtoO.ctl
set t 1
d v    (看1958年6月的图)
c
set t 6
d v   (看1959年6月的图)


发现这两图是一样的。如果没有C，直接就叠加上去了，就是第二章白色的图

爱知天气

蓝莲花 发表于 2014-4-26 15:55
你把gs文件也贴出来吧

如果   set t 1 20
d v
会发现图不会变化

言深深

别的地方没有细看，不错目测标记“ ! 提取6-10月的资料”那一行里面的【n】应该改成【yr】，毕竟yr才是循环变量。n一直不循环，下面你出来的自然出错了···

爱知天气

言深深 发表于 2014-4-26 16:40
别的地方没有细看，不错目测标记“ ! 提取6-10月的资料”那一行里面的【n】应该改成【yr】，毕竟yr才是循环 ...

已修正，果然对了！十分感谢！原来经常用N改变，现在换了竟然一时不知道改过来了。

言深深

爱知天气 发表于 2014-4-26 17:02
已修正，果然对了！十分感谢！原来经常用N改变，现在换了竟然一时不知道改过来了。

解决了就好···

爱知天气

言深深 发表于 2014-4-26 21:33
解决了就好···

自己粗心大意了，真不好意思。再次感谢

默戏

我也遇到类似问题，感谢楼主！