GitLab

- removed ind as a second output and set itera to -1 when point has not converged 

2. added Monin-Obukhov length method switch (L)
2. fixed bug with L for method ERA5

2. added break if method not in options
3. added break if qmeth not in options
4. added break if gust not a 3x1 array
5. added break if tol option is not in options
6. added restriction if abs(L) < 1

AirSeaFluxCode.py with various switch changes for gustiness, tolerance limits, output to missing value for no convergence, humidity options.

flux_subs.py minor changes to improve computations

All u10n functions corrected  

adjusted get_height subroutine to read in multiple (3xn) height input. Added psiq as output. Included all switches jcool, gust, qmeth.