Plot netcdf variable over south pole

I have a netcdf file. https://drive.google.com/file/d/0B0wQbicZOF7xa053WUJ5Q19YWEU/view

I want to plot the tco variable of a specific day over the southern hemisphere between the latitudes -90 to -10 using a stereographic projection. I want to include the coastlines in the plot as well.

Is it possible to do this and how?

Answer 1

[KSSV]

load the variable from ncfile using ncread and plot your data...it is a straight forward task.

[David du Preez]

Sure. I get this far.

tco07=netcdf.getVar(ncid1,3,[0 0 0],[288 180 1]); lon1=netcdf.getVar(ncid1,1,0,288); lat1=netcdf.getVar(ncid1,2,0,180);

How do you create the contour of the variable?

[KSSV]

MAke the respective spatial data...I mean read lon and lat also and call contour

[David du Preez]

This is what I have:

ncid1=netcdf.open('NIWA-BSCombinedOzoneV3.2_2007_Unpatched.nc','NC_NOWRITE')  
  tco07=netcdf.getVar(ncid1,3,[0 0 0],[288 80 1]);
  lon1=netcdf.getVar(ncid1,1,0,288);
  lat1=netcdf.getVar(ncid1,2,0,80);
  tco07=tco07.';
  contour(lon1,lat1,tco07)
  hold on
  load coast
  axesm('stereo','Origin',[-89.5 -120],'MapLatLimit',[-89.5 -10])

It does not change the projection or show the coastlines.