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

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

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

revision 1.3 by dimitri, Thu Dec 5 08:43:03 2002 UTC revision 1.4 by dimitri, Fri May 23 20:19:16 2003 UTC
# Line 8  exp2_path='../input/exp2/'; Line 8  exp2_path='../input/exp2/';
8  load SSMI  load SSMI
9    
10  % B-grid latitude for uice and vice  % B-grid latitude for uice and vice
11  Blat=lat+1; Blon=lon+1;  Blat=lat-1; Blon=lon-1;
12    
13  % load model output  % load model output
14  area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);  area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
# Line 35  vdir2(in)=vice2(in)./icespeed2(in); Line 35  vdir2(in)=vice2(in)./icespeed2(in);
35  % plot comparison figures  % plot comparison figures
36  clf, subplot(321)  clf, subplot(321)
37  mypcolor(lon,lat,area2'); caxis([0 1]), colorbar  mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
38  title('Sea-ice concentration using LSR solver')  title('Sea-ice concentration using ADI solver')
39  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
40    
41  subplot(322)  subplot(322)
42  mypcolor(lon,lat,area2'-area1'); colorbar  mypcolor(lon,lat,area2'-area1'); colorbar
43  title('Difference with ADI solver')  title('Difference with LSR solver')
44  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
45    
46  subplot(323)  subplot(323)
# Line 50  set(gca,'xticklabel',[]) Line 50  set(gca,'xticklabel',[])
50    
51  subplot(324)  subplot(324)
52  mypcolor(lon,lat,heff2'-heff1'); colorbar  mypcolor(lon,lat,heff2'-heff1'); colorbar
53  title('Difference with ADI solver')  title('Difference with LSR solver')
54  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
55    
56  subplot(325)  subplot(325)
# Line 63  mypcolor(Blon,Blat,icespeed2'-icespeed1' Line 63  mypcolor(Blon,Blat,icespeed2'-icespeed1'
63  if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))  if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))
64    hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');    hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');
65  end  end
66  title('Difference with ADI solver')  title('Difference with LSR solver')
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22