1 | #define utilitaires_de_block_archeops
|
---|
2 |
|
---|
3 | #include "diabolo.h"
|
---|
4 | #include "bolo.h"
|
---|
5 | #include "acqui.h"
|
---|
6 | #include "carte_acqui.h"
|
---|
7 | #include "carte_pci.h"
|
---|
8 | #include "trace.h"
|
---|
9 | #include "corel.h"
|
---|
10 | #include "montimer.h"
|
---|
11 | #include "synchro.h"
|
---|
12 | #include "fichier.h"
|
---|
13 | #include "onde.h"
|
---|
14 | #include "tm.h"
|
---|
15 | #include "tache.h"
|
---|
16 | #include "compress.h"
|
---|
17 | #include "controle.h"
|
---|
18 | #include "menu.h"
|
---|
19 |
|
---|
20 |
|
---|
21 |
|
---|
22 |
|
---|
23 | //*****************************************************************************************
|
---|
24 | // ***
|
---|
25 | #ifdef _archeops //--------- pour Archeops ------------------------------ ***
|
---|
26 | //--------------- lit block Archeops ------------------------------ ***
|
---|
27 | // ***
|
---|
28 | //*****************************************************************************************
|
---|
29 |
|
---|
30 |
|
---|
31 | // la fonction litblock
|
---|
32 | // 1 elle lit dans la table de telemesures
|
---|
33 | // 2 elle teste le type de block et sa validité (somme des bits et code fin)
|
---|
34 | // 3 elle appelle la fonction de traitement en lui fournissant un pointeur sur le block
|
---|
35 | // 4 elle sauve le block dans le fichier
|
---|
36 | // C'est la fonction acqui qui apelle lit_block qui avance le pointeur de la table de telemesures
|
---|
37 |
|
---|
38 |
|
---|
39 | block_type_modele bktm; // memoire pour block local pour copie lors de la decompression
|
---|
40 |
|
---|
41 | block_type_synchro_sol blk_sync2; // block memoire pour les synchro sol qui seront traitées lors
|
---|
42 | // de la lecture du block bolo
|
---|
43 |
|
---|
44 |
|
---|
45 |
|
---|
46 | #define blk_tm (&tt->tm.btm[tt->tm.pos_lit].tmtrx) // pointeur sur le block télémesure a lire)
|
---|
47 |
|
---|
48 | void lit_block_archeops(int err)
|
---|
49 | {
|
---|
50 | if(verifie_block_printf(blk_tm)) return;
|
---|
51 |
|
---|
52 | {
|
---|
53 | def_nom_block
|
---|
54 | static int mon_num_block;
|
---|
55 | if(mon_num_block!=numero_block(blk_tm))
|
---|
56 | {
|
---|
57 | mon_num_block=numero_block(blk_tm);
|
---|
58 | printf("\n%d ->",mon_num_block);
|
---|
59 | }
|
---|
60 | printf(" %s ",nom_block[type_block(blk_tm)]);
|
---|
61 | }
|
---|
62 |
|
---|
63 | switch(type_block(blk_tm))
|
---|
64 | {
|
---|
65 | case block_param : break;
|
---|
66 |
|
---|
67 | case block_journal : traite_block_journal((block_type_journal*)blk_tm); break;
|
---|
68 |
|
---|
69 | case block_reglage : traite_block_reglage((block_type_reglage*)blk_tm); break;
|
---|
70 |
|
---|
71 | case block_dilution : traite_block_dilution((block_type_dilution*)blk_tm); break;
|
---|
72 |
|
---|
73 | case block_gps : traite_block_gps((block_type_gps*)blk_tm); break;
|
---|
74 |
|
---|
75 | case block_une_periode : traite_block_une_periode((block_type_une_periode*)blk_tm); break;
|
---|
76 |
|
---|
77 | case block_bolo : traite_block_bolo((block_type_bolo*)blk_tm); break;
|
---|
78 | case block_bolo_comprime : traite_block_bolo_comprime((block_type_bolo_comprime*)blk_tm); break;
|
---|
79 |
|
---|
80 | case block_sst : traite_block_sst((block_type_sst*)blk_tm); break;
|
---|
81 | case block_sst_comprime : traite_block_sst_comprime((block_type_sst_comprime*)blk_tm); break;
|
---|
82 |
|
---|
83 | case block_gyro : traite_block_gyro((block_type_gyro*)blk_tm); break;
|
---|
84 | case block_gyro_comprime : traite_block_gyro_comprime((block_type_gyro_comprime*)blk_tm); break;
|
---|
85 |
|
---|
86 | case block_synchro_sol : traite_block_synchro_sol((block_type_synchro_sol*)blk_tm); break;
|
---|
87 |
|
---|
88 | case block_status_flash : printf("******************status flash recut ********************* \n");son(133); break;
|
---|
89 |
|
---|
90 | #ifdef transmet_data_brute
|
---|
91 | case block_data_brute : traite_block_data_brute((block_type_data_brute*)blk_tm); break;
|
---|
92 | #endif
|
---|
93 |
|
---|
94 | default : {def_nom_block printf("block de type %s non traité par le programme \n",nom_block[type_block(blk_tm)]);}
|
---|
95 | break;
|
---|
96 | }
|
---|
97 |
|
---|
98 | ecrit_fichier_ARCH(blk_tm);
|
---|
99 |
|
---|
100 | }
|
---|
101 | #undef blk_tm
|
---|
102 |
|
---|
103 | //*****************************************************************************************
|
---|
104 | // ***
|
---|
105 | // ------------------ traitement des blocks --------------- ***
|
---|
106 | // ***
|
---|
107 | //*****************************************************************************************
|
---|
108 |
|
---|
109 |
|
---|
110 |
|
---|
111 |
|
---|
112 | //-----------------------------------------------------------------------------------------
|
---|
113 | // ------------- block type 1 : lit le journal ------------------ ------------
|
---|
114 | //-----------------------------------------------------------------------------------------
|
---|
115 | void traite_block_journal(block_type_journal* blk)
|
---|
116 | {
|
---|
117 | int i,code,blocknum;
|
---|
118 | if(!fenetre(fenetre_journal)) nouveauT(fenetre_journal,0,"journal");
|
---|
119 |
|
---|
120 | //printf("ecrit une page du journal \n");
|
---|
121 | for(i=0;i<long_page;i++)
|
---|
122 | {
|
---|
123 | code=type_block(&blk->jj[i]);
|
---|
124 | blocknum=numero_block(&blk->jj[i]);
|
---|
125 | if(code) break; // code vaut 1 pour vider une page de journal incomplete
|
---|
126 | ecritT(fenetre_journal,fin_f,"block_num= %d code %d ==> %d %d %d %d %d %d %d %d \n",
|
---|
127 | blocknum,code,
|
---|
128 | blk->jj[i].mot[0],blk->jj[i].mot[1],blk->jj[i].mot[2],blk->jj[i].mot[3],
|
---|
129 | blk->jj[i].mot[0],blk->jj[4].mot[5],blk->jj[i].mot[6],blk->jj[i].mot[7]);
|
---|
130 | }
|
---|
131 | }
|
---|
132 |
|
---|
133 | //-----------------------------------------------------------------------------------------
|
---|
134 | // ------------- block type 2 : block_reglage ------------------ ------------
|
---|
135 | //-----------------------------------------------------------------------------------------
|
---|
136 | void traite_block_reglage(block_type_reglage* blk)
|
---|
137 | {
|
---|
138 |
|
---|
139 | //printf(" recut un block reglage *********************************** \n\n");
|
---|
140 |
|
---|
141 | gg->reglage=blk->reglage;
|
---|
142 | tt->reglage=blk->reglage;
|
---|
143 | rafraichis_fenetres();
|
---|
144 | }
|
---|
145 |
|
---|
146 |
|
---|
147 | //-----------------------------------------------------------------------------------------
|
---|
148 | // --------------------- block type 3 : lit les data dilution ------------------------
|
---|
149 | //-----------------------------------------------------------------------------------------
|
---|
150 | //
|
---|
151 | // fonction void traite_block_dilution(block_type_dilution*blk)
|
---|
152 | //
|
---|
153 | // la fonction est dans le module dilution.c
|
---|
154 |
|
---|
155 |
|
---|
156 |
|
---|
157 | //-----------------------------------------------------------------------------------------
|
---|
158 | // --------------------- block type 4 : lit la chaine GPS ------------------------
|
---|
159 | //-----------------------------------------------------------------------------------------
|
---|
160 | void traite_block_gps(block_type_gps* blk)
|
---|
161 | {
|
---|
162 | int i;
|
---|
163 | char c1;
|
---|
164 | if(!fenetre(fenetre_GPS)) nouveauT(fenetre_GPS,0,"lecture GPS");
|
---|
165 |
|
---|
166 | //printf("ecrit une ligne GPS \n");
|
---|
167 |
|
---|
168 | if(blk->gps[0])
|
---|
169 | {
|
---|
170 | ecritT(fenetre_GPS,fin_f,"gps --> ");
|
---|
171 | for(i=0;i<80;i++)
|
---|
172 | {
|
---|
173 | c1=blk->gps[i];
|
---|
174 | if (c1) ecritT(fenetre_GPS,fin_f,"%c",c1&0x7f);
|
---|
175 | else i=80;
|
---|
176 | if(c1==13) i=80;
|
---|
177 | }
|
---|
178 | }
|
---|
179 | }
|
---|
180 |
|
---|
181 |
|
---|
182 |
|
---|
183 | //-----------------------------------------------------------------------------------------
|
---|
184 | // ------------- block type 5 : lit les courbes de bolos individuelles ------------
|
---|
185 | // nb_per_block est le nombre de periodes entieres dans le block
|
---|
186 | // --> le nombre de courbes lues est egal a nb_per_block qui doit etre <= nb_max_bolo
|
---|
187 | //-----------------------------------------------------------------------------------------
|
---|
188 | void traite_block_une_periode(block_type_une_periode* blk)
|
---|
189 | {
|
---|
190 | int j,i;
|
---|
191 | for(j=0;j<nb_max_bolo;j++) for (i=0;i<gg->reglage.horloge.nb_mesures;i++)
|
---|
192 | gg->bol_per[j][i]=blk->bol_per[j][i]^0x7fff; // 16 bit -> 15 bit utile + phase
|
---|
193 | // change le signe des 15 bits utiles
|
---|
194 | affi_bolo(); // affichage des bolo individuels
|
---|
195 | }
|
---|
196 |
|
---|
197 |
|
---|
198 |
|
---|
199 |
|
---|
200 | //-----------------------------------------------------------------------------------------
|
---|
201 | // ------------- block type 6 : block_synchro_sol ------------------ ------------
|
---|
202 | //-----------------------------------------------------------------------------------------
|
---|
203 | void traite_block_synchro_sol(block_type_synchro_sol* blk)
|
---|
204 | {
|
---|
205 | blk_sync2 = *blk;
|
---|
206 | // ne fait rien: le block synchro precede toujours le block bolo
|
---|
207 | }
|
---|
208 |
|
---|
209 |
|
---|
210 |
|
---|
211 |
|
---|
212 | //-----------------------------------------------------------------------------------------
|
---|
213 | // --------------------- block type 12 : décomprime les données bolo ------------
|
---|
214 | //-----------------------------------------------------------------------------------------
|
---|
215 | void traite_block_bolo_comprime(block_type_bolo_comprime* blk)
|
---|
216 | {
|
---|
217 | int j,jc;
|
---|
218 | #define blk2 ((block_type_bolo*)(&bktm)) // pointeur sur la memoire locale
|
---|
219 | // maintenant, decomprime blk pour le reecrire dans blk2
|
---|
220 | /*printf("comprime: ds1= %x , %x ds2= %x , %x ",
|
---|
221 | ( (blk->data_bolo[0][0] & 0xffff0000)>>11 ) & 0x001fffe0,
|
---|
222 | ( (blk->data_bolo[0][0] & 0x0000ffff)<<5 ) & 0x001fffe0,
|
---|
223 | ( (blk->data_bolo[1][0] & 0xffff0000)>>11 ) & 0x001fffe0,
|
---|
224 | ( (blk->data_bolo[1][0] & 0x0000ffff)<<5 ) & 0x001fffe0);
|
---|
225 | */
|
---|
226 | /*
|
---|
227 | for(j=0;j<nb_bolo_util;j++)
|
---|
228 | {
|
---|
229 | decompress_7_2((long *)blk->data_bolo[j],blk2->data_bolo[j],nb_per_block*2);
|
---|
230 | }
|
---|
231 | */
|
---|
232 | jc=0;for(j=0;j<nb_max_bolo;j++) // jc = bolo_comprime // j=bolo normal
|
---|
233 | {
|
---|
234 | if( (parametr.bolo[j].bolo_code_util==bolo_normal_transmis) && (jc<nb_bolo_util) )
|
---|
235 | {
|
---|
236 | decompress_7_2((long *)blk->data_bolo[jc],blk2->data_bolo[j],nb_per_block*2);
|
---|
237 | jc++;
|
---|
238 | }
|
---|
239 | }
|
---|
240 | /*printf("ds1= %x , %x ds2= %x , %x \n",
|
---|
241 | blk2->data_bolo[0][0],
|
---|
242 | blk2->data_bolo[0][1],
|
---|
243 | blk2->data_bolo[1][0],
|
---|
244 | blk2->data_bolo[1][1]);
|
---|
245 | */
|
---|
246 | valide_block(&bktm,block_bolo,numero_block(blk)); // pour recopier le numero de block
|
---|
247 | traite_block_bolo(blk2);
|
---|
248 | }
|
---|
249 |
|
---|
250 |
|
---|
251 |
|
---|
252 |
|
---|
253 | //-----------------------------------------------------------------------------------------
|
---|
254 | // ------------- block type 8 : block_bolo ------------------ ------------
|
---|
255 | //-----------------------------------------------------------------------------------------
|
---|
256 | void traite_block_bolo(block_type_bolo* blk)
|
---|
257 | {
|
---|
258 | static num_block_local;
|
---|
259 | int i,j,nb_coups,glitch;
|
---|
260 | static ancien_don_bolo[nb_max_bolo];
|
---|
261 | int secondes;
|
---|
262 |
|
---|
263 | glitch=litD(fenetre_acquisition,e_deglitch,0L);
|
---|
264 |
|
---|
265 | num_block_local++;
|
---|
266 | if(num_block_local!=numero_block(blk)) {printf(" erreur numblock=%d\n",numero_block(blk));num_block_local=numero_block(blk);}
|
---|
267 |
|
---|
268 | secondes=(int) (gg->periode_echantillonage*(double)num_block_local*(double)nb_per_block*2.);
|
---|
269 | //printf("numblock=%d sec=%d \n",num_block_local,secondes);
|
---|
270 | ecritD(fenetre_controle,con_time,"%d",secondes);
|
---|
271 |
|
---|
272 | nb_coups= gg->reglage.horloge.nb_mesures/2 - gg->reglage.horloge.temp_mort;
|
---|
273 |
|
---|
274 | /*
|
---|
275 | if( (num_block_local%10)==0) // impression des valeurs un block sur 10
|
---|
276 | {
|
---|
277 | int aa = (nb_coups<<14) + (nb_coups*190) ;
|
---|
278 | printf(" nb coups dans une somme : %d \n",nb_coups);
|
---|
279 | for(j=0;j<20;j+=2)
|
---|
280 | printf("ds1= %x , %x ds2= %x , %x ",
|
---|
281 | blk->data_bolo[0][j],
|
---|
282 | blk->data_bolo[0][j+1],
|
---|
283 | blk->data_bolo[1][j],
|
---|
284 | blk->data_bolo[1][j+1]);
|
---|
285 | printf("ds1= %8.1f , %8.1f ds2= %8.1f , %8.1f \n",
|
---|
286 | (double)(blk->data_bolo[0][j]- aa),
|
---|
287 | (double)(blk->data_bolo[0][j+1]- aa),
|
---|
288 | (double)(blk->data_bolo[1][j]- aa),
|
---|
289 | (double)(blk->data_bolo[1][j+1]- aa));
|
---|
290 | }
|
---|
291 | */
|
---|
292 |
|
---|
293 | // ------------- lit les valeur detection synchrone ------------------
|
---|
294 |
|
---|
295 | #define val_DS(j,i) (blk->data_bolo[j][i]&0x1fffff)
|
---|
296 | //#define synchro_DS(j,i) ((blk->data_bolo[j][i]>>20)&1)
|
---|
297 |
|
---|
298 | for (i=0;i<nb_per_block*2;i++)
|
---|
299 | {
|
---|
300 | int aa,s;
|
---|
301 | aa = (nb_coups<<14) + (nb_coups*190) ;
|
---|
302 | if( i%2) s=1;else s=-1;
|
---|
303 |
|
---|
304 | for(j=0;j<nb_max_bolo;j++)
|
---|
305 | {
|
---|
306 | gg->don.don_bolo[j]=0.5*(ancien_don_bolo[j]+s*(((val_DS(j,i)-aa)<<1)/nb_coups)); // copy bolos
|
---|
307 | ancien_don_bolo[j]=s*(((val_DS(j,i)-aa)<<1)/nb_coups);
|
---|
308 | }
|
---|
309 |
|
---|
310 | for(j=0;j<nb_sync;j++)
|
---|
311 | {
|
---|
312 | gg->don.sync[j]=blk_sync2.synchro[i][j];
|
---|
313 | }
|
---|
314 |
|
---|
315 |
|
---|
316 | // gg->don.sync[j]=synchro_DS(j,i);
|
---|
317 |
|
---|
318 | gg->temps_cntl=num_block_local*nb_per_block*2+i; // temps controlleur en demi-periodes
|
---|
319 | // printf("gg->temps_cntl=%d \n",gg->temps_cntl);
|
---|
320 | // exec_synchro(0,100,err); // en premier pour construire le flag scan_en_cours
|
---|
321 | exec_synchro(0,100,0); // en premier pour construire le flag scan_en_cours
|
---|
322 | // le trace ds et 1p est fait avec les points corelles du coups precedent !!!
|
---|
323 |
|
---|
324 | GetDateTime(&gg->temps_mac);
|
---|
325 | if(gg->fichier_dia2_ouvert) ecrit_fichier_dia2(gg->don,gg->temps_cntl,gg->temps_mac); // ecrit les donnees dans le fichier
|
---|
326 |
|
---|
327 | if(glitch) deglitch();
|
---|
328 | corel();
|
---|
329 |
|
---|
330 | //if(gg->trace_ouvert&1) exec_trace_tout(fenetre_trace_tout,tache_de_fond,0);
|
---|
331 | for(j=0;j<8;j++) if(gg->trace_tout[j]) exec_trace_tout(gg->trace_tout[j],tache_de_fond,0);
|
---|
332 | if( (gg->trace_ouvert&16) && // trace carte et pas de select synchro
|
---|
333 | (! (gg->trace_ouvert&32))) trace_annexe(gg->corel);
|
---|
334 | if(gg->flag_fou) foufou();
|
---|
335 | if(gg->trace_ouvert&8) exec_ondelette(fenetre_ondelette,0,0);
|
---|
336 |
|
---|
337 |
|
---|
338 | }
|
---|
339 |
|
---|
340 | // _____________ affichage numerique de la valeur d'un bolometre _______________
|
---|
341 | j=litD(fenetre_acquisition,e_bolo,0L)-1;
|
---|
342 | if( (j<0) ou (j>=nb_max_bolo) ) ecritD(fenetre_acquisition,e_valeur,"--");
|
---|
343 | else ecritD(fenetre_acquisition,e_valeur,"b%d=%g µV",j+1,xbol(j));
|
---|
344 |
|
---|
345 |
|
---|
346 | }
|
---|
347 |
|
---|
348 |
|
---|
349 |
|
---|
350 |
|
---|
351 | #ifdef transmet_data_brute
|
---|
352 | //-----------------------------------------------------------------------------------------
|
---|
353 | // ------------- block type 15 : lit data brute (block recut de l'EPLD) ------------ ------------
|
---|
354 | // ------------- utilise uniquement pour la mise au point ------------
|
---|
355 | //-----------------------------------------------------------------------------------------
|
---|
356 | void traite_block_data_brute(block_type_data_brute* blk)
|
---|
357 | {
|
---|
358 | int i,k;
|
---|
359 | static int qq=0;
|
---|
360 | qq++;
|
---|
361 | //if(qq>20)
|
---|
362 | {
|
---|
363 | qq=0;
|
---|
364 | printf(" data brutes : ");
|
---|
365 | for(k=0;k<10;k+=5) {printf(" ");for (i=0;i<5;i++) printf(" %x", (unsigned char) blk->data[k+i+1]);}
|
---|
366 | printf("\n");
|
---|
367 | for(k=0;k<nb_per_block;k++) printf("%d",(blk->data[k*tt->reglage.horloge.nb_mesures*nb_byte_mot+5]>>6)&0x3);
|
---|
368 | printf("\n");
|
---|
369 | // for(k=0;k<2;k++)
|
---|
370 | {
|
---|
371 | // printf("%d",(blk->data[k*tt->reglage.horloge.nb_mesures*nb_byte_mot+2]>>6)&0x3);
|
---|
372 | for(i=0;i<tt->reglage.horloge.nb_mesures;i++)
|
---|
373 | {
|
---|
374 | // printf("%d",(blk->data[(k*tt->reglage.horloge.nb_mesures+i)*nb_byte_mot+5]>>6)&0x3);
|
---|
375 | }
|
---|
376 | // printf("\n");
|
---|
377 | }
|
---|
378 |
|
---|
379 | }
|
---|
380 | }
|
---|
381 | #endif
|
---|
382 |
|
---|
383 |
|
---|
384 |
|
---|
385 |
|
---|
386 | #endif
|
---|
387 |
|
---|