--- MITgcm_contrib/eh3/llc/ecco-godae/README	2006/11/14 00:57:02	1.3
+++ MITgcm_contrib/eh3/llc/ecco-godae/README	2007/03/28 18:29:45	1.5
@@ -17,12 +17,16 @@
  
 cd MITgcm/utils/exch2/matlab-topology-generator/
  
-modify matlab script driver_llc.m
+modify matlab script driver.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
+nr = 360; ng = 90; nb = 90; tnx=45;tny=45;
+for compact format, set 
+mapIO=1
+make sure blanklist.txt is correct 
+(or empty, if no tiles to be omitted)
+then run driver.m
  
 This generates routine and header
 w2_e2setup.F
@@ -42,7 +46,6 @@
 
 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/
 -------------------------------
@@ -57,3 +60,64 @@
 (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/
+
+Required input to generate bathy & initial conditions:
+(1)
+(in directory input/ )
+llc_p90_1.nc
+llc_p90_2.nc
+llc_p90_3.nc
+llc_p90_4.nc
+llc_p90_5.nc
+llc_p90_6.nc
+(2)
+ETOPO2.raw.bin
+(3)
+directory CLIMannual/
+
+Generation steps:
+
+1) run script "gen_bathy"
+   Input:
+   -----
+llc_p90_1.nc
+llc_p90_2.nc
+llc_p90_3.nc
+llc_p90_4.nc
+llc_p90_5.nc
+llc_p90_6.nc
+   Output:
+   -----
+llc_p90_bathy.f001.nc
+llc_p90_bathy.f002.nc
+llc_p90_bathy.f003.nc
+llc_p90_bathy.f004.nc
+llc_p90_bathy.f005.nc
+
+2) run script "gen_initial_cond_50_lev"
+   Input:
+   -----
+(a) above llc_p90_*.nc
+(b) orig. climatology fields in CLIMannual/
+   Output:
+   -----
+
+3) run script "t18_to_global_inputs.m"
+
+4) run script "fill_shallow_regions_and_isolated_basins.m"
+