====== Amp-Pl: Holland, Blundy, 1994 ======

For edenite + 4 quartz = tremolite + albite (A) and edenite + albite = richterite + anorthite (B) equilibrium.

=== Details ===

^ Article | <ref:holland1994> |
^ Accuracy | ± 50 °C  |
^ Ranges | 1 < P < 16 kbar,  400 °C < T < 1100 °C |
^ Renew sensors | Blundy and Holland (1990),  Kohn and Spear (1989, 1990).  |
<use:sensors.pt.amp_pl_holland_blundy>
<test:sensors.pt.amp_pl_holland_blundy>

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


=== Source ===

<Code linenums lang-r>
source('sensors_p_amp_pl_holland_blundy.r');
source('sensors_t_amp_pl_holland_blundy.r');

sensors.pt.amp_pl_holland_blundy <- function(Amp, Pl, Temp, usePFU = FALSE, ...){

stopifnot(usePFU == F);
P_A <- 2;
Temp_A <- 700;
P_B <- 2;
Temp_B <- 700;


interations <- 1;

repeat{
    P_A_DF <- sensors.p.amp_pl_holland_blundy(Amp, Temp = Temp_A);
    Temp_A_DF <- sensors.t.amp_pl_holland_blundy(Amp, Pl, P = P_A);

    P_B_DF <- sensors.p.amp_pl_holland_blundy(Amp, Temp = Temp_B);
    Temp_B_DF <- sensors.t.amp_pl_holland_blundy(Amp, Pl, P = P_B);
    
    P_A <- P_A_DF$P_calc_Holland_Blundy;
    Temp_A <- Temp_A_DF$T_calc_Holland_Blundy_A;

    P_B <- P_B_DF$P_calc_Holland_Blundy;
    Temp_B <- Temp_B_DF$T_calc_Holland_Blundy_B;
    
    interations <- interations + 1;
    if (interations > 15) {break;}
}

return (data.frame(Name=paste0(Amp$Name,'-',Pl$Name),P_calc_Holland_Blundy_A=P_A,T_calc_Holland_Blundy_A=Temp_A,P_calc_Holland_Blundy_B=P_B,T_calc_Holland_Blundy_B=Temp_B));

}

sensors.pt.amp_pl_holland_blundy.args <- function() {
    ars <- list();
    ars$minerals = c('Amp','Pl');
    ars$ranges = c();
    return(ars);
}
</Code>