#!/bin/csh
# rough instruction for building and running the test and compare a
# cubed sphere with 96 tiles and no holes to a cubed sphere with 91 tiles
# and 5 missing tiles. Basically you can use this file as a script after
# modifying a the variable myFavoriteOptfile and maybe commenting in a few
# lines.
# The configuration is based on verification/global_ocean.cs32x15 but with
# seaice instead of thsice.
# The test starts from a pickup files of the 2nd time steps to have non-zero
# velocities etc at the very beginnning of the integration. The momentum 
# forcing is specified as stress, but only because I did not have the correct
# winds interpolated to the sphere. In order to use winds, change the 
# appropriate flag in code/EXF_OPTIONS.h (#define ALLOW_ATM_WIND) and use
# the stresses as input (will be wrong but for one timestep it just means too
# little wind).
# 
set myFavoriteOptfile = '../../../../MITgcm/tools/build_options/darwin_ppc_g77'
#cvs co MITgcm_code (or the entire code if you must)
#cvs co MITgcm_contrib/mlosch/cs32test
cd MITgcm_contrib/mlosch/cs32test
mkdir build_s96t
cd build_s96t
../../../../MITgcm/tools/genmake2 -of $myFavoriteOptfile -rootdir ../../../../MITgcm -mods "../s96t ../code"
make depend
make
cd ../
mkdir build_s91t
cd build_s91t
../../../../MITgcm/tools/genmake2 -of $myFavoriteOptfile -rootdir ../../../../MITgcm -mods "../s91t ../code"
make depend
make
cd ../
# 
mkdir run
cd run
ln -s ../input/* .
# this one works and gives a reasonable result
../build_s96t/mitgcmuv > output.s96t
# this one gives NaNs after the first timestep
../build_s91t/mitgcmuv > output.s91t