Custom Query (96 matches)

Création d'une nouvelle testing basée sur la r4725 ​modipsl.20231011.trunk.tar.gz

Nouveautés

• continuing work on rewriting: suppression of preproccesing keys (XIOS, MPI)
• replayisation: wake, thermals, ratqs, cloudth, lscp_mod, vdif_kcay, lmdz_lscp_old, cloud_optics_prop (alias newmicro)
• code from the TKE workshop: prognostic equation for tke
• work on StratAer
• offline mode back to work
• double call of the ECRAD radiative transfer code
• modification to the cdrag routines to include different bulk formulae to calculate cdrag coefficients over ocean
• bug fixes

Tests

Prendre physiq.def_NPv6.3 par défaut

contrôle qualité de base	validé
tests "multi-physiques" (Lionel)	validé
tutorial.tar / tutorial_prod.tar (Adriana)	validé
1D (Etienne)	validé
Config IPSLCM7: LMDZOR / ICOLMDZOR (Laurent)	validé/validé
trusting LMDZOR (Laurent / Lidia)	validé
trusting IPSLCM6.2 (Laurent / Lidia)	validé
ICOLMDZOR_v7_LAM.1 (Etienne / Laurent?)	validé
HighTune (Fredho)	validé
dynzon	OK sur Jean-Zay tutorial std, mais pb en aqua planete
phydev (Fredho ?)	Commission 4729 sur phydev. validé
phyparam (Fredho)	validé
phyex (Fredho)	validé

Remarques / Détails

• 1d validé avec perte de convergence numérique par rapport aux versions précédentes due à la correction de r4650

Ce problème est présent dans la révision 1140. Dans "phylmd/phys_state_var_mod.F90", "zstd" est alloué avec la taille "klon". Dans "dyn3d/etat0_netcdf.F", à la ligne 461, pour remplir "rugsrel", "zstd(i)" apparaît dans une boucle de 1 à "iip1*jjp1". Il y a donc dépassement de la plage d'indices permis pour "zstd". L'erreur se manifeste dans l'exécution de "create_etat0_limit" en séquentiel sur Brodie avec les options de débogage :

• 240 Subscript error array=zstd size=738 subscript=739 eln=2113 PROG=etat0_netcdf ELN=2113(400016914)

99 Execution suspended PROG=etat0_netcdf ELN=2113(400016914)

Called from create_etat0_limit ELN=183(400001828)

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