@은진,
The line
idx = find(strcmp(NAME, major_constituents{i}));
does not work, because
strcmp(NAME, major_constituents{i})
returns a 0, even when there is a match. To fix this, I converted NAME to a cell array, and used strcmp on the cell array:
NAMEc=cellstr(NAME); % convert character array NAME to cell array
...
idx = find(strcmp(NAMEc, major_constituents{i}));
The line
inclination = TIDECON(idx, 5) * pi / 180; % 도를 라디안으로 변환
does not work, because TIDECON has only 4 columns. Column 4 is the inclination. Therefore I change the line to
inclination = TIDECON(idx, 4) * pi / 180; % 도를 라디안으로 변환
I do not have function ellipsedraw(), so I wrote my own function, drawEllipse(). I call it:
drawEllipse(major_axis, minor_axis, x0, y0, inclination, colors{i});
It works. Now the script runs. There is a warning error at line
legend(major_constituents);
because major_constituents has four strings, but the plot has only two sets of data. Therfore I modified the code to make a cell array of strings, legstr , which includes only names of the tidal components plotted. The script, with the changes described above, is attached. It runs without error. It makes the plot shown below.
Only two components are plotted, because you search for components M2, K1, S2, O1. Of these, only M2 and K1 are among the eight components returned by t_tide, for this data set.
