- 积分
- 7124
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-11-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
!how to compile:: pgf90 -byteswapio liu_read.f90 -o liu.out
!execution:: ./liu.out
!
program one
implicit none
!parameter(nx=31,ny=21,np=4,nt=12)
integer nx
integer ny
integer np
integer nt
integer i,j,ip,it
nx=31
ny=21
np=4
nt=12
!i=1
!j=1
!ip=1
!it=1
character*80 fname
real sx(nx,ny,np,nt)
!nx=e-w grids; ny=s-n grids; nt=12 months
!np=4::industry, energy, transportation, residential
fname='liu_SO2_monthly.grd'
open(11,file=trim(fname),status='old', &
form='unformatted',access='direct',recl=nx*ny*4)
kk=0
do it=1,nt
do ip=1,np
kk=kk+1
read(11,rec=kk) ((sx(i,j,ip,it),i=1,nx),j=1,ny)
enddo
enddo
close(11)
!write out
end
|
|