inconsistencies in find_10thetas

[imab4bsh]

In summer of 2020 (23-27 Jul) we (Birgit, Cécile and Kamila) had an e-mail exchange about some inconsistencies found in the find_10thetas.m script. I am opening this issue here so the discussion does not go lost and some day we fix this :)

There are two main problems:

1. Sometimes the selected theta levels are not present in the entire time series and therefore the diagnostic plots have empty spaces. Even if the operator selects time series breaks, the function does not takes them into account.
2. In frontal regions the selected theta levels are present in very different depths which also leads to misleading plots.
   An example (shared by Birgit) can be found here: 7900496_test_OWCall.docx

Moreover, for us it is unclear how this issues affect the further calculations of salinity differences.

@cabanesc @kamwal @BSHbirgit @apswong @gmaze
ps. should we have this discussion in our public QC forum or leave it here?

[gmaze]

ps. should we have this discussion in our public QC forum or leave it here?

I think this is the perfect place to discuss this issue !