Fortran Coder

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

查看数: 9354 | 评论数: 4 | 收藏 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2015-11-11 22:56

正文摘要:

本人刚刚接触Fortran代码,想将图片红框中所有1的那一部分对应E的最小值获得,然后所有2对应E的最小值获得,3对应E的最小值……以此类推,然后将对应最小值那一行输出到文件中,比如:3个红色框的部分,将它们最小值获 ...

回复

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


捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-23 08:22

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表