matlab_class/gcmfaces_IO/rdmds2gcmfaces.m

function [fldOut]=rdmds2gcmfaces(varargin);
%object:    read with rmds then apply convert2gcmfaces
%input:     varargin are the options to pass to rdmds (type help rdmds)
%output:    fldOut is a gcmfaces object
%
%note:      an earlier version was expecting nFaces to be passed
%           as the last argument; this is not the case anymore.

gcmfaces_global;

v0=rdmds(varargin{1:end});

nn=size(v0);
test1=isfield(mygrid,'xtrct');
test1=test1&(prod(mygrid.ioSize)~=prod(nn(1:2)));
if test1;
    if length(nn)==2; nn=[nn 1]; end;
    v0=reshape(v0,[nn(1)*nn(2) nn(3:end)]);
    v0=v0(mygrid.xtrct_inFull,:,:,:);
    v0=reshape(v0,[mygrid.ioSize nn(3:end)]);
end;

fldOut=convert2gcmfaces(v0);
1	gforget	1.3	function [fldOut]=rdmds2gcmfaces(varargin);
2	gforget	1.4	%object: read with rmds then apply convert2gcmfaces
3	gforget	1.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	gforget	1.1
9	gforget	1.4	gcmfaces_global;
10
11			v0=rdmds(varargin{1:end});
12
13			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	gforget	1.3	fldOut=convert2gcmfaces(v0);