1 |
dimitri |
1.1 |
This directory contains the monthly mean Takahashi 2002 CO2 flux climatology |
2 |
|
|
based on the Wanninkhof parameterization of air-sea gas exchange. Units are |
3 |
|
|
mol/m^2/yr. The ocean CO2 flux has been extrapolated to cover land regions in |
4 |
|
|
order to prevent errors in the interpolation of coastal grid boxes from |
5 |
|
|
differing land-sea masks. |
6 |
|
|
|
7 |
|
|
|
8 |
|
|
Instructions for generating taka02_montlhy.bin |
9 |
|
|
============================================== |
10 |
|
|
|
11 |
|
|
1 ===> get taka02_montlhy.cdf from |
12 |
|
|
quercus.igpp.ucla.edu/OceanInversion/inputs/takahashi/ |
13 |
|
|
|
14 |
|
|
2 ===> generate taka02_montlhy.bin and *_taka02.asc |
15 |
|
|
|
16 |
|
|
% matlab code below requires netcdf toolbox and some other tools |
17 |
|
|
% that can be obtained by contacting menemenlis@jpl.nasa.gov |
18 |
|
|
f=netcdf('taka02_montlhy.cdf'); |
19 |
|
|
mo=f{1}(:); lat=f{2}(:); lon=f{3}(:); tmp=f{4}(:); |
20 |
|
|
dx=[0; diff(lon)]; dy=[0; diff(lat)]; |
21 |
|
|
flux=zeros(length(lon),length(lat),length(mo)); |
22 |
|
|
for i=1:12 |
23 |
|
|
flux(:,:,i)=squeeze(shiftdim(tmp(i,:,:),2)); |
24 |
|
|
end |
25 |
|
|
clf,contourf(lon,lat,mean(flux,3)',-5:5) |
26 |
|
|
plotland,caxis([-5 5]),colorbar |
27 |
|
|
writebin('taka02_montlhy.bin',flux,1); |
28 |
|
|
fid = fopen('lon_taka02.asc','w'); |
29 |
|
|
fprintf(fid,'%7.3f %7.3f\n',[lon dx]'); |
30 |
|
|
fclose(fid); |
31 |
|
|
fid = fopen('lat_taka02.asc','w'); |
32 |
|
|
fprintf(fid,'%8.4f %8.4f\n',[lat dy]'); |