--- MITgcm_contrib/ocean_inversion_project/README	2003/12/18 03:17:56	1.19
+++ MITgcm_contrib/ocean_inversion_project/README	2004/09/30 23:45:07	1.21
@@ -20,9 +20,11 @@
 
 1 ===> get MITgcm code from cvs repository
 
- CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
+ CVSROOT=:pserver:cvsanon@mitgcm.org:/u/gcmpack
  cvs login ( CVS password: cvsanon )
- cvs co -r checkpoint51n_branch MITgcm
+ cvs co -r checkpoint51n_branch MITgcm_code
+ cvs co -r checkpoint51n_branch MITgcm/verification/testreport
+ cvs co -r checkpoint51n_branch MITgcm/verification/global_with_exf
 
 2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory
 
@@ -31,6 +33,7 @@
 
 3 ===> compile and link
 
+ mkdir bin exe
  cd bin
  cp ../verification/global_with_exf/code/* .
  cp ../ocean_inversion_project/code/* .
@@ -138,6 +141,7 @@
  cp ../ocean_inversion_project/code/* .
  cp PTRACERS_OPTIONS.h.netcdf PTRACERS_OPTIONS.h
  cp ../ocean_inversion_project/write_netCDF/*  .
+ \rm mk_output.F
 
 ### need customized option file; example below is for nireas
  ../tools/genmake2 -of linux_ia32_g77_netcdf
@@ -218,3 +222,19 @@
  end
 end
 plot(1:30,0*sumtracer,1:30,sumtracer)
+
+
+% check that tracer flux is 1e18 mols / year
+nb_seconds_per_year=31556880;
+lon=2:4:360; lat=-78:4:78;
+mask=readbin('hFacC.data',[90 40 15],1);
+tracer=zeros(90,40,30); sumtracer=zeros(30,1);
+for i=1:30, mydisp(i)
+ fn=['PtrFlux' myint2str(i) '.0000000180.data'];
+ tracer(:,:,i)=readbin(fn,[90 40],1);
+ for j=1:length(lat)
+    sumtracer(i) = sumtracer(i) + sum(tracer(:,j,i).*mask(:,j,1)) * ...
+        nb_seconds_per_year* (4*1.113195e+05)^2 * cos(pi*lat(j)/180);
+  end
+end
+plot(1:30,0*sumtracer,1:30,sumtracer-mmean(sumtracer))