青藏高原夏季逐日资料中提取678月的资料并做平均 程序能运行但提出的数据画的图...

xiahaha

问题描述：楼主需要提取青藏高原夏季（6、7、8月）热源的资料并做平均（做完应该是每年一个数据），但是自己编的程序提出的数据画的图不对，希望有人能帮看一下哪里有问题，楼主有点懵。
图片：第一张图和最后一张图
原始数据格式（逐日）：纬度-87.5~82.5  经度：2.5~357.5 格点：2.5*2.5
程序：
program EX22
implicit none
real q(143,69,365,35),q1(143,69,35)
integer i,j,d,y
!读取数据
open(1,file='d:\lw\data\q1_daily.grd',form='binary')
do y=1,35
do d=1,365
read(1) ((q(i,j,d,y),i=1,143),j=1,69)
enddo
enddo
!提取678月资料
q1(:,:,:)=0.0
do y=1,35
do d=153,244
do j=1,69
do i=1,143
q1(i,j,y)=q1(i,j,y)+q(i,j,d,y)
!q1(i,j,y)=q(i,j,153+12*(d-1),y)+q(i,j,154+12*(d-1),y)+q(i,j,155+12*(d-1),y)
!write(*,*) k
enddo;enddo;enddo;enddo
do y=1,35
do j=1,69
do i=1,143
q1(i,j,y)=q1(i,j,y)/92.0
enddo;enddo;enddo

open(2,file='d:\lw\data\q1.grd',form='binary')
open(22,file='d:\lw\data\q1.dat')
do y=1,35
write(2) ((q1(i,j,y),i=1,143),j=1,69)
write(22,*) ((q1(i,j,y),i=1,143),j=1,69)
enddo
end

刘晶

亲最好读入的时候顺便再写出了一个二进制文件，然后画个图和初始文件对比一下，看看是不是最开始读入就出现了问题

刘晶

亲最好读入的时候顺便再写出了一个二进制文件，然后画个图和初始文件对比一下，看看是不是最开始读入就出现了问题

xiahaha

刘晶 发表于 2017-4-27 12:30
亲最好读入的时候顺便再写出了一个二进制文件，然后画个图和初始文件对比一下，看看是不是最开始读入就出现 ...

好的 我改改试试

自行车大大爷

你是不是经度弄反了

xiahaha

自行车大大爷 发表于 2017-4-27 17:43
你是不是经度弄反了

啥反了0.0  能不能具体一点点

xiahaha

自行车大大爷 发表于 2017-4-27 17:43
你是不是经度弄反了

我..没有吧 导师让我改成这样的

lgz

写入一个二进制文件

xiahaha

lgz 发表于 2017-4-27 21:51
写入一个二进制文件

哪里..

自行车大大爷

xiahaha 发表于 2017-4-27 21:49
我..没有吧 导师让我改成这样的

试试经度逆序