o2.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpcmodel.h"
/*****************************************************************************/
 
/*****************************************************************************/
/* Default variables for oxygen */
double SaO2=0.97;
double p50Hb=3.6;
double p50Mb=0.319;
double nHb=2.7;
double cHb=150.0;
double cMb=4.7;
/*****************************************************************************/
 
/*****************************************************************************/
double mo2k1k2(
  const double OER,
  const double SaO2,
  const double p50Hb,
  const double p50Mb,
  const double nHb,
  const double cHb,
  const double cMb,
  const int verbose
) {
  if(verbose>0) printf("mo2k1k2()\n");
  if(verbose>1) {
    printf("input:\n");
    printf("  OER := %g\n", OER);
    printf("  SaO2 := %g\n", SaO2);
    printf("  p50Hb := %g\n", p50Hb);
    printf("  p50Mb := %g\n", p50Mb);
    printf("  nHb := %g\n", nHb);
    printf("  cHb := %g\n", cHb);
    printf("  cMb := %g\n", cMb);
  }
  double SHb=(1.0-OER)*SaO2;
  if(verbose>2) printf("SHb := %g\n", SHb);
 
  double pO2=p50Hb*pow(SHb/(1.0-SHb), 1.0/nHb);
  if(verbose>2) printf("pO2 := %g\n", pO2);
 
  double SMb=pO2/(pO2+p50Mb);
  if(verbose>2) printf("SMb := %g\n", SMb);
 
  double rO2=cMb/cHb;
  if(verbose>2) printf("rO2 := %g\n", rO2);
 
  double k1k2=rO2*(SMb/SHb);
  if(verbose>1) printf("K1/k2 := %g\n", k1k2);
 
  return(k1k2);
}
/*****************************************************************************/
 
/*****************************************************************************/
double mo2pO2(
  const double OER,
  const double K1k2,
  const double SaO2,
  const double p50Mb,
  const double cHb,
  const double cMb,
  const int verbose
) {
  if(verbose>0) printf("mo2pO2()\n");
  if(verbose>1) {
    printf("input:\n");
    printf("  OER := %g\n", OER);
    printf("  K1/k2 := %g\n", K1k2);
    printf("  SaO2 := %g\n", SaO2);
    printf("  cHb := %g\n", cHb);
    printf("  cMb := %g\n", cMb);
    printf("  p50Mb := %g\n", p50Mb);
  }
  double SHb=(1.0-OER)*SaO2;
  if(verbose>2) printf("SHb := %g\n", SHb);
 
  double rO2=cMb/cHb;
  if(verbose>2) printf("rO2 := %g\n", rO2);
 
  double SMb=K1k2*(1.0-OER)*SaO2/rO2;
  if(verbose>2) printf("SMb := %g\n", SMb);
 
  double pO2=SMb*p50Mb/(1.0-SMb);
  if(verbose>1) printf("pO2 := %g\n", pO2);
 
  return(pO2);
}
/*****************************************************************************/
 
/*****************************************************************************/