#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     EXT1LTEXT2C     *****      "
  echo " "
  echo "The programs transforms an EXTRA file <file1> with one record"
  echo "              in a TEXT file with two columns"
  echo " "
  echo "     Syntax: ext1ltext2c <file1> <file2> <start_date>"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     EXT1LTEXT2C     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "Start date in the TEXT file="$3
 cp $1 fort.1
fi
cat > ext1ltext2c.f <<EOF
      PROGRAM EXTRA2TEXT
      PARAMETER(NCHAN1=1,UNDEF=0.9E+10)
      
      REAL X(10000)
      INTEGER IDAT,IVAR,ILEV,ILEN
      READ(*,*) IINIT
      CALL PARLES (NCHAN1,IDIM1,ILEN)

      OPEN(2,FORM='FORMATTED',STATUS='UNKNOWN')

      READ(NCHAN1) IDAT,IVAR,ILEV,ILEN
      READ(NCHAN1)(X(I),I=1,ILEN)
      DO I=1,ILEN
        WRITE(2,67) I+IINIT-1,X(I)
67    FORMAT(I4,1X, E12.5)
c        WRITE(2,'(1X,E12.5)') X(I)       
      ENDDO

      PRINT*,'REGULAR END!'
      CLOSE(2)

      END
C*********************************************************************
      SUBROUTINE PARLES (NCHAN,IDIM,ILEN)
      IMPLICIT REAL (A-H,O-Z)
      INTEGER I1,I2,I3
      REAL F

      ILEN=0

      REWIND NCHAN
      READ(NCHAN,END=10) ITIME,INAME,ILEV,ISIZE
      GOTO 20
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
   20 REWIND NCHAN

      IDIM=0  
    1 READ (NCHAN,END=100) I1,I2,I3,I4
      IDIM=MAX0(IDIM,I4)
      READ (NCHAN,END=101) F
      ILEN=ILEN+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,IDIM,ILEN
      REWIND NCHAN
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  IDIM=',I5,',  ILEN=',I5)
      END
EOF
f77 ext1ltext2c.f -o ext1ltext2c.x
ext1ltext2c.x $1 $2 << M
$3
M
cp fort.2 $2
rm fort.1 fort.2
rm ext1ltext2c.x ext1ltext2c.f
exit