LAT-LON-Polar-Cap setup
-----------------------
ed@eh3.com, cnh@mit.edu, heimbach@mit.edu

Setup code/, input/
-------------------
Note that the *.mitgrid files must be named using the following pattern
or they will not be read correctly by MITgcm:

  cd run
  for i in 1 2 3 4 5 6; do 
    ln -s ../input/llc_p90_${i}.mitgrid tile00${i}.mitgrid
  done

Change topology:
 
cd MITgcm/utils/exch2/matlab-topology-generator/
 
modify matlab script driver_llc.m
line
nr = 360; ng = 90; nb = 90; tnx=90;tny=90;
to (e.g.)
nr = 360; ng = 90; nb = 90; tnx=30;tny=30;
run script
 
This generates routine and header
w2_e2setup.F
W2_EXCH2_TOPOLOGY.h
which need to go in code/

Accordingly, bathy.0??.0??.data need changing.

N.B. Domain decomposition nPy not active with exch2, i.e.
possible decomps are (for orig. with 17 tiles):
sNx = 90 30 90  30
sNy = 90 30 90  30
nSx = 17 17  1   1
nSy =  1  1  1   1
nPx =  1  9 17 153
nPy =  1  1  1   1

I've generated artificial tile #18 with all zeros
blanklist.txt needs to be generated by hand.
Not clear how many procs to invoke in presence of blank list.

Setup code_15lev/, input_15lev/
-------------------------------

* Real bathymetry, 15 vertical levels,
* monthly climatological forcing via exf package, uses
  4x4 deg. fields with on-the-fly interpolation
* Bathymetry is generated via script ... in ...
* Initial conditions for T, S are generated via script ... in ...

All input fields can be downloaded using script input_15lev/prepare_run
(which also contains link pattern to standard tile00${i}.mitgrid fields),
then link to run directory.

Setup code_50lev/, input_50lev/
-------------------------------

* Based on ETOPO2 bathy, 
* 50 vertical level exactly matching high_res_cube (ECCO2), incl. partial cell params.
* forced by 6-hourly NCEP atmos. state, uses original reduced Gaussian grid
  and exf interpolation on-the-fly
* initial conditions generated via scripts in 
  MITgcm_contrib/eh3/llc/ecco-godae/climatology/
  resulting fields are
  llc90x90x360_50lev_ptemp.bin, llc90x90x360_50lev_salt.bin
  these fields are run are run with a diffusion/conv.adj. setup of MITgcm
  for 10 days to remove some instabilities
* input fields can be obtained at
  http://channel.mit.edu/datasets/MITgcm_input_files/MITgcm_contrib/eh3/llc/ecco-godae/data_50lev/