/[MITgcm]/MITgcm_contrib/high_res_cube/README_ice
ViewVC logotype

Diff of /MITgcm_contrib/high_res_cube/README_ice

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

revision 1.74 by dimitri, Thu Nov 2 06:05:17 2006 UTC revision 1.78 by dimitri, Tue Dec 5 20:46:48 2006 UTC
# Line 3  Instructions for setting up a cube-spher Line 3  Instructions for setting up a cube-spher
3  =================  =================
4  UPDATES  UPDATES
5    
6  November 1, 2006: 216-CPU code for cube43  December 5, 2006: 216-CPU code for cube46
7  Martin's growth.F routine and Patrick's pkg/seaice diagnostics  Changed diagnostics to accomodate JMCs latest and some
8    changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h
9    to speed up code.
10    
11    November 7, 2006: 216-CPU code for cube43
12    Bug fixes in growth.F routine and new set of diagnostics
13    
14  October 15, 2006: 216-CPU code for cube40  October 15, 2006: 216-CPU code for cube40
15    
# Line 109  comparison output is in: Line 114  comparison output is in:
114  to use matlab to look at the output  to use matlab to look at the output
115  matlab  matlab
116    
117  for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
118           'QNETtave','FVtave','FUtave','EmPmRtave'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
119      tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);      tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);
120      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
121  end  end
122    
123  heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);  heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
124  mask=0*heff; mask(find(heff>.001))=1;  mask=0*heff; mask(find(heff>.001))=1;
125  uice=mask.*permute(readbin(['UICEtave.0000000020.data'],[32,6,32]),[1 3 2]);  uice=mask.*permute(readbin(['SIuice.0000000072.data'],[32,6,32]),[1 3 2]);
126  vice=mask.*permute(readbin(['VICEtave.0000000020.data'],[32,6,32]),[1 3 2]);  vice=mask.*permute(readbin(['SIvice.0000000072.data'],[32,6,32]),[1 3 2]);
127  clf  clf
128  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
129  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')  subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')
Legend:
Removed from v.1.74  
changed lines
  Added in v.1.78
  ViewVC Help
Powered by ViewVC 1.1.22