anchor.diff:  
   tBPly = tRatio*(2*t) ! [mm] thickness of the bounding ply
   tTotal = t + tBPly   ! [mm] thickness of the bounding ply
  +
  +elSize = daOvera*a ! [mm] size of element in refined region close to crack tip
   
   ! Create Geometry
   
  
positive.code: tRatio = 1        ! [-]  ratio of bounding ply thickness to main ply
  atRatio = 0.1     ! [-]  ratio of crack length to main ply thickness
  rhon = 0.01       ! [-]  normalized crack density
  
  EL = ! [MPa] UD longitudinal Young's modulus
  ET = ! [MPa] UD transverse Young's modulus
  
negative.code: MSHKEY, 0
  AMESH, 3, 4, 1
  
  FINISH              ! Finish pre-processing
  
  /SOLU               ! Enter the solution processor
  