1 |
dimitri |
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='../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') |