I took a look at the code, and if we consider the orderededDict "output_structure", the algorithm will export ancillary channels (S_A) in this order:
1) ('R_ns1', S_A), # net shortwave radiation reaching the surface at time t1
2) ('R_nl1', S_A), # net longwave radiation reaching the surface at time t1
3) ('delta_R_n1', S_A), # net radiation divergence in the canopy at time t1
4) ('H_C1', S_A), # canopy sensible heat flux (W/m^2) at time t1
5) ('LE_C1', S_A), # canopy latent heat flux (W/m^2) at time t1
6) ('LE_partition', S_A), # Latent Heat Flux Partition (LEc/LE) at time t1
7) ('T_C1', S_A), # canopy temperature at time t1 (deg C)
8) ('T_S1', S_A), # soil temperature at time t1 (deg C)
9) ('R_A1', S_A), # resistance to heat transport in the surface layer (s/m) at time t1
10) ('R_x1', S_A), # resistance to heat transport in the canopy surface layer (s/m) at time t1
11) ('R_S1', S_A), # resistance to heat transport from the soil surface (s/m) at time t1 fluxes
12) ('L', S_A), # Monin Obukhov Length at time t1
13) ('u_friction', S_A), # Friction velocity
14) ('flag', S_A), # Quality flag
Observation: Fields 7 and 8 showing deg C at code documentation
According to code https://github.com/hectornieto/pyTSEB/blob/master/pyTSEB/PyTSEB.py#L950, there is an error in the order of channels.
I took a look at the code, and if we consider the orderededDict "output_structure", the algorithm will export ancillary channels (S_A) in this order:
1) ('R_ns1', S_A), # net shortwave radiation reaching the surface at time t1 2) ('R_nl1', S_A), # net longwave radiation reaching the surface at time t1 3) ('delta_R_n1', S_A), # net radiation divergence in the canopy at time t1 4) ('H_C1', S_A), # canopy sensible heat flux (W/m^2) at time t1 5) ('LE_C1', S_A), # canopy latent heat flux (W/m^2) at time t1 6) ('LE_partition', S_A), # Latent Heat Flux Partition (LEc/LE) at time t1 7) ('T_C1', S_A), # canopy temperature at time t1 (deg C) 8) ('T_S1', S_A), # soil temperature at time t1 (deg C) 9) ('R_A1', S_A), # resistance to heat transport in the surface layer (s/m) at time t1 10) ('R_x1', S_A), # resistance to heat transport in the canopy surface layer (s/m) at time t1 11) ('R_S1', S_A), # resistance to heat transport from the soil surface (s/m) at time t1 fluxes 12) ('L', S_A), # Monin Obukhov Length at time t1 13) ('u_friction', S_A), # Friction velocity 14) ('flag', S_A), # Quality flag
Observation: Fields 7 and 8 showing deg C at code documentation