/[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.12 by gforget, Tue Jan 20 12:01:19 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 96  timUnits='days since 1992-1-1 0:0:0'; Line 96  timUnits='days since 1992-1-1 0:0:0';
96  [rdm]=read_readme(filReadme);  [rdm]=read_readme(filReadme);
97  disp(rdm');  disp(rdm');
98    
99  %convert to MITgcm format (90x1170 array)  %set output directory/file name
100  myFile=[dirOut nameDiag];%first instance is for subdirectory name  myFile=[dirOut nameDiag];%first instance is for subdirectory name
101  if ~isdir(myFile); mkdir(myFile); end;  if ~isdir(myFile); mkdir(myFile); end;
102  myFile=[myFile '/' nameDiag];%second instance is for file name base  myFile=[myFile filesep nameDiag];%second instance is for file name base
103    
104  %get grid params  %get grid params
105  [grid_diag]=set_grid_diag(avail_diag);  [grid_diag]=set_grid_diag(avail_diag);
# 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.12
  ViewVC Help
Powered by ViewVC 1.1.22