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

Annotation of /MITgcm_contrib/ocean_inversion_project/takahashi/README

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


Revision 1.1 - (hide 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 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]');

  ViewVC Help
Powered by ViewVC 1.1.22