/[MITgcm]/MITgcm/verification/lab_sea/matlab/lookat_exp4.m
ViewVC logotype

Diff of /MITgcm/verification/lab_sea/matlab/lookat_exp4.m

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

revision 1.1 by heimbach, Mon Nov 11 22:03:48 2002 UTC revision 1.1.2.1 by heimbach, Mon Nov 11 22:03:48 2002 UTC
# Line 0  Line 1 
1    % need to be in verification/lab_sea/matlab directory
2    % and to specify location of exp1 and exp2 output
3    cd ../../../verification/lab_sea/matlab
4    exp1_path='../../../exe/exp1/';
5    exp2_path='../../../exe/exp4/';
6    
7    % load monthly-mean SMMR-SSM/I data
8    load SSMI
9    
10    % B-grid latitude for uice and vice
11    Blat=lat+1; Blon=lon+1;
12    
13    % load model output
14    area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
15    heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
16    uice1=readbin([exp1_path 'UICEtave.0000000010.data'],[20 16 1],1);
17    vice1=readbin([exp1_path 'VICEtave.0000000010.data'],[20 16 1],1);
18    area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
19    heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20    uice2=readbin([exp2_path 'UICEtave.0000000010.data'],[20 16 1],1);
21    vice2=readbin([exp2_path 'VICEtave.0000000010.data'],[20 16 1],1);
22    
23    % plot comparison figures
24    clf, orient tall, wysiwyg
25    subplot(421), mypcolor(lon,lat,area2'); colorbar
26    title('Ice concentration during first 10 hours')
27    subplot(422), mypcolor(lon,lat,area2'-area1'); colorbar
28    title('2cpu - 1cpu solution')
29    subplot(423), mypcolor(lon,lat,heff2'); colorbar
30    title('Effective ice thickness')
31    subplot(424), mypcolor(lon,lat,heff2'-heff1'); colorbar
32    title('2cpu - 1cpu solution')
33    subplot(425), mypcolor(lon,lat,uice2'); colorbar
34    title('Uice')
35    subplot(426), mypcolor(lon,lat,uice2'-uice1'); colorbar
36    title('2cpu - 1cpu solution')
37    subplot(427), mypcolor(lon,lat,area2'); colorbar
38    title('Vice')
39    subplot(428), mypcolor(lon,lat,vice2'-vice1'); colorbar
40    title('2cpu - 1cpu solution')
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.2.1
  ViewVC Help
Powered by ViewVC 1.1.22