anchor.diff:  
   nContours = 10 ! [-]  number of contours for J-integral evaluation
   
  -EL = ! [MPa] UD longitudinal Young's modulus
  -ET = ! [MPa] UD transverse Young's modulus
  -nuLT = ! [-] UD in-plane Poisson ratio
  -nuTT = ! [-] UD transverse Poisson ratio
  -GL = ! [MPa] UD in-plane shear modulus
  -GT = ! [MPa] UD transverse shear modulus
  +EL = 3500.0! [MPa] UD longitudinal Young's modulus
  +ET = 3500.0! [MPa] UD transverse Young's modulus
  +nuLT = 0.4! [-] UD in-plane Poisson ratio
  +nuTT = 0.4! [-] UD transverse Poisson ratio
  +!GL = ! [MPa] UD in-plane shear modulus
  +!GT = ! [MPa] UD transverse shear modulus
   
   ! ===> END INPUT DATA
   
  
positive.code: ! ===> START INPUT DATA
  
  CrackDispsFile = 'C:/Ansys_WD/Debug/crackDisps'
  
  Vf = ! [-] Fiber volume fraction
  
  
negative.code: 
  LSEL,S,,,8                                ! Crack
  NSLL,S,1                                  !Select nodes associated to this line
  *GET,NNodes,NODE,0,COUNT             !Get the number of nodes in the selected set
  *DIM, CrackDisps, ARRAY, NNodes, 2     
  *VGET, CrackDisps(1,1), NODE, 1, U, X
  *VGET, CrackDisps(1,2), NODE, 1, U, Y
  *CFOPEN, CrackDispsFile, csv
  *VWRITE, CrackDisps(1,1), ',', CrackDisps(1,2)
  
  LSEL,S,,,7                                ! Crack
  NSLL,S,1                                  !Select nodes associated to this line
  *GET,NNodes,NODE,0,COUNT             !Get the number of nodes in the selected set
  *DIM, React, ARRAY, NNodes, 2     
  *VGET, React(1,1), NODE, 1, U, X
  *VGET, React(1,2), NODE, 1, U, Y
  *CFOPEN, ReactFileLow, csv
  *VWRITE, React(1,1), ',', React(1,2)
  
  LSEL,S,,,6                                ! Crack
  NSLL,S,1                                  !Select nodes associated to this line
  *GET,NNodes,NODE,0,COUNT             !Get the number of nodes in the selected set
  *DIM, React, ARRAY, NNodes, 2     
  *VGET, React(1,1), NODE, 1, U, X
  *VGET, React(1,2), NODE, 1, U, Y
  *CFOPEN, ReactFileUp, csv
  *VWRITE, React(1,1), ',', React(1,2)
  