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);
    
    cat <- c('Sb_apt','Fe_apt','As_apt','Zn_apt');    
    fhlcat <- selectNames(FhlA, cat);
    spfcat <- selectNames(SpfA, c('Fe_apt','Zn_apt'));
    

    X_Fhl_As = fhlcat$As_apt / (fhlcat$Sb_apt + fhlcat$Fe_apt + fhlcat$As_apt + fhlcat$Zn_apt);
    X_Fhl_Fe = fhlcat$Fe_apt / (fhlcat$Sb_apt + fhlcat$Fe_apt + fhlcat$As_apt + fhlcat$Zn_apt);
    X_Fhl_Zn = fhlcat$Zn_apt / (fhlcat$Sb_apt + fhlcat$Fe_apt + fhlcat$As_apt + fhlcat$Zn_apt);
    X_Sph_FeS = spfcat$Fe_apt / (spfcat$Fe_apt + spfcat$Zn_apt);
    X_Sph_ZnS = spfcat$Zn_apt / (spfcat$Fe_apt + spfcat$Zn_apt);
    
    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, X_Fhl_Fe =X_Fhl_Fe));
     return(cbind(FhlA,fhlcat));
}

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