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

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

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

Revision 1.6 - (show annotations) (download)
Thu Aug 7 02:31:30 2003 UTC (20 years, 9 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51j_post, checkpoint52e_pre, checkpoint51n_pre, checkpoint52d_pre, branchpoint-genmake2, branch-netcdf, checkpoint51r_post, checkpoint52b_pre, checkpoint51o_pre, checkpoint51i_post, checkpoint51l_pre, checkpoint51l_post, checkpoint51o_post, checkpoint51f_pre, checkpoint51q_post, checkpoint52, checkpoint52d_post, checkpoint52a_post, checkpoint52b_post, checkpoint52f_post, checkpoint52c_post, checkpoint51h_pre, checkpoint51g_post, ecco_c52_e35, checkpoint51f_post, checkpoint52e_post, checkpoint52a_pre, checkpoint51m_post, checkpoint51t_post, checkpoint52i_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.5: +3 -3 lines 
o Added on-the-fly spatial interpolation capability
    "USE_EXF_INTERPOLATION" to pkg/exf.
  - This is a temporary Cartesian-grid hack until
    the super-duper ESMF coupler becomes available.
  - See verification/global_with_exf/README for usage example.
  - Removed obsolete EXFwindOnBgrid and SEAICEwindOnCgrid
    flags and modified pkg/seaice accordingly.
o Bug fix to pkg/ptracers, pkg/generic_advdiff/gad_calc_rhs.F,
    and pkg/kpp/kpp_transport_ptr.F for dealing with tracer
    non-local transport term.
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
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([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
18 heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
19 uice1=readbin([exp1_path 'UICEtave.0000000010.data'],[20 16 1],1);
20 vice1=readbin([exp1_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')
39 set(gca,'xticklabel',[])
40
41 subplot(322)
42 mypcolor(lon,lat,area1'-area'); colorbar
43 title('Difference with checkpoint51f')
44 set(gca,'xticklabel',[])
45
46 subplot(323)
47 mypcolor(lon,lat,heff1'); caxis([0 .3]), colorbar
48 title('Effective sea-ice thickness (m)')
49 set(gca,'xticklabel',[])
50
51 subplot(324)
52 mypcolor(lon,lat,heff1'-heff'); colorbar
53 title('Difference with checkpoint51f')
54 set(gca,'xticklabel',[])
55
56 subplot(325)
57 mypcolor(Blon,Blat,icespeed1'); caxis([0 10]), colorbar
58 hold on, myquiver(Blon,Blat,udir1',vdir1','k');
59 title('Sea-ice velocity (cm/s)')
60
61 subplot(326)
62 mypcolor(Blon,Blat,icespeed1'-icespeed'); colorbar
63 if mmax(abs(udir1-udir)) | mmax(abs(vdir1-vdir))
64 hold on, myquiver(Blon,Blat,udir1'-udir',vdir1'-vdir','k');
65 end
66 title('Difference with checkpoint51f')
  ViewVC Help
Powered by ViewVC 1.1.22