"snow" dans "phyetat0"

[guez@…]

Dans "phyetat0.F", version 1.26, on a a la ligne 350 :

CALL get_field("SNOW",snow(:,nsrf),found)

"nsrf" vaut "nbsrf+1" (sortie d'une boucle précédente sur "nsrf" jusqu'à "nbsrf"). La taille de "snow" dans la seconde dimension est "nbsrf". On a donc une erreur à l'exécution de "gcm" avec les options de débogage :

• 240 Subscript error array=snow size=2952 subscript=2953 eln=1872 PROG=phyetat0 ELN=1872(400e9a488)

Est-ce qu'on ne voulait pas écrire en fait :

CALL get_field("SNOW",snow,found)

Lionel Guez

[anonymous]

En fait je suis tombé sur le même problème le même jour. C'est une malheureuse coquille que j'ai introduite en mettant à jour phyetat0. En fait il faut comprendre :

CALL get_field("SNOW",snow(:,1),found)

Normalement le tableau snow n'est pas utilisé puisque snow n'existe pas dans le fichier de démarrage. On y trouve par contre les valeurs pour les sous-surfaces. Effectivement sans les options de débogage on passe à coté.

Yann

Replying to guez@lmd.ens.fr:

Dans "phyetat0.F", version 1.26, on a à la ligne 350 :

CALL get_field("SNOW",snow(:,nsrf),found)

"nsrf" vaut "nbsrf+1" (sortie d'une boucle précédente sur "nsrf" jusqu'à "nbsrf"). La taille de "snow" dans la seconde dimension est "nbsrf". On a donc une erreur à l'exécution de "gcm" avec les options de débogage :

• 240 Subscript error array=snow size=2952 subscript=2953 eln=1872 PROG=phyetat0 ELN=1872(400e9a488)

Est-ce qu'on ne voulait pas écrire en fait :

CALL get_field("SNOW",snow,found)

Lionel Guez

[lguez]

Corrected. "snow(:,nsrf)" has been replaced by "snow(:,1)".