4 |
cd MITgcm/verification/lab_sea |
cd MITgcm/verification/lab_sea |
5 |
mkdir build |
mkdir build |
6 |
cd build |
cd build |
7 |
cp ../code/*.h ../code/packages.conf . |
cp ../code/*.h ../code/*.F ../code/packages.conf . |
8 |
../../../tools/genmake2 |
../../../tools/genmake2 |
9 |
make depend |
make depend |
10 |
make |
make |
138 |
|
|
139 |
matlab |
matlab |
140 |
|
|
141 |
pn='nokpp_nogm_withseaicedyn/'; |
clear all, clf reset, pn=''; nt=238; |
142 |
|
heff=zeros(20,16,nt,21); |
143 |
heff=zeros(20,16,21,21); |
area=zeros(20,16,nt,21); |
|
area=zeros(20,16,21,21); |
|
144 |
for e=0:20 |
for e=0:20 |
145 |
for t=0:4, disp([e t]) |
for t=0:nt, disp([e t]) |
146 |
if e==0, fn=[pn 'run/HEFF.']; |
if e==0, fn=[pn 'run/HEFF.']; |
147 |
else fn=[pn 'run' int2str(e) '/HEFF.']; end |
else fn=[pn 'run' int2str(e) '/HEFF.']; end |
148 |
fn=[fn myint2str(t*240,10) '.data']; |
fn=[fn myint2str(t*240,10) '.data']; |
153 |
area(:,:,t+1,e+1)=readbin(fn,[20 16],1); |
area(:,:,t+1,e+1)=readbin(fn,[20 16],1); |
154 |
end |
end |
155 |
end |
end |
156 |
|
for t=1:nt |
|
for t=1:21 |
|
157 |
tmp=squeeze(sum(sum(heff(:,:,t,:)))); |
tmp=squeeze(sum(sum(heff(:,:,t,:)))); |
158 |
clf,subplot(211),plot(tmp-mean(tmp)) |
clf,subplot(211),plot(tmp-mean(tmp)) |
159 |
grid,title([int2str((t-1)*10) ' days']) |
grid,title([int2str((t-1)*10) ' days']) |