import nsrunit; unit conversion on; unit s=1 second^1; unit nM=1E-6 meter^(-3)*mole^1; math main { //Warning: the following variables were set 'extern' or given // an initial value of '0' because the model would otherwise be // underdetermined: t, Ca extern real t s; extern real Ca nM; real M nM; M=6000; real Ntot nM; Ntot=1000; real K_mN nM; K_mN=535; // Var below replaced by constant in model eqns to satisfy unit correction // real n dimensionless; // n=2.92; real K_dN nM; K_dN=1760; real act_N dimensionless; // act_N=(Ca^2.92/(Ca^2.92+K_mN^2.92*(1+K_dN/M))); }