/[MITgcm]/MITgcm_contrib/ocean_inversion_project/takahashi/README
ViewVC logotype

Contents of /MITgcm_contrib/ocean_inversion_project/takahashi/README

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Tue Sep 23 04:34:24 2003 UTC (21 years, 10 months ago) by dimitri
Branch: MAIN
CVS Tags: HEAD
o Mods and bug fixes to pkg/cal and pkg/exf needed for computation
  of tracer Green's fucntions for ocean inversion project.

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

  ViewVC Help
Powered by ViewVC 1.1.22