netcdf to contour plot SST

Trying to turn my SST netcdf plot into a contour plot within my region at a certain time and date?

Not sure what to do next? This is the contents

% Choose netcdf file
% nc_file = '../data/ERA/test 1979-2009.nc'; % put the full path to the netcdf file between the inverted commas
nc_file = '../data/SST/SSTA_data.nc'; % put the full path to the netcdf file between the inverted commas


%% Look at netcdf file metadata

ncdisp(nc_file);


%% Open netcdf file

disp(['Loading ' nc_file]);
ncid = netcdf.open(nc_file);
Source:
           \\studata06\home\GG\Ggp18amw\ManWin\My Documents\MATLAB\Thesis Data\netcdf_example_code (1)\netcdf_example_code\data\SST\SSTA_data.nc
Format:
           classic
Dimensions:
           T = 473
           X = 1
           Y = 1
Variables:
    T   
           Size:       473x1
           Dimensions: T
           Datatype:   single
           Attributes:
                       standard_name = 'time'
                       pointwidth    = 1
                       calendar      = '360'
                       gridtype      = 0
                       units         = 'months since 1960-01-01'
    X   
           Size:       1x1
           Dimensions: X
           Datatype:   single
           Attributes:
                       standard_name = 'longitude'
                       pointwidth    = 1
                       modulus       = 360
                       period        = 360
                       gridtype      = 0
                       units         = 'degree_east'
    Y   
           Size:       1x1
           Dimensions: Y
           Datatype:   single
           Attributes:
                       standard_name = 'latitude'
                       pointwidth    = 1
                       gridtype      = 0
                       units         = 'degree_north'
    ssta
           Size:       1x1x473
           Dimensions: X,Y,T
           Datatype:   single
           Attributes:
                       pointwidth     = 1
                       expires        = 1620000000
                       missing_value  = NaN
                       units          = 'degree_Celsius'
                       long_name      = 'Sea Surface Temperature Anomaly'
                       standard_name  = 'sea_surface_temperature'
                       CE             = 4.2
                       scale_min      = -4.2
                       maxncolor      = 254
                       ncolor         = 254
                       CS             = -4.2
                       scale_max      = 4.2
                       colormap       = '[null 12632256 8388608 [16711680 66] [16760576 43] [15453831 16] 13959039 [13959039 7] 64636 [65535 4] [36095 24] [255 31] [128 60] [128 7] 2763429]'
                       colorscalename = 'sstacolorscale'