- uncert:

	geändert:   devel/calYampSE1.R
	neue Datei:   uncert/se1.uncertAds.R
	neue Datei:   uncert/se1.uncertRg.R
	neue Datei:   uncert/se1.uncertT1.R

devel/calYampSE1.R

@@ -31,6 +31,9 @@ doc <- se1.calPcal(doc)
 doc <- se1.uncertPfill(doc)
 doc <- se1.uncertf(doc)
 doc <- se1.uncertdT(doc)
+doc <- se1.uncertT1(doc)
+doc <- se1.uncertRg(doc)
+doc <- se1.uncertAds(doc)
 
 doc <- dispatchResCal( doc )
 

uncert/se1.uncertAds.R

+se1.uncertAds <- function(ccc){
+    msg   <- "calculated by  se1.uncertCorr()"
+    a     <- abbrevList(ccc)
+
+    pcal  <- getConstVal(a$cav, "cal") # in mbar
+    
+    ## -- Adsorption
+    UF2a <- getSubList(a$cms, "F2_u1_a") ## 
+
+    UF2b <- getSubList(a$cms, "F2_u1_b") ##
+    UF2  <- getSubList(a$cms, "F2_u2") ## 
+
+    uF2  <- rep(NA, length(pcal))
+
+    i.u1 <- which(pcal >  as.numeric(UF2a$From) &
+                  pcal <  as.numeric(UF2a$To))
+    i.u2 <- which(pcal >  as.numeric(UF2$From) &
+                  pcal <  as.numeric(UF2$To))
+
+    if(length(i.u1) > 0){
+        uF2[i.u1] <- as.numeric(UF2a$Value) + as.numeric(UF2a$Value) * pcal[i.u1] 
+    }
+
+    if(length(i.u2) > 0){
+        uF2[i.u2] <- as.numeric(UF2$Value)
+    }
+
+    ccc$Calibration$Analysis$Values$Uncertainty <-
+        setCcl(ccc$Calibration$Analysis$Values$Uncertainty,
+               "uncertAds",
+               "1",
+               uF2,
+               msg)
+
+    
+    return(ccc)
+}

uncert/se1.uncertRg.R

+se1.uncertRg <- function(ccc){
+    msg   <- "calculated by  se1.uncertCorr()"
+    a     <- abbrevList(ccc)
+
+
+    rg    <- getConstVal(a$cav, "rg")
+    pfill <- getConstVal(a$cav, "fill") # in mbar
+    pcal  <- getConstVal(a$cav, "cal") # in mbar
+    
+    ## --- Realgascorrektur
+    uf1  <- getConstVal(a$cms, "F1_u1") ## in 1/mbar
+    uF1  <- (1 - rg) * uf1 * pfill ## rel.
+    
+    ccc$Calibration$Analysis$Values$Uncertainty <-
+        setCcl(ccc$Calibration$Analysis$Values$Uncertainty,
+               "uncertRg",
+               "1",
+               uF1,
+               msg)
+    return(ccc)
+}
+

uncert/se1.uncertT1.R

+se1.uncertT1 <- function(ccc){
+    msg   <- "calculated by  se1.uncertT1()"
+    a     <- abbrevList(ccc)
+
+    Tafter       <- getConstVal(a$cav, "after")
+
+    ## Beiträge der Sensor-Kalib.:
+    us1 <- getConstVal(a$cmco, "keithley_u1") ## abs.
+    us2 <- getConstVal(a$cmco, "keithley_u2") ## abs.
+    us3 <- getConstVal(a$cmco, "keithley_u3") ## abs.
+    us4 <- getConstVal(a$cmco, "keithley_u4") ## abs.
+    us5 <- getConstVal(a$cmco, "keithley_u5") ## abs.
+    us6 <- getConstVal(a$cmco, "keithley_u6") ## abs.
+
+    usc <- sqrt(us1^2 + us2^2 + us3^2 + us4^2 + us5^2 + us6^2)
+
+    ## Beiträge des Normals
+    un1 <- getConstVal(a$cms, "T1_u1") ## abs.
+
+    ## total rel.
+    uT1 <- sqrt(un1^2 + usc^2)/Tafter  ## rel.
+    
+    ccc$Calibration$Analysis$Values$Uncertainty <-
+        setCcl(ccc$Calibration$Analysis$Values$Uncertainty,
+               "uncertT1",
+               "1",
+               uT1,
+               msg)
+
+    return(ccc)
+}