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/).


===============================================
First check that pkg/ptracers works OK by using
salinity initial and boundary conditions
===============================================

1 ===> get MITgcm code from cvs repository

 CVSROOT=:pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
 cvs login ( CVS password: cvsanon )
 cvs co MITgcm

2 ===> put MITgcm_contrib/ocean_inversion_project in MITgcm directory

 cd MITgcm
 cvs co -d ocean_inversion_project MITgcm_contrib/ocean_inversion_project

3 ===> compile

 cd bin
 cp ../verification/global_with_exf/code/* .
 cp ../ocean_inversion_project/code/.genmakerc .
 cp ../ocean_inversion_project/code/* .
 cp ../ocean_inversion_project/code/ptracers_forcing.F.test ptracers_forcing.F
 cp ../ocean_inversion_project/code/ptracers_init.F.test ptracers_init.F
 cp ../ocean_inversion_project/code/ptracers_read_mask.F.test ptracers_read_mask.F
 ../tools/genmake
 make depend
 make

4 ===> execute

 cd ../exe
 cp ../verification/global_with_exf/input/* .
 cp ../ocean_inversion_project/input/* .
 cp ../ocean_inversion_project/input/data.ptracers.test data.ptracers
 mitgcmuv > output .txt
  
5 ===> check that PTRACER output for tracer 1 and tracer 30
  ===> is identical to salinity output

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


======================================================
Second generate binary input files.  netcdf files are not OK

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

 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

===> 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