% used NC files—————————————————————————————————————
% tavg.02021.01.01.nc:  model tavg output of experimental run
% tsi.02021.01.01.nc: model tai output of experimental run
% mask_c. nc: coastal sea mask, need to be copied from "\model runs\coast mask" folder 
% control_O_alk.nc: the total alkalinity from control run "RCP_8.5"
% control_O_dic.nc: the dissolved inorganic carbon from control run "RCP_8.5"

% Variables and their meanings——————————————————————————————
%  diss_ratio: the ratio of dissolved olivine to added olivine
%  alk_dis: dissolved olivine in alkalinity units [Pmol Alk yr-1]
%  alk_add: added olivine in alkalinity units [Pmol Alk yr-1]
%  a_co2: mean atmospheric CO2 concentration [ppm]
%  a_sat: mean atmospheric surface temperature concentration [degree-C]
%  omega: coastal-mean aragonite saturation omega
%  ph: coastal-mean pH
%  fdic_coast: coastal-mean CO2 flux [Pmol yr-1]
%  fdic_global: global-mean CO2 flux [Pmol yr-1]
%  o_oaragsur: global-mean aragonite saturation omega
% o_alk: oceanic total alkalinity [mol m-3]
% o_dic: ocean dissolved inorganic carbon [mol m-3]
% mask: coastal sea mask
% CDR: oceanic DIC change relative to control run [mol m-3]
% rato_x: COA induced oceanic CO2 uptake over alkalinity enhancement 

ferret -nojnl

use tavg.02021.01.01.nc 
use mask_c.nc
use tsi.02021.01.01.nc 
let alk_dis=f_olivint_dis[d=3]
let alk_add=f_olivint[d=3]
let diss_ratio=alk_dis[l=@ave]/alk_add[l=@ave]

list/file=olivX_con_X_diss_ratio.txt/nohead/norow diss_ratio
list/file=olivX_con_X_annual_alk_dis.txt/nohead/norow alk_dis
list/file=olivX_con_X_annual_alk_add.txt/nohead/norow alk_add
list/file=olivX_con_X_co2.txt/nohead/norow a_co2[d=3]
list/file=olivX_con_X_sat.txt/nohead/norow a_sat[d=3]

let omega=o_oaragsur[d=1]*mask[d=2]
list/file=olivX_con_X_omega.txt/nohead/norow omega[x=@ave,y=@ave]

let ph=o_phsur[d=1]*mask[d=2]
list/file=olivX_con_X_ph.txt/nohead/norow ph[x=@ave,y=@ave]

let temp=f_dic[d=1]*mask[d=2]
let fdic_coast=temp[x=@ave,y=@ave]*365*24*3600*3.11*1e13/1e15
let fdic_glob=f_dic[d=3]*365*24*3600*3.6*1e14/1e15
list/file=fdic_coast_x.txt/nohead/norow fdic_coast
list/file=fdic_glob_x.txt/nohead/norow fdic_glob


list/file=omega_glob_olivX_con_X.txt/nohead/norow o_oaragsur[x=@ave,y=@ave]


save/file=olivX_con_X_alk_x.nc o_alk[d=1]
save/file=olivX_con_X_dic_x.nc o_dic[d=1]

exit

%-------------------calculate the COA ratio: CO2/TA

ferret -nojnl

use control_O_alk.nc
use control_O_dic.nc
use olivX_con_X_alk_x.nc
use olivX_con_X_dic_x.nc


let deltaTA=o_alk[d=3]-o_alk[d=1]
let deltaDIC=o_dic[d=4]-o_dic[d=2]
let ratio_x=deltaDIC[l=80,x=@din,y=@din,k=@din]/deltaTA[l=80,x=@din,y=@din,k=@din]

let CDR=deltaDIC[x=@din,y=@din,k=@din,l=80]*12/1e15

list/file=cdr_x.txt/nohead/norow CDR
list/file=ratio_x.txt/nohead/norow ratio_x

exit