Pour une question de commodité, nous avons légèrement modifié la sous-routine de lecture des paramètres.  Maintenant, si on indique que l'on désire prolonger le vecteur d'onde (0,0), cela signifie en fait que nous voulons prolonger tous les vecteurs d'onde.  Auparavant, il était nécessaire de faire la liste explicite de tous les vecteurs d'onde à prolonger.  Bien sûr, pour prolonger tous les vecteurs d'onde, il est nécessaire qu'ils aient tous été calculés lors de la simulation MCQ!

Les modifications qui ont été apportées se trouvent dans la sous-routine LECTURE_PARA, tout juste après la boucle1, et sont contenues dans la condition
if(kx.eq.0.and.ky.eq.0) then...  Il s'agit tout simplement de lister dans l'ordre (1,1); (1,2); ...; (1,ny); (2,1); ...; (nx,ny) l'ensemble des vecteurs d'onde dans les tableaux qx(1,:) et qy(1,:).  Les tableaux nv et ipos servent au repérage des données et leur rôle précis est plutôt obscur (morale: toujours commenter ses programmes!).

Notez qu'auparavant, les données étaient présentées dans les fichiers de sortie ...xi.dat ou ...ps.dat étaient regroupées en 4 colonnes:

Fréquence     chi/w     chimax/w     Erreur

les coordonnées du vecteur d'onde prolongé étant indiquées quelques lignes plus haut.  À présent, ces coordonnées sont aussi présentées sous forme de colonnes, de sorte que les données sont plutôt présentées comme suit:

Coord. X     Coord. Y     Fréquence     chi/w     chimax/w     Erreur

L'objectif de cette modification était de simplifier l'écriture du programme de visualisation graf.f90.