Improved input file handling

Added local_heating_option ‘Wang’ to copy Wang,Gerber, and Polvani JC
2012 localized heating. This does not necessarily work yet!
Inputs for Te, tau, and heating are now not used as zonal means, but w/
longitudinal dependence: This is important in the presence of
topography.
Compute HS94/PK02

MATLAB routines to compute Held and Suares 1994 and Polvani and Kushner
2002 Te and tau profiles.
These can then be given to the code as input files temp.nc, and tau.nc
simpler default grid

made default grid simpler for better control over the used grid