//Remember to set CVODE to MAXSTATES-1 length Assume for now, for //now that I can init states and dstates, pass that memory to the //integrators and state 46 if in algebraic mode. -Modification to use seperate vectors? for mito/force only, how to modify? -> linear vectors with state offsets (I.e. subset vectors) -Where do the 'extra' charges from the alg model get changed -VFO is a calculated but unused variable? -> Move to output? (Commented Out) //VAD_AT //VcnsASP VAAT VACO VANT VATPase VCS VFH VHFe Vhleak VHNe Vhu VIDH VKGDH VMDH VnaCa VNO VSDH VSL Vuni double O1_RyR; double 1000_F_Acap_C_m; double Co; double Istim; double E_Na; double E_K; double E_Ks; double E_Ca; double INa; double G_Na; double IKs; double IK1; double G_Nab; double INab; double G_Kp; double IKp; double ICahalf; double ICamax; double PK; double ICaK; double ICa; double INaK; double INaKmax; double KmNai; double KmKo; double Ko; double Ki1AD_NaK; double Km1AT_NaK; double eta; double Nao; double Cao; double KmCa; double ksat; // double Fcik1; double KmNa; double G_Cab; double ICab; double IpCa; double Km2ATP_CaP; double KiADP_CaP; double Km1ATP_CaP; double KmpCa; double IpCamax; double PnsNa; double KmnsCa; double PnsK; double InsCa; double V_AM_scaler; double f_01; double f_12; double f_23; double V_AM; //Mitochondia double ATPm; double Cm; double DmuH; double DpH; double CPN; double NAD; double VCS; double KmOaa; double AcCoA; double KmAcCoA; double EtCS; double KCS; double VACO; double CIK; double KACOeq; double kfACO; double kKGDH; double EtKG; double KmKG; double Mg; double Kmg; double Kca; double KmKGNAD; double nKG; double VKGDH; double CoA; double kfSL; double KSLeq; double VSL; double kSDH; double EtSDH; double KmSucc; double KiFUM; double KiOxaa; double VSDH; double kfFH; double KFHeq; double VFH; double Kh1; double Kh2; double Kh3; double Kh4; double Koff; double H; double Kioaa; double KmmNAD; double Kmal; double EtMD; double kMDH; double VMDH; double kfAAT; double GLU; double kcnsASP; double KAATeq; double VAAT; double kres; double rhoREN; double Dpsio; double g; double ra; double rc1; double r1; double rc2; double rb; double VNO; double VHNe; double kresf; double FADH2; double FAD; double r2; double r3; double VFO; double VHFe; double kf1; double Pi; double pa; double pc1; double pc2; double VATPase; double Vhu; double gh; double VmDT; double hm; double VANT; double Vhleak; //Junction Paramerters & Vars double Kfb; double Krb; double Nfb; double KSR; double vmaxf; double vmaxr; double KmATP_SR; double Ki_SR; double Ki_prime_SR; double v1; double tauxfer; double tautr; double Jup; double Jrel; double Jtr; double Jxfer; //force Parameters double alpha_SL; double zeta; double FN_Ca; double Fnorm; double Force; double kTrop_pn; double LTRPNtot; double Ktrop_half; double Ntrop; double f_01; double g_01_mod; double f_12; double g_12_mod; double g_23_mod; double f_23; double kltrpn_plus; double LTRPNtot; double kltrpn_minus; //Junction Params double CMDNtot; double KmCMDN; double CSQNtot; double KmCSQN; double Jtrpn; double beta_SS; double beta_JSR; double beta_i; //Vuni double Vmuni; double ktrans; double L; double kact; double na; double Vuni; //VnaCa double b; double VmNC; double Kna; double n; double Knca; double VnaCa; //Mitochondia //Derivative Calculation Parameters double InsNa; double Vmyo; double Acap; double InsK; double VSS; double VNSR; double VJSR; double C_m; double kaplus; double kbplus; double kaplus; double kaminus; double kbminus; double kcminus; double fprime; double gprime; //Stimulation modes double time_on_Is1; double time_off_Is1; double time_on_Is2; double time_off_Is2; double pulse_amplitude; double shift; double t1; double t2; double high_freq; double period; double norm_freq; double high_freq_hz; double norm_freq_hz; double pulse_duration; double start_time_shift; double start_time_shift_time_on; double time_vclamp_on; double vclamp_set; double vclamp_hold;