keep main auxiliary routines here now

ND_progs - nondimensionalisation routines
GEN_progs - general routines
- including special functions, orthogonal functions, bessel functions