1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
| temps_courant=;
chemin=;
N_maille=;
Nbre_ms=;
Nom_calcul = ['NOM_' num2str(N_maille) '_mailles'];
nom = [Nom_calcul '_' type '_' num2str(Nbre_ms) 'ms.dat'];
fid = fopen(fullfile(chemin, nom), 'w');
nom_temps =[Nom_calcul '_' type '_t_' num2str(Nbre_ms) 'ms.dat'];
fid_temps = fopen(fullfile(chemin, nom_temps), 'w');
while temps_courant<Nbre_ms*10^-3
for i=1:N_maille
%Calcul sans rapport avec mon soucis
end
for i=1:N_maille
%Calcul et écriture
Data = [ A(i) , B(i) , C(i) , D(i) , E(i) , F(i) , G(i) , H(i) , I(i) , J(i) , K(i) , L(1,i) , M(i), N(i), O(i), Q(i), R(i) ];
Separateur = '%5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f %5.4f \n';
fprintf(fid, Separateur, Data);
end
fprintf(fid, '\n SEP \n');
compteur = 0;
indice_a=indice_a+1;
%On enregistre le temps courant sur la derniere maille
fprintf(fid_temps, '%5.4e \n', temps_courant);
temps_courant = temps_courant+1;
end
fclose(fid);
fclose(fid_temps); |
Partager