Changeset 645 in Sophya for trunk/Poubelle/archediab.old/archediab.sources/c/dilution.c

Timestamp:

Nov 25, 1999, 2:35:30 PM (26 years ago)

Author:

ansari

Message:

archediab 27

File:

• trunk/Poubelle/archediab.old/archediab.sources/c/dilution.c (modified) (3 diffs)

trunk/Poubelle/archediab.old/archediab.sources/c/dilution.c

@@ -208 +208 @@
 
 
+#define c(j,i)          (1e-4*(double)parametr.nom_coef[parametr.bolo[j].numero_nom_coef].coef[i])      
 
 // les temperatures sur les cartes modifiées
 {
-double I,V,R[4];
+double I,V,R[8],llR;
 int j,k;
 def_gains;
 k=0;
+
 for(j=0;(j<nb_max_bolo) && (k<4);j++)
         {
@@ -222 +224 @@
                 V=0.001*bol_micro_volt(blk->temperature[k],(double)parametr.bolo[j].bolo_gain*gain_ampli(gg->reglage.bolo[j]));
                 if(I>0.0000001)         R[k]=V/I;       else    R[k]=0;                         //      R  en ‡
-                k++;
+ //------  calcul de l'etalonnage en temperature des cartes temperature simplifiées  
+ // j=numero de bebo    T[k] = R[k+4] en Kelvin
+                if ((R[k]-c(j,6))>1.) 
+                        {if((log(R[k]-c(j,6))-c(j,0))>0.001) llR= log(log(R[k]-c(j,6))-c(j,0)) ; else llR=0;} 
+                else llR=0;
+                R[k+4]  = exp(  c(j,1) + c(j,2)* llR + c(j,3)* llR* llR + c(j,4)* llR* llR* llR + c(j,5)* llR* llR* llR* llR) ; 
+                if(R[k+4]>9999) R[k+4]=9999;
+                k++;            
                 }
         }
 
-if(R[0]<90.)    ecritD(fenetre_dilution,dil_temp,"germ= %6.2f‡  4K=%5.1f‡   1.6K=%5.1f‡   10K=%5.1f‡",R[0],R[1],R[2],R[3]);
-else    ecritD(fenetre_dilution,dil_temp,"germ= %6.3fk‡  4K=%5.1f‡   1.6K=%5.1f‡   10K=%5.1f‡",R[0]*0.001,R[1],R[2],R[3]);
+
+if(R[0]<90.)    ecritD(fenetre_dilution,dil_temp
+        , "ge=%6.2f‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fK  4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
+                ,R[0],R[1],R[2],R[3],R[4],R[5],R[6],R[7]);
+else    ecritD(fenetre_dilution,dil_temp
+        ,"ge=%6.3fk‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fmK 4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
+        ,R[0]*0.001,R[1],R[2],R[3],R[4]*1000.,R[5],R[6],R[7]);
+
+
 if(litD(fenetre_dilution,dil_trace_tempe,0))// trace les sondes de niveau helium
         {
@@ -240 +256 @@
                         strcpy(graph->xtitre,"minutes");
                         }
-        tracen(fenetre_temperature_dil,4,minutes,R);
+        tracen(fenetre_temperature_dil,8,minutes,R);
         }
 }