diff --git a/README.md b/README.md
index 9b8169d657a7045869803005665a9ece2e0df2ef..2dfaa028cd185a617d6bdc29ea71025307e2aa89 100644
--- a/README.md
+++ b/README.md
@@ -6,7 +6,7 @@ were generated.
 
 The scripts can be found inside `bs_headons`.
 In `bs_headons/postprocess` there are scripts to extrapolate bamps' `Psi4` data to
-null infinity and to compute the radiated GW energy.
+null infinity and to compute the radiated GW energy as discussed in the paper.
 
 ## Installation
 
@@ -16,7 +16,7 @@ null infinity and to compute the radiated GW energy.
 1. Clone this repository.
 2. Clone `GWanalysis.jl` into a separate directory from `https://git.tpi.uni-jena.de/fatteneder/GWanalysis.jl`
    Check out the version `ef2127891642a7d0c348417299bf8d40e085243` which is needed for this repository to work.
-3. Add `GWanalysis.jl` as dependency:
+3. Add `GWanalysis.jl` as a dependency:
    ```
    $ cd <this-dir>
    $ julia --project
@@ -31,12 +31,15 @@ null infinity and to compute the radiated GW energy.
   julia> Pkg.resolve()
   julia> Pkg.instantiate()
   ```
-5. Run a script:
+5. Adjust the variable `ROOTDIR` in `bs_headons/rootdirs.jl` so that it points to the
+   folder that contains all the bamps results with the subfolders `comparison4, comparison8, etc.`.
+   To get a grip of this data please ask Bernd Bruegmann, as I will let him know where I archive
+   the data.
+6. Run a plotting script:
   ```
   julia> include("bs_headons/dEdt_E/plot.jl")
   ```
 
-
 ## Tips & tricks
 
 - Try keeping the interactive REPL alive and rerun the scripts with `include`.
@@ -51,3 +54,24 @@ null infinity and to compute the radiated GW energy.
   $ julia> # edit something in GWanalysis.jl
   $ julia> include("bs_headons/dEdt_E/plot.jl") # changes should be picked up automatically
   ```
+
+---
+
+## Brief description of the plotting scripts
+
+Figure numbers refer to the arxiv version of the paper: https://arxiv.org/pdf/2311.16251
+
+- `amax_nocts-vs-cts`: Generates Fig 6
+- `cmon_bdry_effect`: Generates Fig 7, the heatmap is generated using Paraview
+- `cmon_nocts-vs-cts`: Generates Fig 2 & 3
+- `cmon_rsc_amax_nocts-vs-cts`: Generates Fig 5
+- `dEdt_comparison`: Generates Fig 14
+- `dEdt_E`: Generates Fig 13
+- `longrun`: Generates Fig 10
+- `N_M_Madm_comparison`: Generates Fig 8 & 9
+- `postprocess`: Scripts to process the Psi4 data from bamps output, their output is stored
+  in folders called `analyze` in the respective bamps output directories.
+- `psi4_comparison_cts_noncts`: Generates Fig 11
+- `psi4_dhdt_h_problem`: Generates Fig 16 for each of the bamps outputs we have.
+- `psi4_self_convergence`: Generates Fig 12
+- `psi4_wiener`: Generates Fig 15