anchor.diff:        InitInp%Morison%MGBottom =  999999.0
         do I = 1,InitInp%Morison%NMGDepths
               ! Adjust the depth values based on MSL2SWL
  -         InitInp%Morison%MGDepths(I)%MGDpth = InitInp%Morison%MGDepths(I)%MGDpth - InitInp%MSL2SWL
  +         InitInp%Morison%MGDepths(I)%MGDpth = InitInp%Morison%MGDepths(I)%MGDpth - InitInp%Morison%MSL2SWL
         end do
         DO I = 1,InitInp%Morison%NMGDepths
               ! Store the boundaries of the marine growth zone
  
positive.code:    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
  
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
  
  