====== Electrum-Sp: Shikazono, 1985 ======

Electrum – sphalerite termometer and fugometer.

=== Details ===

^ Article | <ref:shikazono1985> |
^ Collibration | In this study, sphalerite and electrum from Sanru, Chitose, Sado, Takatama, Seigoshi, Yugashima, Toi, Nebazawa, Takeno, and Kushikino were analyzed by an electron microprobe analyzer  |
^ Ranges | N_Ag > 0.3 |
^ Input | atomic wtp |
^ Accuracy | T ±20 °C (for T > 300 °C)  |
<use:sensors.t.electrum_sp_shikazono>
<test:sensors.t.electrum_sp_shikazono>

<WRAP center round todo 60%>
Для улучшения этой статьи следует:
  * Добавить подробное описание
</WRAP>


=== Source ===

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

sensors.t.electrum_sp_shikazono <- function(Elect, Spf)
{
    aElect <- recalc.wtp_at_to_atp(Elect);
    aSpf <- recalc.wtp_at_to_atp(Spf);

    dElect <- selectNames(aElect,c('Au_atp','Ag_atp'));
    dSpf <- selectNames(aSpf,c('Fe_atp','Zn_atp'));
 
    N_Ag = dElect$Ag_atp/(dElect$Ag_atp+dElect$Au_atp);
    X_FeS = dSpf$Fe_atp/(dSpf$Fe_atp+dSpf$Zn_atp);

    N = 28765 + 22600*(1 - N_Ag)^2 - 6400*(1 - N_Ag)^3;
    D = 49.008 - 9.152*log10(X_FeS) + 18.2961*log10(N_Ag)+ 5.5*(1 - N_Ag)^2;
    Tcalc6 = N / D - 273.15;

    N = 25196 + 22600*(1 - N_Ag)^2 - 6400*(1 - N_Ag)^3;
    D = 44.38 - 9.152*log10(X_FeS) + 18.296*log10(N_Ag)+ 5.5*(1 - N_Ag)^2;
    Tcalc7 = N / D - 273.15;

    Temp = Tcalc7 + 273.15;
    
    log_f_S = (1/(4.576*Temp))*
              (-41980 + 16.52*Temp - 18.296*Temp*log10(N_Ag)+
              (4*(1 - N_Ag)^2*(5650 - 1600*(1 - N_Ag) - 1.375*Temp))); 

    return(data.frame(Name=paste0(Elect$Name,'-',Spf$Name), N_Ag=N_Ag, X_FeS =X_FeS,
           T_calc_ElectSp_Shikazono_eq6 = Tcalc6,
           T_calc_ElectSp_Shikazono_eq7 = Tcalc7,
           log_f_S=log_f_S));
}

sensors.t.electrum_sp_shikazono.args <- function() {
    ars <- list();
    ars$minerals = c('Elect','Sp');
    ars$ranges = c();
    return(ars);
}

sensors.t.electrum_sp_shikazono.test <- function () {
  
  Elect <- data.frame(Name='Elect',Au_wtp=2, Ag_wtp=8);
  Spf <- data.frame(Name='Sp',Fe_wtp=21.23,Zn_wtp=21.23);
 
  result <- sensors.t.electrum_sp_shikazono(Elect, Spf);
  
  expect_equal(result$T_calc_ElectSp_Shikazono_eq6, 312, tolerance = 3);
  expect_equal(result$T_calc_ElectSp_Shikazono_eq7, 297, tolerance = 3);
}
</Code>