Fortran Coder

标题: 如何将连续列的数对应的最小值提取出 [打印本页]

作者: 978142355    时间: 2015-11-11 22:56
标题: 如何将连续列的数对应的最小值提取出
本人刚刚接触Fortran代码,想将图片红框中所有1的那一部分对应E的最小值获得,然后所有2对应E的最小值获得,3对应E的最小值……以此类推,然后将对应最小值那一行输出到文件中,比如:3个红色框的部分,将它们最小值获得后写入文件中变成:
Step number  39  scan point     1  E=  -571.251770020
Step number   5  scan point      2  E=  -571.252807617
Step number   4  scan point     3  E=  -571.253234863
……
这个代码是我自己尝试写的,不过问题还没解决,只解决了红框中那一列中有多少个数。还望各位有什么思路或者代码修改不吝惜赐教,谢谢。
[Fortran] 纯文本查看 复制代码
program takeminenergy
    implicit none
        integer i,j,num,num1,num2,lines,k,num4
    integer,allocatable::x(:),y(:)
    real,allocatable::z(:),minz(:),minumz(:)
    open(12,file='result.txt')!输入文件
        open(13,file='a.txt')   !输出文件
    open(14,file='b.txt')  !输出文件
        num1=1
        num2=0
    write(*,*) 'How many lines:' !读取文件有多少行
    read(*,*) lines
    allocate (x(lines))
    allocate (y(lines))
    allocate (minz(lines))
    allocate (minumz(lines))
    do i=1,lines
        read(12,"(12x,I4,12x,I6)") x(i),y(i)
            if (y(i)== num1) then
                num2=num2+1
            else if (y(i) /= num1) then
                y(num1)=num2
                num1=num1+1
                num2=1
        end if
    end do
    write(*,*) num1
   
        do j=1,num1
            write(13,"(I6)") y(j) !获取数字相同那一列有多少个,比如全是数字1的共有39个
        end do
        close(13)
        close(12)
        
        open(13,file='a.txt')
        open(12,file='result.txt')
        open(15,file='c.txt') !输出文件
        
        do j=1,num1
            read(13,"(I6)") num4
            write(*,*) num4
            do k=1,num4
             read(12,"(38x,f16.9)") z(k)
             write(15,*)z(k)
            end do
            
            minz=minval(z)
            write(14,"(f16.9)") minz
        deallocate(z)
        end do
        
end
     

图片1.jpg (198.64 KB, 下载次数: 253)

图片1.jpg

result.txt

546.88 KB, 下载次数: 2


作者: 楚香饭    时间: 2015-11-12 09:33
本帖最后由 楚香饭 于 2015-11-12 09:44 编辑

你想得太复杂了,其实问题可以很简单。
[Fortran] 纯文本查看 复制代码
Program www_fcode_cn
  Implicit None
  real(kind=kind(0.0d0)) :: rE , rMinE
  integer :: ierr , iRepeat , iPoint , pointerlast , iMinRepeat
  character(len=12) :: c(5)
  Open( 12 , File = 'result.txt' )
  Open( 13 , File = 'out.txt' )
  pointerlast = 1
  rMinE  = 1.0e30
  Do
    Read( 12 , * , ioStat=ierr ) c(1) , c(2) , iRepeat , c(3) , c(4) , iPoint , c(5) , rE
    if ( ( ierr== 0 ) .and. (iPoint == pointerlast) ) then !// 同一个Point
      if( rE <= rMinE ) then
        rMinE      = rE
        iMinRepeat = iRepeat
      end if
    else !//不同Point
      write(13,'(a,1x,a,1x,i4,2x,a,1x,a,i6,2x,a,f16.9)') trim(c(1)) , trim(c(2)) , iMinRepeat , &
        trim(c(3)) , trim(c(4)) , pointerlast , trim(c(5)) , rMinE
      pointerlast = iPoint
      iMinRepeat  = iRepeat
      rMinE = rE
    end if   
    if ( ierr /= 0 ) Exit
  End Do
  Close( 12 )
  Close( 13 )
End Program www_fcode_cn



作者: 978142355    时间: 2015-11-12 10:08
群主太无敌了,唉,我该如何才能达到你的这种境界。目前我也是一边看彭老师的书,一边练习编程,也关注论坛和qq群,但是到真正自己编程的时候发现思想总是很复杂,想问问群主,为什么你的编程思想可以这么简单,我看这代码里也没有什么特别特殊的东西。
作者: 楚香饭    时间: 2015-11-12 10:16
这个没有别的方法,多写多练吧。

很多时候同一个问题是可能有不同的思路的,你把这个代码看懂了,就理解他的思路了。以后遇到同类问题也不妨换个思路看看。
作者: 978142355    时间: 2015-11-12 11:32
懂了,多谢多练,思想自然就有了,再次谢谢群主。




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2