# 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.
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).
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 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)
RH = RH0*(1-0.85 x Z/Zt)**(5/4)
...
@@ -16,7 +16,7 @@ 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).
'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.
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.