Ocean Inversion Project
=======================

Instructions for using pkg/ptracers to compute tracer
Green's functions for Gruber's ocean inversion project
(http://quercus.igpp.ucla.edu/OceanInversion/).

=====================================

To run a pkg/ptracers test case using salt initial and boundary conditions:

 
 
 cd MITgcm/bin
 cp ../verification/global_with_exf/code/* .
 cp ../code/.genmakerc .
 cp ../code/* .
 cp ../code/ptracers_forcing.F.test ptracers_forcing.F
 cp ../code/ptracers_init.F.test ptracers_init.F
 cp ../code/ptracers_read_mask.F.test ptracers_read_mask.F
 ../tools/genmake
 make depend
 make
 cd ../exe
 cp ../verification/global_with_exf/input/* .
 cp ../input/* .
 cp ../input/data.ptracers.test data.ptracers
 mitgcmuv > output .txt
 
===> check that PTRACER output is identical to salinity output

 diff PTRACER30.0000000020.001.001.data PTRACER01.0000000020.001.001.data
 diff PTRACER30.0000000020.001.001.data S.0000000020.001.001.data

===> examine interpolation output for 30-region mask

 matlab
 mask=zeros(94,44,30);
 load fort.10
 for i=1:length(fort)
  mask(fort(i,2)+2,fort(i,3)+2,fort(i,1))=fort(i,4);
 end 
 for i=1:30
  mypcolor(mask(3:92,3:42,i)');
  colorbar,title(i),pause
 end
 mypcolor(sum(mask,3)'); colorbar

=====================================


Global ocean experiment with ptracers
=====================================

===> set up and integrate 30-tracer computation

 cd MITgcm/bin
 cp ../verification/global_with_exf/code/* .
 cp ../code/.genmakerc .
 cp ../code/* .
 ../tools/genmake
 make depend
 make
 cd ../exe
 cp ../verification/global_with_exf/input/* .
 cp ../input/* .
 mitgcmuv > output .txt

===> look at output

 matlab
 for i=1:30
  tr=readbin(['PTRACER' myint2str(i) ...
              '.0000000020.001.001.data'],[90 40 15],1);
  mypcolor(tr(:,:,2)'), colorbar, pause(1)
 end