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

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

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

Revision 1.2 - (hide 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: +40 -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 heimbach 1.2 % 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='../../../exe/exp1/';
5     exp2_path='../../../exe/exp4/';
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     % plot comparison figures
24     clf, orient tall, wysiwyg
25     subplot(421), mypcolor(lon,lat,area2'); colorbar
26     title('Ice concentration during first 10 hours')
27     subplot(422), mypcolor(lon,lat,area2'-area1'); colorbar
28     title('2cpu - 1cpu solution')
29     subplot(423), mypcolor(lon,lat,heff2'); colorbar
30     title('Effective ice thickness')
31     subplot(424), mypcolor(lon,lat,heff2'-heff1'); colorbar
32     title('2cpu - 1cpu solution')
33     subplot(425), mypcolor(lon,lat,uice2'); colorbar
34     title('Uice')
35     subplot(426), mypcolor(lon,lat,uice2'-uice1'); colorbar
36     title('2cpu - 1cpu solution')
37     subplot(427), mypcolor(lon,lat,area2'); colorbar
38     title('Vice')
39     subplot(428), mypcolor(lon,lat,vice2'-vice1'); colorbar
40     title('2cpu - 1cpu solution')
  ViewVC Help
Powered by ViewVC 1.1.22