Update Data treatment authored by Marie Brunel's avatar Marie Brunel
...@@ -16,12 +16,11 @@ if file_name == "mrunoff.nc": ...@@ -16,12 +16,11 @@ if file_name == "mrunoff.nc":
- Convert harvest to yield : - Convert harvest to yield :
``` ```
#gC/m2 -> tDM/ha (dry matter) #gC/m2 -> tDM/ha (dry matter)
harvests_tdm = np.multiply(harvests,0.02222222) harvests_tdm = np.multiply(harvests,0.02222222)
#tDM/ha -> tFM/ha (fresh matter) #tDM/ha -> tFM/ha (fresh matter)
harvests_tfm = np.multiply(harvests_tdm,1.14) #tonne / ha harvests_tfm = np.multiply(harvests_tdm,1.14) #tonne / ha
#production = harvest * surface [years 101; latitude 30; longitude 33]
#production = harvest * surface [years 101; latitude 30; longitude 33] yields = np.multiply(harvests_tfm,cft_surface) #tonne
yields = np.multiply(harvests_tfm,cft_surface) #tonne
``` ```
# Convert latitude and longitude into lpjml 30 minutes system # Convert latitude and longitude into lpjml 30 minutes system
```plaintext ```plaintext
... ...
......