anchor.diff:             DEALLOCATE(trapezoidal_w)
          ENDIF
   
  +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(:,:)
  
positive.code:        end if
  
     ENDDO
  
     blade_CG(:) = blade_CG(:) / blade_mass
  
     CALL Cleanup()
  
negative.code:            ENDDO
         ENDIF
  
  WRITE(*,*) 'nelem:',nelem
         CALL BD_ComputeElementMass(Nuu0,NGPpos,EMass0_GL,&
                                    nqp,quadrature,trapezoidal_pos,trapezoidal_w,&
                                    node_elem,dof_node,&
  