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

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

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

revision 1.5 by gforget, Mon Feb 3 19:47:50 2014 UTC revision 1.8 by gforget, Mon Jul 28 20:54:11 2014 UTC
# Line 80  tim=[1992*ones(nn,1) [1:nn]' 15*ones(nn, Line 80  tim=[1992*ones(nn,1) [1:nn]' 15*ones(nn,
80  tim=datenum(tim)-datenum([1992 1 0]);  tim=datenum(tim)-datenum([1992 1 0]);
81  timUnits='days since 1992-1-1 0:0:0';  timUnits='days since 1992-1-1 0:0:0';
82    
83    %get time step axis
84    [listTimes]=diags_list_times({[dirModel 'diags/STATE/']},{'state_2d_set1'});
85    
86  %get units and long name from available_diagnostics.log  %get units and long name from available_diagnostics.log
87  [avail_diag]=read_avail_diag(filAvailDiag,nameDiag);  [avail_diag]=read_avail_diag(filAvailDiag,nameDiag);
88    
# Line 97  myFile=[myFile '/' nameDiag];%second ins Line 100  myFile=[myFile '/' nameDiag];%second ins
100    
101  %apply mask, and convert to land mask  %apply mask, and convert to land mask
102  msk=grid_diag.msk;  msk=grid_diag.msk;
103  msk=repmat(msk,[1 1 size(myDiag{1},3) size(myDiag{1},4)]);  if length(size(myDiag{1}))==3;
104  myDiag=myDiag.*msk;    msk=repmat(msk,[1 1 size(myDiag{1},3)]);
105    end;
106    for kk=1:size(myDiag{1},4);
107      myDiag(:,:,:,kk)=myDiag(:,:,:,kk).*msk;
108    end;
109  clear msk;  clear msk;
110  %  %
111  land=isnan(grid_diag.msk);  land=isnan(grid_diag.msk);
# Line 136  doCreate=0; Line 143  doCreate=0;
143  write2nctiles(myFile,tim,doCreate,{'tileNo',tileNo},...  write2nctiles(myFile,tim,doCreate,{'tileNo',tileNo},...
144    {'fldName','tim'},{'longName','time'},...    {'fldName','tim'},{'longName','time'},...
145    {'units',timUnits},{'dimIn',dim.tim});    {'units',timUnits},{'dimIn',dim.tim});
146    write2nctiles(myFile,listTimes,doCreate,{'tileNo',tileNo},...
147      {'fldName','step'},{'longName','final time step number'},...
148      {'units','1'},{'dimIn',dim.tim});
149  write2nctiles(myFile,grid_diag.lon,doCreate,{'tileNo',tileNo},...  write2nctiles(myFile,grid_diag.lon,doCreate,{'tileNo',tileNo},...
150    {'fldName','lon'},{'units','degrees_east'},{'dimIn',dim.twoD});    {'fldName','lon'},{'units','degrees_east'},{'dimIn',dim.twoD});
151  write2nctiles(myFile,grid_diag.lat,doCreate,{'tileNo',tileNo},...  write2nctiles(myFile,grid_diag.lat,doCreate,{'tileNo',tileNo},...
# Line 157  function [meta]=read_meta(fileName); Line 167  function [meta]=read_meta(fileName);
167    
168  %read meta file  %read meta file
169  tmp1=dir([fileName '*.meta']); tmp1=tmp1(1).name;  tmp1=dir([fileName '*.meta']); tmp1=tmp1(1).name;
170  tmp2=strfind(fileName,'/');  tmp2=strfind(fileName,filesep);
171  if ~isempty(tmp2); tmp2=tmp2(end); else; tmp2=0; end;  if ~isempty(tmp2); tmp2=tmp2(end); else; tmp2=0; end;
172  tmp1=[fileName(1:tmp2) tmp1]; fid=fopen(tmp1);  tmp1=[fileName(1:tmp2) tmp1]; fid=fopen(tmp1);
173  while 1;  while 1;
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22