/[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.3 by dimitri, Mon Nov 6 20:37:26 2006 UTC
# Line 1  Line 1 
1  % need to be in verification/lab_sea/matlab directory  % need to be in verification/lab_sea/matlab directory
2  % and to specify location of exp1 output  % and to specify location of exp1 output
3  cd ../../../verification/lab_sea/matlab  cd ../../../verification/lab_sea/matlab
4  exp_path='../input/';  exp_path='../build/';
5    
6  % load monthly-mean SMMR-SSM/I data  % load monthly-mean SMMR-SSM/I data
7  load SSMI  load SSMI
# 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.3
  ViewVC Help
Powered by ViewVC 1.1.22