Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 5739)
+++ /issm/trunk/src/c/Makefile.am	(revision 5740)
@@ -150,6 +150,4 @@
 					./objects/Inputs/DoubleInput.h\
 					./objects/Inputs/DoubleInput.cpp\
-					./objects/Elements/BeamRef.h\
-					./objects/Elements/BeamRef.cpp\
 					./objects/Elements/Penta.h\
 					./objects/Elements/Penta.cpp\
@@ -703,6 +701,4 @@
 					./objects/Inputs/DoubleInput.h\
 					./objects/Inputs/DoubleInput.cpp\
-					./objects/Elements/BeamRef.h\
-					./objects/Elements/BeamRef.cpp\
 					./objects/Elements/Penta.h\
 					./objects/Elements/Penta.cpp\
Index: sm/trunk/src/c/objects/Elements/BeamRef.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/BeamRef.cpp	(revision 5739)
+++ 	(revision )
@@ -1,92 +1,0 @@
-/*!\file BeamRef.c
- * \brief: implementation of the BeamRef object
- */
-
-/*Headers:*/
-/*{{{1*/
-#ifdef HAVE_CONFIG_H
-	#include "config.h"
-#else
-#error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
-#endif
-
-#include "stdio.h"
-#include <string.h>
-#include "../objects.h"
-#include "../../EnumDefinitions/EnumDefinitions.h"
-#include "../../shared/shared.h"
-#include "../../Container/Container.h"
-#include "../../include/include.h"
-/*}}}*/
-
-/*Object constructors and destructor*/
-/*FUNCTION BeamRef::BeamRef(){{{1*/
-BeamRef::BeamRef(){
-	this->element_type_list=NULL;
-}
-/*}}}*/
-/*FUNCTION BeamRef::BeamRef(int* types,int nummodels){{{1*/
-BeamRef::BeamRef(const int nummodels){
-
-	/*Only allocate pointer*/
-	element_type_list=(int*)xmalloc(nummodels*sizeof(int));
-
-}
-/*}}}*/
-/*FUNCTION BeamRef::~BeamRef(){{{1*/
-BeamRef::~BeamRef(){
-	xfree((void**)&element_type_list);
-}
-/*}}}*/
-
-/*Management*/
-/*FUNCTION BeamRef::SetElementType{{{1*/
-void BeamRef::SetElementType(int type,int type_counter){
-
-	ISSMASSERT(type==P1Enum || type==P1DGEnum);
-
-	/*initialize element type*/
-	this->element_type_list[type_counter]=type;
-}
-/*}}}*/
-
-/*Reference Element numerics*/
-/*FUNCTION BeamRef::GetJacobianDeterminant {{{1*/
-void BeamRef::GetJacobianDeterminant(double*  Jdet, double* z_list,double gauss){
-
-	*Jdet=1.0/2.0*(z_list[1]-z_list[0]);
-	if(*Jdet<0) ISSMERROR(" negative jacobian determinant!");
-
-}
-/*}}}*/
-/*FUNCTION BeamRef::GetJacobianDeterminant2d {{{1*/
-void BeamRef::GetJacobianDeterminant2d(double*  Jdet, double* xyz_list,double gauss){
-
-	/*Jdet = 0.5 * length*/
-	*Jdet=1.0/2.0*(sqrt(pow(xyz_list[3*1+0]-xyz_list[3*0+0],2)+pow(xyz_list[3*1+1]-xyz_list[3*0+1],2)));
-	if(*Jdet<0) ISSMERROR(" negative jacobian determinant!");
-
-}
-/*}}}*/
-/*FUNCTION BeamRef::GetNodalFunctions {{{1*/
-void BeamRef::GetNodalFunctions(double* l1l2, double gauss){
-	/*This routine returns the values of the nodal functions  at the gaussian point.*/
-
-	l1l2[0]=-0.5*gauss+0.5;
-	l1l2[1]=+0.5*gauss+0.5;
-
-}
-/*}}}*/
-/*FUNCTION BeamRef::GetParameterValue{{{1*/
-void BeamRef::GetParameterValue(double* pvalue, double* value_list,double gauss_coord){
-
-	/*Nodal functions*/
-	double l1l2[2];
-
-	/*Get nodal functions*/
-	GetNodalFunctions(&l1l2[0],gauss_coord);
-
-	/*Assign output pointer*/
-	*pvalue=l1l2[0]*value_list[0]+l1l2[1]*value_list[1];
-}
-/*}}}*/
Index: sm/trunk/src/c/objects/Elements/BeamRef.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/BeamRef.h	(revision 5739)
+++ 	(revision )
@@ -1,31 +1,0 @@
-/*!\file:  BeamRef.h
- * \brief abstract class for handling Beam oriented routines, like nodal functions, 
- * strain rate generation, etc ...
- */ 
-
-
-#ifndef _BEAMREF_H_
-#define _BEAMREF_H_
-
-class BeamRef{
-	
-
-	public: 
-		int* element_type_list; //P1CG, P1DG, MINI, P2...
-		int  element_type;
-		
-		BeamRef();
-		BeamRef(const int nummodels);
-		~BeamRef();
-
-		/*Management*/
-		void SetElementType(int type,int type_counter);
-
-		/*Numerics*/
-		void GetJacobianDeterminant(double* Jdet, double* z_list,double gauss);
-		void GetJacobianDeterminant2d(double* Jdet, double* xyz_list,double gauss);
-		void GetNodalFunctions(double* l1l2, double gauss);
-		void GetParameterValue(double* pvalue, double* value_list,double gauss_coord);
-
-};
-#endif
Index: /issm/trunk/src/c/objects/objects.h
===================================================================
--- /issm/trunk/src/c/objects/objects.h	(revision 5739)
+++ /issm/trunk/src/c/objects/objects.h	(revision 5740)
@@ -36,5 +36,4 @@
 
 /*Elements: */
-#include "./Elements/BeamRef.h"
 #include "./Elements/Element.h"
 #include "./Elements/Penta.h"