- 积分
- 1501
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-8-6
- 最后登录
- 1970-1-1

|

楼主 |
发表于 2015-4-28 21:35:03
|
显示全部楼层
!!!!!!!!!!!!!!130801(2).dat是下面这个批处理程序生成的。130801.dat 是上面那个单时次程序生成的。大小是一样的,但是用grads画图的时候一个可以用一个不能用!!!!!!!!!!!!!!!!!
program main
implicit none
integer,parameter:: n=744
character*8 stid
character*12 filename,filena
real lon,lat,dir,v,tim
integer nlev,flag,error,i
character*12 a(n)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
tim=0.0
nlev=1
flag=1
open (100,file='b1308.txt')
do i=1,n
read(100,*) a(i)
end do
close(100)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do i=1,n
write(filename,"(a12)")a(i)
write(filena,"(a8)")a(i)
open (101,file='G:\data\grads_station\1308\'//filename,status='old',iostat=error)
if(error/=0)then
write(*,*) 'open file failed.'
end if
open (102,file='G:\data\grads_station\1308\'//trim(adjustl(filena))//'.dat',form='binary')
10 read(101,*,end=100)stid,lon,lat,dir,v
! write(*,*) stid,lon,lat,dir,v
! pause
write(102)stid,lat,lon,tim,nlev,flag,dir,v
goto 10
100 continue
nlev=0
write(102)stid,lat,lon,tim,nlev,flag
close(101)
close(102)
end do
end |
|