Discussion :

[Altolino]

Bonjour,

J'ai un script en ascii dans lequel j'aimerai changer le fichier d'entrée et la sortie (en orange) :

$settings                                                   
   MODEL="../love/love_load.txt"
   MODEL2="../love/love_tidal.txt"
   MASKFILE="../topo/mask-180.xyz"
  ICEFILE="../../1_input/ADIPAAO001"
  OUTPUT="../../4_SLM_OUTPUT/ADIPAAO001-OUT"
 LM=180                                                      
 LMM=180                                                                                         
 EXTENT=15.d0                                                
 T0=0.d0                                                     
 TREF=0.d0                                                   
 TLGM=1.d0                                                   
 TEOG=0.d0                                                   
 TLC=2.d0                                                    
 NPC=0                                                       
 PCINT=1.d0                                                  
 DLT=1.d0                                                    
 EPS=1.d-2                                                   
 MAXITER=4                                                   
 METHOD=3                                                    
 $end

Je dois répéter cette opération avec 780 fichiers c'est pourquoi je me demande s'il n'est pas possible d'automatiser ça sur Matlab, sachant que la numérotation des fichiers est toujours la même : de ADIPAAO001 à ADIPAAO780.

Merci de vos réponses

[Jerome Briot]

Une idée :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
41
42
fid = fopen('input.txt', 'w');
fprintf(fid, '$settings\n');
fprintf(fid, 'MODEL="../love/love_load.txt"\n');
fprintf(fid, 'MODEL2="../love/love_tidal.txt"\n');
fprintf(fid, 'MASKFILE="../topo/mask-180.xyz"\n');
fprintf(fid, 'ICEFILE="../../1_input/ADIPAAO001"\n');
fprintf(fid, 'OUTPUT="../../4_SLM_OUTPUT/ADIPAAO001-OUT"\n');
fprintf(fid, 'LM=180\n');
fprintf(fid, 'LMM=180\n');
fprintf(fid, 'EXTENT=15.d0\n');
fprintf(fid, 'T0=0.d0\n');
fprintf(fid, 'TREF=0.d0\n');
fprintf(fid, ' TLGM=1.d0\n');
fprintf(fid, 'TEOG=0.d0\n');
fprintf(fid, 'TLC=2.d0\n');
fprintf(fid, 'NPC=0\n');
fprintf(fid, 'PCINT=1.d0\n');
fprintf(fid, 'DLT=1.d0\n');
fprintf(fid, 'EPS=1.d-2\n');
fprintf(fid, 'MAXITER=4\n');
fprintf(fid, 'METHOD=3\n');
fprintf(fid, '$end\n');
fclose(fid);
 
for n = 1:780
 
    fid = fopen('input.txt', 'r');
    X = textscan(fid, '%s', 'delimiter', '\n');
    fclose(fid);
 
    filename = sprintf('ADIPAAO%03d', n);
 
    X{1}{5} = sprintf('ICEFILE="../../1_input/%s"', filename);
    X{1}{6} = sprintf('OUTPUT="../../4_SLM_OUTPUT/%s-OUT"', filename);
 
    fid = fopen('input.txt', 'w');
    fprintf(fid, '%s\n', X{1}{:});
    fclose(fid);
 
    % Utilisation du fichier ici
 
end

[Altolino]

Merci ! ça à l'air de marcher, il faut juste que je trouve le moyen créer un fichier à chaque itération afin que chaque script soit gardé (ce qui est le but de la manœuvre ).

[Jerome Briot]

Voici une version plus simple si tu veux conserver les fichiers d’entrée :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
for n = 1:780
 
    filename = sprintf('input_%03d.txt', n);
 
    fid = fopen(filename, 'w');
    fprintf(fid, '$settings\n');
    fprintf(fid, 'MODEL="../love/love_load.txt"\n');
    fprintf(fid, 'MODEL2="../love/love_tidal.txt"\n');
    fprintf(fid, 'MASKFILE="../topo/mask-180.xyz"\n');
    fprintf(fid, 'ICEFILE="../../1_input/ADIPAAO%03d"\n', n);
    fprintf(fid, 'OUTPUT="../../4_SLM_OUTPUT/ADIPAAO%03d-OUT"\n', n);
    fprintf(fid, 'LM=180\n');
    fprintf(fid, 'LMM=180\n');
    fprintf(fid, 'EXTENT=15.d0\n');
    fprintf(fid, 'T0=0.d0\n');
    fprintf(fid, 'TREF=0.d0\n');
    fprintf(fid, ' TLGM=1.d0\n');
    fprintf(fid, 'TEOG=0.d0\n');
    fprintf(fid, 'TLC=2.d0\n');
    fprintf(fid, 'NPC=0\n');
    fprintf(fid, 'PCINT=1.d0\n');
    fprintf(fid, 'DLT=1.d0\n');
    fprintf(fid, 'EPS=1.d-2\n');
    fprintf(fid, 'MAXITER=4\n');
    fprintf(fid, 'METHOD=3\n');
    fprintf(fid, '$end\n');
    fclose(fid);
 
 
    % Utilisation du fichier ici
 
end

[Altolino]

Merci beaucoup ! C'est impeccable !

[Jerome Briot]

Une version plus compacte et plus rapide :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
for n = 1:780
 
    filename = sprintf('input_%03d.txt', n);
 
    fid = fopen(filename, 'w');
    fprintf(fid, ...
        '$settings\nMODEL="../love/love_load.txt"\nMODEL2="../love/love_tidal.txt"\nMASKFILE="../topo/mask-180.xyz"\nICEFILE="../../1_input/ADIPAAO%03d"\nOUTPUT="../../4_SLM_OUTPUT/ADIPAAO%03d-OUT"\nLM=180\nLMM=180\nEXTENT=15.d0\nT0=0.d0\nTREF=0.d0\nTLGM=1.d0\nTEOG=0.d0\nTLC=2.d0\nNPC=0\nPCINT=1.d0\nDLT=1.d0\nEPS=1.d-2\nMAXITER=4\nMETHOD=3\n$end\n', ...
        n, n);
    fclose(fid);
 
    % Utilisation du fichier ici
 
end