anchor.diff:     IF ( InitInp%Morison%NAxCoefs > 0 ) THEN
      
         DO I = 1,InitInp%Morison%NAxCoefs 
  -         
  -         !IF (  .NOT. EqualRealNos(InitInp%Morison%AxialCoefs(I)%AxCd, 0.0) ) THEN
  -         !   ErrMsg  = ' AxCd must be equal to zero.  Future versions will allow for non-zero axial coefficients.'
  -         !   ErrStat = ErrID_Fatal
  -         !   RETURN
  -         !END IF   
  -         !IF (  .NOT. EqualRealNos(InitInp%Morison%AxialCoefs(I)%AxCa, 0.0) ) THEN
  -         !   ErrMsg  = ' AxCa must be equal to zero.  Future versions will allow for non-zero axial coefficients.'
  -         !   ErrStat = ErrID_Fatal
  -         !   RETURN
  -         !END IF   
  -         
  -         ! TODO: Once Axial Coefs are working remove the above checks and uncomment the checks below.  GJH 9/29/2013
  +
            IF (  InitInp%Morison%AxialCoefs(I)%AxCd < 0 ) THEN
               CALL SetErrStat( ErrID_Fatal,'AxCd must be greater or equal to zero.',ErrStat,ErrMsg,RoutineName)
               RETURN
  
positive.code:    ! Joint Output List Section
     !-------------------------------------------------------------------------------------------------
  
     !IF ( InitInp%Morison%NJOutputs /= 0 ) THEN  ! TODO Remove this check and add back the other checks once Joint Outputs are supported
     !CALL SetErrStat( ErrID_Fatal,'NJOutputs in the Joint output list must be equal to zero.  Future versions of HydroDyn will support values greater or equal to zero and less than 10.'
     !   ErrStat = ErrID_Fatal
     !   RETURN
     !END IF   
     
     
     IF ( ( InitInp%Morison%NJOutputs < 0 ) .OR. ( InitInp%Morison%NMOutputs > 9 ) ) THEN
        CALL SetErrStat( ErrID_Fatal,'NJOutputs in the Joint output list must be greater or equal to zero and less than 10.',ErrStat,ErrMsg,RoutineName)
        RETURN
  
negative.code: 
        ! Check JointOvrlp values
    !NOTE: This is ignored in the current version of Morison.  3/15/2020 GJH
     
     IF ( InitInp%Morison%NJoints > 1 ) THEN
  
  