====== Fhl-Sp (unpublished) ======

Fhl-Spf – termometer. Для блеклых руд на месторождениях представленных малосеребристыми разновидностями, с использованием уравнение Р.О.Сэка.

=== Details ===

<WRAP center round important 60%>
**Неопубликованный инструмент**

Данный инструмент не опубликован в научной литературе.
</WRAP>

^ Article |   |
^ Input | Elements wt%  |
^ Reference |   |
^ Collibration |  |
^ Accuracy |  |
^[[:ppy:mod_sensors#t_fhl_spf_unpubl|Calculate]] ^ ^

=== Source ===

<Code linenums lang-r>
source('recalc_wtp_at_to_atp.r')

sensors.t.fhl_spf_unpubl <- function(Fhl, Spf)
{
    FhlA = recalc.wtp_at_to_atp(Fhl);
    SpfA = recalc.wtp_at_to_atp(Spf);
    
    fhlcat <- selectNames(FhlA, c('Cu_atp', 'Sb_atp','Fe_atp','As_atp','Zn_atp'));
    spfcat <- selectNames(SpfA, c('Fe_atp','Zn_atp'));
    
    X_Fhl_As <- fhlcat$As_atp / (fhlcat$Sb_atp + fhlcat$As_atp);
    X_Fhl_Fe <- fhlcat$Fe_atp / (fhlcat$Fe_atp + fhlcat$Zn_atp);
    X_Fhl_Zn <- fhlcat$Zn_atp / (fhlcat$Fe_atp + fhlcat$Zn_atp);
    X_Sph_FeS <- spfcat$Fe_atp / (spfcat$Fe_atp + spfcat$Zn_atp);
    X_Sph_ZnS <- spfcat$Zn_atp / (spfcat$Fe_atp + spfcat$Zn_atp);
    
    R = 1.987;
    
    Tcalc = (2071.8+1294.9*X_Fhl_As) / (R*(
        (0.7285-0.9186*X_Sph_FeS - 0.5295*X_Sph_FeS^2 - 0.1772*X_Sph_FeS^3)*
        (1-2*X_Sph_FeS)+
        X_Fhl_Fe/X_Fhl_Zn*X_Sph_ZnS/X_Sph_FeS)
        );

    return(data.frame(Name=paste0(Fhl$Name,'-',Spf$Name), T_calc_fhl_spf = Tcalc-273, X_Fhl_As =X_Fhl_As, X_Fhl_Fe =X_Fhl_Fe, 
                                                                                  X_Fhl_Zn =X_Fhl_Zn, X_Sph_FeS=X_Sph_FeS, X_Sph_ZnS=X_Sph_ZnS));

}

sensors.t.fhl_spf_unpubl.args <- function() {
    ars <- list();
    ars$minerals = c('Fhl','Spf');
    ars$ranges = c();
    return(ars);
}

</Code>