Table des matières

1- Pourquoi cette nouvelle version?
2- Comment utiliser cette nouvelle version?
3- Quelles ont été les modifications apportées au programme de prolongement analytique pour obtenir cette nouvelle version?

1- Pourquoi cette nouvelle version?

En principe, on devrait pouvoir effectuer le prolongement analytique d'une simulation pour laquelle les calculs n'ont été faits que pour un ou quelques vecteur(s) d'onde. Cependant, l'ancienne version supposait implicitement, lors de la lecture des données, que les calculs avaient été faits pour tous les vecteurs d'onde. La nouvelle version
du programme de prolongement analytique rectifie cette situation.
 

2- Comment utiliser cette nouvelle version?

La nouvelle version s'utilise exactement comme l'ancienne, à une différence près: il faut indiquer, à la fin du fichier trnf.dat, le nombre de vecteurs d'onde ayant été calculés lors de la simulation MCQ.
 

3- Quelles ont été les modifications apportées au programme de prolongement analytique pour obtenir cette nouvelle version?

A- Fichier trfn.dat:
    * A la fin du fichier, nous avons ajouté la variable NVCAL, qui indique le nombre de vecteurs d'onde ayant été calculés lors de la simulation MCQ.

B- Sous-routine LECTURE_PARA:
    * Ajout de la lecture de la variable NVCAL (variable déclarée dans le module "parametres")
    * Message d'erreur si NVCAL est plus petit que le nombre total de sites et que l'option "Matrice de covariance" a été choisie
    * Message d'erreur si NVCAL est plus petit que le nombre total de sites et que l'option "Densite d'etat totale" a été choisie

C- Sous-routine LECTURE_DATA:
    * Ajout de variables internes: TEST, GTL, ERRTL,TQX,TQY
    * Modification du nombre de lignes lues dans la fichier re.dat (la boucle sur J=1, 22+nx*ny devient une boucle sur J+1, 22+nvcal)
    * Remplacement de 2 boucles imbriquées (J=1,nx et I=1,ny) par une boucle (J=1,nvcal) pour la lecture de données inutiles
    * Au lieu d'enregistrer les données dans les tableaux GT et ERRT, on les enregistre dans les tableaux GTL et ERRTL.  De plus, les indices des vecteurs d'onde sont conservés dans les tableaux TQX et TQY.
    * On transfert les données de GTL à GT et de ERRTL à ERRT.  Pour ce faire, on vérifie si les indices de vecteurs d'onde des tableaux GT et ERRT sont les mêmes que ceux des tableaux TQX et TQY; si c'est le cas, on transfert les résultats contenus dans GTL et ERRTL, sinon, on entre la valeur 0.0 dans GT et ERRT.
    * Il a fallu mettre en option l'addition des différentes contributions pour les réseaux carrés (symétrie), car il n'y a plus nécessairement de données pour tous les sites.

Autrement dit, les tableaux GT et ERRT contiennent maintenant les résultats MCQ pour les vecteurs d'onde qui ont été calculés et 0.0 pour les vecteurs d'onde qui n'ont pas été calculés.  Il n'a pas été nécessaire d'apporter d'autres modifications pour les calculs subséquents.