| 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]'); |