--- MITgcm_contrib/ocean_inversion_project/README	2003/09/25 03:01:59	1.6
+++ MITgcm_contrib/ocean_inversion_project/README	2003/09/26 20:01:41	1.7
@@ -115,31 +115,14 @@
 % tracer uptake is approximately 1e18 mols.
 lon=2:4:360; lat=-78:4:78;
 thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
-tracer=zeros(90,40,15,30); sumtracer1=zeros(30,1);
-for i=1:30, mydisp(i)
- fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];
- tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
- for x=1:length(lon)
-  for y=1:length(lat)
-   for z=1:length(thk)
-    sumtracer1(i) = sumtracer1(i) + tracer(x,y,z,i) * ...
-                    thk(z) * (4*1.113195e+05)^2 * cos(pi*lat(y)/180);
-   end
-  end
- end
-end
-
-% some matlab code for checking that one year's worth of
-% tracer uptake is approximately 1e18 mols.
-lon=2:4:360; lat=-78:4:78;
-thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
+mask=readbin('hFacC.001.001.data',[90 40 15],1);
 tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
 for i=1:30, mydisp(i)
- fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];
+ fn=['PTRACER' myint2str(i) '.0000001800.001.001.data'];
  tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
  for j=1:length(lat)
   for k=1:length(thk)
-   sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i)) * ...
+   sumtracer(i) = sumtracer(i) + sum(tracer(:,j,k,i).*mask(:,j,k)) * ...
                   thk(k) * (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
   end
  end