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:        end if
  
     ENDDO
  STOP
     blade_CG(:) = blade_CG(:) / blade_mass
  
     CALL Cleanup()
  
negative.code:            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,&
  