--- MITgcm_contrib/ocean_inversion_project/README	2003/09/23 04:34:23	1.3
+++ MITgcm_contrib/ocean_inversion_project/README	2003/10/10 16:30:39	1.9
@@ -35,8 +35,7 @@
  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_*.F
  ../tools/genmake
  make depend
  make
@@ -51,10 +50,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
@@ -88,6 +86,36 @@
  mitgcmuv > output .txt &
 
 
+==========================================================
+MPI instructions for carrying out a 3000-year quasi-stationary
+integration using the global_ocean.90x40x15 configuration.
+==========================================================
+
+1 ===> compile, link, and execute
+
+ cd MITgcm/exe
+ rm *
+ cd ../bin
+ rm *
+ cp ../verification/global_with_exf/code/* .
+ cp ../ocean_inversion_project/code/.genmakerc .
+ cp ../ocean_inversion_project/code/* .
+ cp CPP_EEOPTIONS.h_mpi CPP_EEOPTIONS.h
+ cp SIZE.h_mpi SIZE.h
+ ../tools/genmake
+ make depend
+ make
+ cd ../exe
+ cp ../verification/global_with_exf/input/eedata .
+ cp ../verification/global_with_exf/input/data.* .
+ cp ../verification/global_with_exf/input/POLY3.COEFFS .
+ ln -sf ../verification/global_with_exf/input/*.bin .
+ cp ../ocean_inversion_project/input/* .
+ ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
+ ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
+ mpirun -np 36 dplace -s1 -c3-38 mitgcmuv < /dev/null > & ! errlog &
+
+
 ======================================================
 
 % some matlab code for looking at fort.10 debug files
@@ -102,6 +130,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,37 +138,23 @@
  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.
-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) '.0000000180.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
 end
+plot(1:30,0*sumtracer,1:30,sumtracer)