/[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.77 by dimitri, Wed Nov 22 07:21:31 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  November 22, 2006: 216-CPU code for cube44
7  Martin's growth.F routine and Patrick's pkg/seaice diagnostics  Changed diagnostics, yet again.
8    
9    November 7, 2006: 216-CPU code for cube43
10    Bug fixes in growth.F routine and new set of diagnostics
11    
12  October 15, 2006: 216-CPU code for cube40  October 15, 2006: 216-CPU code for cube40
13    
# Line 109  comparison output is in: Line 112  comparison output is in:
112  to use matlab to look at the output  to use matlab to look at the output
113  matlab  matlab
114    
115  for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...  for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
116           'QNETtave','FVtave','FUtave','EmPmRtave'}           'surForcT','oceTAUX','oceTAUY','surForcS'}
117      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]);
118      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause      cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
119  end  end
120    
121  heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);  heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
122  mask=0*heff; mask(find(heff>.001))=1;  mask=0*heff; mask(find(heff>.001))=1;
123  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]);
124  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]);
125  clf  clf
126  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')  subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
127  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.77
  ViewVC Help
Powered by ViewVC 1.1.22