how to turn netcdf into contour map

I have T2M data plot and want to put it into a contour plot at lat and lon 80.5N and -58W showing T2M over time how do I do this?

I've never used this before

Thanks

Answer 1

[Chad Greene]

Hi Abigail,

The first thing I do with a netcdf file is type

ncdisp('myfile.nc') 

to see what's inside it. If you have some variables named latitude, longitude, and T2M, read them in like this:

lon = ncread('myfile.nc','longitude'); 
lat = ncread('myfile.nc','latitude'); 
T = ncread('myfile.nc','T2M'); 

Most of the time with NetCDF's of climate data, you'll have to swap the first two dimensions of any gridded data. That probably looks like this:

T = permute(T,[2 1]); 

for 2D data or

T = permute(T,[2 1 3]); 

if T is a 3D data cube.

Assuming T is 3D where the third dimension corresponds to time, you could contour the first time slice like this:

contour(lon,lat,T) 

But you say you want to show T over time. What exactly do you mean by that?

[Abigail Waring]

Thank you

I managed to get it working but had to rewrite the code