处理降水资料如何把按月分的日降水资料处理成按站点分的资料

高小一

资料是这样的
D:\杂七杂八\1.png
1961年1月，2月。。。。12月一直到2015年的1月。。。。12月共55年*12月共660个文件，每个文件里有所有站点在这一个月的所有天数的数据
文件格式大概是这样， 59981 1650 11220     49 1961  1  6  32700      1      1 0 0 0
后面六项是降水数据，32700为缺测
D:\杂七杂八\2.png
现在老师要求把这些资料处理成按站点分开的数据，就是xxxxx（站号）.txt，然后按数据说明怎样把数据读取出来并变成以站号分类，还有如何处理经纬度数据单位？拜托大神帮助
批处理读取所有文件的按照清风大神的帖子看了，读取数据也会，但是读取以后如何按不同站分开并输出成不同文件不懂，跪求帮助

等待中幸福

是有点复杂，但还是可以做的。其实就是读和写的问题，按照你的想法写就行了。

高小一

请问能帮忙看看这个程序吗，万分感谢
INTEGER,PARAMETER :: NS=839,NY=55,ND=31,NM=12,WJS=660
INTEGER year,mon,day
real    a(10)
character*5 sta(NS),filename(WJS)
real,allocatable :: staID(:),b(:),lat(:),lon(:),hgt(:),rain1(:,:,:,:),rain2(:,:,:,:),rain3(:,:,:,:)
allocate (staID(NS),b(WJS),lat(NS),lon(NS),hgt(NS),rain1(NS,NY,NM,ND),rain2(NS,NY,NM,ND),rain3(NS,NY,NM,ND))
!!!!!!日值说明，20-8时降水量、8-20时降水量、20-20时累计降水量、
!!!!!!     
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  READ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
OPEN(11,FILE='D:\data\data\station.txt')
do i=1,NS
read(11,*)staID(i)
enddo
CLOSE(11)

!===========1961-2015:
open(10,file='D:\data\data\pre\filename.txt')  
do i=1,WJS
  read(10,*) b(i)
  enddo  
close(10)   
do i=1,WJS
write(filename,'(b)')(b(i))   
  OPEN(11,FILE='D:\data\data\pre\'//filename//'.txt',err=8000)（提示这行有问题）

!!!!!!pre文件夹里含12*55个txt资料，文件名都是“SURF_CLI_CHN_MUL_DAY-pre-13011-xxxx(年）xx（月）”的形式
            read(11,*)
!!!!!!read(11,*)表示从文件11中读出一个数据，但不写入任何变量
500  read(11,*,err=600)a(1:10)
  DO is=1,ns
  if(a(1).eq.staID(is))then   
  lat(is)=a(2)
  lon(is)=a(3)
  hgt(is)=a(4)
  year=a(5)
  mon=a(6)
  day=a(7)
  rain1(is,year,mon,day)=a(8)
  rain2(is,year,mon,day)=a(9)
  rain3(is,year,mon,day)=a(10)
  endif
  enddo
   goto 500
CLOSE(11)
600 continue
8000 continue
ENDDO

DO IS=1,NS
  lat(is)=INT(lat(is)/100)+MOD(lat(is),100.)/60.
  lon(is)=INT(lon(is)/100)+MOD(lon(is),100.)/60.
  hgt(is)=hgt(is)/10
  rain1(is,year,mon,day)= rain1(is,year,mon,day)/10
  rain2(is,year,mon,day)= rain1(is,year,mon,day)/10
  rain3(is,year,mon,day)= rain1(is,year,mon,day)/10
  enddo
!按要求处理数据
do i=1,ns
WRITE(sta(is),'(i5)')INT(staid(is))

  OPEN(11,FILE='D:\data\clh\'//sta(is)//'.txt/',err=8000)
    do IY=1,NY
        DO IM=1,NM
        DO ID=1,ND
        WRITE(11,*)Lat(is),lon(is),hgt(IS),year,mon,day,rain1(is,iy,im,id),rain2(is,iy,im,id),rain3(is,iy,im,id)

  enddo
  enddo
  enddo
  CLOSE(11)
  enddo
!把数据按不同站分开的形式输出到不同以站号命名的文件里
deallocate (staID,lat,lon,hgt,rain1,rain2,rain3)
END