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

Diff of /MITgcm/verification/lab_sea/matlab/lookat_exp6.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/exp1/';
5    exp2_path='../../../exe/exp6/';
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 'FUtave.0000000010.data'],[20 16 1],1);
17    vice1=readbin([exp1_path 'FVtave.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 'FUtave.0000000010.data'],[20 16 1],1);
21    vice2=readbin([exp2_path 'FVtave.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
37    subplot(321)
38    mypcolor(lon,lat,area1'); caxis([0 1]), colorbar
39    title('Sea-ice concentration with dynamics')
40    subplot(322)
41    mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
42    title('Without dynamics')
43    subplot(323)
44    mypcolor(lon,lat,heff1'); caxis([0 .3]), colorbar
45    title('Effective sea-ice thickness with dynamics (m)')
46    subplot(324)
47    mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
48    title('Without dynamics')
49    subplot(325)
50    mypcolor(Blon,Blat,icespeed1'); caxis([0 10]), colorbar
51    hold on,myquiver(Blon,Blat,udir1',vdir1','k');
52    title('Wind stress with dynamics (N/m^2)')
53    subplot(326)
54    mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar
55    hold on,myquiver(Blon,Blat,udir2',vdir2','k');
56    title('Without dynamics')
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22