Create photon_sct_gauss.c

Exceedingly fast, gcc version of Cashwell-Carter / Prahl combined with Kirk ATTEN5 program in c
	 
	 Has been used to follow 10^12 individual photons in the turbid ocean/atomosphere 
	 and stores their trajectory lengths and exit angles fficiently.   
	 Based on the Cashwell-Carter LLNL neutron scattering code highly modified and adapted
	 
	 ***
	 This is intended to be used with a water cooled, i7 to avoid CPU over heating. 
	   Use at your own risk otherwise! you may very well fry you CPU otherwise
***
	 ****
  Uses a Gaussian laser beam initial size wo and includes diffraction
	 SSE2 Meresenne twister random number generator
	 and the very fast, -pi to pi sine and cosine polonomial algorithm
	 by Charles K Garrett, Feb 2012

   standard arrays with option base 0
   calloc is needed for arrays > 100000
   Petzold data must not be CSV

    compile with:
    gcc -Ofast -msse2 -frename-registers -fno-strict-aliasing -DHAVE_SSE2=1 -DDSFMT_MEXP=19937 -o photon_scat_gauss dSFMT.c photon_scat_gauss.c -lrt -lm -lgsl -lgslcblas
    
 Uses: dsfmt.c and dsfmt.h frp, Saito and Matsumoto  
 contained in dsfmt-ser-2.2.3
    Copyright (C) 2007, 2008 Mutsuo Saito, Makoto Matsumoto and
    Hiroshima University. All rights reserved.
    Copyright (C) 2012 Mutsuo Saito, Makoto Matsumoto,
    Hiroshima University and The University of Tokyo.
    All rights reserved. see  BSD License in dsfmt.h