diff --git a/README.md b/README.md
index 1d4cd57b03feea5852e41e0c4056f8883bef1858..60ea4fcd7bc093a3330b93e4437e0e614abc9b55 100644
--- a/README.md
+++ b/README.md
@@ -1,4 +1,6 @@
-# Keshtgar et al. (2022). Cloud-radiative impact on the dynamics and predictability of an idealized extratropical cyclone
+# Scripts for "Cloud-radiative impact on the dynamics and predictability of an idealized extratropical cyclone"
+
+Keshtgar et al. (2023), https://doi.org/10.5194/wcd-2022-35
 
 Code repository for the ICON simulation run scripts, scripts for deriving baroclinic life cycle initial conditions, postprocessing of model output files, and the analysis scripts.
-The post-processed data used in the analysis along with a copy of the Git repository are published at the LMU open data server (https://doi.org/10.57970/h1y02-bjv70).
+The post-processed data used in the analysis along with a copy of the scripts here are also available at the LMU open data server (https://doi.org/10.57970/h1y02-bjv70).
diff --git a/blc_initial_conditions/README.md b/blc_initial_conditions/README.md
index efde33da9897c362b89959e7cc6a97d7de8952cf..93ffb179ffe49e4cbcae21172ecc4f2dd908bd1f 100644
--- a/blc_initial_conditions/README.md
+++ b/blc_initial_conditions/README.md
@@ -1,8 +1,8 @@
 **Initialization Procedure**
 
-The initial conditions for baroclinic life cycle simulations are the same as in Schäfer andVoigt, 2018 (https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017GL076726) and follow the life cycle type 1 configuration of Polvani and Esler, 2007 (https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2007JD008555) with the zonally uniform zonal wind in thermal wind balance with a zonally uniform temperature field.
+The initial conditions for baroclinic life cycle simulations are the same as in Schäfer an dVoigt, 2018 (https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017GL076726) and follow the life cycle type 1 configuration of Polvani and Esler, 2007 (https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2007JD008555) with the zonally uniform zonal wind in thermal wind balance with a zonally uniform temperature field.
 
-The maximum speed of the jet is set to 45 m/s and is centered at 45°N. Meridional and vertical velocities are also set to 0 m/s. The initial surface temperature is 300 . The baroclinic life cycle is triggered by adding a sinusoidal thermal wave with a 1 K amplitude at all levels and with a wavelength equal to the domain length. The surface relative humidity (RH0), is set to 80%, and the moisture scale height (Zt) is set to 12 km. The relative humidity is specified according to:
+The maximum speed of the jet is set to 45 m/s and is centered at 45°N. Meridional and vertical velocities are also set to 0 m/s. The initial surface temperature is 30 K. The baroclinic life cycle is triggered by adding a sinusoidal thermal wave with a 1 K amplitude at all levels and with a wavelength equal to the domain length. The surface relative humidity (RH0), is set to 80%, and the moisture scale height (Zt) is set to 12 km. The relative humidity is specified according to:
 
 RH = RH0*(1-0.85 x Z/Zt)**(5/4)
 
@@ -16,7 +16,7 @@ RH = RH0*(1-0.85 x Z/Zt)**(5/4)
 
 'Channel_4000x9000_2500m_with_boundary.nc' that is the channel grid generated by MPI grid_generator (grid.create.channels.run).
 
-These files will be uploaded to the KITOpen.
 
 2- The bash script '02_limited_channel_setup_2km_extpar_and_ozone.sh' prepares external files for the planar channel simulations: external fields and Ozone from a standard aquaplanet simulation.
 
+The input files are available at the LMU open data server.
diff --git a/icon_runscripts/README.md b/icon_runscripts/README.md
index eee23faeed4d937d13d3fbd446341bef5c43d0d8..d9a7abde971b5335efec229e57f2f2ddab179d40 100644
--- a/icon_runscripts/README.md
+++ b/icon_runscripts/README.md
@@ -1,25 +1,25 @@
 **ICON simulations**
 
-LC1-channel-4000x9000km-2km-0002 : no radiation
+LC1-channel-4000x9000km-2km-0002 : Simulation with no radiation
 
-LC1-channel-4000x9000km-2km-0003 : cloud radiaton
+LC1-channel-4000x9000km-2km-0003 : Simulation with only cloud radiation
 
-LC1-channel-4000x9000km-2km-0004 : 2x cloud radiation
+LC1-channel-4000x9000km-2km-0004 : Simulation with only cloud radiation scaled by a factor of 2
 
-LC1-channel-4000x9000km-2km-0005 : restart at day 3
+LC1-channel-4000x9000km-2km-0005 : Simulation with only cloud radiation until day 3
 
-LC1-channel-4000x9000km-2km-0006 : restart at day 4
+LC1-channel-4000x9000km-2km-0006 : Simulation with only cloud radiation until day 4
 
-LC1-channel-4000x9000km-2km-0007 : restart at day 5
+LC1-channel-4000x9000km-2km-0007 : Simulation with only cloud radiation until day 5
 
-LC1-channel-4000x9000km-2km-0008 : restart at day 6
+LC1-channel-4000x9000km-2km-0008 : Simulation with only cloud radiation until day 6
 
-LC1-channel-4000x9000km-2km-0009 : no radiation, two-moment microphysics
+LC1-channel-4000x9000km-2km-0009 : Simulation with no radiation, and with 2-moment microphysics
 
-LC1-channel-4000x9000km-2km-0010 : cloud radiaton, two-moment microphysics
+LC1-channel-4000x9000km-2km-0010 : Simulation with only cloud radiation, and with 2-moment microphysics
 
-LC1-channel-4000x9000km-80km-0001 : full radiation
+LC1-channel-4000x9000km-80km-0001 : Simulation with all-sky radiation
 
-LC1-channel-4000x9000km-80km-0002 : no radiation
+LC1-channel-4000x9000km-80km-0002 : Simulation with no radiation
 
-LC1-channel-4000x9000km-80km-0003 : cloud radiaton
+LC1-channel-4000x9000km-80km-0003 : Simulation with only cloud radiation