| 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/* . |
| 141 |
cp ../ocean_inversion_project/code/* . |
cp ../ocean_inversion_project/code/* . |
| 142 |
cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h |
cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h |
| 143 |
cp ../ocean_inversion_project/write_netCDF/* . |
cp ../ocean_inversion_project/write_netCDF/* . |
| 144 |
|
\rm mk_output.F |
| 145 |
|
|
| 146 |
### need customized option file; example below is for nireas |
### need customized option file; example below is for nireas |
| 147 |
../tools/genmake2 -of linux_ia32_g77_netcdf |
../tools/genmake2 -of linux_ia32_g77_netcdf |
| 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)) |