Matlab--TS Plots

SO503, Spring 2012

Matlab TS Plots

You are to produce two figures for the following data set, which is a single profile from the NODC site below.

A TS plot with the station profile, and contours of sigma-tee.
A map showing the location of the station. You should look up the commands geoshow and scatterm

You will have writeup which correctly references the source of the data, correctly refers to the two figures, and briefly interprets what the TS diagram shows. This should be short and concise.

Your program should be able to handle a second data set with the same format, but a different number of points, solely by changing the name of the file to read. You will be given that data file at a later time.

Resources:

Matlab sigma-tee function: http://globec.whoi.edu/globec-dir/sigmat-calc-matlab.html. You should know what to do with a function.
CTD data file source: http://www.nodc.noaa.gov/GTSPP/gtspp-home.html
File format documentation: http://www.nodc.noaa.gov/GTSPP/document/datafmt/medsfmt.html

Other web resources include http://web.cecs.pdx.edu/~gerry/MATLAB/plotting/loadingPlotData.html
Insure that you understand what a line like "y(r) = str2num(readin(22:30));" does. You need to use code like this because the file has characters in it, and not just numbers.

This will be due Tuesday 24 April at 1330.

Hints:

File manipulation:

id=fopen('C:\mydocs\Courses\SO503\Labs\matlab_sigma_t\profile2.txt');
while ~feof(fid)
readin = fgetl(id);
end;
fclose(id);

This will work to read into a string:

param1=readin(133:136);

Trying to make param(1) is very hard. I would use an if statement:

if i==1

param1=readin(133:136);

else

param2=readin(133+17:136+17);

end

As a starting point for the map:

worldmap world

geoshow('landareas.shp', 'FaceColor', [0.5 1.0 0.5]);

axesm('MapProjection','robinson','Frame','off','Grid','off')

scatterm(lat,long);