Equilibrium integrations for dissolved inorganic
carbon (DIC) and for radiocarbon (C14).
================================================

Instructions for using pkg/ptracers for equilibrium DIC
and C14 integrations as per OCMIP 2:
http://www.cgd.ucar.edu/oce/OCMIP/HOWTO-Abiotic.pdf

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

1 ===> get MITgcm code from cvs repository

 cvs co -r checkpoint57w_post MITgcm_code
 cvs co -r checkpoint57w_post MITgcm/verification/testreport
 cvs co -r checkpoint57w_post MITgcm/verification/global_with_exf

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 and link

 mkdir bin exe
 cd bin
 \cp ../verification/global_with_exf/code/* .
 \cp ../ocean_inversion_project/code_c14/* .
 ../tools/genmake2
 make depend
 make -j
 \cp mitgcmuv ../exe

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 ../ocean_inversion_project/input_c14/* .
 ./mitgcmuv >! output.txt

5 ===> check that PTRACER output and salinity output are identical.

 diff PTRACER01.0000000020.data S.0000000020.data
 diff PTRACER02.0000000020.data S.0000000020.data


==========================================================
Instructions for carrying out 3000-year quasi-stationary
integration using the ecco1x1 configuration 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 MITgcm/../iter69
==========================================================

 module purge
 module load modules scsl.1.5.0.0 intel-comp.8.1.024 mpt.1.12.0.0
 cd MITgcm
 \rm exe/* bin/*
 cd bin
 \cp ../ocean_inversion_project/code_c14/* .
 \cp SIZE.h.ecco1x1 SIZE.h
 \cp ptracers_forcing_surf.F.ecco1x1 ptracers_forcing_surf.F
 ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas
 make depend
 make -j
 \cp mitgcmuv ../exe
 cd ../exe
 \cp ../ocean_inversion_project/input_c14/* .
 \mv data.ecco1x1 data
 \mv data.ptracers.ecco1x1 data.ptracers
 \ln -sf ../../iter69/* .
 qsub job_c14_ecco1x1