XBTTSG
XBTTSG est une suite de programme permettant de traiter des XBT et de les comparer aux donnees TSG.
Alexandre Ganachaud (www.ird.nc/ECOP) est l'auteur pour tous commentaires.
Resultats
1er partie: tri des XBT et comparaison avec le point TSG proche
-
1-Verifie que les tirs XBT sont dans les 5 milles du point TSG. Ecrit des fichiers d'erreur par bateau: ALIS, ATAL, COIS, EXPL, LOND, NORO, PAIS, PROV, RIMB, ROME, RONS, VOYA, WASH: au format
Deviation RONS Drop 341 5 nm t-118.17 x-118.10 t-11.92 x-11.86 67755
tir milles d'erreur long. TSG long. XBT lat TSG lat XBT tir ID -
2-affiche globalement toutes les deviations a chaque profondeur XBT: deviation1, deviation2, deviation3, deviation4.
-
3-sauve l'ensemble (avec positions corrigees) dans thxbt_corrected.mat.
-
statistiques sur les XBT pour eliminer d'eventuels mauvais points: pour l'instant par bandes de latitudes. A ameliorer pour avoir plutot des carres de lat-long.
-
sauve danshxbt_filtered.mat -sans les "mauvais" tirs.
sauve aussi l'ID des mauvais tirs "badshots" pour corrections et des "notsobadshots" qui sont hors statistiques mais en accord avec le TSG (donc gardes).
-
xbttsg_bias.m: biais par voyage entre les XBT moyens de 0 a 10m et les TSG.
-
xbttsg_pltbpoints.m: plot des carte pour identifier les voyages malades: exemple pour les points avec ecarts superieurs a 4 degres: xbtdev4.0, 2 degres xbtdev2.0, 1 degre xbtdev1.0.
-
xbttsg_biaspertrip.m: calcul des biais par voyage pour chaque bateau.
-
Avant elimination des "mauvais" voyages: ALIS, ATAL, COIS, EXPL, LOND, NORO, PAIS, PROV, RIMB, ROME, RONS, VOYA, WASH.
-
Apres elimination on peu lire en haut des graphiques le biais general du bateau: ALIS, ATAL, COIS, EXPL, LOND, NORO, PAIS, PROV, RIMB, ROME, RONS, VOYA, WASH.On peut lire les voyages elimines dans le fichier xbttsg_bias.m.
Programmes
-
read_tagxbt_tsv.m: lit le fichier .tsv extrait de la base SQL. Output=Data1/thxbt.mat contenant une structure avec toutes les donnees
-
xbttsg_checkposplot.m: Verifie que les tirs XBT sont dans les 5 milles du point TSG. Ecrit des fichiers d'erreur (XBT_badpos_ship.txt) pour chaque bateau. Corrige les positions XBT en employant les positions TSG. Output thxbt_corrected.mat (meme variable thxbt).
-
xbttsg_filter.m: Filtre les XBT en eliminant les XBT hors de 2 deviations standard par bandes de latitudes; garde quand meme ceux qui sont proches (1 degre) du TSG. Output= Data1/hxbt_filtered.mat incluant les badshots (elimines de thxbt) et les notsobadshots en accord avec le tsg.
-
xbttsg_bias.m: calcul le bias + l'ecart type des deviations entre XBT et TSG par voyage, pour chaque bateau. Faire tourner manuellement et iterativement en eliminant les voyages suspects (structure badtrips)
2ème partie: visualisation des séries TSG avec comparaison XBT
- CORAL ISLANDER
- cois199905 peu de biais mais un XBT outlier
- cois199906 biais bruite de l'ordre de 1degre. Qu'en penser ? Bon ? Garder la salinite ?
- cois200001 quelques XBT hors serie TSG. Que faire ? Lesquels eliminer ?
- cois200002
- cois200003 tres bon
- cois200004
- cois200005 3XBT. pas de point de comparaison. Ajoute-t-on Reynolds ?
- cois200102
- cois200104
- cois200105
- cois200106
- cois200107
- lond199801
- lond199802
- lond199901 probleme de données
- lond199902
- pais200103 interessants biais et non-biais: on garde tout ? Dans quelles limites ?
- pais200104
- pais200106
- pais200201
- pais200202
- pais200203
- pais200204
- prov199701 de bien curieux biais
- prov199702 limite parfait !
- rimb199601
- rimb199602
- rimb199603 2 degres de biais sur la fin ! Que faire ?
- rimb199701
- rimb199702 5 a 6 degres sur la fin ! A jeter mais la salinite semble OK.
- rimb199703