关于提取文件数据

奇迹一哥

我想把这个气象文件中date,tmax,tmin,rain四列提出，并且根据第一列判断年份然后按年份分别输出到多个文件（第一列81001的意思是81年第一天，一个文件有十几年）。现在我已经实现整列输提出，但是刚学编程不久，按年份输出还不会，不知谁有类似的代码，能不能分享一下 学习学习。

Penny_Kuan

你可以，提取date的前两位字符作为新的变量或用取整的方法保留年份，再根据这个变量来输出。

奇迹一哥

Penny_Kuan 发表于 2017-1-26 10:18
你可以，提取date的前两位字符作为新的变量或用取整的方法保留年份，再根据这个变量来输出。

你的意思我明白了，但是刚学编程不久，想问个比较基础的问题，就是我现在是把date设为动态数组integer, allocatable :: days(:)
根据年份要输出到多个文件，怎样根据年份自动生成多个文件并写入数据，再者文件名该怎末弄。

jstzysq

  character*8::filename
  integer::yr,dy,date,stat=0
  real::srad,tmax,tmin,rain
  do yr=1981,1999    !假设数据到99年
    write(filename,'(i4,".txt")'),yr     !技巧，年份变为字符串，并加上.txt，变成要输出的文件名1981.txt，......
    open(yr,file=filename)             !打开所有要输出的文件
  end do
  open(7,file=该气象数据)
  read(7,*);read(7,*);read(7,*);read(7,*)    !跳过前4行没用的行
  do while(stat==0)
    read(7,*,iostat=stat),date,srad,tmax,tmin,rain
    yr=date/1000       !参照2楼，把日期中的年份提出
    dy=date-yr*1000  !提取日期中的天
    yr=yr+1900         !年份变成一九几几年
    write(yr,'(i3,4f6.1)'),dy,srad,tmax,tmin,rain    !输出到对应文件
  end do
  close(7)
  do yr=1981,1999
    close(yr)               !关掉所有文件
  end do
  附赠一条在linux shell下只要一行就搞定的命令
  sed '1,4d' 该气象数据 | awk '{yr=substr($1,1,2)+1900;$1=substr($1,3,3);print > yr".txt"}'

奇迹一哥

jstzysq 发表于 2017-1-26 14:10
character*8::filename
  integer::yr,dy,date,stat=0
  real::srad,tmax,tmin,rain

太谢谢了