#! /bin/sh
#
#     Last change:  Mihai Dima 2001
#
if [ $# -lt 2 ]
then
  echo " "
  echo "        *****     FORMINV     *****    "
  echo " "
  echo "          The program transforms "
  echo " a TEXT file (file1) in an  EXTRA file (file2)"
  echo " "
  echo "      Syntax: forminv <file1> <file2>"
  echo " "
  exit 0
else
 echo " "
 echo "        *****     FORMINV     *****    "
 echo " "
 echo "file1="$1
 echo "file2="$2
 cp $1 fort.1
fi
cat > forminv.f <<EOF
      PROGRAM PFORMIN
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'0**** FORMINV(FORMATTED,BINARY)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLESF (NCHAN,NDIM,NSET)
      CALL FORMIN (NDIM)     
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLESF (NCHAN,NDIM,NSET)
      IMPLICIT INTEGER (I-N)
      IMPLICIT REAL (A-H,O-Z)
C
C  SUBROUTINE FOR READING THE DIMENSION PARAMETER OF THE FIELDS ON
C
      NSET=0

      REWIND NCHAN
      READ(UNIT=NCHAN,FMT='(4I10)',END=10) I1,I2,I3,I4
      GOTO 20
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
   20 REWIND NCHAN
      NDIM=0
    1 READ(NCHAN,'(4I10)',END=100) I1,I2,I3,I4
      print*,'parles',NSET,I1,I2,I3,I4
      NDIM=MAX0(NDIM,I4)
      IVIER=I4/6
      ICONT=IVIER*6
      IF(ICONT.NE.I4) THEN
       IVIER=IVIER+1
      END IF
      DO LOOP=1,IVIER
       READ (NCHAN,'(1X,E12.5)',END=101) F
      END DO
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      REWIND NCHAN
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C*****************************************************************
      SUBROUTINE FORMIN(IE)
      IMPLICIT REAL (A-H,O-Z)
C
C  PURPOSE: transformation of formatted (DFORM) output into
C           T21 standard data 
C
      PARAMETER (UNDEF=0.9E+10,UNDEFE=.9*UNDEF)
      DIMENSION ID(4)
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C FOLLOWING LINES
C      PARAMETER (IE77=9E12)
C      DIMENSION  X(IE)
C      POINTER (NX,X)
C      ALLOCATE(NX,IE)
C AND ACTIVATE THE FOLLOWING LINES
      PARAMETER (IE77=20000)
      DIMENSION  X(IE77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      N = 0
      REWIND 1
    1 READ(1,'(4I10)',END=2) ID
        IF(N.EQ.0) WRITE(6,6010) ' FIRST RECORD', ID
        READ(1,'(6(1X,E12.5))',END=3) (X(I),I=1,ID(4))
        N = N+1
        DO II=1,ID(4)
          IF (X(II).GT.0.8E+10) THEN
            X(II)=UNDEF
          ENDIF
        ENDDO
        WRITE(2) ID
        WRITE(2) (X(I),I=1,ID(4))
        GOTO 1
C
    2 WRITE(6,6010) ' LAST RECORD ', ID 
      WRITE(6,6010) ' NUMBER OF RECORDS: ', N
      WRITE(6,6010) ' REGULAR END '
      STOP
C
    3 WRITE(6,6010) 'UNEXPECTED *EOF*: ', ID
      PRINT *, 'DFORM: fatal end ' 
      STOP 
C
 6010 FORMAT(' *FORMINV *',A,5I10)
      END 
EOF
f77 forminv.f -o forminv.x
forminv.x $1 $2
cp fort.2 $2
rm fort.1 fort.2
rm forminv.x forminv.f
exit