"Final update and updated paper reference"

Author: messiem

README.md

@@ -1,6 +1,6 @@
 # sourcecode_GAexport
 
-This toolbox contains the programs used to reproduce results and figures from Messié et al. (submitted).   
+This toolbox contains the programs used to reproduce results and figures from Messié et al. (2025).   
 The programs are written for Matlab. The Statistics toolbox is required for multilinear regression in `ga_validate_gridded_export_vs_SedTraps` (l. 93).
 
 * * *
@@ -16,14 +16,14 @@ The code generates figures, saved in `figures/`, reproducing results from Fig. 1
 ### Description of functions ### 
 
 `ga_analyze_SedTraps_trajectories`: reproduces results using backward Lagrangian trajectories originating from the CCE-LTER sediment trap positions. Trajectories are available in `data/CCE_SedTrap_trajectories.mat`.  
-`ga_validate_gridded_export_vs_SedTraps`: reproduces the validation of the GA gridded product published on Zenodo, against the CCE-LTER dataset. The program also displays results using the EF-OC dataset (available online at https://spg-satdata.ucsd.edu/wc_productivity/wc_productivity.htm). 
+`ga_validate_gridded_export_vs_SedTraps`: reproduces the validation of the GA gridded product published on Zenodo (https://doi.org/10.5281/zenodo.14758969), against the CCE-LTER dataset. The program also displays results using the EF-OC satellite export dataset (available online at https://spg-satdata.ucsd.edu/wc_productivity/wc_productivity.htm). 
 
 ### Description of data ###
 
 **Data to download and save in the `data/` folder:** 
 
 `Nsupply_California_CCMP_REP.nc`: download from https://doi.org/10.5281/zenodo.14641977 (nutrient supply used to force the plankton model)   
-`export_CCMP3km_GlobCurrent_daily_REP.nc`: download from https://doi.org/10.5281/zenodo.14084208 (output of the GA model, forced by Nsupply and GlobCurrent).
+`export_CCMP3km_GlobCurrent_daily_REP.nc`: download from https://doi.org/10.5281/zenodo.14758969 (output of the GA model, forced by Nsupply and GlobCurrent).
 
 **CCE-LTER datasets:**  
 
@@ -43,7 +43,7 @@ backward Lagrangian trajectories originating from the CCE-LTER sediment trap pos
 
 Please refer this paper when using these scripts:  
 
-Messié, M., C.L. Huffard, M. Stukel, and H.A. Ruhl (2024). **Spatial and temporal interplay between oceanic circulation and biological production in shaping carbon export off the California coast**.  *Geophysical Research Letters*, submitted. Preprint available at https://doi.org/10.22541/essoar.173272956.60470361/v1
+Messié, M., C.L. Huffard, M. Stukel, and H.A. Ruhl (2025). **Spatial and temporal interplay between oceanic circulation and biological production in shaping carbon export off the California coast**.  *Geophysical Research Letters*, https://doi.org/10.1029/2024GL113707
 
 * * *
 

figures/Fig1.jpg

@@ -2,10 +2,13 @@
 
 
 %% GA_ANALYZE_SEDTRAPS_TRAJECTORIES: look at sediment trap source water characteristics 
-% This function recalculates water_age, coastal Nsupply, and along-trajectory modeled export.
-% Outputs are compared against CCE_SedTrap_Data.csv (differences are due to precision in the csv file).
+% This function recalculates water_age, coastal Nsupply, and along-trajectory modeled export,
+% and reproduces Fig. 1 & 2a in Messié et al. (2025) as well as results from Table 1.
+% Outputs are compared against CCE_SedTrap_Data.csv.
 %
 % ga_analyze_SedTraps_trajectories
+%
+% Monique Messié, 2024
 
 
 distcoast_lim=10;		% distance to coast used as "source data". This is a trade-off between too far from the coast, and unreliable/no current data near shore.
@@ -278,7 +281,7 @@
 % Scatter plot (Fig. 2a)
 figure, set(gcf,'Units','Centimeters'); pos=get(gcf,'Position');
 pos(3)=6; pos(4)=7; set(gcf,'Position',pos)
-axes('Position',[0.2 0.28 0.75 0.7],'FontSize',7), hold on
+axes('Position',[0.17 0.28 0.8 0.7],'FontSize',7), hold on
 for ipts=find(iok)'
 	plot([SedTrap.Carbon_flux_corr_mgm2day(ipts)-SedTrap.Carbon_flux_standard_error_corr_mgm2day(ipts)/2 ...
 		SedTrap.Carbon_flux_corr_mgm2day(ipts)+SedTrap.Carbon_flux_standard_error_corr_mgm2day(ipts)/2],...

figures/Fig2a.jpg

@@ -2,10 +2,13 @@
 
 
 %% GA_VALIDATE_GRIDDED_EXPORT_VS_SEDTRAPS: match gridded export product with the CCE-LTER sediment trap dataset
-% This function matches gridded GA products (K20 and GA) with the CCE-LTER dataset.
-% Outputs are compared against CCE_SedTrap_Data.csv (differences are due to precision in the csv file).
+% This function matches gridded export products (EF-OC and GA) with the CCE-LTER dataset, 
+% and reproduces Fig. 2b & c in Messié et al. (2025) as well as results from Table 1.
+% Outputs are compared against CCE_SedTrap_Data.csv.
 %
 % ga_validate_gridded_export_vs_SedTraps
+% 
+% Monique Messié, 2024
 
 
 
@@ -131,7 +134,7 @@
 	pos(3)=12; pos(4)=7; set(gcf,'Position',pos)
 	set(gcf,'PaperPositionMode','Auto','PaperUnits','Centimeters','PaperSize',[pos(3), pos(4)])
 	fontsize=7;
-axes('Position',[0.1 0.27 0.37 0.7],'FontSize',fontsize), hold on
+axes('Position',[0.08 0.27 0.4 0.7],'FontSize',fontsize), hold on
 	for ipts=find(iok)'
 		plot([SedTrap.Carbon_flux_corr_mgm2day(ipts) SedTrap.Carbon_flux_corr_mgm2day(ipts)],...
 		[SedTrap.export_EFOC_mgm2day(ipts) SedTrap.GAgridded_CZeu_mgm2day(ipts)],'Color',[0.7 0.7 0.7])
@@ -142,7 +145,7 @@
 	plot(xlim,xlim,'k')
 	xlabel({'{\it In situ} Zeu export [mg m^{-2} d^{-1}]'}), ylabel('Gridded export products [mg m^{-2} d^{-1}]')
 	legend([h1,h2],{'EF-OC','GA C_{Zeu}'},'FontSize',fontsize,'Location','northwest')
-axes('Position',[0.6 0.27 0.37 0.7],'FontSize',fontsize), hold on
+axes('Position',[0.58 0.27 0.4 0.7],'FontSize',fontsize), hold on
 	ibelow30m=SedTrap.Zeu_m>=30;
 	errSat=abs(SedTrap.export_EFOC_mgm2day*b(1)-SedTrap.GAgridded_CZeu_mgm2day*b(2))*sqrt(2);
 	for ipts=find(iok)'
@@ -160,7 +163,7 @@
 	plot(xlim,xlim,'k')
 	xlabel('{\it In situ} Zeu export [mgC m^{-2} d^{-1}]'), ylabel('Combined satellite export [mgC m^{-2} d^{-1}]')
 	hbar4=colorbar('SouthOutside'); hbar4.Title.String='Drift [km]';
-	hbar4.Position=[0.6 0.06 0.37 0.03];
+	hbar4.Position=[0.58 0.06 0.4 0.03];
 print('-djpeg','-r300','figures/Fig2bc.jpg')