Data-treatment.md

@@ -36,11 +36,18 @@ for pt in pts :
     print(new_pt)  
 ```
 # Cell fraction to area
+## 30 minutes resolution
 ```
 gridhectares = np.zeros((len(latitudes),len(longitudes)))
 for lat,latitude in enumerate(latitudes) :
      gridhectares[lat,:] = 110 * 110 * math.cos(math.pi*latitude/180) * 1/2 * 1/2 * 100
 ```
+## 5 minutes resolution
+```
+surfaces = np.zeros_like(cftfracs)
+for lat,latitude in enumerate(latitudes) :
+    surfaces[:,:,lat,:] = 110 * 110 * math.cos(math.pi*latitude/180) * 0.083 * 0.083 * 100
+```
 # Moving average
 ```
 ma_harvest = []