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

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

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

revision 1.1 by dimitri, Sat Mar 13 14:57:54 2004 UTC revision 1.2 by dimitri, Tue Mar 16 04:59:01 2004 UTC
# Line 35  vdir1(in)=vice1(in)./icespeed1(in); Line 35  vdir1(in)=vice1(in)./icespeed1(in);
35    
36  clf, subplot(321)  clf, subplot(321)
37  mypcolor(lon,lat,area1'); caxis([0 1]), colorbar  mypcolor(lon,lat,area1'); caxis([0 1]), colorbar
38  title('Sea-ice concentration')  title('Sea-ice concentration, this run')
39  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
40    
41  subplot(322)  subplot(322)
42  mypcolor(lon,lat,area1'-area'); colorbar  mypcolor(lon,lat,area'); caxis([0 1]), colorbar
43  title('Difference with checkpoint51f')  title('Sea-ice concentration, checkpoint52l')
44  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
45    
46  subplot(323)  subplot(323)
47  mypcolor(lon,lat,heff1'); caxis([0 .3]), colorbar  mypcolor(lon,lat,heff1'); caxis([0 .6]), colorbar
48  title('Effective sea-ice thickness (m)')  title('Effective sea-ice thickness (m), this run')
49  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
50    
51  subplot(324)  subplot(324)
52  mypcolor(lon,lat,heff1'-heff'); colorbar  mypcolor(lon,lat,heff'); caxis([0 .6]), colorbar
53  title('Difference with checkpoint51f')  title('Effective sea-ice thickness (m), checkpoint52l')
54  set(gca,'xticklabel',[])  set(gca,'xticklabel',[])
55    
56  subplot(325)  subplot(325)
57  mypcolor(Blon,Blat,icespeed1'); caxis([0 10]), colorbar  mypcolor(Blon,Blat,icespeed1'); caxis([0 12]), colorbar
58  hold on, myquiver(Blon,Blat,udir1',vdir1','k');  hold on, myquiver(Blon,Blat,udir1',vdir1','k');
59  title('Sea-ice velocity (cm/s)')  title('Sea-ice velocity (cm/s), this run')
60    
61  subplot(326)  subplot(326)
62  mypcolor(Blon,Blat,icespeed1'-icespeed'); colorbar  mypcolor(Blon,Blat,icespeed'); caxis([0 12]), colorbar
63  if mmax(abs(udir1-udir)) | mmax(abs(vdir1-vdir))  hold on, myquiver(Blon,Blat,udir',vdir','k');
64    hold on, myquiver(Blon,Blat,udir1'-udir',vdir1'-vdir','k');  title('Sea-ice velocity (cm/s), checkpoint52l')
 end  
 title('Difference with checkpoint51f')  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22