http://www.wiki.mohid.com/index.php?action=history&feed=atom&title=Matlab_samples 2026-04-04T20:47:48Z Revision history for this page on the wiki MediaWiki 1.28.0 http://www.wiki.mohid.com/index.php?title=Matlab_samples&diff=485&oldid=prev 2008-12-03T10:27:16Z

<p>1 revision</p> <table class="diff diff-contentalign-left" data-mw="interface"> <tr style='vertical-align: top;' lang='en'> <td colspan='1' style="background-color: white; color:black; text-align: center;">← Older revision</td> <td colspan='1' style="background-color: white; color:black; text-align: center;">Revision as of 10:27, 3 December 2008</td> </tr><tr><td colspan='2' style='text-align: center;' lang='en'><div class="mw-diff-empty">(No difference)</div> </td></tr></table>

Guillaume http://www.wiki.mohid.com/index.php?title=Matlab_samples&diff=484&oldid=prev 2008-07-08T09:14:53Z

<p>‎<span dir="auto"><span class="autocomment">Read and write netcdf files</span></span></p> <p><b>New page</b></p><div>This article compiles some useful [[matlab]] samples to use with the MOHID converted [[netcdf]] files.<br /> <br /> ==Netcdf MOHID vector plot==<br /> Sample code for the Portuguese coast application. It integrates several layers of the horizontal velocity. And plots and saves it.<br /> <br /> file='F:\Results200606\20060601_Hydrodynamic_1.nc';<br /> nc_dump(file)<br /> <br /> lon = nc_varget(file,'lon')<br /> lat = nc_varget(file,'lat')<br /> depth = nc_varget(file,'depth')<br /> <br /> siz = [length(lon)-1 length(lat)-1];<br /> lon2d = zeros(siz); <br /> lat2d = zeros(siz);<br /> <br /> for i=1:siz(1) <br /> for j=1:siz(2)<br /> lon2d(i,j) = lon(i); <br /> lat2d(i,j) = lat(j);<br /> end<br /> end<br /> <br /> SumVelX = 0; <br /> SumVelY = 0;<br /> <br /> for k=35:42 <br /> <br /> Vel_xk = nc_varget(file,'Vel_X',[k 1 1],[1 116 176]); <br /> SumVelX = SumVelX + Vel_xk * (depth(k+1) - depth(k)) / (depth(43) - depth(36)); <br /> <br /> Vel_yk = nc_varget(file,'Vel_Y',[k 1 1],[1 116 176]); <br /> SumVelY = SumVelY + Vel_yk * (depth(k+1) - depth(k)) / (depth(43) - depth(36));<br /> <br /> end<br /> <br /> %plot arrows<br /> m_proj('miller', 'lon', [min(lon) max(lon)], 'lat', [min(lat) max(lat)]);<br /> m_quiver(lon2d(1:3:end), lat2d(1:3:end), Vel_xk(1:3:end), Vel_yk(1:3:end), 0, 'k');<br /> hold on;<br /> m_grid('box','fancy','linestyle','none','fontsize', 7);<br /> m_usercoast('[http://www.mohid.com/HydroGroup/data_sources/coastline.mat coastline.mat]','patch', [.5 .5 .5]); <br /> %plot scale<br /> m_quiver(-8.5, 37.5, .1, 0.0, 0, 'k');<br /> htv5 = text(.015, .69, '10 cm s^{-1}');<br /> set(htv5,'FontSize',8);<br /> <br /> %save figure as a png file<br /> saveas(gcf, 'MyArrowsFigure.png', 'png');<br /> <br /> ==Read and write netcdf files==<br /> <br /> %Read level2 residual fields<br /> %%%%%%%%%%%% 3D results %%%%%%%%%%%%%%%%%%%%%<br /> file = 'BiscayResidual.nc';<br /> u = nc_varget(file,'Vel_X');<br /> v = nc_varget(file,'Vel_Y');<br /> lon = nc_varget(file,'lon');<br /> lat = nc_varget(file, 'lat');<br /> depth = nc_varget(file, 'depth');<br /> level = nc_varget(file, 'Lev_Z');<br /> <br /> %Create level1 residual fields<br /> out = 'MOHIDresidual3D.nc';<br /> nc_create_empty(out, nc_write_mode);<br /> nc_add_dimension(out, 'lon', length(lon));<br /> nc_add_dimension(out, 'lat', length(lat));<br /> nc_add_dimension(out, 'deptht', length(depth));<br /> <br /> s_ncvars = struct ( 'Name', ...<br /> { 'lon',...<br /> 'lat',...<br /> 'deptht',...<br /> 'u',...<br /> 'v',...<br /> 'elevation'}, ...<br /> 'Nctype', ...<br /> { 'float',...<br /> 'float',...<br /> 'float',...<br /> 'float',...<br /> 'float',...<br /> 'float'}, ...<br /> 'Dimension', ...<br /> { {'lon'},...<br /> {'lat'},...<br /> {'deptht'},...<br /> {'deptht' 'lat' 'lon'},...<br /> {'deptht' 'lat' 'lon'},...<br /> {'lat' 'lon'}}, ...<br /> 'Attribute', ...<br /> { struct( 'Name','units',...<br /> 'Value','degrees_east' ), ...<br /> struct( 'Name','units',...<br /> 'Value','degrees_north' ), ...<br /> struct( 'Name','units',...<br /> 'Value','m' ), ...<br /> struct( 'Name','units',...<br /> 'Value','m s-1' ), ...<br /> struct( 'Name','units',...<br /> 'Value','m s-1' ), ...<br /> struct( 'Name','units',...<br /> 'Value','m' ), ...<br /> })<br /> <br /> for i=1:length(s_ncvars) <br /> nc_addvar(out, s_ncvars(i));<br /> end<br /> <br /> nc_dump(out);<br /> <br /> ==See also==<br /> *[[Matlab]]<br /> <br /> [[Category:Matlab]]</div>

192.168.20.177