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

Diff of /MITgcm/verification/lab_sea/matlab/lookat_exp5.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.2 by heimbach, Tue Nov 12 20:54:29 2002 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    exp1_path='../../../exe/';
5    
6    % load monthly-mean SMMR-SSM/I data
7    load SSMI
8    
9    % B-grid latitude for uice and vice
10    Blat=lat+1; Blon=lon+1;
11    
12    % load model output
13    area=readbin([exp1_path 'AREAtave.0000010950.data'],[20 16 1],1);
14    heff=readbin([exp1_path 'HEFFtave.0000010950.data'],[20 16 1],1);
15    uice=readbin([exp1_path 'UICEtave.0000010950.data'],[20 16 1],1);
16    vice=readbin([exp1_path 'VICEtave.0000010950.data'],[20 16 1],1);
17    
18    % compute ice speed and direction
19    udir=0*uice; vdir=0*vice;
20    icespeed=100*sqrt(uice.^2+vice.^2);
21    in=find(icespeed);
22    udir(in)=uice(in)./icespeed(in);
23    vdir(in)=vice(in)./icespeed(in);
24    
25    % plot comparison figures
26    clf
27    subplot(221)
28    mypcolor(lon,lat,SSMI(:,:,4)'/100); caxis([0 1]), colorbar
29    title('Observed March sea-ice concentration')
30    subplot(222), mypcolor(lon,lat,area'); caxis([0 1]), colorbar
31    title('Modeled March sea-ice concentration')
32    subplot(223), mypcolor(lon,lat,heff'); caxis([0 1]), colorbar
33    title('Effective thickness in m')
34    subplot(224)
35    mypcolor(Blon,Blat,icespeed'); caxis([0 10]), colorbar
36    hold on,myquiver(Blon,Blat,udir',vdir','k');
37    title('Sea-ice velocity (cm/s)')
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22