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:        temp_id1 = (igp-1)*dof_node
         temp_id2 = (igp-1)*dof_node/2
         Stif(1:6,1:6) = EStif0_GL(1:6,1:6,igp)
         CALL BD_GaussPointData(hhx(:,igp),hpx(:,igp),Nuuu,Nrrr,&
               uu0(temp_id1+1:temp_id1+6),E10(temp_id2+1:temp_id2+3),node_elem,dof_node,&
               uuu,uup,E1,RR0,kapa,Stif,cet,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code: 
             Jacobian(j,i) = SQRT(DOT_PRODUCT(Gup0,Gup0))
  
         ENDDO
     ENDDO
  
  