Merge pull request #23 from lopezvoliver/tests

Fix bug in Td
kaust-halo · Mar 6, 2024 · ed8d683 · ed8d683
2 parents 8bb7881 + f9b1fa6
commit ed8d683
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/geeet/tseb.py b/geeet/tseb.py
@@ -381,7 +381,7 @@ def tseb_series_iteration(img, list):
             Td1 = Tclin.multiply(Image(1).add(rso.divide(rao)))
             Td2 = (Hco.multiply(rxo).divide(rho.multiply(cp))).multiply(Image(1).add(rso.divide(rxo)).add(rso.divide(rao)))
             Td3 = Ta.multiply(rso.divide(rao))
-            Td = Td1.add(Td2).subtract(Td3)
+            Td = Td1.subtract(Td2).subtract(Td3)
             # N95 equation A.11 (deltaTc), i.e. correction to the linear approximation of Tc:
             dTc_num = (Tr.pow(4)).subtract(f_theta.multiply(Tclin.pow(4))).subtract((Image(1).subtract(f_theta)).multiply(Td.pow(4)))
             dTc_denom1 = (Td.pow(3)).multiply(rso.divide(rao).add(1)).multiply(4).multiply(Image(1).subtract(f_theta))