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

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

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

revision 1.1 by dimitri, Thu Dec 5 08:06:18 2002 UTC revision 1.2 by dimitri, Thu Dec 5 08:43:03 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='../input/exp6/';
5    exp2_path='../input/exp7/';
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 'uVeltave.0000000010.data'],[20 16 23],1);
17    vice1=readbin([exp1_path 'vVeltave.0000000010.data'],[20 16 23],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 'uVeltave.0000000010.data'],[20 16 23],1);
21    vice2=readbin([exp2_path 'vVeltave.0000000010.data'],[20 16 23],1);
22    eta1=readbin([exp1_path 'ETAtave.0000000010.data'],[20 16 1],1);
23    eta2=readbin([exp2_path 'ETAtave.0000000010.data'],[20 16 1],1);
24    t1=readbin([exp1_path 'Ttave.0000000010.data'],[20 16 23],1);
25    t2=readbin([exp2_path 'Ttave.0000000010.data'],[20 16 23],1);
26    s1=readbin([exp1_path 'Stave.0000000010.data'],[20 16 23],1);
27    s2=readbin([exp2_path 'Stave.0000000010.data'],[20 16 23],1);
28    
29    % plot comparison figures
30    clf, orient tall, wysiwyg, k=1;
31    subplot(721), mypcolor(lon,lat,area2'); colorbar
32    title('Ice concentration during first 10 hours'),set(gca,'xticklabel',[])
33    subplot(722), mypcolor(lon,lat,area2'-area1'); colorbar
34    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
35    subplot(723), mypcolor(lon,lat,heff2'); colorbar
36    title('Effective ice thickness'),set(gca,'xticklabel',[])
37    subplot(724), mypcolor(lon,lat,heff2'-heff1'); colorbar
38    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
39    subplot(725), mypcolor(lon,lat,uice2(:,:,k)'); colorbar
40    title('uVel at level 1'),set(gca,'xticklabel',[])
41    subplot(726), mypcolor(lon,lat,uice2(:,:,k)'-uice1(:,:,k)'); colorbar
42    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
43    subplot(727), mypcolor(lon,lat,vice2(:,:,k)'); colorbar
44    title('vVel at level 1'),set(gca,'xticklabel',[])
45    subplot(728), mypcolor(lon,lat,vice2(:,:,k)'-vice1(:,:,k)'); colorbar
46    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
47    subplot(7,2,9), mypcolor(lon,lat,eta2'); colorbar
48    title('SSH'),set(gca,'xticklabel',[])
49    subplot(7,2,10), mypcolor(lon,lat,eta2'-eta1'); colorbar
50    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
51    subplot(7,2,11), mypcolor(lon,lat,t2(:,:,k)'); colorbar
52    title('theta at level 1'),set(gca,'xticklabel',[])
53    subplot(7,2,12), mypcolor(lon,lat,t2(:,:,k)'-t1(:,:,k)'); colorbar
54    title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
55    subplot(7,2,13), mypcolor(lon,lat,s2(:,:,k)'); colorbar
56    title('salt at level 1')
57    subplot(7,2,14), mypcolor(lon,lat,s2(:,:,k)'-s1(:,:,k)'); colorbar
58    title('2cpu - 1cpu solution')
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22