README file for the M-AMOC1.0 AMOC-emulator described by Bakker and Schmittner (2016).
Last edited by Pepijn Bakker at May 12th 2016.
 
This file will give a short overview of the included M-AMOC1.0 Matlab scripts and their functioning. The main function to perform the AMOC-emulator fitting procedure is: M-AMOC1.0


Included files:
create_constantsmat.m
density.m
M-AMOC1.0.m
runmodel_amoc.m
runmodel_coupledsystem.m
runmodel_ebm.m
Total_atmos_effect.m:


create_constantsmat.m:
Function that lists and produces all necesarry constants. In this file all values are set.
Called by: M-AMOC1.0.m

density.m:
Function that used temperature and salininty to calculate and output density using a non-linear equation of state.
Called by: runmodel_amoc.m

M-AMOC1.0.m:
This is the main function that controls the fitting procedure of the AMOC-emulator following Simullated Annealing. It calculates random initial parameter values, runs the AMOC-emulator to obtain an AMOC evolution. That is used to calculate a fit between the AMOC-emulator and UVic AMOC evolution. Than an iterative procedure starts that starts looking for an optimal parameter set (that is with best fit) by applying random perturbation to the free parameters until a sufficient result is obtained.
Called by: none, main executable

runmodel_amoc.m:
Function that runs the actual 4-box AMOC model.
Called by: runmodel_coupledsystem.m

runmodel_coupledsystem.m:
Function that loops over runmodel_ebm.m and runmodel_amoc.m for all the different RCP scenarios.
Called by: M-AMOC1.0.m

runmodel_ebm.m:
Function that runs the actual 1D atmospheric energy balance model.
Called by: runmodel_coupledsystem.m

Total_atmos_effect.m:
Function that transferres GCM-based regional temperature evolution into 'total atmosphere effect' time series.
Called by: runmodel_coupledsystem.m

Forcing data:
The folder Forcing_data contains example forcing data based on UVic simulations that can be used to test the model. The folder also contains a README file specifically for the forcing data.

