Changes between Version 15 and Version 16 of INCAbugs


Ignore:
Timestamp:
06/10/11 15:50:40 (13 years ago)
Author:
acosce
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • INCAbugs

    v15 v16  
    1717}}} 
    1818 
    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)__ == 
     20Les 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../).  
    2121  
    2222 
    2323== __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.  
     24l'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.  
    2525 
    2626== __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]] 
     27Depuis le passage à une parallélisation mixte mpi/openMP le calcul de la variable ptrop dans set_ub_vals était à zéro. [[BR]] 
    2828 ptrop=ptrop_glo(nbbeg_loc:nbend_loc) [[BR]] 
    2929commit 197 sur INCA3 
    3030 
    3131== __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 é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  
     32Dans 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  
     33O3I 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  
    3434commenter les lignes qui suivent. [[BR]] 
    3535Commit 200 sur INCA3 
    3636 
    3737 
    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)__ ==  
     39Problème apparu lors de l'optimisation du code pour la machine SX9. Cela entraine une erreur sur le diagnostique ec550aer. [[BR]] 
    4040{{{ 
    4141Il faut remplacer la ligne  
     
    5555 
    5656== __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éclaration  
     57Sur 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  
    5858de la variable '''mapping1''' 
    5959{{{ 
     
    7777 
    7878== __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 routine  
     79Dans 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.  
     80Il faut donc rajouter en début de routine  
    8181{{{ 
    8282 sunon(:)  = 0.  
     
    8989[[BR]] 
    9090[[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]] 
     91== __Config NMHC : EXP_SLV & IMP_SLV (05 décembre 2008)__ == 
     92Dans 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  
     93quelques erreurs se sont glissées dans leur ordre. [[BR]] 
    9494 
    9595__EXP_SLV.F90__[[BR]] 
    96 Les réactions c2h4 et c3h6 sont inversées.  
     96Les réactions c2h4 et c3h6 sont inversées.  
    9797{{{ 
    9898      + reaction_rates(:,177+compt)*base_sol(:,id_c3h6)             & 
     
    114114[[BR]] 
    115115__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és 
     116Pour 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 
    117117entre c2h4 et c3h6 
    118118{{{ 
     
    149149}}} 
    150150 
    151 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 
     151Pour la '''production de hno3''' il manque une équation (onitr + oh -> M.CO3 + 0.75.HNO3 + 0.25.NO2 +0.25.H2O) et 
     152un facteur 0.4 pour l'équation onitr + no3 -> M.CO3 +0.4.HNO3 + 0.8.NO2 + 0.5.NO 
    153153{{{ 
    154154             + reaction_rates(k,279+compt) * base_sol(k,id_onitu) * base_sol(k,id_oh) * 0.5   & 
     
    181181Ce calcul a disparu lors du changement de version : lecture de BC et POM dans un fichier texte et lecture dans le fichier sflx.nc 
    182182 
    183 Ce bug peut impliquer des différences de 30% et plus sur les variables AIN et ASN (et les varibles en dépendant)  
     183Ce bug peut impliquer des différences de 30% et plus sur les variables AIN et ASN (et les varibles en dépendant)  
    184184 
    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éfaut  
     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éfaut  
    186186[[BR]] 
    187187[[BR]]