anchor.diff:             ENDDO
          ENDIF
   
  +WRITE(*,*) 'nelem:',nelem
          CALL BD_ComputeElementMass(Nuu0,NGPpos,EMass0_GL,&
                                     nqp,quadrature,trapezoidal_pos,trapezoidal_w,&
                                     node_elem,dof_node,&
  
positive.code:            DEALLOCATE(trapezoidal_w)
         ENDIF
  
         blade_mass = blade_mass + elem_mass
         blade_CG(:) = blade_CG(:) + elem_CG(:)
         blade_IN(:,:) = blade_IN(:,:) + elem_IN(:,:)
  
negative.code:        end if
  
     ENDDO
  STOP
     blade_CG(:) = blade_CG(:) / blade_mass
  
     CALL Cleanup()
  