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	function [fldOut]=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	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	fldOut=convert2gcmfaces(v0);