帮我找找代码错误吧~~~~

crystalyang44

采用SODA2.0.2的海温资料计算垂直温度梯度最大处所在的深度，资料是0.5*0.5的分辨率，x方向720，y方向330，z有40个，t有600个，是月平均的，想法是先计算每一个格点上的垂直温度梯度，再比较得到垂直方向上最大的温度梯度所在深度，代码如下：
program main
implicit none
real dat(720,330,40,600),dt(720,330,40,600),dtmax(720,330,600),lev(720,330,600)
integer x,y,z,t
real level(40)
data level/5.01,15.07,25.28,35.76,46.61,57.98,70.02,82.92,96.92,112.32,129.49,148.96,171.4,197.79,229.48,268.46,317.65,381.39,465.91,579.31,729.35,918.37,1139.15,1378.57,1625.7,1875.11,2125.01,2375.0,2625.0,2875.0,3125.0,3375.0,3625.0,3875.0,4125.0,4375.0,4625.0,4875.0,5125.0,5375.0/                                                                                                  //垂直方向不是等距的格点，所以给出垂直方向上各格点所在深度

open(1,file='d:\shuju\yuanshi\yuanshi.dat',form='binary',status='old')
open(2,file='d:\shuju\yuanshi\dtmax.dat',form='binary',status='unknown')
open(3,file='d:\shuju\yuanshi\lev.dat',form='binary',status='unknown')
read(1)((((dat(x,y,z,t),x=1,720),y=1,330),z=1,40),t=1,600)


t=1                                                                                //先以第一个时间做实验
do z=1,40
  do y=1,330
   do x=1,720

lqouc

http://bbs.06climate.com/forum.php?mod=viewthread&tid=4571

风往北吹

up~~~

subtropical

楼主解决了吗