Fhl-Sp (unpublished)

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

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('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$Fe_atp + fhlcat$As_atp + fhlcat$Zn_atp);
    X_Fhl_Fe <- fhlcat$Fe_atp / (fhlcat$Sb_atp + fhlcat$Fe_atp + fhlcat$As_atp + fhlcat$Zn_atp);
    X_Fhl_Zn <- fhlcat$Zn_atp / (fhlcat$Sb_atp + fhlcat$Fe_atp + fhlcat$As_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, X_Fhl_As =X_Fhl_As));
 
}
 
sensors.t.fhl_spf_unpubl.args <- function() {
    ars <- list();
    ars$minerals = c('Fhl','Spf');
    ars$ranges = c();
    return(ars);
}