Changes between Version 15 and Version 16 of INCAbugs
- Timestamp:
- 06/10/11 15:50:40 (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
INCAbugs
v15 v16 17 17 }}} 18 18 19 == __fichiers de vents Ã10m en 320x160 (19 janvier)__ ==20 Les fichiers de vents à 10m 165_...nc, 166_....nc ou u10...nc, v10...nc pour les années 2007 2008 2009 2010 n'étaient pas au bon format (short au lieu de double ou float). Cela nâentraînait pas un bug dans la simulation, mais celle-ci lisait des valeurs absurdes. Ces fichiers ont été modifié. Les anciens sont systématiquement sauvegardés dans des répertoires 165_166_short ou u10_v10_short (sous p24data/ECMWF320x160/AN../).19 == __fichiers de vents à 10m en 320x160 (19 janvier)__ == 20 Les fichiers de vents à 10m 165_...nc, 166_....nc ou u10...nc, v10...nc pour les années 2007 2008 2009 2010 n'étaient pas au bon format (short au lieu de double ou float). Cela n'entraînait pas un bug dans la simulation, mais celle-ci lisait des valeurs absurdes. Ces fichiers ont été modifié. Les anciens sont systématiquement sauvegardés dans des répertoires 165_166_short ou u10_v10_short (sous p24data/ECMWF320x160/AN../). 21 21 22 22 23 23 == __utilisation de Adjust (17 janvier 2011)__ == 24 l'option adjust introduit des erreurs dans le code. Il faut lancer une pr é-simulation avec l'option adjust et ensuite la mettre à zéro et pointer vers le fichier !BandResol de la pré-simulation.24 l'option adjust introduit des erreurs dans le code. Il faut lancer une pré-simulation avec l'option adjust et ensuite la mettre à zéro et pointer vers le fichier !BandResol de la pré-simulation. 25 25 26 26 == __ptrop (17 janvier 2011)__ == 27 Depuis le passage à une parallélisation mixte mpi/openMP le calcul de la variable ptrop dans set_ub_vals était à zéro. [[BR]]27 Depuis le passage à une parallélisation mixte mpi/openMP le calcul de la variable ptrop dans set_ub_vals était à zéro. [[BR]] 28 28 ptrop=ptrop_glo(nbbeg_loc:nbend_loc) [[BR]] 29 29 commit 197 sur INCA3 30 30 31 31 == __Exp_slv (17 janvier 2011)__ == 32 Dans toutes les versions de exp_slv on arr ête le calcul de base_sol à 2 espÚces de la fin. Les deux derniÚres sont systématiquement considérées comme étant33 O3I et O3S. Or ces esp Úces n'existent pas dans les configurations AER et GES. Dans ce cas lÃil vous faut modifier le calcul de base_sol (enlever le "-2") et32 Dans toutes les versions de exp_slv on arrête le calcul de base_sol à 2 espèces de la fin. Les deux dernières sont systématiquement considérées comme étant 33 O3I et O3S. Or ces espèces n'existent pas dans les configurations AER et GES. Dans ce cas là il vous faut modifier le calcul de base_sol (enlever le "-2") et 34 34 commenter les lignes qui suivent. [[BR]] 35 35 Commit 200 sur INCA3 36 36 37 37 38 == __aeropt_5wv.F90 (16 d écembre 2010)__ ==39 Probl Úme apparu lors de l'optimisation du code pour la machine SX9. Cela entraine une erreur sur le diagnostique ec550aer. [[BR]]38 == __aeropt_5wv.F90 (16 décembre 2010)__ == 39 Problème apparu lors de l'optimisation du code pour la machine SX9. Cela entraine une erreur sur le diagnostique ec550aer. [[BR]] 40 40 {{{ 41 41 Il faut remplacer la ligne … … 55 55 56 56 == __SETHET (30 janvier 2009)__ == 57 Sur la machine platine du CCRT il y a un '''probl Úme de compilation''' dans la routine sethet du modÚle. Il s'agit d'enlever un espace lors de la déclaration57 Sur la machine platine du CCRT il y a un '''problème de compilation''' dans la routine sethet du modèle. Il s'agit d'enlever un espace lors de la déclaration 58 58 de la variable '''mapping1''' 59 59 {{{ … … 77 77 78 78 == __DIURNAL_GEOM (14 janvier 2009)__ == 79 Dans diurnal_geom les variables sunon et sunoff ne sont pas initialis ées ce qui peut poser problÚme car leurs calculs dépendent de conditions IF qui ne sont pas toujours vérifiées.80 Il faut donc rajouter en d ébut de routine79 Dans diurnal_geom les variables sunon et sunoff ne sont pas initialisées ce qui peut poser problème car leurs calculs dépendent de conditions IF qui ne sont pas toujours vérifiées. 80 Il faut donc rajouter en début de routine 81 81 {{{ 82 82 sunon(:) = 0. … … 89 89 [[BR]] 90 90 [[BR]] 91 == __Config NMHC : EXP_SLV & IMP_SLV (05 d écembre 2008)__ ==92 Dans les routines exp_sol, imp_sol les r éactions sont représentées par leurs numéros. Cependant ces numéros sont indiqués en dur dans le code et93 quelques erreurs se sont gliss ées dans leur ordre. [[BR]]91 == __Config NMHC : EXP_SLV & IMP_SLV (05 décembre 2008)__ == 92 Dans les routines exp_sol, imp_sol les réactions sont représentées par leurs numéros. Cependant ces numéros sont indiqués en dur dans le code et 93 quelques erreurs se sont glissées dans leur ordre. [[BR]] 94 94 95 95 __EXP_SLV.F90__[[BR]] 96 Les r éactions c2h4 et c3h6 sont inversées.96 Les réactions c2h4 et c3h6 sont inversées. 97 97 {{{ 98 98 + reaction_rates(:,177+compt)*base_sol(:,id_c3h6) & … … 114 114 [[BR]] 115 115 __IMP_SLV.F90__ [[BR]] 116 Pour les '''pertes de O3''' la r éaction O3+c2h2 n'existe pas, il faut la modifier par o3+c2h4. De plus de nouveau les numéros des réactions sont inversés116 Pour les '''pertes de O3''' la réaction O3+c2h2 n'existe pas, il faut la modifier par o3+c2h4. De plus de nouveau les numéros des réactions sont inversés 117 117 entre c2h4 et c3h6 118 118 {{{ … … 149 149 }}} 150 150 151 Pour la '''production de hno3''' il manque une équation (onitr + oh -> M.CO3 + 0.75.HNO3 + 0.25.NO2 +0.25.H2O) et152 un facteur 0.4 pour l' équation onitr + no3 -> M.CO3 +0.4.HNO3 + 0.8.NO2 + 0.5.NO151 Pour la '''production de hno3''' il manque une équation (onitr + oh -> M.CO3 + 0.75.HNO3 + 0.25.NO2 +0.25.H2O) et 152 un facteur 0.4 pour l'équation onitr + no3 -> M.CO3 +0.4.HNO3 + 0.8.NO2 + 0.5.NO 153 153 {{{ 154 154 + reaction_rates(k,279+compt) * base_sol(k,id_onitu) * base_sol(k,id_oh) * 0.5 & … … 181 181 Ce calcul a disparu lors du changement de version : lecture de BC et POM dans un fichier texte et lecture dans le fichier sflx.nc 182 182 183 Ce bug peut impliquer des diff érences de 30% et plus sur les variables AIN et ASN (et les varibles en dépendant)183 Ce bug peut impliquer des différences de 30% et plus sur les variables AIN et ASN (et les varibles en dépendant) 184 184 185 '''ATTENTION''' : dans certaines versions ces variables sont calcul ées à l'intérieur d'une boucle "IF (.not. bcpom_through_sflx) " dans laquelle on n'entre pas par défaut185 '''ATTENTION''' : dans certaines versions ces variables sont calculées à l'intérieur d'une boucle "IF (.not. bcpom_through_sflx) " dans laquelle on n'entre pas par défaut 186 186 [[BR]] 187 187 [[BR]]