Plot the meridional extent of fronts associated with the Antarctic Circumpolar Current as described by Alex Orsi and others in their landmark work, "On the meridional extent and fronts of the Antarctic Circumpolar Current, Deep-Sea Res. I, 42, 641-673, 1995".
This script downloads data from the World Ocean Circulation Experiments Southern Ocean Atlas, which can be found online at http://woceatlas.tamu.edu/Sites/html/atlas/SOA_DATABASE_DOWNLOAD.html.
This function is now a plugin for Antarctic Mapping Tools.
Chad Greene (2019). Plot Antarctic Circumpolar Current Fronts (https://www.mathworks.com/matlabcentral/fileexchange/40627-plot-antarctic-circumpolar-current-fronts), MATLAB Central File Exchange. Retrieved .
1.1.0.0 | A cleaner rewrite as a plugin for Antarctic Mapping Tools. Example file now included. |
Inspired by: ncpolarm, Antarctic Mapping Tools
Inspired: SODB Southern Ocean Database Toolbox, Argo Toolbox
Create scripts with code, output, and formatted text in a single executable document.
Lavnish Gupta (view profile)
Hello Chad,
Could you please tell me how can superimpose this on my pcolor plot:
mymap=pcolor(longitude,latitude,G');
mymap.EdgeAlpha=0;
load coast;
hold on;
plot(long,lat,'k');
plot(long+360,lat,'k');
colorbar;
xlabel('Longitude');
ylabel('Latitude');
title('Sea Surface Temperature Tendency');
I want to kind of overlay the acc lines on my map but it is erasing my map. Any ideas?
Chad Greene (view profile)
Oh dang, yes, I've heard from a number of folks that the TAMU site no longer seems to be supporting the SODB. Feel free to email me directly and I'll send you the zip file containing the ACC front data.
Vince Clementi (view profile)
Hi Chad,
Thanks for the assistance. However, I am getting a notification that the datasets (pf.txt, etc) cannot be found. When I go to the tamu link above, I find that the link is dead. Any idea how I can download the front data? Thanks!
Chad Greene (view profile)
Hi Vincent,
If you have Matlab's Mapping Toolbox, this'll do it:
worldmap([-60 -25],[155 190])
geoshow('landareas.shp')
pf = importdata('pf.txt');
plotm(pf.data(:,2),pf.data(:,1),'r-');
saf = importdata('saf.txt');
plotm(saf.data(:,2),saf.data(:,1),'b-');
stf = importdata('stf.txt');
plotm(stf.data(:,2),stf.data(:,1),'g-');
If not, you can use my borders function (also on file exchange) like this:
borders('new zealand','facecolor',rgb('tan'),'nomappingtoolbox')
axis([155 190 -60 -25])
pf = importdata('pf.txt');
plot(pf.data(:,1),pf.data(:,2),'r-');
saf = importdata('saf.txt');
plot(saf.data(:,1),saf.data(:,2),'b-');
stf = importdata('stf.txt');
plot(stf.data(:,1),stf.data(:,2),'g-');
Let me know if you have any other questions.
Vince Clementi (view profile)
I should elaborate--looking to plot the PF, SAF, STF on my map of New Zealand. Thanks.
Vince Clementi (view profile)
Brilliant package here, Chad (as well as your Antarctic Mapping Tools toolbox). Wondering if there is a way to overlay the antarctic fronts on a map I've created using this package without the Antarctic Mapping Tools kit? Cheers.