Hi Signe,
You can force the surface layer of md.thermal.spctemperature. Use timestamps in the last row,
the same way you would force smb. Each column will represent temperature at a specific time.
All non-constrained values should be set always to NaN. Where you wish to impose a single
point constraint (in your case, the surface layer), you will set a temperature value.
Good luck!
Nicole