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); |
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},... |