## goal:
-To generate "figure 1" from model output "tavg.xx.xx.xx.nc" and "tai.xx.xx.xx.nc" for each individual model run

## related contents:
# sat.mat: global surface temperature matrix
# dis.mat: dissolved olivine (alkalinity units) matrix
# input.mat: added olivine (alkalinity units) matrix
# co2.mat: atmospheric CO2 matrix
# ph_global.mat: global-mean pH matrix
# ph.mat: coastal-mean pH matrix
# c_dic.mat: coastal-mean carbon flux matrix
# c_omega: coastal-mean omega matrix
# g_dic.mat: global-mean carbon flux matrix
# g_omega: coastal-mean omega matrix
# time_series.xlsx: above matrixes in Excel
# alk_dis.txt: matlab script to plot coastal-mean dissolved olivine time series (Alk units)
# alk_input.txt: matlab script to plot added olivine time series (Alk units)
# co2.txt: matlab script to plot atmospheric CO2 time series
# c_omega.txt matlab script to plot coastal-mean aragonite omega
# legend.txt: a legend example for time series.


## How to use it:
- (1) Use ferret script to generate the ASCII output for key variables. (script name: time_series.txt)
- (2) Convert results from (1) into the Excel and Matlab matrixes (*.mat and *.xlsx)
- (3) run matlab scripts (all the other *.txt files other than "time_series.txt") to generate the plots

## Note:
- (1) In Matlab matrix, the y axis indicates time from year 2020 to year 2099 with one year interval; x axis indicates model runs, namely "olivInf_con","oliv10_con","oliv100_con","oliv1000_con","olivInf_con_t","oliv10_con_t","oliv100_con_t","oliv1000_con_t","olivInf_omega3.4","oliv10_omega3.4","oliv100_omega3.4","oliv1000_omega3.4","olivInf_omega9","oliv10_omega9","oliv100_omega9","oliv1000_omega9","RCP_8.5" from left to right
- (2) the generated plots are the elements to create Indesign-based figure for publications. The refinery work in Indeisgn is not explained here.