--- MITgcm/verification/lab_sea/README	2002/12/28 10:11:11	1.4
+++ MITgcm/verification/lab_sea/README	2003/04/30 07:04:08	1.6
@@ -24,8 +24,8 @@
   vwindFile      = 'v10m.labsea79'    # 10-m meridional wind
   atempFile      = 'tair.labsea1979'  # 2-m air temperature
   aqhFile        = 'qa.labsea1979'    # 2-m specific humidity
-  lwfluxFile     = 'flo.labsea1979'   # longwave radiation
-  swfluxFile     = 'fsh.labsea1979'   # shortwave radiation
+  lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
+  swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
   precipFile     = 'prate.labsea1979' # precipitation
   evapFile       = 'evap.labsea1979'  # evaporation
 
@@ -84,6 +84,21 @@
   make
   mv mitgcmuv mitgcmuv_exf_forcing
 
+1-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP
+  cd ../../../verification/lab_sea/input
+  ln -sf ../code/SIZE.h .
+  ln -sf ../code/CPP_EEOPTIONS.h .
+  ln -sf ../code/CPP_OPTIONS.h .
+  ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h
+  ln -sf ../code/SEAICE_OPTIONS.h .
+  ../../../tools/genmake -makefile
+    ==> on alhena use:
+    ../../../tools/genmake -platform=o2k_noopt -makefile
+  make clean
+  make depend
+  make
+  mv mitgcmuv mitgcmuv_compute_evap
+
 Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES:
   cd ../../../verification/lab_sea/input
   ln -sf ../code/SIZE.h .
@@ -180,16 +195,14 @@
   mkdir exp3b
   mv *tave.0000000001.data exp3b
 
-Use matlab script lookat_exp3.m to compare
-the output of exp2 to that of exp1:
+Use matlab script lookat_exp3.m to look at
+results from exp3a and exp3b:
   cd ../../../verification/lab_sea/matlab
   matlab
   lookat_exp3
 
 If the dynamic ice solvers are correct, they should converge to a
-constant solution.  The tests above show that both the ADI and the LSR
-solver are unable to handle periodic domains and therefore that they
-cannot be "correctly" parallelized.
+solution that is constant at each latitude.
 
 
 Instructions for running Experiment 4
@@ -217,7 +230,8 @@
 Also as sea-ice ages (thickens) the difference at the tiles becomes
 increasingly small.  For solver pkg/seaice/adi.F a value of NPSEUDO=10
 appears adequate for forward integrations but cannot be used for sea-ice
-adjoint model.  Work is underway to fix this problem.
+adjoint model.  With NPSEUDO=300 or greater, the tile disappears, but the
+computational cost is very large (11 times that of teh forward model).
 
 
 Instructions for running Experiment 5
@@ -337,3 +351,24 @@
   cd ../../../verification/lab_sea/matlab
   matlab
   lookat_exp9
+
+
+Instructions for running Experiment 10
+======================================
+
+This is a 1-cpu, 10-hour integration similar to exp8, but
+with bulk formula computation of evaporation fields.
+
+To run Experiment 10:
+  cd ../../../verification/lab_sea/input
+  ln -sf DATA.10hours_exf data
+  ln -sf DATA.SEAICE.adi_exf data.seaice
+  mitgcmuv_compute_evap >&! output.txt
+  mkdir exp10
+  mv *tave.0000000010.data exp10
+
+Use matlab script lookat_exp10.m to compare
+the output of exp10 to that of exp8:
+  cd ../../../verification/lab_sea/matlab
+  matlab
+  lookat_exp10