| 74 | |
| 75 | === 25/03/2011 Utilisation du script de correction des variables IPCC nbp et fLuc === |
| 76 | |
| 77 | Ce script est à installer dans modipsl/libIGCM et nécessite de faire un ins_job sur la machine de post-traitement. |
| 78 | |
| 79 | * Ce qui a été fait sur cesium pour la simulation piControl2 : |
| 80 | {{{ |
| 81 | cd modipsl ; svn co http://forge.ipsl.jussieu.fr/igcmg/svn/TOOLS TOOLS ; cp AA_stomate_flux_correction libIGCM/. |
| 82 | cd util ; ./ins_job ; |
| 83 | cd ../libIGCM ; |
| 84 | vi stomate_flux_correction.job (pour lire la partie préciser le repertoire SUBMIT_DIR dans lequel se trouve le config.card correspondant à la simulation à traiter) |
| 85 | ccc_msub stomate_flux_correction.job |
| 86 | }}} |
| 87 | * En cas de plantage il faut ressoumettre après les modifications suivantes : |
| 88 | {{{ |
| 89 | vi stomate_flux_correction.job |
| 90 | # Mise à jour des variables CURRENT_LOOP, imin, imax. Prendre les valeurs de ces variables dans le dernier output SBGcorr.* |
| 91 | # Par défaut, le job retraitera tous les fichiers depuis le début. Plus d'infos sur ces variables dans le job lui-même. |
| 92 | }}} |
| 93 | * A noter : le job sauve au fur et à mesure les fichiers stomate_history_ipcc.nc originaux dans des repertoires SAVE* dans SBG/Output/MO. Une fois vérifié que tout s'est bien passé (voir ci dessous), pensez à supprimer ces fichiers sauvés. |
| 94 | * Comment verifier que tout s'est bien passé ? |
| 95 | Par exemple avec les deux fichiers suivants : piControl2_26201201_26201231_1M_stomate_ipcc_history.nc (fichier après correction) et SAVE_9900/piControl2_26201201_26201231_1M_stomate_ipcc_history.nc (fichier avant correction) |
| 96 | {{{ |
| 97 | cdo diffv piControl2_26201201_26201231_1M_stomate_ipcc_history.nc SAVE_9900/piControl2_26201201_26201231_1M_stomate_ipcc_history.nc |
| 98 | 12 : 2620-12-16 12:00:00 fLuc 0 9216 6364 : F T 3.8952e-06 0.0000 |
| 99 | 13 : 2620-12-16 12:00:00 nbp 0 9216 6364 : T T 3.8952e-06 1.00000 |
| 100 | cdo infov piControl2_26201201_26201231_1M_stomate_ipcc_history.nc |
| 101 | 12 : 2620-12-16 12:00:00 fLuc 0 9216 6364 : 0.0000 0.0000 0.0000 |
| 102 | 13 : 2620-12-16 12:00:00 nbp 0 9216 6364 : -7.1456e-08 -2.0327e-09 5.9766e-08 |
| 103 | cdo infov SAVE_9900/piControl2_26201201_26201231_1M_stomate_ipcc_history.nc |
| 104 | 12 : 2620-12-16 12:00:00 fLuc 0 9216 6364 : 0.0000 1.2280e-06 3.8952e-06 |
| 105 | 13 : 2620-12-16 12:00:00 nbp 0 9216 6364 : -3.9240e-06 -1.2301e-06 5.9766e-08 |
| 106 | }}} |
| 107 | On peut aussi vérifier la valeur de la variable nbp dans le nouveau fichier à l'aide du monitoring. |
| 108 | Par exemple pour picontrol2 : [[BR]] |
| 109 | [[BR]] |
| 110 | [[Image(nbp_piControl2.png, 120px)]] |
95 | | Ce script est à faire installé dans modipsl/libIGCM pour faire un ins_job et à faire tourner dans une copie de la configuration sur une frontale, |
96 | | sans argument. Il modifiera alors l'ensemble des '''stomate_ipcc_history.nc'''. Il faudra ensuite recréer les TS IPCC et les MONITORINGs si des variables IPCC sont monitorées. |
97 | | Ce n'est pas le cas pour le couplé par défaut. |
98 | | Il ne transfère que les stomate_history.nc et les stomate_history_ipcc.nc pour ne corriger que les deux variables nbp et fLuc de ces derniers. |
99 | | Son temps d'exécution est de 30 minutes au minimum pour 10 ans (120x2 fichiers à démigrer), en fonction des capacités de démigration du dmnfs. |
100 | | Pour l'instant (rev 1310), se relance est problématique car il reprend la liste complète des fichiers. |
| 132 | |
| 133 | Ce script modifie l'ensemble des '''stomate_ipcc_history.nc''' d'une simulation. Il faut ensuite recréer les TS IPCC et les MONITORINGs (si des variables du fichier stomate_ipcc_history.nc sont monitorées). Ce script necessite seulement le transfert des fichiers stomate_history.nc et stomate_history_ipcc.nc pour ne corriger que les deux variables nbp et fLuc de ces derniers. Son temps d'exécution est de 30 minutes au minimum pour 10 ans (120x2 fichiers à démigrer), en fonction des capacités de démigration du dmnfs. |
| 134 | Plus d'infos sur son utilisation là [wiki:ModifIPSLCM5A#a25/03/2011] |