Contour maps from Netcdf files

I want to make a contour map of the TCO variable for 10 January 2007 from a Netcdf file. I also only want to display southern hemisphere latitudes.

Link to Netcdf file: https://drive.google.com/file/d/0B0wQbicZOF7xcWx0RmtoXzZHWHM/view

I have tried the code below but it doesn't work.

ncdisp('BodekerScientificCombinedOzoneV3.0_2007_Unpatched.nc')
        ncid1 = netcdf.open('BodekerScientificCombinedOzoneV3.0_2007_Unpatched.nc','NC_NOWRITE');
        TCO1 = netcdf.getVar(ncid1,3,[0 0 0],[288 180 1]);
        lon1 = netcdf.getVar(ncid1,1,0,180);
        lat1 = netcdf.getVar(ncid1,0,0,288);
        for p = 1:180
        for q = 1:288
        map1(q,p) = TCO1(p,q);
        end
        end
        contour(lon1,lat1,map1)

Answer 1

[Star Strider]

You can eliminate the loop with the transpose function or operator ('):

map1 = TCO1';

You either need to reverse the first two arguments to use the transposed matrix:

contourf(lat1,lon1,map1)

or not transpose it to use your original code.

[David du Preez]

Thank you. The first part makes sense.

Which lines are you referring to by saying "the first two arguments"

[Star Strider]

My pleasure.

Those refer to the contourf call. Generically:

contourf(Arg1, Arg2, Matrix)

so the first two arguments (inputs) are the vectors (or matrices) that create the x- and y-axes scales.