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:            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(:,:)
  
negative.code:        mmm  = 0.0D0
         mmm  = EMass0_GL(1,1,igp)
  
         elem_mass = elem_mass + gw(igp) * Jacobian * mmm
         elem_CG(1:3) = elem_CG(1:3) + gw(igp) * Jacobian * mmm * NGPpos(1:3,igp)
         elem_IN(1:3,1:3) = elem_IN(1:3,1:3) + gw(igp) * Jacobian * mmm * &
  