anchor.diff:    hammerHeight = initialHammerHeight
     hammerVelocity = initialHammerVelocity
     hammerForce = 0
  +  stringForce = 0
     hammerDone = .FALSE.
   
     !!! Constants
  
positive.code:     end do
  
      if (hammerDone.eqv..FALSE.) then
        call rigid_hammer_strike(y(:,localt),hammerHeight,hammerMass,hammerVelocity,hammerForce, &
           lastHammerVelocity,compression,deltaT,g,stringForce,hammerDone,tension,bendingPrefactor, &
           k,b,rho,deltaX)
      end if
      
  
  
negative.code: !!!!!!!!!!!!!!!!! localt having a constant value of 3 (which
  !!!!!!!!!!!!!!!!! helps continuity as well.)
    F = 0
    vy = 0
  
      y(1,localt) = bL1*initY(1)+bL2*initY(2)+bL3*initY(3)+bL4*initY(1)+bLF*0
      y(2,localt) = a1*(initY(4)-initY(2)+2*initY(1)) + a2*(initY(3)+initY(1)) &
  