Data-treatment.md

@@ -16,4 +16,19 @@ if file_name == "mrunoff.nc":
 - Convert harvest to yield :  
 ```
 yield = np.multiply(np.multiply(harvest,0.02222222),1.14)
+```
+# Convert latitude and longitude into lpjml 30 minutes system
+```plaintext
+pts = [ 
+(-43.42368753718026, -12.682038634410734), (-42.71628164452376, -10.559820956441273), (-41.70570179787164, -9.549241109789147), (-42.81733962918898, -15.71377817436711)]
+for pt in pts :
+    new_pt = []
+    for i in [1,0]:
+        p = pt[i]
+        if round(p) > p :
+            p = round(p)-0.25
+        else : 
+            p = round(p)+0.25
+        new_pt.append(p)
+    print(new_pt)  
 ```
\ No newline at end of file