Warning, /verification/lab_sea/matlab/lookat_ice.m is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 83079591 on 2006-11-06 20:37:26 UTC
e7c33124ce Dimi*0001 % need to be in verification/lab_sea/matlab directory
0002 % and to specify location of exp1 output
0003 cd ../../../verification/lab_sea/matlab
8307959141 Dimi*0004 exp_path='../build/';
e7c33124ce Dimi*0005
0006 % load monthly-mean SMMR-SSM/I data
0007 load SSMI
0008
0009 % B-grid latitude for uice and vice
0010 Blat=lat-1; Blon=lon-1;
0011
0012 % load model output
0013 area =readbin('../results/AREAtave.0000000010.data' ,[20 16 1],1);
0014 heff =readbin('../results/HEFFtave.0000000010.data' ,[20 16 1],1);
0015 uice =readbin('../results/UICEtave.0000000010.data' ,[20 16 1],1);
0016 vice =readbin('../results/VICEtave.0000000010.data' ,[20 16 1],1);
0017 area1=readbin([exp_path 'AREAtave.0000000010.data'],[20 16 1],1);
0018 heff1=readbin([exp_path 'HEFFtave.0000000010.data'],[20 16 1],1);
0019 uice1=readbin([exp_path 'UICEtave.0000000010.data'],[20 16 1],1);
0020 vice1=readbin([exp_path 'VICEtave.0000000010.data'],[20 16 1],1);
0021
0022 % compute ice speed and direction
0023 udir=0*uice; vdir=0*vice;
0024 udir1=0*uice; vdir1=0*vice;
0025 icespeed=100*sqrt(uice.^2+vice.^2);
0026 in=find(icespeed);
0027 udir(in)=uice(in)./icespeed(in);
0028 vdir(in)=vice(in)./icespeed(in);
0029 icespeed1=100*sqrt(uice1.^2+vice1.^2);
0030 in=find(icespeed1);
0031 udir1(in)=uice1(in)./icespeed1(in);
0032 vdir1(in)=vice1(in)./icespeed1(in);
0033
0034 % plot comparison figures
0035
0036 clf, subplot(321)
0037 mypcolor(lon,lat,area1'); caxis([0 1]), colorbar
0e4895cdf4 Dimi*0038 title('Sea-ice concentration, this run')
e7c33124ce Dimi*0039 set(gca,'xticklabel',[])
0040
0041 subplot(322)
0e4895cdf4 Dimi*0042 mypcolor(lon,lat,area'); caxis([0 1]), colorbar
0043 title('Sea-ice concentration, checkpoint52l')
e7c33124ce Dimi*0044 set(gca,'xticklabel',[])
0045
0046 subplot(323)
0e4895cdf4 Dimi*0047 mypcolor(lon,lat,heff1'); caxis([0 .6]), colorbar
0048 title('Effective sea-ice thickness (m), this run')
e7c33124ce Dimi*0049 set(gca,'xticklabel',[])
0050
0051 subplot(324)
0e4895cdf4 Dimi*0052 mypcolor(lon,lat,heff'); caxis([0 .6]), colorbar
0053 title('Effective sea-ice thickness (m), checkpoint52l')
e7c33124ce Dimi*0054 set(gca,'xticklabel',[])
0055
0056 subplot(325)
0e4895cdf4 Dimi*0057 mypcolor(Blon,Blat,icespeed1'); caxis([0 12]), colorbar
e7c33124ce Dimi*0058 hold on, myquiver(Blon,Blat,udir1',vdir1','k');
0e4895cdf4 Dimi*0059 title('Sea-ice velocity (cm/s), this run')
e7c33124ce Dimi*0060
0061 subplot(326)
0e4895cdf4 Dimi*0062 mypcolor(Blon,Blat,icespeed'); caxis([0 12]), colorbar
0063 hold on, myquiver(Blon,Blat,udir',vdir','k');
0064 title('Sea-ice velocity (cm/s), checkpoint52l')