All simulations are to be compared with present-day simulations forced by AMIP SSTs and sea-ice

1) limit conditions for LGM from IPSL model:

* SST and sea-ice:
LGM_IPSL_minPI_plusAMIP_histmth_SST_SIC.nc

* land ice fraction: use fract_lic (between 0 and 1) from:
composite_R96x71_controle_AMIP.204-204.nc
for higher resolution, you can also use:
landiceref.nc

* topography: use phis (in m) from:
composite_R96x71_controle_AMIP.204-204.nc
This contains also the land fraction: fract_ter
for higher resolution, you can also use: 
Relief_ICE5G_3_camille.nc

To make sure that everythink is consistent with your control simulation, you can actually add the LGM-PD difference of SST, sea-ice, topography, land and land-ice fraction:
controle_AMIP_histmth_SST_SIC.nc
landiceref.nc
Relief.nc

* gas concentration and orbital parameters for LGM:
excentricity = 0.018994
precession = 114.42
obliquity = 22.949
solar_constant = 1365.
co2_ppm = 185.
CH4_ppb = 350.
N2O_ppb = 200.
CFC11_ppt = 0.
CFC12_ppt = 0.

* as a reminder, gas concentration and orbital parameters that I used for present-day:
excentricity = 0.016715
precession = 102.7
obliquity = 23.441
solar_constant = 1365.
co2_ppm = 348.
CH4_ppb = 1650.
N2O_ppb = 306.
CFC11_ppt = 280.
CFC12_ppt = 484.

2) limit conditions for LGM from CLIMAP reconstruction:

Same as 1) but use:
LGM_climap_histmth_SST_SIC.nc
for SST and sea-ice.

3) limit conditions for MH from IPSL model:

* just change SSTs and sea-ice:
MH_IPSL_minPI_plusAMIP_histmth_SST_SIC.nc

* gas concentration and orbital parameters 
excentricity = 0.018682
precession = 0.87
obliquity = 24.105
solar_constant = 1365.
co2_ppm = 280.
CH4_ppb = 650.
N2O_ppb = 200.
CFC11_ppt = 0.
CFC12_ppt = 0.

4) limit conditions for 4xCO2 from ISPL model:

* just change SSTs and sea-ice:
4xCO2_IPSL_minPI_plusAMIP_histmth_SST_SIC.nc

* gas concentration and orbital parameters 
excentricity = 0.016715
precession = 102.7
obliquity = 23.441
solar_constant = 1365.
co2_ppm = 1400.
CH4_ppb = 1650.
N2O_ppb = 306.
CFC11_ppt = 280.
CFC12_ppt = 484.