added header to scripts

plus some clarification in the code

README.md

@@ -19,14 +19,21 @@ Humboldt-Universität zu Berlin.
 
 ### Ressources
 This study builds on several resources:
-- [Farinotti et al. (2019): A consensus estimate for the ice thickness distribution of all glaciers on Earth, In Nat. Geosci. 12 (3), pp. 168–173. ](https://doi.org/10.1038/s41561-019-0300-3)
-- [RGI Consortium (2017): Randolph Glacier Inventory - A Dataset of Global Glacier Outlines: Version 6.0, Global Land Ice Measurements from Space Initiative.](https://doi.org/10.7265/N5-RGI-60)
-- [Takaku et al. (2018): Quality Improvements of ‘AW3D’ Global Dsm Derived from Alos Prism. In : Proc. IGARSS 2018](10.1109/IGARSS.2018.8518360)
+- [Farinotti et al. (2019): A consensus estimate for the 
+  ice thickness distribution of all glaciers on Earth, 
+  In Nat. Geosci. 12 (3), pp. 168–173. 
+  ](https://doi.org/10.1038/s41561-019-0300-3)
+- [RGI Consortium (2017): Randolph Glacier Inventory - 
+  A Dataset of Global Glacier Outlines: Version 6.0, 
+  Global Land Ice Measurements from Space Initiative.
+  ](https://doi.org/10.7265/N5-RGI-60)
+- [Takaku et al. (2018): Quality Improvements of ‘AW3D’ 
+  Global Dsm Derived from Alos Prism. In : Proc. IGARSS 
+  2018](10.1109/IGARSS.2018.8518360)
 
 ### Improving the dataset
-
 Please use the [Issues](https://github.com/cryotools/subglacial-overdeepenings/issues)
-section to report any errors. In order to share your own improvements to the 
+section to report any errors. If you would like to share your own improvements to the 
 code or the shapefiles, please fork the repository, commit your changes, 
 and file a merge request.
 

code/README.md

@@ -1,14 +1,20 @@
 ![GLAMoR](https://cryo-tools.org/wp-content/uploads/2020/07/GLAMoR-LOGO-400px.png)
 ## Scripts
-The scripts in this repository can be used to calculate subglacial overdeepenings.
-The main work is done by *calculations.py*, which consists of two subsequent routines
- for all glaciers. As the `region group` routine in ArcGIS is still buggy for
- large datasets, we outsourced this step to R. Therefore, the routine needs to be stopped
- and the results have to be run through the R script in order to provide the
- necessary data for the second python routine.
+The scripts in this repository can be used to calculate subglacial overdeepenings
+and estimate the hazard of mass movement impacts into potential future lakes.
+The main work is done by the python script *calculations.py*, 
+which consists of two subsequent routines for all glaciers. 
+
+As a key function in ArcGIS (`region group`) sometimes breaks
+when handling large datasets, we recommend to outsource this step 
+to the `raster` package of R. With this, R steps in to provide the 
+necessary overdeepening identification for the second routine of 
+the python script. 
+For more details, please see the header of *calculations.py*. 
 
 ### Required data
-In order to run properly, the scripts in this repository require different datasets to be downloaded.
+In order to run properly, the scripts in this repository require 
+different datasets to be downloaded.
 For the glacier ice thickness estimates, we used the dataset of 
 [Farinotti et al. (2019)](https://doi.org/10.1038/s41561-019-0300-3)
 The glacier outlines can be downloaded from the 
@@ -17,9 +23,11 @@ While theoretically a number of DEMs are available for HMA, we recommend using t
 [ALOS World 3D - 30m (AW3D30)](https://www.eorc.jaxa.jp/ALOS/en/aw3d30/index.htm) due to its superior accuracy.
 
 ### Required software
-Valid licences for several extensions of the ArcGIS software are needed, 
-i.e. the *ArcGIS 3D Analyst*, the *ArcGIS Spatial Analyst*, and the *ArcGIS Geostatistical Analyst*.
-While most of the code is written in Python, we use R for some steps as well.
+The python script was developed using the `arcpy` package af ArcGIS 10.7. 
+Therefore, valid licences for several extensions of the ArcGIS software are needed, 
+i.e. the *3D Analyst*, the *Spatial Analyst*, 
+and the *Geostatistical Analyst*.
+While most of the code is written in Python 2.7, we use R 3.6.1 for some steps as well.
 
 ### Preprocessing
 Some preprocessing steps are not included in the scripts as they heavily depend on the employed