<p><b>ringler@lanl.gov</b> 2012-02-09 10:55:17 -0700 (Thu, 09 Feb 2012)</p><p><br>
add file to plot test case output data<br>
</p><hr noshade><pre><font color="gray">Added: branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/plot.m
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/plot.m                                (rev 0)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/plot.m        2012-02-09 17:55:17 UTC (rev 1490)
@@ -0,0 +1,47 @@
+%clear all
+
+nVertLevels = 20
+nx = 16
+ny = 50
+iLevel = 1
+var = 'temperature'
+
+ncid_hex = netcdf.open('10000m_20levs/output.nc','nc_nowrite');
+
+[dimname, dimlen] = netcdf.inqDim(ncid_hex,0);
+
+ID_hex = netcdf.inqVarID(ncid_hex,var);
+work = netcdf.getVar(ncid_hex,ID_hex);
+size(work)
+
+workq = work(1,:);
+size(workq)
+
+r = min(workq);
+zmin = min(r)
+r = max(workq);
+zmax = max(r)
+
+inc = (zmax - zmin) / 20;
+zlevs = zmin:inc:zmax
+
+hex = reshape(workq, nx, ny);
+size(hex)
+
+workavg=hex;
+for j=2:2:ny
+    for i=2:nx-1
+        hex(i,j,:) = (workavg(i+1,j,:)+workavg(i,j,:))/2.0;
+    end
+end
+clear workavg
+
+hexp = hex(1:nx,1:ny);
+b=hexp;
+figure(1)
+b = transpose(b);
+contourf(b,zlevs)
+caxis([zmin zmax])
+colormap(jet)
+colorbar
+

</font>
</pre>