Very interesting, had no idea about "flushing"
Could be that many of our assumptions about ice melting and refreeze need confirming.
For the time being, instead of having a set value of -1.8 I just set ice bottom as the thermistor below which the temp is 0.2degC less than whatever the ocean temp is at that measurement time.
''
ice_thermistors=ice_or_water_thermistors.where(ice_or_water_thermistors<mean_ocean_heat-0.2)
''
Looks good if my ice temperature doubts turn out to be wrong.
17 older tbuoys deployment report update
I have been spending most of my time working on the HEAT data as I agree there are probably issues in the TEMP data. The former is in some ways more powerful for estimating thickness,
and i suspect when they eventually publish their own method it will centre on HEAT not TEMP, or maybe a combination.
In my opinion its more powerful because you will never see a flat profile through the chain in HEAT, but you do sometimes see it in TEMP.
Its good they are providing their top ice thermistor now, that make things easier.
If i get anything solid ill send you a message with the code and you can take a look.