anchor.diff:  
   O%Beddoes%FP   = O%Beddoes%FSP(J,IBLADE) - O%Beddoes%DF(J,IBLADE)
   O%Beddoes%FPC  = O%Beddoes%FSPC(J,IBLADE) - O%Beddoes%DFC(J,IBLADE)
  -SRFP = SQRT( ABS(O%Beddoes%FP) )  * SIGN( 1., O%Beddoes%FP ) + 1.
  -SRFPC= SQRT( ABS(O%Beddoes%FPC) ) * SIGN( 1.,O%Beddoes%FPC )
  +SRFP = SQRT( ABS(O%Beddoes%FP) )  * SIGN( 1.0_ReKi, O%Beddoes%FP ) + 1.
  +SRFPC= SQRT( ABS(O%Beddoes%FPC) ) * SIGN( 1.0_ReKi,O%Beddoes%FPC )
   
   O%Beddoes%FK   = 0.25 * SRFP * SRFP
   O%Beddoes%CN   = O%Beddoes%CNCP * O%Beddoes%FK + O%Beddoes%CNIQ
  
positive.code:       O%Beddoes%FSPC(J,IBLADE) = 1.0
     ELSE
        TEMP = O%Beddoes%FSPC(J,IBLADE)/((O%Beddoes%AFE(J,IBLADE)-AOL1)*O%Beddoes%AFE(J,IBLADE))
        O%Beddoes%FSPC(J,IBLADE) = TEMP * TEMP * SIGN ( 1., TEMP )
        IF ( O%Beddoes%FSPC(J,IBLADE) >  1.0 ) O%Beddoes%FSPC(J,IBLADE) =  1.0
        IF ( O%Beddoes%FSPC(J,IBLADE) < -1.0 ) O%Beddoes%FSPC(J,IBLADE) = -1.0
     ENDIF
  
  ENDIF
  
  SRFP = SQRT( ABS( O%Beddoes%FSP(J,IBLADE) ) ) * SIGN( 1., O%Beddoes%FSP(J,IBLADE) ) + 1.
  O%Beddoes%FK   = 0.25 * SRFP * SRFP
  O%Beddoes%CVN(J,IBLADE) = O%Beddoes%CNPOT(J,IBLADE) * ( 1. - O%Beddoes%FK )
  
  
negative.code:    O%Beddoes%FSPC(J,IBLADE) = 1.0
  ELSE
     TEMP = 2.*SQRT(ABS(O%Beddoes%FSP(J,IBLADE)/(O%Beddoes%AFE(J,IBLADE)-AOL1)))-1.
     O%Beddoes%FSP(J,IBLADE) = TEMP * TEMP * SIGN ( 1.0_ReKi, TEMP )
     IF ( O%Beddoes%FSP(J,IBLADE) >  1.0 ) O%Beddoes%FSP(J,IBLADE) =  1.0
     IF ( O%Beddoes%FSP(J,IBLADE) < -1.0 ) O%Beddoes%FSP(J,IBLADE) = -1.0
  
  