利用FORTRAN生成二进制数据之后编辑CTL文件，但是总是出错，代码如下，求大神指教！

weiyuntao

Fortran代码如下：
SUBROUTINE POST
IMPLICIT NONE
INTEGER,PARAMETER :: M = 73, N = 25, TT = 61
REAL*8 :: U(M,N,TT), V(M,N,TT), T(M,N,TT), Q(M,N,TT), WA(M,N,TT), H(M,N,TT), SST(M,N,TT)
REAL*8 :: SAN(M,N,TT), WEE(M,N,TT)
REAL*8 :: U1(M*N*tt), V1(M*N*tt), T1(M*N*tt), Q1(M*N*tt), WA1(M*N*tt), H1(M*N*tt), SST1(M*N*tt)
REAL*8 :: SAN1(M*N*tt), WEE1(M*N*tt)
integer :: it,i,j
! open original data files
OPEN(100, FILE = 'D:\work\CNOPI\cnopi\post_analysis\u', FORM = 'BINARY')
OPEN(200, FILE = 'D:\work\CNOPI\cnopi\post_analysis\v', FORM = 'BINARY')
OPEN(300, FILE = 'D:\work\CNOPI\cnopi\post_analysis\t', FORM = 'BINARY')
OPEN(400, FILE = 'D:\work\CNOPI\cnopi\post_analysis\q', FORM = 'BINARY')
OPEN(500, FILE = 'D:\work\CNOPI\cnopi\post_analysis\wa', FORM = 'BINARY')
OPEN(600, FILE = 'D:\work\CNOPI\cnopi\post_analysis\h', FORM = 'BINARY')
OPEN(700, FILE = 'D:\work\CNOPI\cnopi\post_analysis\sst', FORM = 'BINARY')
OPEN(800, FILE = 'D:\work\CNOPI\cnopi\post_analysis\san', FORM = 'BINARY')
OPEN(900, FILE = 'D:\work\CNOPI\cnopi\post_analysis\wee', FORM = 'BINARY')
READ(100) u1 ;READ(200) v1 ;READ(300) t1 ;READ(400) q1
READ(500) wa1 ;READ(600) h1 ;READ(700) sst1
READ(800) san1 ;READ(900) wee1
CLOSE(100); CLOSE(200); CLOSE(300); CLOSE(400); CLOSE(500)
CLOSE(600); CLOSE(700); CLOSE(800); CLOSE(900)
u=reshape(u1,(/m,n,tt/))
v=reshape(v1,(/m,n,tt/))
t=reshape(t1,(/m,n,tt/))
q=reshape(q1,(/m,n,tt/))
wa=reshape(wa1,(/m,n,tt/))
h=reshape(h1,(/m,n,tt/))
sst=reshape(sst1,(/m,n,tt/))
san=reshape(san1,(/m,n,tt/))
wee=reshape(wee1,(/m,n,tt/))

! open object data file
open(1000,file='ddd.dat',form='binary')
do it=1,tt
write(1000) ( (u(i,j,it),i=1,m),j=1,n)
write(1000) ( (v(i,j,it),i=1,m),j=1,n)
write(1000) ( (t(i,j,it),i=1,m),j=1,n)
write(1000) ( (q(i,j,it),i=1,m),j=1,n)
write(1000) ( (wa(i,j,it),i=1,m),j=1,n)
write(1000) ( (h(i,j,it),i=1,m),j=1,n)
write(1000) ( (sst(i,j,it),i=1,m),j=1,n)
write(1000) ( (san(i,j,it),i=1,m),j=1,n)
write(1000) ( (wee(i,j,it),i=1,m),j=1,n)
enddo
close(1000)

RETURN
END SUBROUTINE POST

CTL文件：
DSET  ^ddd.dat
TITLE Theorical Model Output Variables
UNDEF -9.99E33
XDEF 73 linear   0.000   5.000
YDEF 25 linear -60.000   5.000
ZDEF 1  levels 850
tdef 61 linear 12Z01jun2011 1dy
VARS 9
u 0 99 zonal wind
v 0 99 meridional wind
t 0 99 air temperature
q 0 99 specific humidity
wa 0 99  convective precipitation
h 0 99 mixed layer depth
sst 0 99  sea surface temperature
san 0 99  diverange
wee 0 99  Ekman Pumping
ENDVARS

用GrADS打开该CTL文件时，会出现下面情况：
虽然能打开，但是，时间维只有一个值，我定义的变量时间长度是61啊，
为什么打开之后，序列长度是1呢。而且，set t 1, d t之后，显示的都是0啊。
不知道哪里出问题了，有没有大神指教一下！！！多谢

543267158

写数据的时候，分开写试试，不要9个变量写到一个do里面，写成9个do

river

为什么时间维只有一？grads没有任何提示告诉你时间维只有一啊

lqouc

read的时候加循环

stone

应该是fortran写入的问题，

julie_7908

主要还是循环和recl要对应好