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

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

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

revision 1.1 by dimitri, Fri Dec 27 08:01:56 2002 UTC revision 1.1.2.2 by dimitri, Thu May 22 20:49:15 2003 UTC
# Line 0  Line 1 
1    % need to be in verification/lab_sea/matlab directory
2    % and to specify location of exp1 output
3    cd ../../../verification/lab_sea/matlab
4    exp8_path='../input/exp8/';
5    exp9_path='../input/exp9/';
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([exp8_path 'AREAtave.0000000010.data'],[20 16 1],1);
15    heff1=readbin([exp8_path 'HEFFtave.0000000010.data'],[20 16 1],1);
16    uice1=readbin([exp8_path 'UICEtave.0000000010.data'],[20 16 1],1);
17    vice1=readbin([exp8_path 'VICEtave.0000000010.data'],[20 16 1],1);
18    area2=readbin([exp9_path 'AREAtave.0000000010.data'],[20 16 1],1);
19    heff2=readbin([exp9_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20    uice2=readbin([exp9_path 'UICEtave.0000000010.data'],[20 16 1],1);
21    vice2=readbin([exp9_path 'VICEtave.0000000010.data'],[20 16 1],1);
22    
23    % compute ice speed and direction
24    udir1=0*uice1; vdir1=0*vice1;
25    udir2=0*uice2; vdir2=0*vice2;
26    icespeed1=100*sqrt(uice1.^2+vice1.^2);
27    in=find(icespeed1);
28    udir1(in)=uice1(in)./icespeed1(in);
29    vdir1(in)=vice1(in)./icespeed1(in);
30    icespeed2=100*sqrt(uice2.^2+vice2.^2);
31    in=find(icespeed2);
32    udir2(in)=uice2(in)./icespeed2(in);
33    vdir2(in)=vice2(in)./icespeed2(in);
34    
35    % plot comparison figures
36    clf, subplot(321)
37    mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
38    title('Sea-ice concentration with pkg/exf fluxes')
39    set(gca,'xticklabel',[])
40    
41    subplot(322)
42    mypcolor(lon,lat,area2'-area1'); colorbar
43    title('Difference with Experiment 8')
44    set(gca,'xticklabel',[])
45    
46    subplot(323)
47    mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
48    title('Effective sea-ice thickness (m)')
49    set(gca,'xticklabel',[])
50    
51    subplot(324)
52    mypcolor(lon,lat,heff2'-heff1'); colorbar
53    title('Difference with Experiment 8')
54    set(gca,'xticklabel',[])
55    
56    subplot(325)
57    mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar
58    hold on, myquiver(Blon,Blat,udir2',vdir2','k');
59    title('Sea-ice velocity (cm/s)')
60    
61    subplot(326)
62    mypcolor(Blon,Blat,icespeed2'-icespeed1'); colorbar
63    if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))
64      hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');
65    end
66    title('Difference with Experiment 8')
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.2.2
  ViewVC Help
Powered by ViewVC 1.1.22