//Created by libAntimony v2.4
model heldt_circ_pbpk_2002__Circ_Environment(time_, HR, hrf)

  // Assignment Rules:
  hrf := HR / 0.06;

  // Variable initializations:
  time_ = ;
  HR = 77;
end

model heldt_circ_pbpk_2002__LATiming(realtime, HR, beattime, hrf)

  // Assignment Rules:
  beattime := (realtime - floor(realtime / hrf) * hrf) / 1000;

  // Variable initializations:
  realtime = ;
  HR = ;
  hrf = ;
end

model heldt_circ_pbpk_2002__LAElastanceFunction(time_, hrf, E_LA)

  // Assignment Rules:
  Ts := TsK * root(hrf / 1000);
  E_LA := piecewise( Edia + ((Esys - Edia) * (1 - cos((pi * time_) / Ts))) / 2 , ( geq(time_, 0)) && (time_ <= Ts ), Edia + ((Esys - Edia) * (1 + cos((2 * pi * (time_ - Ts)) / Ts))) / 2 , ( time_ < 1.5 * Ts) && (geq(time_, Ts) ), Edia );

  // Variable initializations:
  time_ = ;
  Esys = 0.145;
  Edia = 0.131;
  TsK = 0.2;
  hrf = ;
end

model heldt_circ_pbpk_2002__LeftAtrium(time_, E_LA, P_PV, F_LV, V_LA, P_LA, F_LA, conc_X_LA, conc_X_PV)

  // Assignment Rules:
  P_LA := E_LA * (V_LA - V_LA_0) + P_LA_ext;
  F_LA := (P_PV - P_LA) / R_LA;
  conc_X_LA := sub_X_LA / V_LA;

  // Rate Rules:
  V_LA' = (F_LA - F_LV) / 1000;
  sub_X_LA' = (conc_X_PV * F_LA - conc_X_LA * F_LV) / 1000;

  // Variable initializations:
  time_ = ;
  E_LA = ;
  P_PV = ;
  F_LV = ;
  P_LA_ext = 0;
  V_LA = 144.038898121062;
  V_LA_0 = 15;
  R_LA = 0.01;
  conc_X_PV = ;
  sub_X_LA = 0;
end

model heldt_circ_pbpk_2002__LVTiming(realtime, HR, beattime, hrf)

  // Assignment Rules:
  beattime := (realtime - floor(realtime / hrf) * hrf) / 1000 - PRint;

  // Variable initializations:
  realtime = ;
  HR = ;
  PRint = 0.12;
  hrf = ;
end

model heldt_circ_pbpk_2002__LVElastanceFunction(time_, hrf, E_LV)

  // Assignment Rules:
  Ts := TsK * root(hrf / 1000);
  E_LV := piecewise( Edia + ((Esys - Edia) * (1 - cos((pi * time_) / Ts))) / 2 , ( geq(time_, 0)) && (time_ <= Ts ), Edia + ((Esys - Edia) * (1 + cos((2 * pi * (time_ - Ts)) / Ts))) / 2 , ( time_ < 1.5 * Ts) && (geq(time_, Ts) ), Edia );

  // Variable initializations:
  time_ = ;
  Esys = 5.6;
  Edia = 0.19;
  TsK = 0.35;
  hrf = ;
end

model heldt_circ_pbpk_2002__LeftVentricle(time_, E_LV, F_AO, F_CO, P_LA, V_LV, P_LV, F_LV, conc_X_LV, conc_X_LA)

  // Assignment Rules:
  P_LV := E_LV * (V_LV - V_LV_0) + pLVPext;
  F_LV := piecewise( (P_LA - P_LV) / R_LV , P_LA > P_LV , 0 );
  conc_X_LV := sub_X_LV / V_LV;
  fiber_stress := P_LV * (1 + (3 * V_LV) / 200);

  // Rate Rules:
  V_LV' = ((F_LV - F_AO) - F_CO) / 1000;
  sub_X_LV' = (conc_X_LA * F_LV - conc_X_LV * (F_AO + F_CO)) / 1000;

  // Variable initializations:
  time_ = ;
  E_LV = ;
  F_AO = ;
  F_CO = ;
  P_LA = ;
  pLVPext = 0;
  V_LV = 118.75315693072;
  V_LV_0 = 30;
  R_LV = 0.001;
  conc_X_LA = ;
  sub_X_LV = 0;
end

model heldt_circ_pbpk_2002__Aorta(time_, F_CR, P_LV, F_AR, V_AO, F_AO, P_AO, conc_X_AO, conc_X_LV)

  // Assignment Rules:
  F_AO := piecewise( F_AOalways , P_LV > P_AO , 0 );
  P_AO := (1 / C_AO) * (V_AO - V_AO_0);
  conc_X_AO := sub_X_AO / V_AO;

  // Rate Rules:
  V_AO' = ((F_AO - F_AR) - F_CR) / 1000;
  F_AOalways' = ((1 / L_AO) * (P_LV - P_AO) - (R_AO / L_AO) * F_AOalways) / 1000;
  sub_X_AO' = (conc_X_LV * F_AO - conc_X_AO * (F_CR + F_AR)) / 1000;

  // Variable initializations:
  time_ = ;
  F_CR = ;
  P_LV = ;
  F_AR = ;
  C_AO = 1.2;
  L_AO = 5e-5;
  R_AO = 0.01;
  V_AO_0 = 200;
  V_AO = 302.714548324312;
  F_AOalways = -6898.97886842467;
  conc_X_LV = ;
  sub_X_AO = 0;
end

model heldt_circ_pbpk_2002__CoronaryCirc(time_, P_LV, P_RA, F_CO, F_COV, V_CO, P_CO, conc_X_CO, conc_X_LV)

  // Assignment Rules:
  F_CO := piecewise( (P_LV - P_CO) / R_CO , P_LV > P_CO , 0 );
  F_COV := (P_CO - P_RA) / R_COV;
  P_CO := (1 / C_CO) * (V_CO - V_CO_0);
  conc_X_CO := sub_X_CO / V_CO;

  // Rate Rules:
  V_CO' = (F_CO - F_COV) / 1000;
  sub_X_CO' = (conc_X_LV * F_CO - conc_X_CO * F_COV) / 1000;

  // Variable initializations:
  time_ = ;
  P_LV = ;
  P_RA = ;
  R_CO = 1;
  R_COV = 0.2;
  V_CO_0 = 10;
  C_CO = 0.1;
  V_CO = 11.4012164759732;
  conc_X_LV = ;
  sub_X_CO = 0;
end

model heldt_circ_pbpk_2002__CarotidCirc(time_, P_AO, P_VC, F_CR, F_CRV, V_CR, P_CR, conc_X_CR, conc_X_AO)

  // Assignment Rules:
  F_CR := (P_AO - P_CR) / R_CR;
  F_CRV := (P_CR - P_VC) / R_CRV;
  P_CR := (1 / C_CR) * (V_CR - V_CR_0);
  conc_X_CR := sub_X_CR / V_CR;
  conc_X_CRtissue := sub_X_CRtissue / V_CRtissue;

  // Rate Rules:
  V_CR' = (F_CR - F_CRV) / 1000;
  sub_X_CR' = ((conc_X_AO * F_CR - conc_X_CR * F_CRV) - (conc_X_CR - conc_X_CRtissue) * D_CRtissue) / 1000;
  sub_X_CRtissue' = ((conc_X_CR - conc_X_CRtissue) * D_CRtissue) / 1000;

  // Variable initializations:
  time_ = ;
  P_AO = ;
  P_VC = ;
  R_CR = 5;
  R_CRV = 0.02;
  V_CR_0 = 10;
  C_CR = 1;
  V_CR = 27.2147903053333;
  conc_X_AO = ;
  sub_X_CR = 0;
  D_CRtissue = 1;
  V_CRtissue = 1;
  sub_X_CRtissue = 0;
end

model heldt_circ_pbpk_2002__Arteries(time_, F_AD, F_MU, F_GI, F_LI, F_KI, F_OT, F_SK, P_AO, V_AR, F_AR, P_AR, conc_X_AR, conc_X_AO)

  // Assignment Rules:
  F_AR := (P_AO - P_AR) / R_AR;
  P_AR := (1 / C_AR) * (V_AR - V_AR_0);
  conc_X_AR := sub_X_AR / V_AR;

  // Rate Rules:
  V_AR' = (((((((F_AR - F_AD) - F_MU) - F_GI) - F_LI) - F_KI) - F_OT) - F_SK) / 1000;
  sub_X_AR' = (conc_X_AO * F_AR - conc_X_AR * (F_SK + F_AD + F_MU + F_GI + F_LI + F_KI + F_OT)) / 1000;

  // Variable initializations:
  time_ = ;
  F_AD = ;
  F_MU = ;
  F_GI = ;
  F_LI = ;
  F_KI = ;
  F_OT = ;
  F_SK = ;
  P_AO = ;
  R_AR = 1.5;
  C_AR = 0.5;
  V_AR_0 = 370;
  V_AR = 383.557125288364;
  conc_X_AO = ;
  sub_X_AR = 0;
end

model heldt_circ_pbpk_2002__Adipose(time_, P_AR, P_VE, V_AD, F_AD, F_ADV, P_AD, conc_X_AD, conc_X_AR)

  // Assignment Rules:
  F_AD := (P_AR - P_AD) / R_AD;
  F_ADV := (P_AD - P_VE) / R_ADV;
  P_AD := (1 / C_AD) * (V_AD - V_AD_0);
  conc_X_AD := sub_X_AD / V_AD;

  // Rate Rules:
  V_AD' = (F_AD - F_ADV) / 1000;
  sub_X_AD' = (conc_X_AR * F_AD - conc_X_AD * F_ADV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_AD = 1.81;
  R_AD = 1;
  R_ADV = 0.33;
  V_AD_0 = 75;
  V_AD = 111.639088723426;
  conc_X_AR = ;
  sub_X_AD = 0;
end

model heldt_circ_pbpk_2002__Muscle(time_, P_AR, P_VE, V_MU, F_MU, F_MUV, P_MU, conc_X_MU, conc_X_AR)

  // Assignment Rules:
  F_MU := (P_AR - P_MU) / R_MU;
  F_MUV := (P_MU - P_VE) / R_MUV;
  P_MU := (1 / C_MU) * (V_MU - V_MU_0);
  conc_X_MU := sub_X_MU / V_MU;

  // Rate Rules:
  V_MU' = (F_MU - F_MUV) / 1000;
  sub_X_MU' = (conc_X_AR * F_MU - conc_X_MU * F_MUV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_MU = 1.81;
  R_MU = 1;
  R_MUV = 0.33;
  V_MU_0 = 75;
  V_MU = 111.639088723426;
  conc_X_AR = ;
  sub_X_MU = 0;
end

model heldt_circ_pbpk_2002__GutIntestine(time_, P_AR, P_LI, V_GI, F_GI, F_GIV, P_GI, conc_X_GI, conc_X_AR)

  // Assignment Rules:
  F_GI := (P_AR - P_GI) / R_GI;
  F_GIV := (P_GI - P_LI) / R_GIV;
  P_GI := (1 / C_GI) * (V_GI - V_GI_0);
  conc_X_GI := sub_X_GI / V_GI;

  // Rate Rules:
  V_GI' = (F_GI - F_GIV) / 1000;
  sub_X_GI' = (conc_X_AR * F_GI - conc_X_GI * F_GIV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_LI = ;
  C_GI = 1.81;
  R_GI = 1;
  R_GIV = 0.33;
  V_GI_0 = 75;
  V_GI = 111.639088723426;
  conc_X_AR = ;
  sub_X_GI = 0;
end

model heldt_circ_pbpk_2002__Liver(time_, P_AR, P_VE, V_LI, F_LI, F_GIV, F_LIV, P_LI, conc_X_LI, conc_X_AR, conc_X_GI)

  // Assignment Rules:
  F_LI := (P_AR - P_LI) / R_LI;
  F_LIV := (P_LI - P_VE) / R_LIV;
  P_LI := (1 / C_LI) * (V_LI - V_LI_0);
  conc_X_LI := sub_X_LI / V_LI;

  // Rate Rules:
  V_LI' = ((F_GIV + F_LI) - F_LIV) / 1000;
  sub_X_LI' = (((conc_X_GI * F_GIV + conc_X_AR * F_LI) - conc_X_LI * F_LIV) + conc_X_LI * metabolism) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_LI = 1.81;
  R_LI = 1;
  R_LIV = 0.33;
  V_LI_0 = 75;
  V_LI = 111.639088723426;
  F_GIV = ;
  conc_X_AR = ;
  conc_X_GI = ;
  sub_X_LI = 0;
  metabolism = 0;
end

model heldt_circ_pbpk_2002__Kidney(time_, P_AR, P_VE, V_KI, F_KI, F_KIV, P_KI, conc_X_KI, conc_X_AR)

  // Assignment Rules:
  F_KI := (P_AR - P_KI) / R_KI;
  F_KIV := (P_KI - P_VE) / R_KIV;
  P_KI := (1 / C_KI) * (V_KI - V_KI_0);
  conc_X_KI := sub_X_KI / V_KI;

  // Rate Rules:
  V_KI' = (F_KI - F_KIV) / 1000;
  sub_X_KI' = (conc_X_AR * F_KI - conc_X_KI * F_KIV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_KI = 1.81;
  R_KI = 1;
  V_KI_0 = 75;
  R_KIV = 0.33;
  V_KI = 111.639088723426;
  conc_X_AR = ;
  sub_X_KI = 0;
end

model heldt_circ_pbpk_2002__Skin(time_, P_AR, P_VE, V_SK, F_SK, F_SKV, P_SK, conc_X_SK, conc_X_AR)

  // Assignment Rules:
  F_SK := (P_AR - P_SK) / R_SK;
  F_SKV := (P_SK - P_VE) / R_SKV;
  P_SK := (1 / C_SK) * (V_SK - V_SK_0);
  conc_X_SK := sub_X_SK / V_SK;

  // Rate Rules:
  V_SK' = (F_SK - F_SKV) / 1000;
  sub_X_SK' = (conc_X_AR * F_SK - conc_X_SK * F_SKV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_SK = 1.81;
  R_SK = 1;
  R_SKV = 0.33;
  V_SK_0 = 75;
  V_SK = 111.639088723426;
  conc_X_AR = ;
  sub_X_SK = 0;
end

model heldt_circ_pbpk_2002__OtherTissue(time_, P_AR, P_VE, V_OT, F_OT, F_OTV, P_OT, conc_X_OT, conc_X_AR)

  // Assignment Rules:
  F_OT := (P_AR - P_OT) / R_OT;
  F_OTV := (P_OT - P_VE) / R_OTV;
  P_OT := (1 / C_OT) * (V_OT - V_OT_0);
  conc_X_OT := sub_X_OT / V_OT;

  // Rate Rules:
  V_OT' = (F_OT - F_OTV) / 1000;
  sub_X_OT' = (conc_X_AR * F_OT - conc_X_OT * F_OTV) / 1000;

  // Variable initializations:
  time_ = ;
  P_AR = ;
  P_VE = ;
  C_OT = 1.81;
  R_OT = 1;
  R_OTV = 0.33;
  V_OT_0 = 75;
  V_OT = 111.639088723426;
  conc_X_AR = ;
  sub_X_OT = 0;
end

model heldt_circ_pbpk_2002__Veins(time_, F_ADV, F_MUV, F_LIV, F_KIV, F_OTV, F_SKV, P_VC, V_VE, F_VE, F_VEV, P_VE, conc_X_VE, conc_X_AD, conc_X_MU, conc_X_LI, conc_X_KI, conc_X_SK, conc_X_OT)

  // Assignment Rules:
  F_VE := F_SKV + F_ADV + F_MUV + F_LIV + F_KIV + F_OTV;
  F_VEV := (P_VE - P_VC) / R_VEV;
  P_VE := (1 / C_VE) * (V_VE - V_VE_0);
  conc_X_VE := sub_X_VE / V_VE;

  // Rate Rules:
  V_VE' = (F_VE - F_VEV) / 1000;
  sub_X_VE' = ((conc_X_SK * F_SKV + conc_X_AD * F_ADV + conc_X_MU * F_MUV + conc_X_LI * F_LIV + conc_X_KI * F_KIV + conc_X_OT * F_OTV) - conc_X_VE * F_VEV) / 1000;

  // Variable initializations:
  time_ = ;
  F_ADV = ;
  F_MUV = ;
  F_LIV = ;
  F_KIV = ;
  F_OTV = ;
  F_SKV = ;
  P_VC = ;
  R_VEV = 0.0223;
  V_VE = 833.291684587657;
  C_VE = 13.24;
  V_VE_0 = 596;
  conc_X_AD = ;
  conc_X_MU = ;
  conc_X_LI = ;
  conc_X_KI = ;
  conc_X_SK = ;
  conc_X_OT = ;
  sub_X_VE = 5;
end

model heldt_circ_pbpk_2002__VenaCava(time_, F_VEV, P_RA, F_CRV, V_VC, F_VC, F_VCV, P_VC, conc_X_VC, conc_X_VE, conc_X_CR)

  // Assignment Rules:
  F_VC := F_CRV + F_VEV;
  F_VCV := (P_VC - P_RA) / R_VCV;
  P_VC := (1 / C_VC) * (V_VC - V_VC_0);
  conc_X_VC := sub_X_VC / V_VC;

  // Rate Rules:
  V_VC' = (F_VC - F_VCV) / 1000;
  sub_X_VC' = ((conc_X_VE * F_VEV + conc_X_CR * F_CRV) - conc_X_VC * F_VCV) / 1000;

  // Variable initializations:
  time_ = ;
  F_VEV = ;
  P_RA = ;
  F_CRV = ;
  R_VCV = 0.0267;
  V_VC = 1849.93788284599;
  V_VC_0 = 600;
  C_VC = 73.88;
  conc_X_VE = ;
  conc_X_CR = ;
  sub_X_VC = 0;
end

model heldt_circ_pbpk_2002__RATiming(realtime, HR, beattime, hrf)

  // Assignment Rules:
  beattime := (realtime - floor(realtime / hrf) * hrf) / 1000;

  // Variable initializations:
  realtime = ;
  HR = ;
  hrf = ;
end

model heldt_circ_pbpk_2002__RAElastanceFunction(time_, hrf, E_RA)

  // Assignment Rules:
  Ts := TsK * root(hrf / 1000);
  E_RA := piecewise( Edia + ((Esys - Edia) * (1 - cos((pi * time_) / Ts))) / 2 , ( geq(time_, 0)) && (time_ <= Ts ), Edia + ((Esys - Edia) * (1 + cos((2 * pi * (time_ - Ts)) / Ts))) / 2 , ( time_ < 1.5 * Ts) && (geq(time_, Ts) ), Edia );

  // Variable initializations:
  time_ = ;
  Esys = 0.11;
  Edia = 0.099;
  TsK = 0.2;
  hrf = ;
end

model heldt_circ_pbpk_2002__RightAtrium(time_, E_RA, F_RV, F_VCV, F_COV, V_RA, P_RA, F_RA, conc_X_RA, conc_X_VC, conc_X_CO)

  // Assignment Rules:
  P_RA := E_RA * (V_RA - V_RA_0) + P_RA_ext;
  F_RA := F_VCV + F_COV;
  conc_X_RA := sub_X_RA / V_RA;

  // Rate Rules:
  V_RA' = (F_RA - F_RV) / 1000;
  sub_X_RA' = ((conc_X_VC * F_VCV + conc_X_CO * F_COV) - conc_X_RA * F_RV) / 1000;

  // Variable initializations:
  time_ = ;
  E_RA = ;
  F_RV = ;
  F_VCV = ;
  F_COV = ;
  P_RA_ext = 0;
  V_RA = 152.150765432928;
  V_RA_0 = 15;
  conc_X_VC = ;
  conc_X_CO = ;
  sub_X_RA = 0;
end

model heldt_circ_pbpk_2002__RVTiming(realtime, HR, beattime, hrf)

  // Assignment Rules:
  beattime := (realtime - floor(realtime / hrf) * hrf) / 1000 - PRinterval;

  // Variable initializations:
  realtime = ;
  PRinterval = 0.12;
  HR = ;
  hrf = ;
end

model heldt_circ_pbpk_2002__RVElastanceFunction(time_, hrf, E_RV)

  // Assignment Rules:
  Ts := TsK * root(hrf / 1000);
  E_RV := piecewise( Edia + ((Esys - Edia) * (1 - cos((pi * time_) / Ts))) / 2 , ( geq(time_, 0)) && (time_ <= Ts ), Edia + ((Esys - Edia) * (1 + cos((2 * pi * (time_ - Ts)) / Ts))) / 2 , ( time_ < 1.5 * Ts) && (geq(time_, Ts) ), Edia );

  // Variable initializations:
  time_ = ;
  Esys = 0.67;
  Edia = 0.104;
  TsK = 0.35;
  hrf = ;
end

model heldt_circ_pbpk_2002__RightVentricle(time_, E_RV, P_RA, F_PA, V_RV, P_RV, F_RV, conc_X_RV, conc_X_RA)

  // Assignment Rules:
  P_RV := E_RV * (V_RV - V_RV_0) + P_RV_ext;
  F_RV := piecewise( (P_RA - P_RV) / R_RV , P_RA > P_RV , 0 );
  conc_X_RV := sub_X_RV / V_RV;

  // Rate Rules:
  V_RV' = (F_RV - F_PA) / 1000;
  sub_X_RV' = (conc_X_RA * F_RV - conc_X_RV * F_PA) / 1000;

  // Variable initializations:
  time_ = ;
  E_RV = ;
  P_RA = ;
  F_PA = ;
  P_RV_ext = 0;
  V_RV = 154.337916417774;
  V_RV_0 = 25;
  R_RV = 0.002;
  conc_X_RA = ;
  sub_X_RV = 0;
end

model heldt_circ_pbpk_2002__PulmonaryArtery(time_, F_Pa, P_RV, F_PA, V_PA, P_PA, conc_X_PA, conc_X_RV)

  // Assignment Rules:
  F_PA := piecewise( (P_RV - P_PA) / R_PA , P_RV > P_PA , 0 );
  P_PA := (1 / C_PA) * (V_PA - V_PA_0);
  conc_X_PA := sub_X_PA / V_PA;

  // Rate Rules:
  V_PA' = (F_PA - F_Pa) / 1000;
  sub_X_PA' = (conc_X_RV * F_PA - conc_X_PA * F_Pa) / 1000;

  // Variable initializations:
  time_ = ;
  F_Pa = ;
  P_RV = ;
  R_PA = 0.0227;
  C_PA = 2.222;
  V_PA_0 = 50;
  V_PA = 112.059036699839;
  conc_X_RV = ;
  sub_X_PA = 0;
end

model heldt_circ_pbpk_2002__PulmonaryArteries(time_, P_PA, F_PC, F_SH, V_Pa, F_Pa, P_Pa, conc_X_Pa, conc_X_PA)

  // Assignment Rules:
  F_Pa := (P_PA - P_Pa) / R_Pa;
  P_Pa := (1 / C_Pa) * (V_Pa - V_Pa_0);
  conc_X_Pa := sub_X_Pa / V_Pa;

  // Rate Rules:
  V_Pa' = ((F_Pa - F_PC) - F_SH) / 1000;
  sub_X_Pa' = (conc_X_PA * F_Pa - conc_X_Pa * (F_PC + F_SH)) / 1000;

  // Variable initializations:
  time_ = ;
  P_PA = ;
  F_PC = ;
  F_SH = ;
  V_Pa = 66.9910537203609;
  C_Pa = 1.481;
  R_Pa = 0.053;
  V_Pa_0 = 30;
  conc_X_PA = ;
  sub_X_Pa = 0;
end

model heldt_circ_pbpk_2002__Shunt(time_, F_SH, P_Pa, P_PC)

  // Assignment Rules:
  F_SH := (P_Pa - P_PC) / R_SH;

  // Variable initializations:
  time_ = ;
  R_SH = 2;
  P_Pa = ;
  P_PC = ;
end

model heldt_circ_pbpk_2002__PulmonaryCapillaries(time_, P_Pa, P_PV, V_PC, F_PC, F_PCV, P_PC, conc_X_PC, conc_X_Pa)

  // Assignment Rules:
  F_PC := (P_Pa - P_PC) / R_PC;
  F_PCV := (P_PC - P_PV) / R_PCV;
  P_PC := (1 / C_PC) * (V_PC - V_PC_0);
  conc_X_PC := sub_X_PC / V_PC;

  // Rate Rules:
  V_PC' = (F_PC - F_PCV) / 1000;
  sub_X_PC' = (conc_X_Pa * F_PC - conc_X_PC * F_PCV) / 1000;

  // Variable initializations:
  time_ = ;
  P_Pa = ;
  P_PV = ;
  C_PC = 1.778;
  R_PC = 0.0379;
  R_PCV = 0.0379;
  V_PC_0 = 53;
  V_PC = 91.7784213531534;
  conc_X_Pa = ;
  sub_X_PC = 0;
end

model heldt_circ_pbpk_2002__PulmonaryVein(time_, F_PCV, F_SH, F_LA, V_PV, F_PV, P_PV, conc_X_PV, conc_X_PC, conc_X_Pa)

  // Assignment Rules:
  F_PV := F_PCV + F_SH;
  P_PV := (1 / C_PV) * (V_PV - V_PV_0);
  conc_X_PV := sub_X_PV / V_PV;

  // Rate Rules:
  V_PV' = (F_PV - F_LA) / 1000;
  sub_X_PV' = ((conc_X_PC * F_PCV + conc_X_Pa * F_SH) - conc_X_PV * F_LA) / 1000;

  // Variable initializations:
  time_ = ;
  F_PCV = ;
  F_SH = ;
  F_LA = ;
  C_PV = 5;
  V_PV_0 = 150;
  V_PV = 239.25897115602;
  conc_X_PC = ;
  conc_X_Pa = ;
  sub_X_PV = 0;
end

model *heldt_circ_pbpk_2002____main()

  // Sub-modules, and any changes to those submodules:
  Circ_Environment: heldt_circ_pbpk_2002__Circ_Environment(time_, HR, hrf);
  LATiming: heldt_circ_pbpk_2002__LATiming(time_, HR, beattime0, hrf);
  LAElastanceFunction: heldt_circ_pbpk_2002__LAElastanceFunction(beattime0, hrf, E_LA);
  LeftAtrium: heldt_circ_pbpk_2002__LeftAtrium(time_, E_LA, P_PV, F_LV, V_LA, P_LA, F_LA, conc_X_LA, conc_X_PV);
  LVTiming: heldt_circ_pbpk_2002__LVTiming(time_, HR, beattime, hrf);
  LVElastanceFunction: heldt_circ_pbpk_2002__LVElastanceFunction(beattime, hrf, E_LV);
  LeftVentricle: heldt_circ_pbpk_2002__LeftVentricle(time_, E_LV, F_AO, F_CO, P_LA, V_LV, P_LV, F_LV, conc_X_LV, conc_X_LA);
  Aorta: heldt_circ_pbpk_2002__Aorta(time_, F_CR, P_LV, F_AR, V_AO, F_AO, P_AO, conc_X_AO, conc_X_LV);
  CoronaryCirc: heldt_circ_pbpk_2002__CoronaryCirc(time_, P_LV, P_RA, F_CO, F_COV, V_CO, P_CO, conc_X_CO, conc_X_LV);
  CarotidCirc: heldt_circ_pbpk_2002__CarotidCirc(time_, P_AO, P_VC, F_CR, F_CRV, V_CR, P_CR, conc_X_CR, conc_X_AO);
  Arteries: heldt_circ_pbpk_2002__Arteries(time_, F_AD, F_MU, F_GI, F_LI, F_KI, F_OT, F_SK, P_AO, V_AR, F_AR, P_AR, conc_X_AR, conc_X_AO);
  Adipose: heldt_circ_pbpk_2002__Adipose(time_, P_AR, P_VE, V_AD, F_AD, F_ADV, P_AD, conc_X_AD, conc_X_AR);
  Muscle: heldt_circ_pbpk_2002__Muscle(time_, P_AR, P_VE, V_MU, F_MU, F_MUV, P_MU, conc_X_MU, conc_X_AR);
  GutIntestine: heldt_circ_pbpk_2002__GutIntestine(time_, P_AR, P_LI, V_GI, F_GI, F_GIV, P_GI, conc_X_GI, conc_X_AR);
  Liver: heldt_circ_pbpk_2002__Liver(time_, P_AR, P_VE, V_LI, F_LI, F_GIV, F_LIV, P_LI, conc_X_LI, conc_X_AR, conc_X_GI);
  Kidney: heldt_circ_pbpk_2002__Kidney(time_, P_AR, P_VE, V_KI, F_KI, F_KIV, P_KI, conc_X_KI, conc_X_AR);
  Skin: heldt_circ_pbpk_2002__Skin(time_, P_AR, P_VE, V_SK, F_SK, F_SKV, P_SK, conc_X_SK, conc_X_AR);
  OtherTissue: heldt_circ_pbpk_2002__OtherTissue(time_, P_AR, P_VE, V_OT, F_OT, F_OTV, P_OT, conc_X_OT, conc_X_AR);
  Veins: heldt_circ_pbpk_2002__Veins(time_, F_ADV, F_MUV, F_LIV, F_KIV, F_OTV, F_SKV, P_VC, V_VE, F_VE, F_VEV, P_VE, conc_X_VE, conc_X_AD, conc_X_MU, conc_X_LI, conc_X_KI, conc_X_SK, conc_X_OT);
  VenaCava: heldt_circ_pbpk_2002__VenaCava(time_, F_VEV, P_RA, F_CRV, V_VC, F_VC, F_VCV, P_VC, conc_X_VC, conc_X_VE, conc_X_CR);
  RATiming: heldt_circ_pbpk_2002__RATiming(time_, HR, time_0, hrf);
  RAElastanceFunction: heldt_circ_pbpk_2002__RAElastanceFunction(time_0, hrf, E_RA);
  RightAtrium: heldt_circ_pbpk_2002__RightAtrium(time_, E_RA, F_RV, F_VCV, F_COV, V_RA, P_RA, F_RA, conc_X_RA, conc_X_VC, conc_X_CO);
  RVTiming: heldt_circ_pbpk_2002__RVTiming(time_, HR, time_1, hrf);
  RVElastanceFunction: heldt_circ_pbpk_2002__RVElastanceFunction(time_1, hrf, E_RV);
  RightVentricle: heldt_circ_pbpk_2002__RightVentricle(time_, E_RV, P_RA, F_PA, V_RV, P_RV, F_RV, conc_X_RV, conc_X_RA);
  PulmonaryArtery: heldt_circ_pbpk_2002__PulmonaryArtery(time_, F_Pa, P_RV, F_PA, V_PA, P_PA, conc_X_PA, conc_X_RV);
  PulmonaryArteries: heldt_circ_pbpk_2002__PulmonaryArteries(time_, P_PA, F_PC, F_SH, V_Pa, F_Pa, P_Pa, conc_X_Pa, conc_X_PA);
  Shunt: heldt_circ_pbpk_2002__Shunt(time_, F_SH, P_Pa, P_PC);
  PulmonaryCapillaries: heldt_circ_pbpk_2002__PulmonaryCapillaries(time_, P_Pa, P_PV, V_PC, F_PC, F_PCV, P_PC, conc_X_PC, conc_X_Pa);
  PulmonaryVein: heldt_circ_pbpk_2002__PulmonaryVein(time_, F_PCV, F_SH, F_LA, V_PV, F_PV, P_PV, conc_X_PV, conc_X_PC, conc_X_Pa);
end