Cpx: Putirka nomenative procedure

source('recalc_px_putirka.r');

interpmin.cpx_putirka <- function (wtp_data, ...)
{
p <- recalc.px_putirka(wtp_data);
comps <- data.frame(Name=p$Name);

comps$X_Jd <- ifelse(p$Na_pfu<p$Al_c6_pfu, p$Na_pfu, p$Al_c6_pfu);
comps$X_CaTs <- ifelse(p$Al_c6_pfu>p$Na_pfu, p$Al_c6_pfu - p$Na_pfu, 0);
comps$X_CaTi <- ifelse(p$Al_c4_pfu>comps$X_CaTs, (p$Al_c4_pfu-comps$X_CaTs)/2, 0);
comps$X_CrCaTs <- p$Cr_pfu /2;
comps$X_DiHd <- p$Ca_pfu - comps$X_CaTs - comps$X_CaTi - comps$X_CrCaTs;
comps$X_EnFs <- (p$Fe_pfu+p$Mg_pfu-comps$X_DiHd)/2;
comps$X_En <- comps$X_EnFs*(p$Mg_pfu/(p$Mg_pfu+p$Mn_pfu+p$Fe_pfu));
comps$X_Di <- comps$X_DiHd*(p$Mg_pfu/(p$Mg_pfu+p$Mn_pfu+p$Fe_pfu));

return(comps);
}