SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,LAYER,KSPT)
C
      INCLUDE 'ABA_PARAM.INC' 
C
      DIMENSION STATEV(NSTATV),COORDS(NCRDS)
C
      STATEV(1)=1.D0
      STATEV(2)=2.D0
      STATEV(3)=3.D0
      STATEV(4)=4.D0
      STATEV(5)=5.D0
      STATEV(6)=6.D0
C      WRITE(6,*)'SDVINI-STATEV',(STATEV(IXX),IXX=1,6)
      RETURN
      END 
      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) 
C
      INCLUDE 'ABA_PARAM.INC' 
C
      CHARACTER*80 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
     1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
     2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
     3 DFGRD0(3,3),DFGRD1(3,3)
C
      IF(NDI.EQ.1) THEN
         DDSDDE(1,1)=PROPS(1)
         STRESS(1)=DDSDDE(1,1)*(STRAN(1)+DSTRAN(1))
      ELSE IF(NDI.EQ.2) THEN
         DO 10 K1=1,3
         DO 10 K2=1,3
         DDSDDE(K1,K2)=0.0D0
   10    CONTINUE
         DDSDDE(1,1)=PROPS(1)
         DDSDDE(2,2)=PROPS(1)
         DDSDDE(3,3)=0.5*PROPS(1)
         STRESS(1)=DDSDDE(1,1)*(STRAN(1)+DSTRAN(1))
         STRESS(2)=DDSDDE(2,2)*(STRAN(2)+DSTRAN(2))
         STRESS(3)=DDSDDE(3,3)*(STRAN(3)+DSTRAN(3))
      END IF
      STATEV(1)=STATEV(1)+1.D0*KINC  
      STATEV(2)=STATEV(2)+2.D0*KINC  
      STATEV(3)=STATEV(3)+3.D0*KINC  
      STATEV(4)=STATEV(4)+4.D0*KINC  
      STATEV(5)=STATEV(5)+5.D0*KINC  
      STATEV(6)=STATEV(6)+6.D0*KINC  
C      WRITE(6,*)'UMAT-KINC,STATEV',KINC,(STATEV(IXX),IXX=1,6)
      RETURN
      END