problème de monitoring

[sdipsl]

Pour les sorties orchidee spinup, le monitoring avec ferret/6.9.3 ne fonctionne pas correctement. A creuser.

###### les mails de Josefine qui détaillent le souci ######

Salut,

le probleme vient du version de ferret. Avec ferret/6.8.4 ca marche mais pas avec ferret/6.9.3.

A+
Josefine

Le 12/11/2015 15:38, Josefine Ghattas a écrit :

Salut

voici mon problème monitoring avec un run spinup. Vois le résultat ici :
​http://dods.extra.cea.fr/work/p86ghatt/OL2/PROD/ref3013/FG1low.CWRR.3013/MONITORING

Le problème se voit deja dans les fichiers dans MONITORING. Ca a l'air de venir du fait que l'attribut units (units = "seconds since 1901-01-01 00:00:00") est faux, le TIME_COUNTER est en "year since" ou uniquement un axe qui increment par 1 (1, 2, 3,..) :
/ccc/work/cont003/dsm/p86ghatt/IGCM_OUT/OL2/PROD/ref3013/FG1low.CWRR.3013/MONITORING/files/SBG_HET_RESP_lands_ave.nc

Le TS correspondant a l'air mieux. L'attribut units est le meme (en seconds since) mais ca correspond a l'axe :
/ccc/store/cont003/dsm/p86ghatt/IGCM_OUT/OL2/PROD/ref3013/FG1low.CWRR.3013/SBG/Analyse/TS_YE/FG1low.CWRR.3013_19010101_22201231_1Y_HET_RESP.nc

La mise en place de la simul est ici:
/ccc/work/cont003/dsm/p86ghatt/SIMULATION_ORCHIDEE/ORCr3013_LMDZr2314/modipsl/config/ORCHIDEE_OL/FG1low.CWRR.3013
et les sorties des job ici :
/ccc/scratch/cont003/dsm/p86ghatt/IGCM_OUT/OL2/PROD/ref3013/FG1low.CWRR.3013/Out

Avant, dans la simul ref2724/FG1low.CWRR.2724 ca a marché. Vois monitoring :
​http://dods.extra.cea.fr/work/p86ghatt/OL2/PROD/ref2724/FG1low.CWRR.2724/MONITORING/

Le fichier du monitoring a l'axe de temps annuelle 1, 2, 3... mais les attributs n'y était pas:
/ccc/work/cont003/dsm/p86ghatt/IGCM_OUT/OL2/PROD/ref2724/FG1low.CWRR.2724/MONITORING/files/SBG_HET_RESP_lands_ave.nc

Les TS ont l'air pareil avec les meme attributs que maintenant:
/ccc/store/cont003/dsm/p86ghatt/IGCM_OUT/OL2/PROD/ref2724/FG1low.CWRR.2724/SBG/Analyse/TS_YE/FG1low.CWRR.2724_19010101_22201231_1Y_HET_RESP.nc

Donc pour créer les fichiers de monitoring, on a toujours changé l'axe pour un axe "year since". Avant on enlevait tout les attribue sur l'axe de temps et ca marchait. Maintenant on a gardé les attributes du fichier TS et ca ne correspond plus a l'axe. Ca ne marche donc pas. Si tu sais d'où vient cette changement ou quoi faire ? Merci d'avance!!

A+
Josefine

[jgipsl]

Le problème avec l'axe de temps vient du fait qu'il n'est pas monotonic. Avec l'ancien ferret/6.8.4 l'axe etait transformé en indices et les attribues enlevés. Maintenant avec ferret/6.9.3, l'axe est transformé en indices mais les attribues sont gardés. Ce sont les attribues qui posent le problème.

Apres discussion avec Patrick, il faudrait tester pour le cas IgnoreNonMonotonic=TRUE (option deja existant dans config.card), d'enlever au moins l'attribut units du time_counter lors de la création des time-serises. Ensuite le monitoring pourrait fonctionner comme avant.

[jgipsl]

J'ai testé le suivant sur curie a la fin du fonction IGCM_sys_ncrcat et ca a l'air de fonctionné bien. Test fait avec libIGCM_v2.7 et ferret/6.9.3.

 # Special case for spinup, cyclic time axis: change attribute calendar to none
  if [ X${config_Post_IgnoreNonMonotonic} = XTRUE ] ; then
      eval lastArg=\${$#}
      echo "IGCM_sys_ncrcat : change attribute calendar to none in file:" ${lastArg}
      IGCM_sys_ncatted -a calendar,time_counter,m,c,none ${lastArg}
  fi

A commiter ? Quel tests a faire avant ?

[jgipsl]

Fait in rev [1278]