/[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.9 by gforget, Sun Aug 3 19:49:58 2014 UTC revision 1.11 by gforget, Sun Jan 18 14:07:15 2015 UTC
# Line 24  gcmfaces_global; Line 24  gcmfaces_global;
24  doClim=0;  doClim=0;
25    
26  %directory names  %directory names
27  listDirs={'STATE/','TRSP/'};  listDirs={'STATE/','TRSP/'};%BUDG?
28  filAvailDiag=[dirModel 'available_diagnostics.log'];  filAvailDiag=[dirModel 'available_diagnostics.log'];
29  filReadme=[dirModel 'README'];  filReadme=[dirModel 'README'];
30  dirOut=[dirModel 'nctiles_tmp/'];  dirOut=[dirModel 'nctiles_tmp/'];
# Line 87  tim=datenum(tim)-datenum([1992 1 0]); Line 87  tim=datenum(tim)-datenum([1992 1 0]);
87  timUnits='days since 1992-1-1 0:0:0';  timUnits='days since 1992-1-1 0:0:0';
88    
89  %get time step axis  %get time step axis
90  [listTimes]=diags_list_times({[dirModel 'diags/STATE/']},{'state_2d_set1'});  [listTimes]=diags_list_times({dirIn},{fileModel});
91    
92  %get units and long name from available_diagnostics.log  %get units and long name from available_diagnostics.log
93  [avail_diag]=read_avail_diag(filAvailDiag,nameDiag);  [avail_diag]=read_avail_diag(filAvailDiag,nameDiag);
# Line 107  myFile=[myFile '/' nameDiag];%second ins Line 107  myFile=[myFile '/' nameDiag];%second ins
107  %apply mask, and convert to land mask  %apply mask, and convert to land mask
108  msk=grid_diag.msk;  msk=grid_diag.msk;
109  if length(size(myDiag{1}))==3;  if length(size(myDiag{1}))==3;
110    msk=repmat(msk,[1 1 size(myDiag{1},3)]);    msk=repmat(msk(:,:,1),[1 1 size(myDiag{1},3)]);
111  end;  end;
112  for kk=1:size(myDiag{1},4);  for kk=1:size(myDiag{1},4);
113    myDiag(:,:,:,kk)=myDiag(:,:,:,kk).*msk;    myDiag(:,:,:,kk)=myDiag(:,:,:,kk).*msk;
# Line 130  if doClim; Line 130  if doClim;
130    avail_diag.longNameDiag=[avail_diag.longNameDiag ' (climatology) '];    avail_diag.longNameDiag=[avail_diag.longNameDiag ' (climatology) '];
131  end;  end;
132    
133  %create netcdf file using write2nctiles (works only with old matlab, thus far ...)  %create netcdf file using write2nctiles
134  doCreate=1;  doCreate=1;
135  dimlist=write2nctiles(myFile,myDiag,doCreate,{'tileNo',tileNo},...  dimlist=write2nctiles(myFile,myDiag,doCreate,{'tileNo',tileNo},...
136      {'fldName',nameDiag},{'longName',avail_diag.longNameDiag},...      {'fldName',nameDiag},{'longName',avail_diag.longNameDiag},...
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22