minor changes

cal/ce3.calDeltaVDeltat.R

@@ -76,13 +76,13 @@ ce3.calDeltaVDeltat <- function(ccc){
             ci         <- mp - slope *  mt
             t0         <- (mean(mp, na.rm =TRUE) - ci) / slope
             
-            out        <- which(correlCoef < 0.8)
+            out        <- which(correlCoef < 0.9)
             if(length(out) > 0){
                 t0[out] <- NA
                 msg <- paste(msg
                            , "exclude saw tooth:"
                            , toString(out)
-                           , "because of bad (<0.8) correlation")    
+                           , "because of bad (<0.9) correlation")    
             }
             
 
@@ -124,8 +124,8 @@ ce3.calDeltaVDeltat <- function(ccc){
 
             ## Leitwert = dV/dt *(1- m.D/m.SZ)
 
-            dVdt       <- mean(deltaV, na.rm =TRUE) / mean(deltat, na.rm =TRUE) * (1 - dcorr[j])
-
+#            dVdt       <- mean(deltaV, na.rm =TRUE) / mean(deltat, na.rm =TRUE) * (1 - dcorr[j])
+            dVdt       <- mean(deltaV, na.rm =TRUE) / mean(deltat, na.rm =TRUE) 
             L[j]       <- mean(dVdt, na.rm =TRUE)       
            
             sddV[j]    <- sd(deltaV, na.rm =TRUE)/ mean(deltaV, na.rm =TRUE)

devel/calYampCe3.R

@@ -31,8 +31,9 @@ source("load.R")
 ## cdb$id     <- "ac36e53c6c34d30886c210a6f60aa61b"
 ## cdb$id     <- "ac36e53c6c34d30886c210a6f60a9091" # daresb H2
 ## cdb$id     <- "ac36e53c6c34d30886c210a6f608ad31" ## lwc, lw1, lw2
-cdb$id     <- "ac36e53c6c34d30886c210a6f6089e5a"
+## cdb$id     <- "ac36e53c6c34d30886c210a6f6089e5a"
 ## cdb$id     <- "ac36e53c6c34d30886c210a6f60a9091" # daresb H2
+cdb$id     <- "ac36e53c6c34d30886c210a6f61bd27e" # LEY
 doc <- cdbGetDoc(cdb)$res
 
 if(TRUE){    

devel/srg_vg.R

@@ -21,10 +21,11 @@ source("load.R")
 
 pUnit  <- "mbar"
 
-#cdb$id <- "89e85b1c42ee6a9b70864a7612251a7e" # 14. & 15.10.2014 EDW Kalib 75093_0002
-#cdb$id <- "4a3f871adca4554d648e3a75fa1e05fe" # 03.12.2014 EDW Kalib 75079_0002
-cdb$id <- "1c096311bab3b614ebe28117af1a1766" # 17.2.2015 EDW Kalib 75012_0002
-sigma.F13 <- 1.096
+## cdb$id <- "89e85b1c42ee6a9b70864a7612251a7e" # 14. & 15.10.2014 EDW Kalib 75093_0002
+## cdb$id <- "4a3f871adca4554d648e3a75fa1e05fe" # 03.12.2014 EDW Kalib 75079_0002
+## cdb$id <- "1c096311bab3b614ebe28117af1a1766" # 17.2.2015 EDW Kalib 75012_0002
+cdb$id <- "98b1736d2711ece64469fa99ea74bcc9" # 12.5.2015 EDW Kalib 75030_0002
+sigma.F13 <- 1.095
 doc    <- cdbGetDoc(cdb)$res
 doc    <- refreshAnalysis(cdb, doc)
 doc    <- refreshResult(cdb, doc)
@@ -119,6 +120,6 @@ doc$Calibration$Analysis$Values$Uncertainty <-
 doc <- ce3.uncertTotal(doc)
 doc <- dispatchResSum(doc)
 
-cdb$DBName     <- "vaclab_work" ## DBName
-cdb$dataList <- doc
-res <- cdbUpdateDoc(cdb)$res
\ No newline at end of file
+## cdb$DBName     <- "vaclab_db" ## DBName
+## cdb$dataList <- doc
+## res <- cdbUpdateDoc(cdb)$res
\ No newline at end of file

scripts/analyse_CE3_YAMP.R

@@ -17,7 +17,7 @@ if(length(doc$Calibration) > 0){
     doc <- ce3.writePind(doc)
     ##
     doc <- dispatchResCal(doc)
-
+    if(FALSE){
     ## uncertainty ...
     ## ... fm3 related
     doc <- fm3.uncertPfill(doc)
@@ -56,6 +56,6 @@ if(length(doc$Calibration) > 0){
     doc <- dispatchResSum(doc)
 #
     doc <- ce3.compareCDGs(doc)
-    if(FALSE){
+
     }
 }

scripts/cal_CE3_SZ-PARAM.R

@@ -24,8 +24,8 @@ if((a$cmscg == "N2" || a$cmscg == "Ar" || a$cmscg == "D2"|| a$cmscg == "H2") & i
     if(pcal < 1e-4){
         lw <- "lw1"
     }
-    #if(pcal < 9e-7){
-    if(pcal < 1e-5){ ## H2   
+    if(pcal < 9e-7){
+    #if(pcal < 1e-5){ ## H2   
         lw <- "lw0"
     }
     if(pcal < 1e-10){

uncert/getConductIndex.R

@@ -16,17 +16,12 @@ getConductIndex <- function(ccc){
     ## --- Lw2 --- kl. Lw ---
     lw2List  <- getSubList(a$cms, "useLw2")
     
-    iLw2     <- which((cnom > as.double(lw2List$From)) &
-                          (cnom < 8e-6)#as.double(lw2List$To))
-                      )
-    
+    iLw2     <- which((cnom > as.double(lw2List$From)) & as.double(lw2List$To))#    (cnom < 8e-6)
         
     ## --- Lw1 --- gr. Lw ---
     lw1List  <- getSubList(a$cms, "useLw1")
     
-    iLw1     <- which((cnom > as.double(lw1List$From)) &
-                          (cnom < 1e-4)#as.double(lw1List$To))
-                      )
+    iLw1     <- which((cnom > as.double(lw1List$From)) & as.double(lw1List$To)) # (cnom < 1e-4)
     
     ## --- LwC --- constLw ---
     lwCList    <- getSubList(a$cms, "useLwC")

utils/resSrgError.R

@@ -2,13 +2,16 @@ resSrgError <- function(ccc){
     msg <- "calculated by resSrgError"
     ## - data reduction as given by ToDo
     ## - round
-
+    ## 10% max. Abw. vom Zieldruck als default
+    maxdev <- 0.1 
     a      <- abbrevList(ccc)
+  
     p.target <- as.numeric(a$cpt$Values$Pressure$Value)
-    maxdev   <- as.numeric(a$cpt$MaxDev)
+    if(length(a$cpt$Values$Pressure$MaxDev) == 1){
+        maxdev   <- as.numeric(a$cpt$Values$Pressure$MaxDev)
+    }
+    
     k        <- 2 # Erweiterungsfaktor
-    ## 10% max. Abw. vom Zieldruck als default
-    if(length(maxdev) == 0) maxdev <- 0.30
     
     ##-------------##
     ## pcal
@@ -45,11 +48,12 @@ resSrgError <- function(ccc){
         i.out  <- NULL
         i.take <- which(pcal >  (p.target[i] *(1- maxdev)) &
                         pcal <  (p.target[i] *(1+ maxdev)))
-        print(itake)
+        print(pcal)
+        print(p.target)
         msg <- paste(msg,"; For target pressure:",p.target[i],
                      "I take the points:", toString(i.take))
 
-        if(length(i.take) > 1){
+        if(length(i.take) > 2){
             ## ut ist schon k=2, eswerden alle Punkte genommen,
             ## bei dem e.delta kleiner als 3-sigma ist
             ## wobei e.delte die Abweichung vom Reverenzwert ist