geomview_out.f

来自「Multi Dimensional Scaling」· F 代码 · 共 44 行

      subroutine geomview_out(rval,notus1,nd1)
c
      implicit none
c
      integer notus1,nd1,ndc,ndc1
c
      parameter (ndc=96,ndc1=ndc-1)
c
      real rval(0:nd1,0:notus1),atgc(0:ndc1,0:notus1)
      real totatgc(notus1+1)
c
      integer iperm(notus1+1)
c
      integer iotus,id,kflag,ier,notus
c
      notus=notus1+1
c
      rewind (51)
      do 100 iotus=0,notus1
        read(51,*) (atgc(id,iotus),id=0,ndc1)
 100  continue
c
      do 110 iotus=1,notus
        totatgc(iotus)=0
        do 110 id=0,ndc1
          totatgc(iotus)=totatgc(iotus)+atgc(id,iotus-1)
 110    continue
c
        kflag=1
        call spsort(totatgc,notus1+1,iperm,kflag,ier)
c
        write(36,*) 'VECT'
        write(36,*) '1 ',notus,' 1 ',notus,' 1'
        do 120 iotus=1,notus1+1
          write(36,3000) (rval(id,iperm(iotus)-1),id=0,nd1)
 3000     format(3f10.5)
 120    continue
        write(36,*) '1 0 0 1'
c
        return
c
        end