o2.c File Reference

Default parameters and helper functions for oxygen metabolism. More...

#include "libtpcmodel.h"

Go to the source code of this file.

Functions
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)
	Calculates K1/k2 ratio for [O-15]O2 in muscle, based on OER.
double	mo2pO2 (const double OER, const double K1k2, const double SaO2, const double p50Mb, const double cHb, const double cMb, const int verbose)

Variables
double	SaO2 =0.97
double	p50Hb =3.6
double	p50Mb =0.319
double	nHb =2.7
double	cHb =150.0
double	cMb =4.7

Detailed Description

Default parameters and helper functions for oxygen metabolism.

Author

Vesa Oikonen

Definition in file o2.c.

Function Documentation

mo2k1k2()

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 )

Calculates K1/k2 ratio for [O-15]O2 in muscle, based on OER.

Returns

Returns K1/k2.

Parameters

OER	Oxygen extraction ratio; assumptions hold only in the range OER=0.2-0.9
SaO2	Arterial oxygen saturation fraction; you can enter DEFAULT_SAO2
p50Hb	Half-saturation pressure p50 (kPa) for hemoglobin; you can enter DEFAULT_P50HB
p50Mb	Half-saturation pressure p50 (kPa) p50 for myoglobin; you can enter DEFAULT_P50MB
nHb	Hill coefficient n for hemoglobin; you can enter DEFAULT_NHB
cHb	Hemoglobin concentration in blood (mg/g); you can enter DEFAULT_CHB
cMb	Myoglobin concentration in muscle (mg/g); you can enter DEFAULT_CMB
verbose	Verbose level; if zero, then nothing is printed into stdout or stderr

Definition at line 29 of file o2.c.

  {
  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);
}

mo2pO2()

double mo2pO2	(	const double	OER,
		const double	K1k2,
		const double	SaO2,
		const double	p50Mb,
		const double	cHb,
		const double	cMb,
		const int	verbose )

Calculates the partial pressure of oxygen in muscle, based on OER and K1/k2.

Returns

Returns pO2 (kPa).

Parameters

OER	Oxygen extraction ratio; assumptions hold only in the range OER=0.2-0.9
K1k2	K1/k2 ratio for [O-15]O2
SaO2	Arterial oxygen saturation fraction; you can enter DEFAULT_SAO2
p50Mb	Half-saturation pressure p50 (kPa) p50 for myoglobin; you can enter DEFAULT_P50MB
cHb	Hemoglobin concentration in blood (mg/g); you can enter DEFAULT_CHB
cMb	Myoglobin concentration in muscle (mg/g); you can enter DEFAULT_CMB
verbose	Verbose level; if zero, then nothing is printed into stdout or stderr

Definition at line 83 of file o2.c.

  {
  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);
}

Variable Documentation

cHb

double cHb =150.0

Hemoglobin concentration in blood (mg/g)

Definition at line 20 of file o2.c.

Referenced by mo2k1k2(), and mo2pO2().

cMb

double cMb =4.7

Myoglobin concentration in muscle (mg/g)

Definition at line 22 of file o2.c.

Referenced by mo2k1k2(), and mo2pO2().

nHb

double nHb =2.7

Hill coefficient n for hemoglobin

Definition at line 18 of file o2.c.

Referenced by mo2k1k2().

p50Hb

double p50Hb =3.6

Half-saturation pressure p50 (kPa) for hemoglobin

Definition at line 14 of file o2.c.

Referenced by mo2k1k2().

p50Mb

double p50Mb =0.319

Half-saturation pressure p50 (kPa) p50 for myoglobin

Definition at line 16 of file o2.c.

Referenced by mo2k1k2(), and mo2pO2().

SaO2

double SaO2 =0.97

Arterial oxygen saturation fraction

Definition at line 12 of file o2.c.

Referenced by mo2k1k2(), and mo2pO2().