%% Calulate HNLC POC and NPP close all; clearvars addpath data/ addpath utils/ % Load model surface NO3 load('fec_CESM2opt_var_annual_surface.mat'); lat = output.lat; lon = output.lon; tarea = output.TAREA; % variable NPP_var = output.NPP; POC_flux_var = output.POC_flux; Chl_var = output.Chl; load('fec_CESM2opt_var_lim_2004_annual_surface.mat'); NPP_var_lim = output.NPP; POC_flux_var_lim = output.POC_flux; Chl_var_lim = output.Chl; %fixed 10 load('fec_CESM2opt_fixed_10_annual_surface.mat') NPP_10 = output.NPP; POC_flux_10 = output.POC_flux; Chl_10 = output.Chl; %fixed 7 load('fec_CESM2opt_fixed_7_annual_surface.mat') NPP_7 = output.NPP; POC_flux_7 = output.POC_flux; Chl_7 = output.Chl; %fixed 5 load('fec_CESM2opt_fixed_5_annual_surface.mat') NPP_5 = output.NPP; POC_flux_5 = output.POC_flux; Chl_5 = output.Chl; %fixed 3 load('fec_CESM2opt_fixed_3_annual_surface.mat') NPP_3 = output.NPP; POC_flux_3 = output.POC_flux; Chl_3 = output.Chl; load('data/HNLC_grid.mat') HNLCmask = HNLC.HNLCmask_gx3v7; POC_flux_HNLC_var = POC_flux_var.*HNLCmask; NPP_HNLC_var = NPP_var.*HNLCmask; Chl_HNLC_var = Chl_var.*HNLCmask; POC_flux_HNLC_var_lim = POC_flux_var_lim.*HNLCmask; NPP_HNLC_var_lim = NPP_var_lim.*HNLCmask; Chl_HNLC_var_lim = Chl_var_lim.*HNLCmask; POC_flux_HNLC_10 = POC_flux_10.*HNLCmask; NPP_HNLC_10 = NPP_10.*HNLCmask; Chl_HNLC_10 = Chl_10.*HNLCmask; POC_flux_HNLC_7 = POC_flux_7.*HNLCmask; NPP_HNLC_7 = NPP_7.*HNLCmask; Chl_HNLC_7 = Chl_7.*HNLCmask; POC_flux_HNLC_5 = POC_flux_5.*HNLCmask; NPP_HNLC_5 = NPP_5.*HNLCmask; Chl_HNLC_5 = Chl_5.*HNLCmask; POC_flux_HNLC_3 = POC_flux_3.*HNLCmask; NPP_HNLC_3 = NPP_3.*HNLCmask; Chl_HNLC_3 = Chl_3.*HNLCmask; POC_tot_HNLC_var= sum(tarea.*POC_flux_HNLC_var,'all','omitnan')*1e-15; POC_tot_HNLC_var_lim= sum(tarea.*POC_flux_HNLC_var_lim,'all','omitnan')*1e-15; POC_tot_HNLC_10= sum(tarea.*POC_flux_HNLC_10,'all','omitnan')*1e-15; POC_tot_HNLC_7= sum(tarea.*POC_flux_HNLC_7,'all','omitnan')*1e-15; POC_tot_HNLC_5= sum(tarea.*POC_flux_HNLC_5,'all','omitnan')*1e-15; POC_tot_HNLC_3= sum(tarea.*POC_flux_HNLC_3,'all','omitnan')*1e-15; NPP_tot_HNLC_var= sum(tarea.*NPP_HNLC_var,'all','omitnan')*1e-15; NPP_tot_HNLC_var_lim= sum(tarea.*NPP_HNLC_var_lim,'all','omitnan')*1e-15; NPP_tot_HNLC_10= sum(tarea.*NPP_HNLC_10,'all','omitnan')*1e-15; NPP_tot_HNLC_7= sum(tarea.*NPP_HNLC_7,'all','omitnan')*1e-15; NPP_tot_HNLC_5= sum(tarea.*NPP_HNLC_5,'all','omitnan')*1e-15; NPP_tot_HNLC_3= sum(tarea.*NPP_HNLC_3,'all','omitnan')*1e-15;