/[MITgcm]/MITgcm_contrib/enderton/Diagnostics/DiagUtility/rdmnc_mod2.m
ViewVC logotype

Diff of /MITgcm_contrib/enderton/Diagnostics/DiagUtility/rdmnc_mod2.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by enderton, Mon Sep 12 20:33:26 2005 UTC revision 1.3 by enderton, Tue Oct 18 19:22:15 2005 UTC
# Line 1  Line 1 
1  function [S] = rdmnc_mod(varargin)  function [S] = rdmnc_mod2(varargin)
2    
3  % Usage:  % Usage:
4  %   S=RDMNC(FILE1,FILE2,...)  %   S=RDMNC(FILE1,FILE2,...)
# Line 14  function [S] = rdmnc_mod(varargin) Line 14  function [S] = rdmnc_mod(varargin)
14  %   ITER    Vector of iterations in the MNC/netcdf files, not model time.  %   ITER    Vector of iterations in the MNC/netcdf files, not model time.
15  %  %
16  % Output:  % Output:
17  %   S       Structure with mutliple fields  %   S       Structure with fields corresponding to 'VAR1', 'VAR2', ...
18  %  %
19  % Description:  % Description:
20  %   This function rudimentary wrapper for joining and reading netcdf files  %   This function is a rudimentary wrapper for joining and reading netcdf
21  %   produced by MITgcm.  It does not give the same flexibility as opening  %   files produced by MITgcm.  It does not give the same flexibility as
22  %   the netcdf files directly using netcdf(). It is useful for quick  %   opening the netcdf files directly using netcdf(), but is useful for
23  %   loading of entire model fields which are distributed in multiple netcdf  %   quick loading of entire model fields which are distributed in multiple
24  %   files.  %   netcdf files.
25  %  %
26  % Example:  % Example:
27  %   >> S=rdmnd('state.*','XC','YC','RC','T');  %   >> S=rdmnd('state.*','XC','YC','T');
28  %   >> imagesc( S.XC, S.YC, S.T(:,:,1)' );  %   >> imagesc( S.XC, S.YC, S.T(:,:,1)' );
29    %
30    %  Author:  Alistair Adcroft
31    %  Modifications:  Daniel Enderton
32    
33  % Initializations  % Initializations
34  file={};  file={};
# Line 100  end Line 103  end
103    
104  function [A] = read_att(nc);  function [A] = read_att(nc);
105      allatt=ncnames(att(nc));      allatt=ncnames(att(nc));
106      A='none';      if ~isempty(allatt)
107      for attr=allatt;          for attr=allatt;
108          A.(char(attr))=nc.(char(attr))(:);              A.(char(attr))=nc.(char(attr))(:);
109            end
110        else
111            A = 'none';
112      end      end
113    
114  function [i0,j0,fn] = findTileOffset(S);  function [i0,j0,fn] = findTileOffset(S);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22