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

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

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

Revision 1.2 - (show annotations) (download)
Tue Nov 12 20:54:29 2002 UTC (21 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47a_post, checkpoint47
Changes since 1.1: +58 -0 lines 
Merging from release1_p8:
o New verification/lab_sea exp. for seaice package
  (cf. README file).
   N.B.: fairly large differences in accuracy occur across
        different platforms/compilers; comparison between
        g77 (Linux) and F77 (SGI Origin 2000/MIPSpro)
        generated output gives:
   T           S           U           V
  C D M    c        m  s        m  s        m  s        m  s
  n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
  f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
  g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

  Y Y Y Y  5  5  7  7  7  8 10  9  6  6  6  6  7  5  7  5  7 FAIL  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 exp1_path='../../../exe/exp1/';
5 exp2_path='../../../exe/exp2/';
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([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 'UICEtave.0000000010.data'],[20 16 1],1);
21 vice2=readbin([exp2_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
37 subplot(321)
38 mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
39 title('Sea-ice concentration using LSR solver')
40 subplot(322)
41 mypcolor(lon,lat,area2'-area1'); colorbar
42 title('Difference with ADI solver')
43 subplot(323)
44 mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
45 title('Effective sea-ice thickness (m)')
46 subplot(324)
47 mypcolor(lon,lat,heff2'-heff1'); colorbar
48 title('Difference with ADI solver')
49 subplot(325)
50 mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar
51 hold on,myquiver(Blon,Blat,udir2',vdir2','k');
52 title('Sea-ice velocity (cm/s)')
53 subplot(326)
54 mypcolor(Blon,Blat,icespeed2'-icespeed1'); colorbar
55 if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))
56 hold on,myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');
57 end
58 title('Difference with ADI solver')
  ViewVC Help
Powered by ViewVC 1.1.22