Discussion :

[Butterfly83]

Bonjour,
alors j aimerais avoir un peu d aide, sachant que je suis novice en Fortran. Alors je raisous un système et suivant que ma méthode converge ou pas je continu ou je veux diviser mon pas de temps pour adapter, en gros, mon maillage.
J ai pensé à faire un "goto" mais a priori ca ne marche pas, sachant que mon pas de temps et toutes les interpolation qui en découle est définit dans un autre fichier.

Si quelqu un a une petite idée ce serrait sympa.
Merci.

[genteur slayer]

fait montrer comment sont imbriqué (en gros) tes boucles, sinon tu a la boucle DO avec condition final et la commande EXIT:

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
 
Boucletemporelle:do while(t<TempsFinal)
  Boucle_pas_de_temps:Do
    <.................>
    <tout ton calcul + calcul de l'erreur (servant à adapter le pas)
    <.................>
    if (erreur<seuil)then 
      t=t+dt
      c=(Constante/erreur)**(1.d0/5.d0)   !facteur d"adaptation du pas          
      if (c>2.d0) c=2.d0
      if (c<5.d-1) c=5.d-1 
      dt=c*dt
      EXIT
    else
      dt=dt/2.d0
    endif
  end do Boucle_pas_de_temps
end do boucle temporelle

[Butterfly83]

En gros ce que j ai c est:

Boucle de temps:
subroutine(1) qui fait tout les calculs
subroutine(2) qui dit si ca converge
Fin boucle de temps

et je voudrais que quand (2) me dit que ca converge pas retourner au début de la boucle temps.

[Butterfly83]

?le facteur d adaptation c est quoi?

[Butterfly83]

la subroutine (2) est dans la (1), c est très simplifié mais en gros c est ce qui ce passe.

[Butterfly83]

Mais le schéma que tu m as donné ressemble à ce que j ai dans la tête et ce que je veux programmer.