Grt-Bt: Kleemann, Reinhardt, 1994

Использовать Проверить

source("utilities_tpf.r");

sensors.t.grt_bt_kleemann_reinhardt_1994 <- function(Grt, Bt, P, ...){
    Bt <- tpf.recalc_minaral("BT",Bt);
    Grt <- tpf.recalc_minaral("GRT",Grt);

    XMF <-  Grt$Mg_pfu  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu);
    XFM <-  (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu)  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu);
    XMG <-  Grt$Mg_pfu  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu  +  Grt$Mn_pfu  +  Grt$Ca_pfu);
    XMN <-  Grt$Mn_pfu  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu  +  Grt$Mn_pfu  +  Grt$Ca_pfu);
    XFE <-  (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu)  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu  +  Grt$Mn_pfu  +  Grt$Ca_pfu);
    XCA <-  Grt$Ca_pfu  /  (Grt$Mg_pfu  +  Grt$Fe_p2_pfu + Grt$Fe_p3_pfu  +  Grt$Mn_pfu  +  Grt$Ca_pfu);
    MF <- Grt$Mg_pfu / (Grt$Fe_p2_pfu + Grt$Fe_p3_pfu);
    Z <- Bt$Al_c6_pfu;
    XAL <-  Bt$Al_c6_pfu  / (Bt$Mg_pfu  +  Bt$Fe_p2_pfu + Bt$Fe_p3_pfu  +  Bt$Mn_pfu  +  Bt$Al_c6_pfu  +  Bt$Ti_pfu);
    XBG <-  Bt$Mg_pfu  / (Bt$Mg_pfu  +  Bt$Fe_p2_pfu + Bt$Fe_p3_pfu  +  Bt$Mn_pfu  +  Bt$Al_c6_pfu  +  Bt$Ti_pfu);
    XTI <-  Bt$Ti_pfu  / (Bt$Mg_pfu  +  Bt$Fe_p2_pfu + Bt$Fe_p3_pfu  +  Bt$Mn_pfu  +  Bt$Al_c6_pfu  +  Bt$Ti_pfu);
    XBE <-  (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu)  / (Bt$Mg_pfu  +  Bt$Fe_p2_pfu + Bt$Fe_p3_pfu  +  Bt$Mn_pfu  +  Bt$Al_c6_pfu  +  Bt$Ti_pfu);
    BMF <- Bt$Mg_pfu / (Bt$Fe_p2_pfu + Bt$Fe_p3_pfu);
    VANN <-  (514.2  -  (14.3  *  Z))  * 0.3011;
    VPHL <-  (514.2  -  17.3  -  (14.3  *  Z)  +  (7.3  *  Z))  * 0.3011;
    DVCM <-  ((3  -  Z)  *  113.2)  +  (3  *  VANN)  -  ((3  -  Z)  *  115.11)  -  (3  *  VPHL);
    DVHW <-  (DVCM  *  0.1)  / 4.184 /  (3  *  (3  -  Z));
    FAK20 <-  (XCA ^ 2)  -  (2  *  XCA  *  XCA  *  XMG);
    FAK21 <-  (2  *  XCA  *  XMG)  -  (2  *  XCA  *  XMG  *  XMG);
    FAK22 <-   - (2  *  XCA  *  XCA  *  XFE);
    FAK23 <-   - (2  *  XCA  *  XFE  *  XFE);
    FAK24 <-  (2  *  XMG  *  XFE)  -  (2  *  XMG  *  XMG  *  XFE);
    FAK25 <-  (XFE ^ 2)  -  (2  *  XMG  *  XFE  *  XFE);
    FAK26 <-  (XCA  *  XFE)  -  (2  *  XCA  *  XMG  *  XFE);
    FAK27 <-  (XCA  *  XMN)  -  (2  *  XCA  *  XMG  *  XMN);
    FAK28 <-   - (2  *  XCA  *  XFE  *  XMN);
    FAK29 <-  (XFE  *  XMN)  -  (2  *  XMG  *  XFE  *  XMN);
    WHPSUM <-  (5153  *  FAK20)  +  (16539  *  FAK21)  +  (4857  *  FAK22)   +  (626  *  FAK23)  +  (55  *  FAK24)  +  (889  *  FAK25)   +  (14060  *  FAK26)   +  (10857  *  FAK27)  +  (2741  *  FAK28)  +  (472  *  FAK29);
    WSPSUM <-  (4.49  *  FAK20)  +  (4.49  *  FAK21)  +  (1.21  *  FAK22)   +  (1.21  *  FAK23)  +  (5.71  *  FAK26)  +  (4.49  *  FAK27)   +  (1.21  *  FAK28);
    WVPSUM <-  (0.024  *  FAK20)  +  (0.024  *  FAK21)  +  (0.041  *  FAK22)   +  (0.022  *  FAK23)  +  (0.002  *  FAK24)  +  (0.014  *  FAK25)   +  (0.063  *  FAK26)   +  (0.024  *  FAK27)  +  (0.031  *  FAK28)  +  (0.008000001  *  FAK29);
    FAK30 <-   - (2  *  XCA  *  XCA  *  XMG);
    FAK31 <-   - (2  *  XCA  *  XMG  *  XMG);
    FAK32 <-  (XCA ^ 2)  -  (2  *  XCA  *  XCA  *  XFE);
    FAK33 <-  (2  *  XCA  *  XFE)  -  (2  *  XCA  *  XFE  *  XFE);
    FAK34 <-  (XMG ^ 2)  -  (2  *  XMG  *  XMG  *  XFE);
    FAK35 <-  (2  *  XMG  *  XFE)  -  (2  *  XMG  *  XFE  *  XFE);
    FAK36 <-  (XCA  *  XMG)  -  (2  *  XCA  *  XMG  *  XFE);
    FAK37 <-   - (2  *  XCA  *  XMG  *  XMN);
    FAK38 <-  (XCA  *  XMN)  -  (2  *  XCA  *  XFE  *  XMN);
    FAK39 <-  (XMG  *  XMN)  -  (2  *  XMG  *  XFE  *  XMN);
    WHASUM <-  (5153  *  FAK30)  +  (16539  *  FAK31)  +  (4857  *  FAK32)   +  (626  *  FAK33)  +  (55  *  FAK34)  +  (889  *  FAK35)   +  (14060  *  FAK36)   +  (10857  *  FAK37)  +  (2741  *  FAK38)  +  (472  *  FAK39);
    WSASUM <-  (4.49  *  FAK30)  +  (4.49  *  FAK31)  +  (1.21  *  FAK32)   +  (1.21  *  FAK33)  +  (5.71  *  FAK36)  +  (4.49  *  FAK37)   +  (1.21  *  FAK38);
    WVASUM <-  (0.024  *  FAK30)  +  (0.024  *  FAK31)  +  (0.041  *  FAK32)   +  (0.022  *  FAK33)  +  (0.002  *  FAK34)  +  (0.014  *  FAK35)   +  (0.063  *  FAK36)   +  (0.024  *  FAK37)  +  (0.031  *  FAK38)  +  (0.008000001  *  FAK39);
    NENFAK <-  (1.9872  * log(MF  /  BMF))  -  ((22.5  *  XAL)  -  (2.8  *  XTI))            +  ((WSASUM  -  WSPSUM)  /  3);
    T <-  (4840.5  +  ((18590   *  XAL)  -  (4335   *  XTI))  -      ((WHASUM  -  WHPSUM)  /  3))  / (2.548  -  NENFAK)  +      ((DVHW  -  ((WVASUM  -  WVPSUM)  /  3))  / (2.548  -  NENFAK) *  P * 1000)  -  273.15;
    return(data.frame(Name=paste0(Grt$Name,"-", Bt$Name),P=P,T_calc=T));
}

sensors.t.grt_bt_kleemann_reinhardt_1994.args <- function(){
    ars <- list();
    ars$minerals <- c("Grt", "Bt");
    ars$ranges <- c("P");
    return(ars);
}

#sensors.t.grt_bt_kleemann_reinhardt_1994.test <- function(){
#
#}