- 积分
- 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 |
|