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

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

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

revision 1.6 by gforget, Sun Jan 18 14:07:15 2015 UTC revision 1.7 by gforget, Wed Mar 16 15:02:19 2016 UTC
# Line 24  end; Line 24  end;
24    
25  doCheck=0;%set to one to print stuff to screen  doCheck=0;%set to one to print stuff to screen
26    
27    fldInIsaGcmfaces=isa(fldIn,'gcmfaces');
28    
29  %set more optional paramaters to default values  %set more optional paramaters to default values
30  descr='';  descr='';
31  rdm='';  rdm='';
# Line 71  tileList=unique(convert2vector(tileNo)); Line 73  tileList=unique(convert2vector(tileNo));
73  tileList=tileList(~isnan(tileList));  tileList=tileList(~isnan(tileList));
74  ntile=length(tileList);  ntile=length(tileList);
75  %  %
76  if isa(fldIn,'gcmfaces');  if fldInIsaGcmfaces;
77   for ff=1:ntile;   for ff=1:ntile;
78     tmp1=[];     tmp1=[];
79     for gg=1:mygrid.nFaces;     for gg=1:mygrid.nFaces;
# Line 84  if isa(fldIn,'gcmfaces'); Line 86  if isa(fldIn,'gcmfaces');
86     end;     end;
87     fldTiles{ff}=tmp1;     fldTiles{ff}=tmp1;
88   end;   end;
89     clear fldIn;
90  end;  end;
91    
92  %start processing loop  %start processing loop
93  for ff=1:ntile;  for ff=1:ntile;
94    
95  if isa(fldIn,'gcmfaces');  if fldInIsaGcmfaces;
96    fldTile=fldTiles{ff};    fldTile=fldTiles{ff};
97    %reverse order of dimensions    %reverse order of dimensions
98    nn=length(size(fldTile));    nn=length(size(fldTile));
99    fldTile=permute(fldTile,[nn:-1:1]);    fldTile=permute(fldTile,[nn:-1:1]);
100      if ntile==1; clear fldTiles; end;
101  else;  else;
102    fldTile=fldIn;    fldTile=fldIn;
103  end;  end;
# Line 138  dimOut{ff}=dimlist; Line 142  dimOut{ff}=dimlist;
142  if doCreate;  if doCreate;
143    %create netcdf file:    %create netcdf file:
144    %-------------------    %-------------------
145  % ncid=nccreate(fileTile,'NETCDF4');%to allow for big files    if prod(size(fldTile))*4/1e9<1.5;%use (always available) basic ncetcf:
146    ncid=nccreate(fileTile,'clobber');      mode='clobber';
147      else;%to allow for large file:
148        mode='NETCDF4';
149      end;
150      ncid=nccreate(fileTile,mode);
151    nc_global=netcdf.getConstant('NC_GLOBAL');    nc_global=netcdf.getConstant('NC_GLOBAL');
152    
153    if ~isempty(rdm);    if ~isempty(rdm);
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22