anchor.diff:                                    elem_mass,elem_CG,elem_IN,ErrStat2,ErrMsg2)
             CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
   
  -WRITE(*,*) 'nelem:',nelem
  -WRITE(*,*) 'elem_mass:',elem_mass
  -WRITE(*,*) 'elem_CG:',elem_CG
  -DO j=1,3
  -WRITE(*,*) 'elem_IN:',elem_IN(j,:)
  -ENDDO
          blade_mass = blade_mass + elem_mass
          blade_CG(:) = blade_CG(:) + elem_CG(:)
          blade_IN(:,:) = blade_IN(:,:) + elem_IN(:,:)
   
      ENDDO
  -STOP
  +
      if (ErrStat >= AbortErrLev) then
          call Cleanup()
          return
  
positive.code:        mEta(2)      =-EMass0_GL(1,6,igp)
         mEta(3)      = EMass0_GL(1,5,igp)
         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
         CALL BD_GaussPointDataMass(hhx(:,igp),hpx(:,igp),temp_Nvvv,temp_Naaa,RR0,&
               node_elem,dof_node,vvv,aaa,vvp,mmm,mEta,rho,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         CALL BD_GravityForce(mmm,mEta,gravity,Fg,ErrStat2,ErrMsg2)
  
negative.code: !DO i=1,p%elem_total
  !WRITE(*,*) p%Jacobian(:,i)
  !ENDDO
  !STOP
     !CALL WrScr( "Finished reading input" )
     ! Allocate continuous states
     CALL AllocAry(x%q,p%dof_total,'x%q',ErrStat2,ErrMsg2)
  