anchor.diff:     CALL EigenSolve( Init%K, Init%M, Init%TDOF, NOmega, .True., Init, p, Phi, Omega, Init%RootName, ErrStat, ErrMsg )
      IF ( ErrStat /= ErrID_None ) RETURN
        
  +      ! Clean up Omega and Phi arrays because they are not needed
  +   DEALLOCATE(Omega)
  +   DEALLOCATE(Phi)
      
            ! Craig-Bampton reduction
   
  
positive.code:    
     ! normalize PhiM
     MU = MATMUL ( MATMUL( TRANSPOSE(PhiM), MLL ), PhiM )
     
     DO I = 1, DOFM
        MU(I, I) = 1./SQRT( MU(I, I) )  !RRD this is was wrong and I fixed it 6/10/2013
        OmegaM(I) = SQRT( OmegaM(I) )
     ENDDO
     
     PhiM = MATMUL( PhiM, MU )
  
     MBB = MRR + MATMUL(MRL, PhiR) + TRANSPOSE( MATMUL(MRL, PhiR) ) &
               + MATMUL( MATMUL(TRANSPOSE(PhiR), MLL), PhiR )
  
negative.code: !------------------------------------------------------------------------------------------------------
  !------------------------------------------------------------------------------------------------------
  
  End Module SubDyn
  