- 积分
- 106
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-1-20
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
fortran程序: parameter(ns=18,nm=1)
integer sta(ns)
real t(ns,nm),lon(ns),lat(ns)
print*, 'read sta inf'
open(11,file='c:\users\tt\desktop\spring.txt')
do i=1,ns
read(11,*) lat(i),lon(i),t(i,1)
enddo
close(11)
OPEN(22,FILE='c:\users\tt\desktop\spring.dat',FORM='binary')
call staout(t,ns,lon,lat,nm)
close(22)
stop
end
subroutine staout(coryp,ny,lon,lat,kvt)
real coryp(ny,kvt)
real lon(ny),lat(ny)
character*8 sta
irec=0
do it=1,kvt
do i=1,ny
irec=irec+1
tim=0.0
ilev=1
nflag=1
sta=char(i)
rr=float(i)
write(22) sta,lat(i),lon(i),tim,ilev,nflag,coryp(i,it)
enddo
irec=irec+1
ilev=0
write(22) sta,lat(ny),lon(ny),tim,ilev,nflag
enddo
return
end
txt文件格式:
149.50 28.40 0001
139.00 28.20 0002
169.10 38.30 0203
169.50 34.80 0303
153.00 29.00 0401
145.80 37.40 0402
148.70 27.70 0503
124.90 31.60 0601
160.00 27.00 0701
156.00 30.80 0702
141.90 31.80 0802
142.50 31.90 0804
147.30 30.20 0901
143.30 36.40 1101
134.10 32.80 1102
151.90 30.70 1202
155.90 23.00 1403
145.00 24.90 1405
运行之后得不出dat文件,请各位帮忙看看,感谢感谢感谢~~~
|
|