anchor.diff:             ENDDO
          ENDIF
   
  +WRITE(*,*) 'Trapezoidal_pos'
  +WRITE(*,*) trapezoidal_pos
  +WRITE(*,*) 'Trapezoidal_w'
  +WRITE(*,*) trapezoidal_w
  +STOP
  +
          CALL BD_ElementMatrixAcc(Nuu0,Nuuu,Nrr0,Nrrr,Nvvv,&
                                   EStif0_GL,EMass0_GL,gravity,DistrLoad_GL,&
                                   nqp,quadrature,trapezoidal_pos,trapezoidal_w,&
  
positive.code:            ENDDO
         ENDIF
  
  WRITE(*,*) 'nelem:',nelem
         CALL BD_ComputeElementMass(Nuu0,NGPpos,EMass0_GL,&
                                    nqp,quadrature,trapezoidal_pos,trapezoidal_w,&
                                    node_elem,dof_node,&
  
negative.code:            DEALLOCATE(trapezoidal_w)
         ENDIF
  
         blade_mass = blade_mass + elem_mass
         blade_CG(:) = blade_CG(:) + elem_CG(:)
         blade_IN(:,:) = blade_IN(:,:) + elem_IN(:,:)
  