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:       IF ( ( ABS( SPhi ) > 0.01 ) .AND. ( ABS( CPhi ) > 0.01 ) )  THEN
           A2P = SOLFACT*( CLA*SPhi - CDA*CPhi )*( 1.0 + A2P )*VNROTOR2/( 4.0*O%LOSS*SPhi*CPhi )
        ELSEIF ( ABS( SPhi ) > 0.01 )  THEN   ! Tangential velocity near zero, phi near 90 degrees.
           A2P = SOLFACT*( CLA*SPhi - CDA*SIGN( 0.01, CPhi ) )*( 1.0 + A2P )*VNROTOR2/( 4.0*O%LOSS*SPhi*SIGN( 0.01, CPhi ) )
        ELSE   ! Normal velocity near zero, phi near 0 degrees.
           A2P = SOLFACT*( CLA*SIGN( 0.01, SPhi ) - CDA*CPhi )*( 1.0 + A2P )*VNROTOR2/( 4.0*O%LOSS*SIGN( 0.01, SPhi )*CPhi )
        ENDIF
     ELSE
        SWRLARG = 1.0 + 4.0*O%LOSS*A2*VNW*VNA*VT2_Inv
  
negative.code: 
  O%Beddoes%CMQ = -.25 * CNQ - (XKA*CO/3.) * (O%Beddoes%QX(J,IBLADE) - O%Beddoes%DQP(J,IBLADE))
  
  O%Beddoes%CNIQ = MIN( ABS( CNI+CNQ ), 1. ) * SIGN( 1.0_ReKi, CNI+CNQ )
  
  O%Beddoes%XN(J,IBLADE) = O%Beddoes%OLDXN(J,IBLADE)*EXP(-.14*BS) + .3*DA*EXP(-.07*BS)
  O%Beddoes%YN(J,IBLADE) = O%Beddoes%OLDYN(J,IBLADE)*EXP(-.53*BS) + .7*DA*EXP(-.265*BS)
  