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

Contents of /MITgcm/verification/lab_sea/matlab/lookat_exp11.m

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

Revision 1.2 - (show annotations) (download)
Fri May 23 20:19:16 2003 UTC (20 years, 11 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51f_post, checkpoint51l_post, checkpoint51k_post, checkpoint51o_post, branchpoint-genmake2, checkpoint51j_post, checkpoint51q_post, checkpoint52a_pre, branch-netcdf, checkpoint52e_pre, checkpoint51d_post, checkpoint51, checkpoint51r_post, checkpoint52, checkpoint52b_pre, checkpoint52d_post, checkpoint51o_pre, checkpoint51t_post, checkpoint51b_pre, checkpoint52a_post, checkpoint51i_post, checkpoint52e_post, checkpoint52b_post, checkpoint52i_post, checkpoint51a_post, checkpoint51n_pre, checkpoint52f_post, checkpoint51c_post, checkpoint52c_post, checkpoint51p_post, checkpoint51n_post, checkpoint52d_pre, checkpoint51e_post, checkpoint51b_post, checkpoint51h_pre, checkpoint51i_pre, checkpoint51l_pre, checkpoint50h_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint51g_post, checkpoint52f_pre, ecco_c52_e35, checkpoint50g_post, checkpoint51f_pre, checkpoint51m_post, checkpoint51s_post, checkpoint50i_post
Branch point for: checkpoint51n_branch, netcdf-sm0, branch-nonh, tg2-branch, branch-genmake2
Changes since 1.1: +66 -0 lines 
checkpoint50g_post
o merged with release1_p17 (pkg/seaice and verification/lab_sea)
  - added SEAICE_MULTILEVEL for 8-category sea-ice thermodynamics
  - LSR sea-ice dynamic solver moved to SouthWest B-grid location and
    made the default because of faster convergence than ADI
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='../input/exp1/';
5 exp11_path='../input/exp11/';
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 'UICEtave.0000000010.data'],[20 16 1],1);
17 vice1=readbin([exp1_path 'VICEtave.0000000010.data'],[20 16 1],1);
18 area2=readbin([exp11_path 'AREAtave.0000000010.data'],[20 16 1],1);
19 heff2=readbin([exp11_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20 uice2=readbin([exp11_path 'UICEtave.0000000010.data'],[20 16 1],1);
21 vice2=readbin([exp11_path 'VICEtave.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, subplot(321)
37 mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
38 title('Sea-ice concentration with multilevel sea-ice')
39 set(gca,'xticklabel',[])
40
41 subplot(322)
42 mypcolor(lon,lat,area2'-area1'); colorbar
43 title('Difference with Experiment 1')
44 set(gca,'xticklabel',[])
45
46 subplot(323)
47 mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
48 title('Effective sea-ice thickness (m)')
49 set(gca,'xticklabel',[])
50
51 subplot(324)
52 mypcolor(lon,lat,heff2'-heff1'); colorbar
53 title('Difference with Experiment 1')
54 set(gca,'xticklabel',[])
55
56 subplot(325)
57 mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar
58 hold on, myquiver(Blon,Blat,udir2',vdir2','k');
59 title('Sea-ice velocity (cm/s)')
60
61 subplot(326)
62 mypcolor(Blon,Blat,icespeed2'-icespeed1'); colorbar
63 if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))
64 hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');
65 end
66 title('Difference with Experiment 1')
  ViewVC Help
Powered by ViewVC 1.1.22