--- MITgcm_contrib/ocean_inversion_project/README	2003/09/23 04:34:23	1.3
+++ MITgcm_contrib/ocean_inversion_project/README	2003/09/24 06:50:26	1.5
@@ -35,8 +35,8 @@
  cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/.genmakerc .
  cp ../ocean_inversion_project/code/* .
- cp ptracers_forcing.F.test ptracers_forcing.F
- rm ptracers_init.F ptracers_forcing.F
+ rm ptracers_init.F ptracers_forcing_surf.F PTRACERS.h
+ rm ptracers_read_mask.F ptracers_read_takahashi.F
  ../tools/genmake
  make depend
  make
@@ -51,10 +51,9 @@
  cp ../ocean_inversion_project/input/* .
  cp data.test data
  cp data.ptracers.test data.ptracers
- mitgcmuv > output .txt
+ mitgcmuv > output.txt
   
-5 ===> check that PTRACER output for tracer 1 and tracer 30
-  ===> is identical to salinity output
+5 ===> check that PTRACER output and salinity output are identical.
 
  diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
  diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data
@@ -102,6 +101,7 @@
 caxis([-6 6]), colorbar, plotland
 
 % some matlab code for looking at PTRACER output files
+salt=readbin('S.0000000020.001.001.data',[90 40 15],1);
 tracer=zeros(90,40,15,30);
 for i=1:30
  fn=['PTRACER' myint2str(i) '.0000000020.001.001.data'];
@@ -109,6 +109,8 @@
  clf, mypcolor(tracer(:,:,1,i)'); pause(1)
 end
 tmp=sum(tracer,4); clf, mypcolor(tmp(:,:,1)')
+mypcolor(
+
 
 % some matlab code for checking that one year's worth of
 % tracer uptake is approximately 1e18 mols.
@@ -134,7 +136,7 @@
 thk=[50 70 100 140 190 240 290 340 390 440 490 540 590 640 690];
 tracer=zeros(90,40,15,30); sumtracer=zeros(30,1);
 for i=1:30, mydisp(i)
- fn=['PTRACER' myint2str(i) '.0000000180.001.001.data'];
+ fn=['PTRACER' myint2str(i) '.0000003600.001.001.data'];
  tracer(:,:,:,i)=readbin(fn,[90 40 15],1);
  for j=1:length(lat)
   for k=1:length(thk)
@@ -143,3 +145,4 @@
   end
  end
 end
+plot(1:30,0*sumtracer,1:30,sumtracer)