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

Annotation of /MITgcm/verification/lab_sea/matlab/lookat_exp6.m

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

Revision 1.4 - (hide annotations) (download)
Fri May 23 20:19:16 2003 UTC (21 years, 1 month 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.3: +17 -38 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 heimbach 1.2 % need to be in verification/lab_sea/matlab directory
2     % and to specify location of exp1 output
3     cd ../../../verification/lab_sea/matlab
4 dimitri 1.3 exp1_path='../input/exp1/';
5     exp2_path='../input/exp6/';
6 heimbach 1.2
7     % load monthly-mean SMMR-SSM/I data
8     load SSMI
9    
10     % load model output
11     area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
12     heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
13     area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
14     heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
15    
16     % plot comparison figures
17 dimitri 1.4 clf, subplot(221)
18 heimbach 1.2 mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
19 dimitri 1.4 title('Sea-ice concentration, no dynamics')
20     set(gca,'xticklabel',[])
21    
22     subplot(222)
23     mypcolor(lon,lat,area2'-area1'); colorbar
24     title('Difference with LSR solver')
25     set(gca,'xticklabel',[])
26    
27     subplot(223)
28 heimbach 1.2 mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
29 dimitri 1.4 title('Effective sea-ice thickness (m)')
30     set(gca,'xticklabel',[])
31    
32     subplot(224)
33     mypcolor(lon,lat,heff2'-heff1'); colorbar
34     title('Difference with LSR solver')
35     set(gca,'xticklabel',[])
  ViewVC Help
Powered by ViewVC 1.1.22