| 
 
	积分254贡献 精华在线时间 小时注册时间2017-2-7最后登录1970-1-1 
 | 
 
 
 楼主|
发表于 2017-2-27 10:51:07
|
显示全部楼层 
| 这是详细代码: 
 program turn
 implicit none
 integer nx,ny,nz,nt
 parameter(nx=37,ny=17,nz=2,nt=48)
 integer it,iz,i,j
 integer u(nx,ny,nz,nt),v(nx,ny,nz,nt)
 open(1,file='d:/grads/u.dat',form='formatted',status='old')
 open(2,file='d:/grads/v.dat',form='formatted',status='old')
 
 open(13,file='d:/gards/turn.grd',form='binary')
 
 do it=1,nt
 do  iz=1,nz
 read(1,*)((u(i,j,iz,it),i=1,nx),j=1,ny)
 read(2,*)((v(i,j,iz,it),i=1,nx),j=1,ny)
 end do
 end do
 close(1)
 close(2)
 
 do  it=1,nt
 do  iz=1,nz
 write(13)((u(i,j,iz,it),i=1,nx),j=1,ny)
 write(13)((v(i,j,iz,it),i=1,nx),j=1,ny)
 end do
 end do
 close(13)
 stop
 end
 | 
 |