integra.for

来自「德国人开发的地震处理分析软件」· FOR 代码 · 共 43 行

$DEBUG

      PROGRAM INTEGRA

      DIMENSION ACEL(40000)

      WRITE(*,*)
      WRITE(*,*)
      WRITE(*,*)
      WRITE(*,*) '           SE ESTA INTEGRANDO EL REGISTRO...'
      WRITE(*,*)
      WRITE(*,*) '              CUANDO TERMINE EL CALCULO'
      WRITE(*,*) '            CIERRE LA VENTANA DOS Window'
      WRITE(*,*)
      WRITE(*,*)
      WRITE(*,*)

      OPEN(5,FILE='ACEL.DAT',STATUS='OLD')
      OPEN(6,FILE='SDAT.DAT',STATUS='OLD')
      OPEN(7,FILE='Avel.SAL',STATUS='UNKNOWN')
 
      READ(6,'(I8,F8.3)')NACEL,DT

      DO I=1,NACEL
         READ(5,'(E12.4)') ACEL(I)
      END DO

      A0=0.0
      VEL=0.0
      DO I=1,NACEL
         VEL=(ACEL(I)+A0)/2.*DT+VEL
         WRITE(7,'(E12.5)') VEL
         A0=ACEL(I)
      END DO

      CLOSE (5)
      CLOSE (6)
      CLOSE (7)

      STOP
      END