1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
---|
2 | % |
---|
3 | % Make a plot from the results of the VORTEX test case |
---|
4 | % |
---|
5 | % Further Information: |
---|
6 | % http://www.brest.ird.fr/Roms_tools/ |
---|
7 | % |
---|
8 | % This file is part of ROMSTOOLS |
---|
9 | % |
---|
10 | % ROMSTOOLS is free software; you can redistribute it and/or modify |
---|
11 | % it under the terms of the GNU General Public License as published |
---|
12 | % by the Free Software Foundation; either version 2 of the License, |
---|
13 | % or (at your option) any later version. |
---|
14 | % |
---|
15 | % ROMSTOOLS is distributed in the hope that it will be useful, but |
---|
16 | % WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
17 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
18 | % GNU General Public License for more details. |
---|
19 | % |
---|
20 | % You should have received a copy of the GNU General Public License |
---|
21 | % along with this program; if not, write to the Free Software |
---|
22 | % Foundation, Inc., 59 Temple Place, Suite 330, Boston, |
---|
23 | % MA 02111-1307 USA |
---|
24 | % |
---|
25 | % Copyright (c) 2005-2006 by Pierrick Penven |
---|
26 | % e-mail:Pierrick.Penven@ird.fr |
---|
27 | % |
---|
28 | % Ref: Penven, P., L. Debreu, P. Marchesiello and J.C. McWilliams, |
---|
29 | % Application of the ROMS embedding procedure for the Central |
---|
30 | % California Upwelling System, Ocean Modelling, 2006. |
---|
31 | % |
---|
32 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
---|
33 | clear all |
---|
34 | close all |
---|
35 | % |
---|
36 | % User defined parameters |
---|
37 | % |
---|
38 | vname='temp'; |
---|
39 | tindex=11; |
---|
40 | % |
---|
41 | % Caxis depending of the variable name |
---|
42 | % |
---|
43 | if length(vname)>1 |
---|
44 | if vname(1:2)=='te' |
---|
45 | type='r'; |
---|
46 | ddd=1; |
---|
47 | cmin=17; |
---|
48 | dc=0.1; |
---|
49 | cmax=22; |
---|
50 | cff=1; |
---|
51 | elseif vname(1:2)=='ze' |
---|
52 | type='r'; |
---|
53 | ddd=0; |
---|
54 | cmin=-100; |
---|
55 | dc=5; |
---|
56 | cmax=100; |
---|
57 | cff=100; |
---|
58 | elseif vname(1:2)=='ub' |
---|
59 | type='u'; |
---|
60 | ddd=0; |
---|
61 | cmin=-50; |
---|
62 | dc=5; |
---|
63 | cmax=50; |
---|
64 | cff=100; |
---|
65 | elseif vname(1:2)=='vb' |
---|
66 | type='v'; |
---|
67 | ddd=0; |
---|
68 | cmin=-50; |
---|
69 | dc=5; |
---|
70 | cmax=50; |
---|
71 | cff=100; |
---|
72 | end |
---|
73 | else |
---|
74 | if vname(1)=='u' |
---|
75 | type='u'; |
---|
76 | ddd=1; |
---|
77 | cmin=-100; |
---|
78 | dc=20; |
---|
79 | cmax=100; |
---|
80 | cff=100; |
---|
81 | elseif vname(1)=='v' |
---|
82 | type='v'; |
---|
83 | ddd=1; |
---|
84 | cmin=-100; |
---|
85 | dc=20; |
---|
86 | cmax=100; |
---|
87 | cff=100; |
---|
88 | end |
---|
89 | end |
---|
90 | % |
---|
91 | % Parent |
---|
92 | % |
---|
93 | nc=netcdf('vortex_his.nc'); |
---|
94 | N=length(nc('s_rho')); |
---|
95 | time=round(nc{'scrum_time'}(tindex)/(24*3600)); |
---|
96 | disp(['Day : ',num2str(time)]) |
---|
97 | X=1e-3*nc{'x_rho'}(:); |
---|
98 | Y=1e-3*nc{'y_rho'}(:); |
---|
99 | if ddd==1 |
---|
100 | t1=cff*squeeze(nc{vname}(tindex,N,:,:)); |
---|
101 | else |
---|
102 | t1=cff*squeeze(nc{vname}(tindex,:,:)); |
---|
103 | end |
---|
104 | close(nc) |
---|
105 | [Xu,Xv,Xrp]=rho2uvp(X); |
---|
106 | [Yu,Yv,Yrp]=rho2uvp(Y); |
---|
107 | if type=='r' |
---|
108 | X1=X; |
---|
109 | Y1=Y; |
---|
110 | elseif type=='u' |
---|
111 | X1=Xu; |
---|
112 | Y1=Yu; |
---|
113 | elseif type=='v' |
---|
114 | X1=Xv; |
---|
115 | Y1=Yv; |
---|
116 | end |
---|
117 | % |
---|
118 | % Child |
---|
119 | % |
---|
120 | nc=netcdf('vortex_his.nc.1'); |
---|
121 | nestvortex=0; |
---|
122 | if ~isempty(nc) |
---|
123 | nestvortex=1; |
---|
124 | X=1e-3*nc{'x_rho'}(:); |
---|
125 | Y=1e-3*nc{'y_rho'}(:); |
---|
126 | if ddd==1 |
---|
127 | t=cff*squeeze(nc{vname}(tindex,N,:,:)); |
---|
128 | else |
---|
129 | t=cff*squeeze(nc{vname}(tindex,:,:)); |
---|
130 | end |
---|
131 | close(nc) |
---|
132 | [Xu,Xv,Xp]=rho2uvp(X); |
---|
133 | [Yu,Yv,Yp]=rho2uvp(Y); |
---|
134 | if type=='r' |
---|
135 | X2=X; |
---|
136 | Y2=Y; |
---|
137 | elseif type=='u' |
---|
138 | X2=Xu; |
---|
139 | Y2=Yu; |
---|
140 | elseif type=='v' |
---|
141 | X2=Xv; |
---|
142 | Y2=Yv; |
---|
143 | end |
---|
144 | Xbox=cat(1,Xp(1:end,1), ... |
---|
145 | Xp(end,1:end)' ,... |
---|
146 | Xp(end:-1:1,end),... |
---|
147 | Xp(1,end:-1:1)'); |
---|
148 | Ybox=cat(1,Yp(1:end,1), ... |
---|
149 | Yp(end,1:end)' ,... |
---|
150 | Yp(end:-1:1,end),... |
---|
151 | Yp(1,end:-1:1)'); |
---|
152 | end |
---|
153 | % |
---|
154 | % Plots |
---|
155 | % |
---|
156 | contour(X1,Y1,t1,[cmin:dc:cmax],'r') |
---|
157 | xlabel('X [km]') |
---|
158 | ylabel('Y [km]') |
---|
159 | title([vname,' - Day ',num2str(time)]) |
---|
160 | if nestvortex==1 |
---|
161 | hold on |
---|
162 | contour(X2,Y2,t,[cmin:dc:cmax],'k--') |
---|
163 | plot(Xbox,Ybox,'k'); |
---|
164 | hold off |
---|
165 | figure |
---|
166 | tint=interp2(X1,Y1,t1,X2,Y2,'cubic'); |
---|
167 | tdiff=tint-t; |
---|
168 | disp(['max difference = ',num2str(max(max(abs(tdiff))),2)]) |
---|
169 | disp(['relative difference = ',... |
---|
170 | num2str(100*max(max(abs(tdiff)))/max(max(abs(tint))),2),' %']) |
---|
171 | imagesc(flipud(tdiff)) |
---|
172 | colorbar |
---|
173 | title(['Parent - Child : ',vname,' - Day ',num2str(time)]) |
---|
174 | end |
---|