%
% Carga los archivos de resultados del modelo QG_barotrop.f
% y los parametros de QG_param.dat


dir_salida=[cd '/out_tmp/']

[param]=textread('QG_param.dat','%s%*[^\n]','headerlines',7);
for j=1:length(param), eval(param{j}); end


archivos=dir(dir_salida); % nombre de los archivos en el directorio actual + \output

k1=0; % indice de matriz funcion corriente
k2=0; % indice de matriz vorticidad

for i=3:length(archivos)
    if archivos(i).name(1:3)=='psi'
        k1=k1+1;
        psi_tmp(:,:,k1)=load([dir_salida archivos(i).name]);
    end
    if archivos(i).name(1:3)=='vor'
        k2=k2+1;
        vort_tmp(:,:,k2)=load([dir_salida archivos(i).name]);
    end
    if archivos(i).name(1:7)=='QG_diag'
        QG_diag=load([dir_salida archivos(i).name]);
    end
    if archivos(i).name(1:7)=='QG_wind'
        QG_curlw=load([dir_salida archivos(i).name]);
    end
end

%condicion de contorno:
vort_tmp(:,1,:)=vort_tmp(:,3,:);
vort_tmp(:,2,:)=vort_tmp(:,3,:);

vort_tmp(:,end,:)=vort_tmp(:,end-2,:);
vort_tmp(:,end-1,:)=vort_tmp(:,end-2,:);
%%%

vor_adim=vort_tmp;
psi_adim=psi_tmp;
vor_adimF=vor_adim(:,:,end); % salva ultimo time step en matriz vorF
psi_adimF=psi_adim(:,:,end);

clear archivos dir_salida i k1 k2 vort_tmp psi_tmp param