 
     // Now run the MCMC sampler
 
-    std::ofstream chifile("chis.dat");
+    //std::ofstream chifile("chis.dat");
     std::ofstream thetafile("thetas.dat");
     
-    int mcmc_iter = 1000;
+    int mcmc_iter = 10000;
     int naccept_theta = 0;
+    std::cout << "Running MCMC Sampler...." << std::endl;
+    
     for (int i=0; i<mcmc_iter; i++) {
         // Now grab the pointers to the vectors, needed to run the kernel since it doesn't understand Thrust
         // We do this here because the thrust vector are smart, and we want to make sure they don't reassign

static __constant__ int c_m = 5;
    // observational frequencies, corresponding to the Herschel PACS and SPIRES bands
static const double nu[m] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
static __constant__ c_nu[m] = {4.286e12, 1.765e12, 1.200e12, 8.571e11, 6.000e11};
static const double hplanck = 6.6260755e-27;  // Planck's constant, in cgs
static const double clight = 2.997925e10;
static const double kboltz = 1.380658e-16;

     collection of the mean and variances of the normal distributions.
    */
    
    int n = 10000; // # of data points
    
    double mu_norm = 8.5 * log(10);  // Average value of the natural logarithm of the SED normalization
    double sig_norm = 0.5 * log(10);  // Standard deviation in the SED normalization
