报文文件名循环

小孩别要太多 · 发表于 2019-3-21 11:55:08

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由小孩别要太多于 2019-3-21 11:55 编辑

我编写了一段读取报文的程序，但是只能读取一天的报文，做了文件名循环后我想得出一个月的数据，但是程序出了一些问题，想请教一下大家该如何改正这个程序，谢谢大家！程序如下：
program main
      real jingdu(1:78),weidu(1:78),tz(1:201),j1,j2,D(01:30)
      character sta(80)*5,shici(1:7,1:4)*3
      real::tmin(7),tmax(7),t1,t2,E
   real::ttmin(7),ttmax(7),tt1,tt2
      character(len=100) ::lujing,luu
   CHARACTER(10)::t
   Character(8)::str
      Character(83)::str1,str2,str3,str4,str5,str6


   CHARACTER(5)::z
   CHARACTER(4)::Y,iii1
      CHARACTER(2)::Y1,M,DD,jjj,kkk,iii
      integer::yy,yy1,mm,dd1,st
      integer:: ALIVE
      CHARACTER(4)::zifu1
   CHARACTER(2)::D1
   CHARACTER(2)::E1
   CHARACTER(18)::zifu2
      CHARACTER(34)::zifu3
      CHARACTER(17)::zifu4
      CHARACTER(3)::zifu5
      CHARACTER(32)::zifu6
      CHARACTER(129)::zifu7
   CHARACTER(34)::zifu10
      CHARACTER(64)::zifuq
      CHARACTER(54)::zifuh
   CHARACTER(5)::zifu11

      CHARACTER(5)::zhanhao1,zhanhao2,zhanhao3,zhanhao4,
   &zhanhao5,zhanhao6

      zhifu5=' 19'

   !    读出系统时间
   !CALL DATE_AND_TIME( DATE =str, TIME = t, ZONE = z)
      !Y=str(1:4)

   !Y1=str(3:4)
      !M=str(5:6)
      !D=str(7:8)
   !print*,Y,'年',M,'月',D,'日'


      do D=01,30

   E=D+1
   D1='D'
      E1='E'

   str1='Z_SEVP_C_BFLZ_201901'//D1//'213002_P_RFFC-SLMOF-201901'//E1/
   $/'0000-16803.txt'

   str2='Z_SEVP_C_BFLX_201901'//D1//'224500_P_RFFC-SPCC-201901'//E1//
   $'0000-16812.txt'

   !str1='Z_SEVP_C_BCLZ_20190312213002_P_RFFC-SLMOF-201903130000-16803
   !&.txt'

   !str2='Z_SEVP_C_BCLX_20190312083000_P_RFFC-SPCC-201903121200-16812.
   !&txt'

      open(00,file='F:\baowen\quanzh\'//str1)
      open(01,file='F:\baowen\zhidaoyubao\zhongxintai\'//str1)
      open(02,file='F:\baowen\zhidaoyubao\zhouju\zaoshang\'//str2)
   do i=1,1772
      read(01,*)
      enddo
      do n=1,51
   read(02,*)
      enddo
   read(01,600)zifu6
      write(00,600)zifu6
      do i=1,7
         do j=1,3
      read(01,*)
         enddo
         read(01,800)t1max,t1min
      do j=1,3
      read(01,*)
         enddo
         read(01,800)t2max,t2min
      if(t1max>t2max) then
         tmax(i)=t1max
         else
         tmax(i)=t2max
         end if
      if(t1min>t2min) then
         tmin(i)=t2min
         else
         tmin(i)=t1min
         end if
      read(02,*)
         read(02,800)ttmax(i),ttmin(i)
         write(00,889)tmin(i),tmax(i),ttmin(i),ttmax(i)
   end do
      END DO
      close(01)
      close(02)


      write(00,998)'NNNN'
      close(00)



100 format(a4)
200 format(a18)
300 format(a34)
400 format(a17)
500 format(a3)
600 format(a32)
700 format(a129)
800 format(64x,f5.1,1x,f5.1)
900 format(70x,f5.1)
901 format(59x,f5.1)
998 format(a4)
999 format(a31)
888 format(f5.1,1x,f5.1)
889 format(f5.1,1x,f5.1,1x,f5.1,1x,f5.1)
      end

jingli41 · 发表于 2019-3-21 12:29:55

do循环里，天是数字型，但在文件路径名称里，天是字符型，如果天数小于10，例如：2，在字符型里显示的是2，这样文件就会找不到。虽然你循环里是02.

小孩别要太多 · 发表于 2019-3-21 13:59:47

jingli41 发表于 2019-3-21 12:29
do循环里，天是数字型，但在文件路径名称里，天是字符型，如果天数小于10，例如：2，在字符型里显示的是2， ...

谢谢您！我理解您的意思了，但是字符型的文件名该如何循环呢？

小孩别要太多 · 发表于 2019-3-21 17:58:10

修改后的编码将日数改成了字符型，但还是不对，循环出现了问题，望大佬指正，谢谢！

program main
real jingdu(1:78),weidu(1:78),tz(1:201),j1,j2,D(01:30)
character sta(80)*5,shici(1:7,1:4)*3
real::tmin(7),tmax(7),t1,t2
   real::ttmin(7),ttmax(7),tt1,tt2
character(len=100) ::lujing,luu
   CHARACTER(10)::t
   Character(8)::str
Character(83)::str1,str2,str3,str4,str5,str6
Character(30)::D(ii),E(jj)
Character(2)::ii(01,30),jj(02,31)

   CHARACTER(5)::z
   CHARACTER(4)::Y,iii1
CHARACTER(2)::Y1,M,DD,jjj,kkk,iii
integer::yy,yy1,mm,dd1,st
integer:: ALIVE
CHARACTER(4)::zifu1
   CHARACTER(18)::zifu2
CHARACTER(34)::zifu3
CHARACTER(17)::zifu4
CHARACTER(3)::zifu5
CHARACTER(32)::zifu6
CHARACTER(129)::zifu7
   CHARACTER(34)::zifu10
CHARACTER(64)::zifuq
CHARACTER(54)::zifuh
   CHARACTER(5)::zifu11

CHARACTER(5)::zhanhao1,zhanhao2,zhanhao3,zhanhao4,
   &zhanhao5,zhanhao6



zhifu5=' 19'

   !    读出系统时间
   !CALL DATE_AND_TIME( DATE =str, TIME = t, ZONE = z)
!Y=str(1:4)

   !Y1=str(3:4)
!M=str(5:6)
!D=str(7:8)
   !print*,Y,'年',M,'月',D,'日'




   str1='Z_SEVP_C_BFLZ_201901'//D(ii)//'213002_P_RFFC-SLMOF-201901'//
   $E(jj)//'0000-16803.txt'

   str2='Z_SEVP_C_BFLX_201901'//D(ii)//'224500_P_RFFC-SPCC-201901'//
   $E(jj)//'0000-16812.txt'

   !str1='Z_SEVP_C_BCLZ_20190312213002_P_RFFC-SLMOF-201903130000-16803
   !&.txt'

   !str2='Z_SEVP_C_BCLX_20190312083000_P_RFFC-SPCC-201903121200-16812.
   !&txt'

open(00,file='F:\baowen\quanzh\'//str1)
open(01,file='F:\baowen\zhidaoyubao\zhongxintai\'//str1)
open(02,file='F:\baowen\zhidaoyubao\zhouju\zaoshang\'//str2)
   do i=1,1772
read(01,*)
enddo
do n=1,51
   read(02,*)
enddo
   read(01,600)zifu6
write(00,600)zifu6
do i=1,7
   do j=1,3
      read(01,*)
   enddo
   read(01,800)t1max,t1min
      do j=1,3
      read(01,*)
   enddo
   read(01,800)t2max,t2min
      if(t1max>t2max) then
   tmax(i)=t1max
   else
   tmax(i)=t2max
   end if
      if(t1min>t2min) then
   tmin(i)=t2min
   else
   tmin(i)=t1min
   end if
      read(02,*)
   read(02,800)ttmax(i),ttmin(i)
   write(00,889)tmin(i),tmax(i),ttmin(i),ttmax(i)
   end do

close(01)
close(02)


write(00,998)'NNNN'
close(00)

100 format(a4)
200 format(a18)
300 format(a34)
400 format(a17)
500 format(a3)
600 format(a32)
700 format(a129)
800 format(64x,f5.1,1x,f5.1)
900 format(70x,f5.1)
901 format(59x,f5.1)
998 format(a4)
999 format(a31)
888 format(f5.1,1x,f5.1)
889 format(f5.1,1x,f5.1,1x,f5.1,1x,f5.1)
end

jingli41 · 发表于 2019-3-21 20:08:57

举个例子：
char*2 a
do i=1,30
  write(a,"(I2)")i    *将整数转换成字符
  if(i.lt.10)then
   str1='2019030'//trim(adjustl（a）)
  else
   str1='201903'//trim(adjustl（a）)
  endif
enddo
上面的程序你理解了，文件名循环应该就没问题了。

小孩别要太多 · 发表于 2019-11-18 15:47:46

jingli41 发表于 2019-3-21 12:29
do循环里，天是数字型，但在文件路径名称里，天是字符型，如果天数小于10，例如：2，在字符型里显示的是2， ...

谢谢您，现在知道了！

yijianmingliu · 发表于 2019-11-24 16:16:20

为什么不用python呢

		自动登录	找回密码
密码			立即注册

[求助] 报文文件名循环

登录后查看更多精彩内容~

回帖奖励 +1 金钱