如何将连续列的数对应的最小值提取出

978142355 · 发表于 2015-11-11 22:56:32

本人刚刚接触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] 纯文本查看 复制代码

01program takeminenergy
02    implicit none
03  integer i,j,num,num1,num2,lines,k,num4
04    integer,allocatable::x(:),y(:)
05    real,allocatable::z(:),minz(:),minumz(:)
06    open(12,file='result.txt')!输入文件
07  open(13,file='a.txt')   !输出文件
08    open(14,file='b.txt')  !输出文件
09  num1=1
10  num2=0
11    write(*,*) 'How many lines:' !读取文件有多少行
12    read(*,*) lines
13    allocate (x(lines))
14    allocate (y(lines))
15    allocate (minz(lines))
16    allocate (minumz(lines))
17    do i=1,lines
18        read(12,"(12x,I4,12x,I6)") x(i),y(i)
19      if (y(i)== num1) then
20          num2=num2+1
21      else if (y(i) /= num1) then
22          y(num1)=num2
23          num1=num1+1
24          num2=1
25        end if
26    end do
27    write(*,*) num1
28     
29        do j=1,num1
30            write(13,"(I6)") y(j) !获取数字相同那一列有多少个，比如全是数字1的共有39个
31        end do
32        close(13)
33        close(12)
34         
35        open(13,file='a.txt')
36        open(12,file='result.txt')
37        open(15,file='c.txt') !输出文件
38         
39        do j=1,num1
40            read(13,"(I6)") num4
41            write(*,*) num4
42            do k=1,num4
43             read(12,"(38x,f16.9)") z(k)
44             write(15,*)z(k)
45            end do
46             
47            minz=minval(z)
48            write(14,"(f16.9)") minz
49        deallocate(z)
50        end do
51         
52end
53     

978142355 · 发表于 2015-11-12 11:32:29

懂了，多谢多练，思想自然就有了，再次谢谢群主。

楚香饭 · 发表于 2015-11-12 10:16:15

这个没有别的方法，多写多练吧。

很多时候同一个问题是可能有不同的思路的，你把这个代码看懂了，就理解他的思路了。以后遇到同类问题也不妨换个思路看看。

978142355 · 发表于 2015-11-12 10:08:28

群主太无敌了，唉，我该如何才能达到你的这种境界。目前我也是一边看彭老师的书，一边练习编程，也关注论坛和qq群，但是到真正自己编程的时候发现思想总是很复杂，想问问群主，为什么你的编程思想可以这么简单，我看这代码里也没有什么特别特殊的东西。

楚香饭 · 发表于 2015-11-12 09:33:14

本帖最后由楚香饭于 2015-11-12 09:44 编辑

你想得太复杂了，其实问题可以很简单。

[Fortran] 纯文本查看 复制代码

01Program www_fcode_cn
02  Implicit None
03  real(kind=kind(0.0d0)) :: rE , rMinE
04  integer :: ierr , iRepeat , iPoint , pointerlast , iMinRepeat
05  character(len=12) :: c(5)
06  Open( 12 , File = 'result.txt' )
07  Open( 13 , File = 'out.txt' )
08  pointerlast = 1
09  rMinE  = 1.0e30
10  Do
11    Read( 12 , * , ioStat=ierr ) c(1) , c(2) , iRepeat , c(3) , c(4) , iPoint , c(5) , rE
12    if ( ( ierr== 0 ) .and. (iPoint == pointerlast) ) then !// 同一个Point
13      if( rE <= rMinE ) then
14        rMinE      = rE
15        iMinRepeat = iRepeat
16      end if
17    else !//不同Point
18      write(13,'(a,1x,a,1x,i4,2x,a,1x,a,i6,2x,a,f16.9)') trim(c(1)) , trim(c(2)) , iMinRepeat , &
19        trim(c(3)) , trim(c(4)) , pointerlast , trim(c(5)) , rMinE
20      pointerlast = iPoint
21      iMinRepeat  = iRepeat
22      rMinE = rE
23    end if   
24    if ( ierr /= 0 ) Exit
25  End Do
26  Close( 12 )
27  Close( 13 )
28End Program www_fcode_cn

[求助] 如何将连续列的数对应的最小值提取出

规矩勋章

新人勋章

爱心勋章

水王勋章

元老勋章

热心勋章