   	{
   		y_dummy[i] /= sum;
   	}
-  	Real T = 1200;
-  	y_dummy[0] = T;
+  	y_dummy[0] = LOW_T;
   	Real P = 101325;
   	for (int j = 0; j < NUM; j++)
   	{
   		pres_host[j] = P;
-  		T_host[j] = T;
+  		if (j % 2 == 0)
+  			T_host[j] = LOW_T;
+  		else
+  			T_host[j] = HI_T;
   	}
 
   	populate(NUM, P, y_dummy, conc_host);


// load same initial conditions for all threads
#define SAME_IC
#define REPEATS 1

// shuffle initial conditions randomly
//#define SHUFFLE

  	//come up with a random state
  	for (int i = 1; i < NN; i++)
  	{
  		Real r = 1;
  		sum += r * r;
  		y_dummy[i] = r;
  	}
