* profiles_gen reorganized to be a top level code, and compiled after TGYRO
* CGYRO continues to be developed
* revamp of gyro and cgyro python plotting
* clean up some platforms for GYRO
* remove VENUS platform
* fix size of GYRO exchange term
* make use of newer fortran standard that defines bessel functions
* made fftw3 default for GYRO and CGYRO
* revised reading/writing of restart files for MPI-IO capable platforms
* new parallel NEO capability
* toroidal NEO current is exactly converted from the parallel current
* new platforms: CORE_IFORT, CORI_INTEL, CORI_KNL_HT2_IFORT, CORI_KNL_IFORT, EDISON_IFORT, GATHETA
* profiles_gen uses ne psi_grid unless ne psi_grid has psi_n>1, in which case use a uniform grid
* remove definitions of GAMMA and Bessel functions
* write the TGLF field spectrum
* consolidate TGLF quasilinear weight variables to a single module
* new logic for handling `IBRANCH`
* new TGYRO flag: `TGYRO_ZERO_DENS_GRAD_FLAG` an index to set that species' density gradient to 0 to easily get the D and V of NEO/TGLF.
* fix TGYRO output when using 5 (or more) species