%--------------------------------------------------------------------------
%Archivo de datos de las propiedades dinámicas de la estructura,
%provenientes de Spice
%--------------------------------------------------------------------------
%m=5;%número de modos
mat=[];


opcion=menu('Cargar archivos con las propiedades dinámicas de la barra SIN daño','Buscar archivo con las frecuencia y amortiguamiento',...
 'Salir');

%Cargar archivo con los vectores de frecuencia y amortiguación
if opcion==1
uiload
else %Salida del programa

disp('Gracias por usar "DIbEMA".')
break
end


opcion=menu('Cargar archivos con las propiedades dinámicas de la barra SIN daño','Buscar archivo con los modos reales',...
 'Salir');

%Cargar archivo con la matriz de modos
if opcion==1
uiload
else %Salida del programa

disp('Gracias por usar "DIbEMA".')
break
end




%Ordenar datos en vector columna para introducirlo en Dibema
for i=1:m 
mat(:,i)=[ i;freq_sel(i);damp_sel(i);mod_real(:,i)];   %nº modo; Frecuencia natural (Hz);Amortiguación,Desplazamiento de cada grado de libertad
end

[q,m]=size(mat);
matriz(1:q,1)=mat(:,1);

if m>1
for i=2:m
  matriz((i-1)*q+1:i*q,1)=mat(:,i);
end
end



matriz(:,2)=zeros(length(matriz(:,1)),1);

for i=1:m %m es el nº de modos
    matriz(((i-1)*n+3*i+1):i*(n+3),1)=matriz(((i-1)*n+3*i+1):i*(n+3),1)*...
        (1/sqrt((matriz(((i-1)*n+3*i+1):i*(n+3),1))'*(matriz(((i-1)*n+3*i+1):i*(n+3),1))));
end

