anchor.diff:  
   L = t/rhon        ! [mm] length of the RVE
   a = aLRatio*L     ! [mm] 2a = crack length
  +
  +refHlow = (1-reftRatio)*t 
  +refHup = (1+reftRatio)*t
   
   tBPly = tRatio*(2*t) ! [mm] thickness of the bounding ply
   tTOT = t + tBPly   ! [mm] thickness of the bounding ply
  
positive.code: daOvera = 0.05    ! [-]  ratio of crack increment (i.e. crack tip element size) to crack length
  epsx = 0.01       ! [-]  applied strain
  uniP = 0.0        ! [-]  uniform pressure applied to crack face
  
  nContours = 10 ! [-]  number of contours for J-integral evaluation
  
  
negative.code: ! Create Geometry
  
  ! Points
  
  K, 1, 0.0, 0.0
  K, 2, L, 0.0
  K, 3, L, tTOT
  K, 4, 0.0, tTOT
  
  K, 5, 0.0, refHlow
  K, 6, L, refHlow
  K, 7, L, refHup
  K, 8, 0.0, refHup
  