20 |
|
|
21 |
1 ===> get MITgcm code from cvs repository |
1 ===> get MITgcm code from cvs repository |
22 |
|
|
23 |
CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack |
CVSROOT=:pserver:cvsanon@mitgcm.org:/u/gcmpack |
24 |
cvs login ( CVS password: cvsanon ) |
cvs login ( CVS password: cvsanon ) |
25 |
cvs co -r checkpoint51n_branch MITgcm |
cvs co -r checkpoint51n_branch MITgcm_code |
26 |
|
cvs co -r checkpoint51n_branch MITgcm/verification/testreport |
27 |
|
cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf |
28 |
|
|
29 |
2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory |
2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory |
30 |
|
|
33 |
|
|
34 |
3 ===> compile and link |
3 ===> compile and link |
35 |
|
|
36 |
|
mkdir bin exe |
37 |
cd bin |
cd bin |
38 |
cp ../verification/global_with_exf/code/* . |
cp ../verification/global_with_exf/code/* . |
39 |
cp ../ocean_inversion_project/code/* . |
cp ../ocean_inversion_project/code/* . |
222 |
end |
end |
223 |
end |
end |
224 |
plot(1:30,0*sumtracer,1:30,sumtracer) |
plot(1:30,0*sumtracer,1:30,sumtracer) |
225 |
|
|
226 |
|
|
227 |
|
% check that tracer flux is 1e18 mols / year |
228 |
|
nb_seconds_per_year=31556880; |
229 |
|
lon=2:4:360; lat=-78:4:78; |
230 |
|
mask=readbin('hFacC.data',[90 40 15],1); |
231 |
|
tracer=zeros(90,40,30); sumtracer=zeros(30,1); |
232 |
|
for i=1:30, mydisp(i) |
233 |
|
fn=['PtrFlux' myint2str(i) '.0000000180.data']; |
234 |
|
tracer(:,:,i)=readbin(fn,[90 40],1); |
235 |
|
for j=1:length(lat) |
236 |
|
sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ... |
237 |
|
nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180); |
238 |
|
end |
239 |
|
end |
240 |
|
plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer)) |