|
1 %====================================================================== |
|
2 % C H E C K I N V . M |
|
3 % doc: Tue Jun 13 14:29:15 2017 |
|
4 % dlm: Tue Jun 13 14:33:27 2017 |
|
5 % (c) 2017 A.M. Thurnherr |
|
6 % uE-Info: 10 77 NIL 0 0 72 0 2 4 NIL ofnI |
|
7 %====================================================================== |
|
8 |
|
9 % CHANGES BY ANT: |
|
10 % Jun 13, 2017: Matlab incompatibility (legend had extraneous -1 final arg) |
|
11 |
1 function p=checkinv(dr,di,de,der,p,ps) |
12 function p=checkinv(dr,di,de,der,p,ps) |
2 % function p=checkinv(dr,di,de,der,p,ps) |
13 % function p=checkinv(dr,di,de,der,p,ps) |
3 % check inversion for consistency |
14 % check inversion for consistency |
4 % Martin Visbeck 2004 |
15 % Martin Visbeck 2004 |
5 |
16 |
82 ic3=[ic3(1:2:end),fliplr(ic3(2:2:end))]; |
93 ic3=[ic3(1:2:end),fliplr(ic3(2:2:end))]; |
83 colormap(col(ic3,:)); |
94 colormap(col(ic3,:)); |
84 |
95 |
85 subplot(211) |
96 subplot(211) |
86 bar(dr.z,de.ocean_constraints(ic,:)','stack') |
97 bar(dr.z,de.ocean_constraints(ic,:)','stack') |
87 legend(de.type_constraints(ic,:),-1) |
98 % legend(de.type_constraints(ic,:),-1) |
|
99 legend(de.type_constraints(ic,:)) |
88 ylabel('sum of weights') |
100 ylabel('sum of weights') |
89 title('ocean velocity constraints') |
101 title('ocean velocity constraints') |
90 xlabel('depth [m]') |
102 xlabel('depth [m]') |
91 axis tight |
103 axis tight |
92 |
104 |
93 subplot(212) |
105 subplot(212) |
94 bar(de.ctd_constraints(ic,:)','stack') |
106 bar(de.ctd_constraints(ic,:)','stack') |
95 legend(de.type_constraints(ic,:),-1) |
107 % legend(de.type_constraints(ic,:),-1) |
|
108 legend(de.type_constraints(ic,:)) |
96 title('CTD velocity constraints') |
109 title('CTD velocity constraints') |
97 ylabel('sum of weights') |
110 ylabel('sum of weights') |
98 xlabel('super ensemble') |
111 xlabel('super ensemble') |
99 axis tight |
112 axis tight |
100 |
113 |