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

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

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

Revision 1.3 - (show annotations) (download)
Fri May 23 20:19:16 2003 UTC (21 years ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51j_post, checkpoint52e_pre, checkpoint51n_pre, checkpoint52d_pre, checkpoint50g_post, branchpoint-genmake2, branch-netcdf, checkpoint51r_post, checkpoint52b_pre, checkpoint51o_pre, checkpoint51i_post, checkpoint51e_post, checkpoint51b_post, checkpoint51l_pre, checkpoint51l_post, checkpoint51o_post, checkpoint51f_pre, checkpoint51q_post, checkpoint51, checkpoint52, checkpoint52d_post, checkpoint51b_pre, checkpoint52a_post, checkpoint52b_post, checkpoint52f_post, checkpoint52c_post, checkpoint51h_pre, checkpoint51g_post, ecco_c52_e35, checkpoint51f_post, checkpoint52e_post, checkpoint51c_post, checkpoint52a_pre, checkpoint51d_post, checkpoint51m_post, checkpoint51t_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51i_pre, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint51s_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch
Changes since 1.2: +1 -1 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 and exp2 output
3 cd ../../../verification/lab_sea/matlab
4 exp1_path='../input/exp6/';
5 exp2_path='../input/exp7/';
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 'uVeltave.0000000010.data'],[20 16 23],1);
17 vice1=readbin([exp1_path 'vVeltave.0000000010.data'],[20 16 23],1);
18 area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
19 heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20 uice2=readbin([exp2_path 'uVeltave.0000000010.data'],[20 16 23],1);
21 vice2=readbin([exp2_path 'vVeltave.0000000010.data'],[20 16 23],1);
22 eta1=readbin([exp1_path 'ETAtave.0000000010.data'],[20 16 1],1);
23 eta2=readbin([exp2_path 'ETAtave.0000000010.data'],[20 16 1],1);
24 t1=readbin([exp1_path 'Ttave.0000000010.data'],[20 16 23],1);
25 t2=readbin([exp2_path 'Ttave.0000000010.data'],[20 16 23],1);
26 s1=readbin([exp1_path 'Stave.0000000010.data'],[20 16 23],1);
27 s2=readbin([exp2_path 'Stave.0000000010.data'],[20 16 23],1);
28
29 % plot comparison figures
30 clf, orient tall, wysiwyg, k=1;
31 subplot(721), mypcolor(lon,lat,area2'); colorbar
32 title('Ice concentration during first 10 hours'),set(gca,'xticklabel',[])
33 subplot(722), mypcolor(lon,lat,area2'-area1'); colorbar
34 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
35 subplot(723), mypcolor(lon,lat,heff2'); colorbar
36 title('Effective ice thickness'),set(gca,'xticklabel',[])
37 subplot(724), mypcolor(lon,lat,heff2'-heff1'); colorbar
38 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
39 subplot(725), mypcolor(lon,lat,uice2(:,:,k)'); colorbar
40 title('uVel at level 1'),set(gca,'xticklabel',[])
41 subplot(726), mypcolor(lon,lat,uice2(:,:,k)'-uice1(:,:,k)'); colorbar
42 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
43 subplot(727), mypcolor(lon,lat,vice2(:,:,k)'); colorbar
44 title('vVel at level 1'),set(gca,'xticklabel',[])
45 subplot(728), mypcolor(lon,lat,vice2(:,:,k)'-vice1(:,:,k)'); colorbar
46 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
47 subplot(7,2,9), mypcolor(lon,lat,eta2'); colorbar
48 title('SSH'),set(gca,'xticklabel',[])
49 subplot(7,2,10), mypcolor(lon,lat,eta2'-eta1'); colorbar
50 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
51 subplot(7,2,11), mypcolor(lon,lat,t2(:,:,k)'); colorbar
52 title('theta at level 1'),set(gca,'xticklabel',[])
53 subplot(7,2,12), mypcolor(lon,lat,t2(:,:,k)'-t1(:,:,k)'); colorbar
54 title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
55 subplot(7,2,13), mypcolor(lon,lat,s2(:,:,k)'); colorbar
56 title('salt at level 1')
57 subplot(7,2,14), mypcolor(lon,lat,s2(:,:,k)'-s1(:,:,k)'); colorbar
58 title('2cpu - 1cpu solution')
  ViewVC Help
Powered by ViewVC 1.1.22