TITLE IND.MOD - An Indicator Dilution Model of the Circulation
: From "Mathematical and Computer Modeling of Physiological Systems"
:  By: Dr. Vincent Rideout
:  Adapted from ACSL to SCoP by Dr. Wesley M. Granger

PARAMETER
{
Fin = 0.5  (mg/sec)
PW = 0.5   (sec)
F = 100.   (ml/sec)
FVSD = 0.0 (ml/sec)
TP = 5.0   (sec)
TS = 8.0   (sec)
T0 = 1.25  (sec)

T1 = 2.50  (sec)
T2 = 1.25  (sec)
T3 = 1.25  (sec)
T4 = 7.50  (sec)
T5 = 2.00  (sec)
T6 = 10.0  (sec)
T7 = 1.25  (sec)
}

STATE
{
C1 FROM 0 TO 0.001    (mg/ml)
C2 FROM 0 TO 0.001    (mg/ml)
C3 FROM 0 TO 0.001    (mg/ml)
C4 FROM 0 TO 0.001    (mg/ml)

C5 FROM 0 TO 0.0002    (mg/ml)
C6 FROM 0 TO 0.0002    (mg/ml)
C7 FROM 0 TO 0.0002    (mg/ml)
C0 FROM 0 TO 0.0002    (mg/ml)
}

ASSIGNED
{
FD FROM 0 TO 0.005     (mg/sec)
FS FROM 0 TO 100       (ml/sec)
K7 FROM 0 TO 1       : Fraction of flow through systemic circ
K3 FROM 0 TO 1       : Fraction of flow through VSD


CP FROM 0 TO 0.001
CS FROM 0 TO 0.0002
}

INDEPENDENT
{
Time FROM 0 TO 60 WITH 200  (sec)
}

PLOT C1, CP, C5, CS VS Time

DERIVATIVE indicator
{
LAG C1 BY TP
LAG C5 BY TS
FS = F - FVSD
K7 = FS/F
K3 = FVSD/F
FD = pulse(Time, 0.0, Fin, PW)
C0' = ((K7*C7) + (K3*C3) - C0)/T0   : Right Ventricle
C1' = (C0 + (FD/F) - C1)/T1         : Pulmonary Artery
CP = lag_C1_TP
C2' = ((CP - C2)/T2)                : Pulmonary Veins, Left atrium
C3' = (C2 - C3)/T3                 : Left Ventricle
C4' = (C3 - C4)/T4                  : Aorta, Large Arteries
C5' = (C4 - C5)/T5                  : Small Arteries
CS = lag_C5_TS
C6' = (CS - C6)/T6                  : Small Veins
C7'  = (C6-C7)/T7                   : Large Veins, Right Atrium
}

BREAKPOINT
{
  SOLVE indicator METHOD runge
}