src/mock_sample_functions.c
double alpha_distribution(long *seed)  
{
	double alpha_avg,alpha_sigma,val;
	
	alpha_avg = -0.87;
        alpha_sigma = 0.33;
        val = alpha_avg + alpha_sigma*gasdev(seed);

	
	//flat distribution
	//val = (2.0-(-2.0))*ran2d(seed)+(-2.0);

	return val;
}
//find alpha; Gaussian distribution based on Sakamoto et al. 2009
find alpha; Gaussian distribution based on Sakamoto et al. 2009
find alpha ; Gaussian distribution based on Sakamoto et al . 2009
double alpha_distribution(long *seed)  
{
	double alpha_avg,alpha_sigma,val;
	alpha_avg = -0.87;
        alpha_sigma = 0.33;
        val = alpha_avg + alpha_sigma*gasdev(seed);
	return val;
}
double alpha_distribution ( long * seed ) { double alpha_avg , alpha_sigma , val ; alpha_avg = -0.87 ; alpha_sigma = 0.33 ; val = alpha_avg + alpha_sigma * gasdev ( seed ) ; return val ; }
find alpha; Gaussian distribution based on Sakamoto et al.
find alpha ; Gaussian distribution based on Sakamoto et al .
//flat distribution //val = (2.0-(-2.0))*ran2d(seed)+(-2.0);