# Verification experiment, initially based on
# MITgcm/verification/tutorial_global_oce_biogeo

# ========
# 1. Get code
 cvs -d :pserver:cvsanon:cvsanon@mitgcm.org:/u/gcmpack co -D "11/28/17" MITgcm_code
 cvs -d :pserver:cvsanon:cvsanon@mitgcm.org:/u/gcmpack co MITgcm_contrib/ecco_darwin/v4_3deg
 cvs -d :pserver:cvsanon:cvsanon@mitgcm.org:/u/gcmpack co MITgcm_contrib/ecco_darwin/v4_llc270
 cvs -d :pserver:cvsanon:cvsanon@mitgcm.org:/u/gcmpack co -D "03/22/18" MITgcm_contrib/darwin/pkg/darwin
 cd MITgcm/pkg
 ln -sf ../../MITgcm_contrib/darwin/pkg/darwin .
 cd ..
 mkdir build run

# ================
# 2. Build executable
#    Prerequisite: 1. Get code
 cd build
 ../tools/genmake2 -ieee -mo \
  '../../MITgcm_contrib/ecco_darwin/v4_3deg/code ../../MITgcm_contrib/ecco_darwin/v4_llc270/code_darwin ../../MITgcm_contrib/ecco_darwin/v4_llc270/code'
 make depend
 make -j 8

# ======================
# 3. Run verification setup
#    Prerequisite: 2. Build executable
 cd ../run
 ln -sf ../build/mitgcmuv .
 cp ../../MITgcm_contrib/ecco_darwin/v4_llc270/input/data* .
 cp ../../MITgcm_contrib/ecco_darwin/v4_llc270/input_darwin/data* .
 cp ../../MITgcm_contrib/ecco_darwin/v4_3deg/input/*data* .
 ln -sf ../../MITgcm_contrib/ecco_darwin/v4_3deg/data/* .
 mkdir diags
 ./mitgcmuv > output.txt
# Compare to verification output
 diff <(grep %MON output.txt) <(grep %MON ../../MITgcm_contrib/ecco_darwin/v4_3deg/results/output.txt)

# ============================
# 4. Build and run MPI executable
#    Prerequisite: 1. Get code
 cd build
 rm *
 cp ../../MITgcm_contrib/ecco_darwin/v4_3deg/code/SIZE.h_mpi SIZE.h
 ../tools/genmake2 -mpi -mo \
  '../../MITgcm_contrib/ecco_darwin/v4_3deg/code ../../MITgcm_contrib/ecco_darwin/v4_llc270/code_darwin ../../MITgcm_contrib/ecco_darwin/v4_llc270/code'
 make depend
 make -j 8
 cd ../run
 mkdir diags
 ln -sf ../build/mitgcmuv .
 cp ../../MITgcm_contrib/ecco_darwin/v4_llc270/input/data* .
 cp ../../MITgcm_contrib/ecco_darwin/v4_llc270/input_darwin/data* .
 cp ../../MITgcm_contrib/ecco_darwin/v4_3deg/input/*data* .
 mv data_mpi data
 ln -sf ../../MITgcm_contrib/ecco_darwin/v4_3deg/data/* .
 mpirun -np 8 ./mitgcmuv &
# Monitor run
 tail -f STDOUT.0000 | grep advcfl_W

# ============================
# 5. MATLAB code for computing volume, salt, salinity, and DIC budgets
#    Prerequisite: 4. Build and run MPI executable
#    Can be executed as soon as 3 or more months of output are available
 cd ../../MITgcm_contrib/ecco_darwin/v4_3deg/matlab
# start matlab
# if using gcmfaces: budget_v4_3deg_with_gcmfaces
# if not using gcmfaces: budget_v4_3deg_without_gcmfaces