/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_IO/rdmds2gcmfaces.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/matlab_class/gcmfaces_IO/rdmds2gcmfaces.m

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

revision 1.1 by gforget, Wed Feb 10 14:46:49 2010 UTC revision 1.5 by gforget, Fri Mar 17 17:41:56 2017 UTC
# Line 1  Line 1 
1  function [fldOut]=rdmds2gcmfaces(fileName,varargin);  function [fldOut,IT,M]=rdmds2gcmfaces(varargin);
2    %object:    read with rmds then apply convert2gcmfaces
3    %input:     varargin are the options to pass to rdmds (type help rdmds)
4    %output:    fldOut is a gcmfaces object
5    %
6    %note:      an earlier version was expecting nFaces to be passed
7    %           as the last argument; this is not the case anymore.
8    
9  if nargin==3; tt=varargin{1}; nFaces=varargin{2};  gcmfaces_global;
10  elseif nargin==2; tt=varargin{1}; nFaces=5;  
11  else; tt=[]; nFaces=5;  [v0,IT,M]=rdmds(varargin{1:end});
 end;  
   
 if ~isempty(tt); v0=rdmds(fileName,tt); else; v0=rdmds(fileName); end;  
12    
13  fldOut=convert2gcmfaces(v0,nFaces);  nn=size(v0);
14    test1=isfield(mygrid,'xtrct');
15    test1=test1&(prod(mygrid.ioSize)~=prod(nn(1:2)));
16    if test1;
17        if length(nn)==2; nn=[nn 1]; end;
18        v0=reshape(v0,[nn(1)*nn(2) nn(3:end)]);
19        v0=v0(mygrid.xtrct_inFull,:,:,:);
20        v0=reshape(v0,[mygrid.ioSize nn(3:end)]);
21    end;
22    
23    fldOut=convert2gcmfaces(v0);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22