PROGRAM SET1D1
C=======================================================================
C                   DATA SET UP PROGRAM FOR ONEDEE2.FOR
C     ONE DIMENSIONAL TIME DEPENDENT CONVECTION-DIFFUSION EQUATION
C                    METHOD: FINITE ELEMENT METHOD
C                         ELEMENT:  LINEAR
C                   EIJI FUKUMORI      SEPTEMBER 1985
C=======================================================================
      PARAMETER ( IW = 1, MXE=2000, MXN=MXE+1 )
      IMPLICIT REAL * 8 ( A-H, O-Z )
      DIMENSION X(MXN), RHO(MXN), V(MXN), DB(MXN), S(MXN), C(MXN)
C=======================================================================
C---------- CREATING INPUTDATA ---------
      NE = 100
      NNODE = NE + 1
      TL = 10.D0
      DT = 0.05D0
      TMAX = 34.D0
      OPEN ( IW, FILE ='ONEDEE.DAT', STATUS='UNKNOWN')
      WRITE (IW,*) NNODE
      WRITE (IW,*) DT
      WRITE (IW,*) TMAX
C=======================================================================
      IB1 = 1
      BV1 = 0.D0
      IBN = 1
      BVN = 0.D0
      WRITE (IW,*) IB1, BV1
      WRITE (IW,*) IBN, BVN
C=======================================================================
      DX = TL / NE
      DENSITY = 1.
      VELOCITY = 0.2
      DIFFUSE = 0.001
      DO  I = 1 , NNODE
      X(I) = (I-1)*DX
      RHO(I) = DENSITY
      V(I) = VELOCITY
      DB(I) = DIFFUSE
      WRITE (IW,*) X(I),RHO(I),V(I),DB(I)
      END DO
C=======================================================================
      SPIKE = 1000.D0
      DO I = 1 , NNODE
      S(I) = 0.D0
      END DO
      S(5) = SPIKE
      WRITE (IW,*) ( S(I) , I = 1 , NNODE )
      CLOSE (IW)
      STOP
      END