回 帖 发 新 帖 刷新版面

主题:求助怎么连续坐标编号

子程序是离散一条线,分成若干段,写出坐标还有标号,n是单位长度x0y0xnyn是端点坐标
主程序是离散两条线,有一个端点坐标是相同的   
  
    
     program ji    
        PARAMETER (NPP=2100000,MEE=2100000)
    real x(1000),y(1000)
    integer   linenum(10),i,n
      open(1,file='line.txt',status='old')
    open(2,file='g4.dat')   
         read(1,*)  n    !   
     read(1,*) lnum
      do 10 i=1,lnum
     
      read(1,*)linenum(i),x(i),y(i)   
   10 continue
       call disline (n,x(1),y(1),x(2),y(2))
       call disline (n,x(2),y(2),x(3),y(3))

       write (*,*)
     stop
     end

         subroutine disline(n,x0,y0,xn,yn)
    real x(100),y(100)

        long=sqrt((x0-xn)*(x0-xn)+(y0-yn)*(y0-yn))!线段长度
        num=int(long/n)                       !节点数
        unitlo=long/num                           !等分线段
        dx=unitlo*(xn-x0)/long
        dy=unitlo*(yn-y0)/long 
    l=0
    do 10 j=1,num+1
        l=l+1
        x(l)=x0+dx*(j-1)
    y(l)=y0+dy*(j-1)
10    continue
        nod=l
        do 81 i=1,nod
    write(2,109) i,x(i),y(i)
 81    continue
109   format(i8,2E15.6)
    end
这个是line文件

3
1   0.0 0.0 
2   3.0 0.0
3   1.5 1.5

回复列表 (共10个回复)

沙发

好像您贴了很几回了,不过确实没看懂。

板凳

我也是猜不出楼主什么意思. 不过楼主些代码的风格不好, 不小心就有bug了.

3 楼

着急用啊,有谁能帮我解答下

4 楼


我是新手,能指导下吗

5 楼

你连问题都描述不清,试问我们该如何帮助?

6 楼

我是学土木的,不懂fortran,怎么能描述明白啊

7 楼

你知道或者有什么, 你想实现什么, 打算通过什么方法实现, 现在如何尝试.

8 楼

[quote]我是学土木的,不懂fortran,怎么能描述明白啊[/quote]

  与Fortran无关,与专业无关,是你对自己的提问不负责任。你连问题都不表达清楚,谁能帮你,谁愿帮你?

9 楼


不能光批评啊,给点具体的意见怎么能说明白

10 楼

[quote]
不能光批评啊,给点具体的意见怎么能说明白[/quote]如7楼。
你在做什么方向的什么问题
遇到了什么样的问题
想要怎么处理

我来回复

您尚未登录,请登录后再回复。点此登录或注册