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

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

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


Revision 1.2 - (show annotations) (download)
Tue Mar 16 04:59:01 2004 UTC (17 years, 10 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint52n_post, checkpoint53d_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint58n_post, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint58h_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint58q_post, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint58j_post, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint53g_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint58k_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.1: +12 -14 lines
update verification/lab_sea

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 exp_path='../input/';
5
6 % load monthly-mean SMMR-SSM/I data
7 load SSMI
8
9 % B-grid latitude for uice and vice
10 Blat=lat-1; Blon=lon-1;
11
12 % load model output
13 area =readbin('../results/AREAtave.0000000010.data' ,[20 16 1],1);
14 heff =readbin('../results/HEFFtave.0000000010.data' ,[20 16 1],1);
15 uice =readbin('../results/UICEtave.0000000010.data' ,[20 16 1],1);
16 vice =readbin('../results/VICEtave.0000000010.data' ,[20 16 1],1);
17 area1=readbin([exp_path 'AREAtave.0000000010.data'],[20 16 1],1);
18 heff1=readbin([exp_path 'HEFFtave.0000000010.data'],[20 16 1],1);
19 uice1=readbin([exp_path 'UICEtave.0000000010.data'],[20 16 1],1);
20 vice1=readbin([exp_path 'VICEtave.0000000010.data'],[20 16 1],1);
21
22 % compute ice speed and direction
23 udir=0*uice; vdir=0*vice;
24 udir1=0*uice; vdir1=0*vice;
25 icespeed=100*sqrt(uice.^2+vice.^2);
26 in=find(icespeed);
27 udir(in)=uice(in)./icespeed(in);
28 vdir(in)=vice(in)./icespeed(in);
29 icespeed1=100*sqrt(uice1.^2+vice1.^2);
30 in=find(icespeed1);
31 udir1(in)=uice1(in)./icespeed1(in);
32 vdir1(in)=vice1(in)./icespeed1(in);
33
34 % plot comparison figures
35
36 clf, subplot(321)
37 mypcolor(lon,lat,area1'); caxis([0 1]), colorbar
38 title('Sea-ice concentration, this run')
39 set(gca,'xticklabel',[])
40
41 subplot(322)
42 mypcolor(lon,lat,area'); caxis([0 1]), colorbar
43 title('Sea-ice concentration, checkpoint52l')
44 set(gca,'xticklabel',[])
45
46 subplot(323)
47 mypcolor(lon,lat,heff1'); caxis([0 .6]), colorbar
48 title('Effective sea-ice thickness (m), this run')
49 set(gca,'xticklabel',[])
50
51 subplot(324)
52 mypcolor(lon,lat,heff'); caxis([0 .6]), colorbar
53 title('Effective sea-ice thickness (m), checkpoint52l')
54 set(gca,'xticklabel',[])
55
56 subplot(325)
57 mypcolor(Blon,Blat,icespeed1'); caxis([0 12]), colorbar
58 hold on, myquiver(Blon,Blat,udir1',vdir1','k');
59 title('Sea-ice velocity (cm/s), this run')
60
61 subplot(326)
62 mypcolor(Blon,Blat,icespeed'); caxis([0 12]), colorbar
63 hold on, myquiver(Blon,Blat,udir',vdir','k');
64 title('Sea-ice velocity (cm/s), checkpoint52l')

  ViewVC Help
Powered by ViewVC 1.1.22