Data-treatment.md

@@ -15,7 +15,13 @@ if file_name == "mrunoff.nc":
 # Harvest
 - Convert harvest to yield :  
 ```
-yield = np.multiply(np.multiply(harvest,0.02222222),1.14)
+#gC/m2 -> tDM/ha (dry matter)
+            harvests_tdm = np.multiply(harvests,0.02222222)
+            #tDM/ha -> tFM/ha (fresh matter)
+            harvests_tfm = np.multiply(harvests_tdm,1.14) #tonne / ha
+        
+            #production = harvest * surface  [years 101; latitude 30; longitude 33]
+            yields = np.multiply(harvests_tfm,cft_surface) #tonne
 ```
 # Convert latitude and longitude into lpjml 30 minutes system
 ```plaintext