--- MITgcm_contrib/ocean_inversion_project/README	2003/12/18 03:17:56	1.19
+++ MITgcm_contrib/ocean_inversion_project/README	2004/10/06 20:26:54	1.22
@@ -20,9 +20,9 @@
 
 1 ===> get MITgcm code from cvs repository
 
- CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/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,27 +31,29 @@
 
 3 ===> compile and link
 
+ mkdir bin exe
  cd bin
- cp ../verification/global_with_exf/code/* .
- cp ../ocean_inversion_project/code/* .
- rm ptracers_*.F PTRACERS_OPTIONS.h
+ \cp ../verification/global_with_exf/code/* .
+ \cp ../ocean_inversion_project/code/* .
+ \rm ptracers_*.F PTRACERS_OPTIONS.h
  ../tools/genmake2
-   (on orion: ../tools/genmake2 -of linux_ia64_efc+mpi)
+   (on orion:    ../tools/genmake2 -of linux_ia64_efc+mpi       )
+   (on columbia: ../tools/genmake2 -of linux_ia64_efc+mpi_altix )
  make depend
- make
+ make -j 16
 
 4 ===> execute
 
  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/* .
- cp data.test data
- cp data.ptracers.test data.ptracers
- mitgcmuv > output.txt
-  
+ \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/* .
+ \cp data.test data
+ \cp data.ptracers.test data.ptracers
+ ./mitgcmuv >! output.txt
+
 5 ===> check that PTRACER output and salinity output are identical.
 
  diff PTRACER01.0000000020.001.001.data S.0000000020.001.001.data
@@ -90,6 +92,7 @@
  mpirun -np 24 dplace -s1 mitgcmuv < /dev/null > & ! errlog &
    (on orion: qsub job.orion)
 
+
 ==========================================================
 MPI instructions for carrying out time-dependent, 1765-2005,
 anthropogenic carbon perturbation tracer experiments using
@@ -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
@@ -157,7 +161,35 @@
  mitgcmuv > output .txt &
 
 
-======================================================
+==========================================================
+Instructions for carrying out time-dependent, 1765-2005,
+anthropogenic carbon perturbation tracer experiments using
+the ecco1x1 configuration and netcdf output on columbia.
+This assumes that the appropriate libnetcdf.a exists
+and that it be specified in the genmake2 option file.
+It also assumes that surface forcing files are available
+in /nobackup2/menemenl/ocmip/iter69
+==========================================================
+
+ cd MITgcm/exe
+ \rm *
+ cd ../bin
+ \rm *
+ \cp ../ocean_inversion_project/code_ecco1x1/* .
+ \cp ../ocean_inversion_project/write_netCDF/*  .
+ \rm mk_output.F
+ \cp SIZE.h.144 SIZE.h
+ ../tools/genmake2 -of linux_ia64_efc+mpi_altix
+ make depend
+ make -j 16
+ cd ../exe
+ \cp ../ocean_inversion_project/input_ecco1x1/* .
+ ln -sf /nobackup2/menemenl/ocmip/iter69/* .
+ ln -sf ../ocean_inversion_project/region_mask/30reg_regionmask.bin .
+ ln -sf ../ocean_inversion_project/takahashi/taka02_montlhy.bin .
+ ln -sf ../ocean_inversion_project/atm_co2/splco2_cis92a.dat .
+ qsub job01
+
 
 ==========================================================
 Instructions for generating netcdf output files on orion
@@ -218,3 +250,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))