PARCDL_rateconstants_main.py

import numpy as np


def PARCDL_rateconstants_main(Vc,Ve,Vm,Vn):

    # going to return (kD,kC,kA,kR,kP, \
        # kRP1,kRP2,kRP3,kRP4,kRP5,kRP6,kRP7,kRP8,kRP9,kRP10,kRP11,kRP12,kRP13,kRP14,kRP15,kRP16,kRP17,kRP18,kRP19,kRP20,kRP21,kRP22,kRP23,kRP24,kRP25,kRP26,kRP27,kRP28,kRP29,kRP30,kRP31,kRP32,kRP33,kRP34, \
        # kDP,kPA)


    mpc2nmcf_Vc=1E9/(Vc*6.023E23)
    mpc2nmcf_Vm=1E9/(Vm*6.023E23)
    mpc2nmcf_Vn=1E9/(Vn*6.023E23)

    kD = np.zeros(shape=(59,1))

    # kD **
    kD[0,0]=0.9*(1000/3600)
    kD[1,0]=5/(1000*3600)
    kD[2,0]=2/3600;
    kD[3,0]=10/3600;
    kD[4,0]=4;
    kD[5,0]=1*1000;
    kD[6,0]=1.4/(1000*3600) ;
    kD[7,0]=0.14/(1000*3600);
    kD[8,0]=0.9/3600;
    kD[9,0]=0.2*(1000/3600);
    kD[10,0]=1/3600;
    kD[11,0]=0.5/(1000*3600);
    kD[12,0]=0.25/3600;
    kD[13,0]=0.7/3600;
    kD[14,0]=0.1/(1000*3600);
    kD[15,0]=50/3600;
    kD[16,0]=4;
    kD[17,0]=0.2*1000;
    kD[18,0]=7.5/3600;
    kD[19,0]=10*(1000/3600);
    kD[20,0]=4.5*(1000/3600);
    kD[21,0]=0.033948*3600;
    kD[22,0]=0.058614*3600;
    kD[23,0]=0.0035778*3600;
    kD[24,0]=0.33705*3600;
    kD[25,0]=0.013496*3600;
    kD[26,0]=0.030939*3600;
    kD[27,0]=0.0042462*3600;
    kD[28,0]=0.081109*3600;
    kD[29,0]=0.10093*3600;
    kD[30,0]=0.13295*3600;
    kD[31,0]=0.080377*3600;
    kD[32,0]=0.12247*3600;
    kD[33,0]=0.16522*3600;
    kD[34,0]=0.099235*3600;
    kD[35,0]=0.13629*3600;
    kD[36,0]=0.28505*3600;
    kD[37,0]=0.10575*3600;
    kD[38,0]=0.066299*3600;
    kD[39,0]=0.11092*3600;
    kD[40,0]=0.042434*3600;
    kD[41,0]=0;
    kD[42,0]=0;
    kD[43,0]=0;
    kD[44,0]=0;
    kD[45,0]=1E-5*5;
    kD[46,0]=1E-5;
    kD[47,0]=1E-5;
    kD[48,0]=0.0005;
    kD[49,0]=0;
    kD[50,0]=(1E-4)*50;
    kD[51,0]=5000;
    kD[52,0]=0;
    kD[53,0]=5;
    kD[54,0]=20;
    kD[55,0]=20;
    kD[56,0]=20;
    kD[57,0]=2;
    kD[58,0]=2;

    # These taus (kD24-kD43) were predicted with a genetic algorithm in matlab with a
    # function tolerance of 1e-20, 1000 generations, attempting to minimize sum
    # of all conditions (DSB only, UV only, and DSB+UV together).

    kC = np.zeros(shape=(173,1))

    kC[0,0]=0.2*100;
    kC[1,0]=0.11*100;
    kC[2,0]=0.25*100;
    kC[3,0]=0.5*100;
    kC[4,0]=2*100;
    kC[5,0]=1.5*100;
    kC[6,0]=0.5*100;
    kC[7,0]=17;
    kC[8,0]=0.75*100;
    kC[9,0]=0.5*100;
    kC[10,0]=0.5*100;
    kC[11,0]=0.5*100;
    kC[12,0]=0.1*100;
    kC[13,0]=0.1*100;
    kC[14,0]=0.1*100;
    kC[15,0]=0.01*100;
    kC[16,0]=0.5*100;
    kC[17,0]=0.1*100;
    kC[18,0]=0.1*100;
    kC[19,0]=5*100;
    kC[20,0]=0.5*100;
    kC[21,0]=0.1*100;
    kC[22,0]=0.1*100;
    kC[23,0]=0.1*100;
    kC[24,0]=0.01*100;
    kC[25,0]=0.5*100;
    kC[26,0]=0.1*100;
    kC[27,0]=0.1*100;
    kC[28,0]=5*100;
    kC[29,0]=0.5*100;
    kC[30,0]=0.1*100;
    kC[31,0]=0.1*100;
    kC[32,0]=0.1*100;
    kC[33,0]=0.1*100;
    kC[34,0]=0.1*100;
    kC[35,0]=0.1*100;
    kC[36,0]=0.1*100;
    kC[37,0]=0.1*100;
    kC[38,0]=0.1*100;
    kC[39,0]=0.1*100;
    kC[40,0]=0.1*100;
    kC[41,0]=0.1*100;
    kC[42,0]=2*100;
    kC[43,0]=0.15/(100*3600);
    kC[44,0]=0.05/3600;
    kC[45,0]=0.2/(100*3600);
    kC[46,0]=0.1/3600;
    kC[47,0]=0.15/(100*3600);
    kC[48,0]=0.125/3600;
    kC[49,0]=0.12/(100*3600);
    kC[50,0]=0.2/3600;
    kC[51,0]=0.15/3600;
    kC[52,0]=0.2/(100*3600);
    kC[53,0]=0.1/3600;
    kC[54,0]=0.15/(100*3600);
    kC[55,0]=0.125/3600;
    kC[56,0]=0.12/(100*3600);
    kC[57,0]=0.2/3600;
    kC[58,0]=0.15/(100*3600);
    kC[59,0]=0.0375/3600;
    kC[60,0]=0.005/3600;
    kC[61,0]=0.4*100;
    kC[62,0]=0.29/3600;
    kC[63,0]=0.175/(100*3600);
    kC[64,0]=0.2/(100*3600);
    kC[65,0]=0.2/(100*3600);
    kC[66,0]=0.25/(100*3600);
    kC[67,0]=1.1*100;
    kC[68,0]=0.005*100;
    kC[69,0]=0.2*100;
    kC[70,0]=0.1*100;
    kC[71,0]=0.05/3600;
    kC[72,0]=0.14/3600;
    kC[73,0]=0.1/3600;
    kC[74,0]=0.2/3600;
    kC[75,0]=2*100;
    kC[76,0]=0.1*100;
    kC[77,0]=0.005/3600;
    kC[78,0]=0.005/3600;
    kC[79,0]=0.005/3600;
    kC[80,0]=0.005/3600;
    kC[81,0]=0.06/3600;
    kC[82,0]=0.06/3600;
    kC[83,0]=0.01/3600;
    kC[84,0]=0.005/3600;
    kC[85,0]=0.1*100;
    kC[86,0]=0.002/3600;
    kC[87,0]=1.1/3600;
    kC[88,0]=0.1/3600;
    kC[89,0]=0.1/3600;
    kC[90,0]=0.1/3600;
    kC[91,0]=0.1/3600;
    kC[92,0]=0.1*100;
    kC[93,0]=0.1*100;
    kC[94,0]=0.075/3600;
    kC[95,0]=0.15/3600;
    kC[96,0]=0.1/3600;
    kC[97,0]=0.2/3600;
    kC[98,0]=0.075/3600;
    kC[99,0]=0.15/3600;
    kC[100,0]=0.01/3600;
    kC[101,0]=0.06/3600;
    kC[102,0]=0.04/3600;
    kC[103,0]=0.5*100;
    kC[104,0]=0.1/3600;
    kC[105,0]=0.2/3600;
    kC[106,0]=2*100;
    kC[107,0]=0.1*100;
    kC[108,0]=2*100;
    kC[109,0]=0.1*100;
    kC[110,0]=2*100;
    kC[111,0]=0.1*100;
    kC[112,0]=2*100;
    kC[113,0]=0.1*100;
    kC[114,0]=0.05/(100*3600);
    kC[115,0]=0.5/3600;
    kC[116,0]=0.025/(100*3600);
    kC[117,0]=0.5/3600;
    kC[118,0]=2.2/3600;
    kC[119,0]=1.25*(100/3600);
    kC[120,0]=4/3600;
    kC[121,0]=4/3600;
    kC[122,0]=100/3600;
    kC[123,0]=2*(100/3600);
    kC[124,0]=8/3600;
    kC[125,0]=0.1*(100/3600);
    kC[126,0]=0.75*(100/3600);
    kC[127,0]=0.1*(100/3600);
    kC[128,0]=1/3600;
    kC[129,0]=2*(100/3600);
    kC[130,0]=1*(100/3600);
    kC[131,0]=1*(100/3600);
    kC[132,0]=0.8*(100/3600);
    kC[133,0]=0.05*(100/3600);
    kC[134,0]=2.5*(100/3600);
    kC[135,0]=0.06*(100/3600);
    kC[136,0]=5*(100/3600);
    kC[137,0]=3*(100/3600);
    kC[138,0]=5*(100/3600);
    kC[139,0]=1.1*(100/3600);
    kC[140,0]=2/3600;
    kC[141,0]=3.9/3600;
    kC[142,0]=1*(100/3600);
    kC[143,0]=2/3600;
    kC[144,0]=1.85/3600;
    kC[145,0]=2.1/3600;
    kC[146,0]=0.2*(100/3600);
    kC[147,0]=1.4/3600;
    kC[148,0]=8/3600;
    kC[149,0]=0.7*(100/3600);
    kC[150,0]=4/3600;
    kC[151,0]=5/3600;
    kC[152,0]=5/3600;
    kC[153,0]=1.2/3600;
    kC[154,0]=1*(100/3600);
    kC[155,0]=0.8*(100/3600);
    kC[156,0]=0.1/3600;
    kC[157,0]=0.1*(100/3600);
    kC[158,0]=0.11*(100/3600);
    kC[159,0]=0.15*(100/3600);
    kC[160,0]=0.105*(100/3600);
    kC[161,0]=0.12*(100/3600);
    kC[162,0]=0.13*(100/3600);
    kC[163,0]=0.8*(100/3600);
    kC[164,0]=0.15*(100/3600);
    kC[165,0]=0.06*(100/3600);
    kC[166,0]=1;
    kC[167,0]=1/100;
    kC[168,0]=1;
    kC[169,0]=1/100;
    kC[170,0]=1;
    kC[171,0]=1/100;
    kC[172,0]=0.4/3600;


    kA = np.zeros(shape=(87,1))


    kA[0,0]=4E-7/mpc2nmcf_Vc;
    kA[1,0]=1E-3;
    kA[2,0]=1E-5;
    kA[3,0]=1E-6/mpc2nmcf_Vc;
    kA[4,0]=1E-3;
    kA[5,0]=1E-6/mpc2nmcf_Vc;
    kA[6,0]=1E-3;
    kA[7,0]=1;
    kA[8,0]=1E-6/mpc2nmcf_Vc;
    kA[9,0]=1E-3;
    kA[10,0]=1E-7/mpc2nmcf_Vc;
    kA[11,0]=1E-3;
    kA[12,0]=1;
    kA[13,0]=1E-6/mpc2nmcf_Vc;
    kA[14,0]=1E-3;
    kA[15,0]=1;
    kA[16,0]=3E-8/mpc2nmcf_Vc;
    kA[17,0]=1E-3;
    kA[18,0]=1;
    kA[19,0]=2E-6/mpc2nmcf_Vc;
    kA[20,0]=1E-3;
    kA[21,0]=0.1;
    kA[22,0]=1E-6/mpc2nmcf_Vc;
    kA[23,0]=1E-2;
    kA[24,0]=1;
    kA[25,0]=1E-7/mpc2nmcf_Vc;
    kA[26,0]=1E-3;
    kA[27,0]=1;
    kA[28,0]=1E-6/mpc2nmcf_Vc;
    kA[29,0]=1E-3;
    kA[30,0]=1E-7/mpc2nmcf_Vc;
    kA[31,0]=1E-3;
    kA[32,0]=1;
    kA[33,0]=.01;
    kA[34,0]=.01;
    kA[35,0]=1E-6/mpc2nmcf_Vc;
    kA[36,0]=1E-3;
    kA[37,0]=1E-6/mpc2nmcf_Vc;
    kA[38,0]=1E-3;
    kA[39,0]=1E-6/mpc2nmcf_Vc;
    kA[40,0]=1E-3;
    kA[41,0]=1E-6/mpc2nmcf_Vc;
    kA[42,0]=1E-3;
    kA[43,0]=1E-6/mpc2nmcf_Vc;
    kA[44,0]=1E-3;
    kA[45,0]=1E-6/mpc2nmcf_Vc;
    kA[46,0]=1E-3;
    kA[47,0]=1;
    kA[48,0]=2E-6/mpc2nmcf_Vc;
    kA[49,0]=1E-3;
    kA[50,0]=10;
    kA[51,0]=2E-6/mpc2nmcf_Vc;
    kA[52,0]=1E-3;
    kA[53,0]=10;
    kA[54,0]=.01;
    kA[55,0]=.01;
    kA[56,0]=5E-7/mpc2nmcf_Vc;
    kA[57,0]=1E-3;
    kA[58,0]=1;
    kA[59,0]=5E-8/mpc2nmcf_Vc;
    kA[60,0]=1E-3;
    kA[61,0]=5E-9/mpc2nmcf_Vc;
    kA[62,0]=1E-3;
    kA[63,0]=1;
    kA[64,0]=.01;
    kA[65,0]=.01;
    kA[66,0]=2E-6/mpc2nmcf_Vc;
    kA[67,0]=1E-3;
    kA[68,0]=7E-6/mpc2nmcf_Vc;
    kA[69,0]=1E-3;
    kA[70,0]=1E-6/mpc2nmcf_Vc;
    kA[71,0]=1E-3;
    kA[72,0]=1E-6/mpc2nmcf_Vc;
    kA[73,0]=1E-3;
    kA[74,0]=1E-6/mpc2nmcf_Vc;
    kA[75,0]=1E-3;
    kA[76,0]=1E-7/mpc2nmcf_Vc;
    kA[77,0]=1E-3;
    kA[78,0]=1;
    kA[79,0]=1E-6/mpc2nmcf_Vc;
    kA[80,0]=1E-3;
    kA[81,0]=.01;
    kA[82,0]=.01;
    kA[83,0]=.01;
    kA[84,0]=0.1;
    kA[85,0]=0.1;
    kA[86,0]=0;


    kR = np.zeros(shape=(158,1))


    kR[0,0]=0.001;
    kR[1,0]=0.001;
    kR[2,0]=0.01;
    kR[3,0]=0.001;
    kR[4,0]=0.01;
    kR[5,0]=0.01;
    kR[6,0]=0.001;
    kR[7,0]=0.01;
    kR[8,0]=0.01;
    kR[9,0]=0.01;
    kR[10,0]=0.01;
    kR[11,0]=0.001;
    kR[12,0]=0.001;
    kR[13,0]=0.01;
    kR[14,0]=0.01;
    kR[15,0]=0.01;
    kR[16,0]=0.01;
    kR[17,0]=0.01;
    kR[18,0]=0.01;
    kR[19,0]=0.01;
    kR[20,0]=0.001;
    kR[21,0]=0.01;
    kR[22,0]=0.01;
    kR[23,0]=0.01;
    kR[24,0]=0.001;
    kR[25,0]=0.001;
    kR[26,0]=0.01;
    kR[27,0]=0.01;
    kR[28,0]=0.01;
    kR[29,0]=0.01;
    kR[30,0]=0.001;
    kR[31,0]=0.01;
    kR[32,0]=0.01;
    kR[33,0]=0.01;
    kR[34,0]=0.001;
    kR[35,0]=0.01;
    kR[36,0]=0.01;
    kR[37,0]=0.01;
    kR[38,0]=0.01;
    kR[39,0]=0.01;
    kR[40,0]=0.01;
    kR[41,0]=0.01;
    kR[42,0]=0.001;
    kR[43,0]=0.01;
    kR[44,0]=0.01;
    kR[45,0]=0.01;
    kR[46,0]=0.001;
    kR[47,0]=0.001;
    kR[48,0]=0.01;
    kR[49,0]=0.01;
    kR[50,0]=0.01;
    kR[51,0]=0.01;
    kR[52,0]=0.001;
    kR[53,0]=0.01;
    kR[54,0]=0.01;
    kR[55,0]=0.01;
    kR[56,0]=0.001;
    kR[57,0]=0.01;
    kR[58,0]=0.01;
    kR[59,0]=0.01;
    kR[60,0]=0.01;
    kR[61,0]=0.01;
    kR[62,0]=0.001;
    kR[63,0]=0.01;
    kR[64,0]=0.01;
    kR[65,0]=0.01;
    kR[66,0]=0.001;
    kR[67,0]=0.1;
    kR[68,0]=0.001;
    kR[69,0]=0.1;
    kR[70,0]=0.001;
    kR[71,0]=0.1;
    kR[72,0]=0.05;
    kR[73,0]=0.0016;
    kR[74,0]=0.01;
    kR[75,0]=0.01;
    kR[76,0]=0.05;
    kR[77,0]=0.0016;
    kR[78,0]=0.01;
    kR[79,0]=0.01;
    kR[80,0]=0.05;
    kR[81,0]=0.05;
    kR[82,0]=0.01;
    kR[83,0]=0.0001;
    kR[84,0]=1;
    kR[85,0]=0.001;
    kR[86,0]=0.001;
    kR[87,0]=0.00002;
    kR[88,0]=0.01;
    kR[89,0]=0.001;
    kR[90,0]=0.002;
    kR[91,0]=0.001;
    kR[92,0]=0.001;
    kR[93,0]=0.01;
    kR[94,0]=0.01;
    kR[95,0]=0;
    kR[96,0]=0.01;
    kR[97,0]=0;
    kR[98,0]=0.1;
    kR[99,0]=1;
    kR[100,0]=0.1;
    kR[101,0]=1;
    kR[102,0]=0.1;
    kR[103,0]=1;
    kR[104,0]=0.001;
    kR[105,0]=0.1*10;
    kR[106,0]=0.001;
    kR[107,0]=0.1*10;
    kR[108,0]=0.001;
    kR[109,0]=0.1*10;
    kR[110,0]=0.001;
    kR[111,0]=0.1*10;
    kR[112,0]=0.001;
    kR[113,0]=0.1*10;
    kR[114,0]=0.001;
    kR[115,0]=0.1*10;
    kR[116,0]=0.001;
    kR[117,0]=0.1*10;
    kR[118,0]=0.001;
    kR[119,0]=0.1*10;
    kR[120,0]=0.001;
    kR[121,0]=0.1*10;
    kR[122,0]=0.001;
    kR[123,0]=0.001;
    kR[124,0]=0.001;
    kR[125,0]=0.001;
    kR[126,0]=0.001;
    kR[127,0]=0.001;
    kR[128,0]=0.001;
    kR[129,0]=0.001;
    kR[130,0]=0.001;
    kR[131,0]=0.001;
    kR[132,0]=0.001;
    kR[133,0]=0.001;
    kR[134,0]=0.001;
    kR[135,0]=0.001;
    kR[136,0]=0.001;
    kR[137,0]=0.001;
    kR[138,0]=0.001;
    kR[139,0]=0.001;
    kR[140,0]=0.001;
    kR[141,0]=0.001;
    kR[142,0]=0.001;
    kR[143,0]=0.1*10;
    kR[144,0]=0.05;
    kR[145,0]=0.0016;
    kR[146,0]=0.001;
    kR[147,0]=0.1*10;
    kR[148,0]=0.001;
    kR[149,0]=0.00002;
    kR[150,0]=3.47E-4;
    kR[151,0]=3.12E-4/10;
    kR[152,0]=3.47E-4;
    kR[153,0]=0.0031;
    kR[154,0]=4.9E-4/5;
    kR[155,0]=0.0031;
    kR[156,0]=2.79E-5/5;
    kR[157,0]=0.01;




    # kRP **
    td=29;

    kRP1 = np.zeros(shape=(td,1))
    kRP1[0:td,0]=1;

    kRP2 = np.zeros(shape=(td,1))
    kRP2[0:td,0]=1;

    kRP3 = np.zeros(shape=(td,1))
    kRP3[0:td,0]=0.001;

    kRP4 = np.zeros(shape=(td,1))
    kRP4[0:td,0]=0.001;

    kRP5 = np.zeros(shape=(29,1))
    kRP5[2:4,0]=0.004667;
    kRP5[[0,1,10,15],0]=0.001667;
    kRP5[4:8,0]=0.001667;


    kRP5[[8,9],0]=1.6667E-4;
    kRP5[11:15,0]=1.6667E-4;
    kRP5[16:25,0]=1.6667E-4;


    kRP5[[25,27],0]=7.9309E-5;
    kRP5[26,0]=1.6667E-4;
    kRP5[28,0]=1.6667E-4;


    kRP6 = np.zeros(shape=(td,1))
    kRP6[0:td,0]=1;

    kRP7 = np.zeros(shape=(td,1))
    kRP7[0:td,0]=1;

    kRP8 = np.zeros(shape=(29,1))
    kRP8[2:4,0]=8.3711E-4;
    kRP8[[0,1,10,15],0]=8.5205E-4;
    kRP8[4:8,0]=8.5205E-4;
    kRP8[[8,9],0]=0.0014;
    kRP8[11:15,0]=0.0014;
    kRP8[16:25,0]=0.0014;

    kRP8[[25,27],0]=7.7E-7;
    kRP8[[26,28],0]=3.1955E-5;

    kRP9 = np.zeros(shape=(td,1))
    kRP9[0:td,0]=0.001;

    kRP10 = np.zeros(shape=(td,1))
    kRP10[0:td,0]=2.6;

    kRP11 = np.zeros(shape=(td,1))
    kRP11[0:td,0]=0.001;

    kRP12 = np.zeros(shape=(td,1))
    kRP12[0:td,0]=2.6;

    kRP13 = np.zeros(shape=(td,1))
    kRP13[0:td,0]=0.001;

    kRP14 = np.zeros(shape=(td,1))
    kRP14[0:td,0]=0.51;

    kRP15 = np.zeros(shape=(td,1))
    kRP15[0:td,0]=0.001;

    kRP16 = np.zeros(shape=(td,1))
    kRP16[0:td,0]=0.07;

    kRP17 = np.zeros(shape=(td,1))
    kRP17[0:td,0]=0.001;

    kRP18 = np.zeros(shape=(td,1))
    kRP18[0:td,0]=0.07;

    kRP19 = np.zeros(shape=(td,1))
    kRP19[0:td,0]=0.001;

    kRP20 = np.zeros(shape=(td,1))
    kRP20[0:td,0]=0.01;

    kRP21 = np.zeros(shape=(td,1))
    kRP21[0:td,0]=1;

    kRP22 = np.zeros(shape=(td,1))
    kRP22[0:td,0]=0.001*1000;

    kRP23 = np.zeros(shape=(td,1))
    kRP23[0:td,0]=0.01*1000;

    kRP24 = np.zeros(shape=(td,1))
    kRP24[0:td,0]=1*10;

    kRP25 = np.zeros(shape=(td,1))
    kRP25[0:td,0]=0.001;

    kRP26 = np.zeros(shape=(td,1))
    kRP26[0:td,0]=0.01*10;

    kRP27 = np.zeros(shape=(td,1))
    kRP27[0:td,0]=1;

    kRP28 = np.zeros(shape=(td,1))
    kRP28[0:td,0]=0.001;

    kRP29 = np.zeros(shape=(td,1))
    kRP29[0:td,0]=0.01;

    kRP30 = np.zeros(shape=(td,1))
    kRP30[0:td,0]=1;

    kRP28[25:29,0]=kRP28[25:29,0]*5;
    kRP29[25:29,0]=kRP29[25:29,0]*5;
    kRP30[25:29,0]=kRP30[25:29,0]*5;

    kRP31 = np.zeros(shape=(td,1))
    kRP31[0:td,0]=0.001;

    kRP32 = np.zeros(shape=(td,1))
    kRP32[0:td,0]=0.01;

    kRP33 = np.zeros(shape=(td,1))
    kRP33[0:td,0]=1;

    kRP34 = np.zeros(shape=(16,1))
    kRP34[0:16:2,0]=0.001;
    kRP34[1:16:2,0]=0.01;


    kP = np.zeros(shape=(190,1))


    kP[0,0]=0.001;
    kP[1,0]=0.01;
    kP[2,0]=0.001;
    kP[3,0]=0.01;
    kP[4,0]=0.001*100;
    kP[5,0]=0.1*10;
    kP[6,0]=0.001*100;
    kP[7,0]=0.1*10;
    kP[8,0]=0.001*100;
    kP[9,0]=0.01*10;
    kP[10,0]=0;
    kP[11,0]=0;
    kP[12,0]=0.001*100;
    kP[13,0]=0.1*10;
    kP[14,0]=0.1*10;
    kP[15,0]=0.001*100;
    kP[16,0]=0.1*10;
    kP[17,0]=0.1*10;
    kP[18,0]=0.001*100;
    kP[19,0]=0.1*10;
    kP[20,0]=1*10;
    kP[21,0]=0.001*100;
    kP[22,0]=0.1*10;
    kP[23,0]=1*10;
    kP[24,0]=0.001*100;
    kP[25,0]=0.1*10;
    kP[26,0]=0.5*10;
    kP[27,0]=0.001*100;
    kP[28,0]=0.1*10;
    kP[29,0]=0.5*10;
    kP[30,0]=0.001*100;
    kP[31,0]=0.001*10;
    kP[32,0]=0.001;
    kP[33,0]=0.1*10;
    kP[34,0]=1;
    kP[35,0]=0.001;
    kP[36,0]=0.1*10;
    kP[37,0]=1;
    kP[38,0]=0.001;
    kP[39,0]=0.1;
    kP[40,0]=1;
    kP[41,0]=0.001;
    kP[42,0]=0.001*100;
    kP[43,0]=0.001*10;
    kP[44,0]=0.013;
    kP[45,0]=0.001;
    kP[46,0]=0.1;
    kP[47,0]=1;
    kP[48,0]=0.01;
    kP[49,0]=0.001;
    kP[50,0]=0.1;
    kP[51,0]=1;
    kP[52,0]=0.01;
    kP[53,0]=0.001;
    kP[54,0]=0.1;
    kP[55,0]=1;
    kP[56,0]=0.001*20;
    kP[57,0]=0.1*20;
    kP[58,0]=1*10;
    kP[59,0]=0.001;
    kP[60,0]=0.1;
    kP[61,0]=1;
    kP[62,0]=1;
    kP[63,0]=0.001;
    kP[64,0]=0.1;
    kP[65,0]=1;
    kP[66,0]=0.001;
    kP[67,0]=0.1;
    kP[68,0]=1;
    kP[69,0]=0.001;
    kP[70,0]=0.1;
    kP[71,0]=1;
    kP[72,0]=0.1*10;
    kP[73,0]=0.011;
    kP[74,0]=0;
    kP[75,0]=0.001;
    kP[76,0]=0.1;
    kP[77,0]=1;
    kP[78,0]=0.018*50;
    kP[79,0]=0.011;
    kP[80,0]=0.001;
    kP[81,0]=0.1;
    kP[82,0]=1;
    kP[83,0]=0.001;
    kP[84,0]=0.1;
    kP[85,0]=1;
    kP[86,0]=0.001;
    kP[87,0]=0.1;
    kP[88,0]=1;
    kP[89,0]=0.001;
    kP[90,0]=0.1;
    kP[91,0]=1;
    kP[92,0]=0.001;
    kP[93,0]=0.1/10;
    kP[94,0]=0.001*1000;
    kP[95,0]=0.1*1000;
    kP[96,0]=0.001;
    kP[97,0]=0.1;
    kP[98,0]=0.001;
    kP[99,0]=0.1;
    kP[100,0]=1;
    kP[101,0]=0.001;
    kP[102,0]=0.001;
    kP[103,0]=0.1;
    kP[104,0]=0.001*10;
    kP[105,0]=0.1*10;
    kP[106,0]=10;
    kP[107,0]=0.001;
    kP[108,0]=0.1;
    kP[109,0]=0.001;
    kP[110,0]=0.1;
    kP[111,0]=0.001;
    kP[112,0]=0.001/10;
    kP[113,0]=0.001*50;
    kP[114,0]=0.1*50;
    kP[115,0]=1*10;
    kP[116,0]=0.001*100;
    kP[117,0]=0.001*50;
    kP[118,0]=0.1*50;
    kP[119,0]=1*10;
    kP[120,0]=0.001*10;
    kP[121,0]=0.01*100;
    kP[122,0]=0.01/100;
    kP[123,0]=0.001;
    kP[124,0]=0.1;
    kP[125,0]=1;
    kP[126,0]=0.01;
    kP[127,0]=60;
    kP[128,0]=0.1;
    kP[129,0]=1;
    kP[130,0]=0.01;
    kP[131,0]=0.012;
    kP[132,0]=0.018;
    kP[133,0]=0.001*10;
    kP[134,0]=0.1*10;
    kP[135,0]=1*1;
    kP[136,0]=0.001/10;
    kP[137,0]=0.1;
    kP[138,0]=1/10;
    kP[139,0]=0.001*100;
    kP[140,0]=0.01;
    kP[141,0]=0;
    kP[142,0]=0;
    kP[143,0]=0;
    kP[144,0]=0;
    kP[145,0]=0.001;
    kP[146,0]=0.1*10;
    kP[147,0]=0.001;
    kP[148,0]=0.001/10;
    kP[149,0]=0.001*100;
    kP[150,0]=0.1*10;
    kP[151,0]=1*10;
    kP[152,0]=0.001*10;
    kP[153,0]=0.001*100;
    kP[154,0]=0.1*10;
    kP[155,0]=1*10;
    kP[156,0]=0.001*30;
    kP[157,0]=0.001;
    kP[158,0]=0.1;
    kP[159,0]=1;
    kP[160,0]=0.001;
    kP[161,0]=0.001;
    kP[162,0]=0.1;
    kP[163,0]=1;
    kP[164,0]=0.001;
    kP[165,0]=0.001;
    kP[166,0]=0.0001;
    kP[167,0]=0.001;
    kP[168,0]=0.001;
    kP[169,0]=0.0001;
    kP[170,0]=0.001;
    kP[171,0]=0.01;
    kP[172,0]=0.0483*10;
    kP[173,0]=0.00001*100;
    kP[174,0]=0.001*100;
    kP[175,0]=0.001*10;
    kP[176,0]=0.001*50;
    kP[177,0]=0.001*100;
    kP[178,0]=0.001*1000;
    kP[179,0]=0.001/100;
    kP[180,0]=0.1;
    kP[181,0]=0.01;
    kP[182,0]=0.00001/2;
    kP[183,0]=0.01;
    kP[184,0]=0.0002*2;
    kP[185,0]=0.01;
    kP[186,0]=0.001;
    kP[187,0]=0.0001;
    kP[188,0]=0.001;
    kP[189,0]=0.0001;



    # kDA **
    # NONE as these are all transcriptional control mechanisms.

    kDP = np.zeros(shape=(4,1))

    # kDP **
    kDP[0,0]=0;
    kDP[1,0]=0;
    kDP[2,0]=0;
    kDP[3,0]=0;


    kPA = np.zeros(shape=(11,1))

    # kPA **
    kPA[0,0]=0.001;
    kPA[1,0]=0.1;
    kPA[2,0]=1;
    kPA[3,0]=0.001;
    kPA[4,0]=0.001;
    kPA[5,0]=0.1;
    kPA[6,0]=1;
    kPA[7,0]=0.001;
    kPA[8,0]=0.1;
    kPA[9,0]=1;
    kPA[10,0]=0.001;

    # kPC **
    # NONE as these are all transcriptional control mechanisms.

    # DC **
    # NONE as these are all transcriptional control mechanisms.





    return (kD,kC,kA,kR,kP, \
        kRP1,kRP2,kRP3,kRP4,kRP5,kRP6,kRP7,kRP8,kRP9,kRP10,kRP11,kRP12,kRP13,kRP14,kRP15,kRP16,kRP17,kRP18,kRP19,kRP20,kRP21,kRP22,kRP23,kRP24,kRP25,kRP26,kRP27,kRP28,kRP29,kRP30,kRP31,kRP32,kRP33,kRP34, \
        kDP,kPA)