Code and data for "Sensitivity of Northern Hemisphere climate to ice-ocean interface heat flux parameterizations"

This repository provides the source code and modeled data for three different ice-ocean heat flux
parameterizations of a 1-D idealized model, as well as 3-D climate models including CICE,
MPIOM and COSMOS, which are used in a GMD manuscript called "Sensitivity of Northern Hemisphere climate to ice-ocean interface heat flux parameterizations".  The NCL-based scripts for plotting the figures are also provided.

The file all.tar.gz consists of 4 folders as the following:

1. 1-D
In the 1-D folder one can find the matlab source code for the 1-D idealized model with main.m being the main script and the others sub-scripts for calculating seasonal changes of different forcings, involving the surface albedo (albedo.m), shortwave fluxes (shortwave.m) and all other kinds of fluxes (otherfluxes.m).

2. code
The code folder provides the source code of the three models used in our study: CICE, MPIOM and COSMOS.

The most important code in terms of the ice-ocean heat flux in CICE can be found at code/cice/source/ice_therm_vertical.F90. The switch of the options of the three parameterizations can be achieved by changing the parameter "oceanic_heat" (1 for icebath, 2 for 2eq and 3 for 3eq) in the namelist when running the model.

The mpiom folder contains the three different set of MPIOM source code for the three ice-ocean heat flux parameterizations respectively.

In cosmos, one could find 4 sub-folders, with the folder echam5 containing the source code for the atmosphere module ECHAM5, and the other 3 folders containing the MPIOM source code incorporating with the three different ice-ocean heat flux parameterizations, similar as the mpiom folder.

3. data
This folder provides the simulated output from the three models: CICE, MPIOM and COSMOS. Each model folder contains three sub-folders called 2eq, 3eq and icebath, representing the modelled data for the 2eq, 3eq and icebath parameterizations respectively. For CICE, we upload the modeled results of the last 10 years. For MPIOM and COSMOS, as the original data set are too large, here we upload the climatology of the data from the last 100 simulation years. Note that in cosmos, there are some additional variables which are listed seperately, namely the AMOC (amoc.nc), the sea surface pressure (slp.nc) as well as the surface temperature (tsurf.nc). reg.nc contains ocean temperature and salinity which have been interpolated onto a 1x1 regular grid.

4. plot_figures
This folder gives the scripts based on NCL to plot the figures in the manuscript, with the *.ncl files being the plotting scripts and the *.eps being the figures.  
All the code, data, scripts can be used by anyone who has interest.