Rough charts based on the data above. Some different curves.
Note that the heat120 method appears to shift the air-snow and snow-ice thermistor no.s upwards, possibly due to heat rising. I haven't corrected for this. Even so, that doesn't explain 0 or negative starting position unless data was sent before (or during) buoy installation. hey ho, it's an estimate.
Im actually working on a "basic method" script at the moment.
Within that script is the algo for detecting Th(snow) and Th(ice) thermistors .They look similar to your estimates, though not identical, ill send you it within 30mins.
New points in the code;
-basic method implemented
-lat and lon now taken from TS buoys, corrected for sampling period (to 1 per day)
Im also working on an additional program that lets you
manually enter Th(snow) and Th(ice) thermistors, when you wish to plot. In some ways I think the human eye might actually be better at this than an algorithm
//////
FYI all the code is pretty horrible at the moment, im reading a book called
"Clean Code, A Handbook of Agile Software Craftmanship"
Once were happy, ill be tightening everything up very strictly