clear
close all

dir_start=cd;

path='/usrtmp/fremy/';
fic='fla_essai1';

load([path fic])

cd(path)
doss=[fic '_avs'];
eval(['!rm -r ' doss])
mkdir(doss)

for im=1:size(depX,3)
    
    cd(dir_start)
    [h,b,d,g]=coup(depX,im);
    cd(path)
    
    depXc=depX(1+h:end-b,1+g:end-d,im);
    depYc=depY(1+h:end-b,1+g:end-d,im);
    depZc=depZ(1+h:end-b,1+g:end-d,im);
    Xc=X(1+h:end-b,1+g:end-d);
    Yc=Y(1+h:end-b,1+g:end-d);
    Zc=Z(1+h:end-b,1+g:end-d);

    largeur=size(Xc,2)-1;
    hauteur=size(Xc,1)-1;
    
    tabcon=sprintf('tabcon_%i.txt',im);
    cd(doss)
    eval(['!touch ' tabcon])
    fid = fopen(tabcon,'r+');
    fseek(fid,1,'bof');
    fprintf(fid,'%s\n',num2str([numel(Xc) largeur*hauteur 0 0 0]));
    for kk=1:numel(Xc)
    fprintf(fid,'%s\n',[num2str([kk Xc(kk) Yc(kk) Zc(kk)]) ' ']);
    end
    elem(1:largeur*hauteur,1)=1:largeur*hauteur;
    ii=0;
    col=1;
    for el=1:numel(elem)
        ii=ii+1;
        if ii>hauteur
        ii=1;
        col=col+1;
        end
        el;
        [(hauteur+1)*(col-1)+ii (hauteur+1)*(col-1)+ii+1 (hauteur+1)*col+ii+1 (hauteur+1)*col+ii];
        if el==133
           % return
        end
        fprintf(fid,'%s\n',[num2str(el) ' 0 quad ' num2str([(hauteur+1)*(col-1)+ii (hauteur+1)*(col-1)+ii+1 (hauteur+1)*col+ii+1 (hauteur+1)*col+ii]) ' ']);
    end
    fclose(fid) ;
    
    nomdep=sprintf('dep_%i.txt',im);
    donnees=[(1:numel(Xc))' depXc(:) depYc(:) depZc(:)];
    dlmwrite(nomdep, donnees, 'delimiter', '\t','precision', 6)
    
    eval(['!mv ' tabcon ' ' tabcon(1:end-4) '.avs'])
    eval(['!mv ' nomdep ' ' nomdep(1:end-4) '.avs'])
    
    cd ..
end

cd(dir_start)