GitLab

renamed L options "S80" -> "tsrv" and "ecmwf" -> "Rb"
changed default L option relative to method used in get_init

changed gamma_moist in hum_subs to gamma for different options of adiabatic lapse rate (dry-constant, definition, or for unsaturated air with water vapor-, moist)
adjusted it in ASFC 

- set default method S88

- for first three iterations u10n<0 is set to 0.5 to allow convergence for small negative values
- added initial restriction for n not to be less than 5; if so it is set to 5

changed flag "l" in AirSeaFluxCode to be set when Rib<-0.5 or Rib>0.2

Deleted Documentation.pdf, data_all_out.csv, data_all_stats.txt, readme.txt files

updated toy_ASFC to run with new version.
updated data_all_out.csv to include psiq in header.
updated data_all_stats.txt and
added readme.txt that gives the input choices in toy_ASFC

AirSeaFluxCode.py updated with more options for the cool skin correction functions and a warm layer correction function

Update AirSeaFluxCode.py, flux_subs.py, get_init.py, hum_subs.py, data_all.nc, Documentation.pdf, era5_r360x180.nc, run_ASFC.py files
Deleted data_Mxtr.nc, data_mod.nc, data_trp.nc, data_xtr.nc files

- get_init (for initial values check and setting defaults)
- hum_subs that includes all subroutines related to humidity calculations
- flux_subs that contains all subroutines related to flux calculations (exchange coefficients, stability functions etc.)
- util_subs that contains utility subroutines as get_height, gc, visc_air etc.
In subroutine get_L the set up of UA option is set up to include Richardson number instead of deriving it from virtual temperature

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