;
zp2=t2_grd.z(:,ii,jj);
tp2=t2(:,ii,jj);
</pre>
<h2>Plot data<a name="4"></a></h2><pre class="codeinput"><span class="comment">% make the comparison plot</span>
plot(t,t_grd.z,tp1,zp1,tp2,zp2);
thicker;grid;
legend(<span class="string">'Data'</span>,<span class="string">'UMAINE'</span>,<span class="string">'WHOI'</span>);
xlabel(<span class="string">'Temp (C)'</span>);
ylabel(<span class="string">'Depth (m)'</span>);
title(sprintf(<span class="string">'Gulf of Maine: Time:%s, Lon:%8.4f, Lat:%8.4f