MITgcm_contrib/jmc_script/check_rdmnc.m


% $Header:  $
% $Name:  $

%- matlab script to test mnc output file: 
%  use rdmnc to read netcdf output files and print dimensions of fields.

%rDir='fizhi-cs-aqualev20/build/';
%rDir='global_ocean.90x40x15/build/';
%rDir='hs94.cs-32x32x5/build/';
%rDir='ideal_2D_oce/build/';
%rDir='isomip/build/';
%rDir='isomip/tr_run.htd/';
%rDir='lab_sea/build/';
%rDir='MLAdjust/build/';
%rDir='MLAdjust/tr_run.0.leithD/';
%rDir='MLAdjust/tr_run.0.leith/';
%rDir='MLAdjust/tr_run.0.smag/';
%rDir='MLAdjust/tr_run.1.leith/';
%rDir='rotating_tank/build/';
%rDir='tutorial_barotropic_gyre/build/';
 rDir='internal_wave/build/';


cd(rDir);
here=pwd;
fprintf(['======== check output from: ',here,' ========\n']);
listF=dir('mnc_test_0001/*.t001.nc');

for n=1:size(listF,1),
  namf=listF(n).name;
  nam2=['mnc_test_0001/',namf(1:end-7),'*.nc'];
  fprintf(['file:>',nam2,'<\n']);
  S=rdmnc(nam2);
  listV=fieldnames(S);
  for i=1:size(listV,1), cvar=char(listV(i)); fprintf([' ',cvar]); end; fprintf('\n');
  j=0; for i=1:size(listV,1), cvar=char(listV(i)); siz=size(S.(cvar)); 
   if prod(siz) ~= siz(1), j=j+1; fprintf([' > ',cvar,' :']); 
       fprintf(' %i',size(S.(cvar))); fprintf(' < |');
       if rem(j,4) == 0, fprintf('\n'); end
   end
  end;
  if rem(j,4) ~= 0, fprintf('\n'); end
end
cd ../..
1	jmc	1.1
2			% $Header: $
3			% $Name: $
4
5			%- matlab script to test mnc output file:
6			% use rdmnc to read netcdf output files and print dimensions of fields.
7
8			%rDir='fizhi-cs-aqualev20/build/';
9			%rDir='global_ocean.90x40x15/build/';
10			%rDir='hs94.cs-32x32x5/build/';
11			%rDir='ideal_2D_oce/build/';
12			%rDir='isomip/build/';
13			%rDir='isomip/tr_run.htd/';
14			%rDir='lab_sea/build/';
15			%rDir='MLAdjust/build/';
16			%rDir='MLAdjust/tr_run.0.leithD/';
17			%rDir='MLAdjust/tr_run.0.leith/';
18			%rDir='MLAdjust/tr_run.0.smag/';
19			%rDir='MLAdjust/tr_run.1.leith/';
20			%rDir='rotating_tank/build/';
21			%rDir='tutorial_barotropic_gyre/build/';
22			rDir='internal_wave/build/';
23
24
25			cd(rDir);
26			here=pwd;
27			fprintf(['======== check output from: ',here,' ========\n']);
28			listF=dir('mnc_test_0001/*.t001.nc');
29
30			for n=1:size(listF,1),
31			namf=listF(n).name;
32			nam2=['mnc_test_0001/',namf(1:end-7),'*.nc'];
33			fprintf(['file:>',nam2,'<\n']);
34			S=rdmnc(nam2);
35			listV=fieldnames(S);
36			for i=1:size(listV,1), cvar=char(listV(i)); fprintf([' ',cvar]); end; fprintf('\n');
37			j=0; for i=1:size(listV,1), cvar=char(listV(i)); siz=size(S.(cvar));
38			if prod(siz) ~= siz(1), j=j+1; fprintf([' > ',cvar,' :']);
39			fprintf(' %i',size(S.(cvar))); fprintf(' < \|');
40			if rem(j,4) == 0, fprintf('\n'); end
41			end
42			end;
43			if rem(j,4) ~= 0, fprintf('\n'); end
44			end
45			cd ../..