Annotation of /MITgcm/verification/lab_sea/README

Example: Labrador Sea Region with Sea-Ice
=========================================

This example sets up a small (20x16x23) Labrador Sea experiment
coupled to a dynamic thermodynamic sea-ice model.
A brief description of the sea-ice model is in "seaice.ps".

The domain of integration spans 280E to 320E and 46N to 78N.
Horizontal grid spacing is 2 degrees.
The 23 vertical levels and the bathymetry file
  bathyFile      = 'bathy.labsea'
are obtained from the the 2-degree ECCO configuration.

Integration is initialized from annual-mean Levitus climatology
 hydrogThetaFile = 'LevCli_temp.labsea'
 hydrogSaltFile  = 'LevCli_salt.labsea'

Surface salinity relaxation is to the monthly mean Levitus climatology
 saltClimFile    = 'SSS.labsea'

Forcing files are a 1979-1999 monthly climatology computed from the
NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)
  uwindFile      = 'u10m.labsea79'    # 10-m zonal wind
  vwindFile      = 'v10m.labsea79'    # 10-m meridional wind
  atempFile      = 'tair.labsea1979'  # 2-m air temperature
  aqhFile        = 'qa.labsea1979'    # 2-m specific humidity
  lwdownFile     = 'flo.labsea1979'   # downward longwave radiation
  swdownFile     = 'fsh.labsea1979'   # downward shortwave radiation
  precipFile     = 'prate.labsea1979' # precipitation

The experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
The test is a 1-cpu, 10-hour integration.   Both the atmospheric
state and the open-water surface fluxes are provided by pkg/exf.


Using testscript to test sea-ice code
=====================================

Running the testscript experiment:
  cd MITgcm/verification
  ./testreport -t lab_sea
The default experiment is Experiment 8, below.
It uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.

Note that fairly large differences in accuracy occur across different
platforms.  For example, testscript comparisons between g77 (Linux)
and f77 (SGI) generated output gives:

                T           S           U           V
C D M    c        m  s        m  s        m  s        m  s
n p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
f n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
g d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

Y Y Y Y  8 10  9 11 10  9 11 13 10  9  8  8  9  8  9  8  9 FAIL  lab_sea
1	heimbach	1.2	Example: Labrador Sea Region with Sea-Ice
2			=========================================
3
4			This example sets up a small (20x16x23) Labrador Sea experiment
5			coupled to a dynamic thermodynamic sea-ice model.
6	dimitri	1.3	A brief description of the sea-ice model is in "seaice.ps".
7	heimbach	1.2
8			The domain of integration spans 280E to 320E and 46N to 78N.
9			Horizontal grid spacing is 2 degrees.
10			The 23 vertical levels and the bathymetry file
11			bathyFile = 'bathy.labsea'
12			are obtained from the the 2-degree ECCO configuration.
13
14			Integration is initialized from annual-mean Levitus climatology
15			hydrogThetaFile = 'LevCli_temp.labsea'
16			hydrogSaltFile = 'LevCli_salt.labsea'
17
18			Surface salinity relaxation is to the monthly mean Levitus climatology
19			saltClimFile = 'SSS.labsea'
20
21			Forcing files are a 1979-1999 monthly climatology computed from the
22			NCEP reanalysis (see pkg/seaice/SEAICE_FFIELDS.h for units and signs)
23	dimitri	1.4	uwindFile = 'u10m.labsea79' # 10-m zonal wind
24			vwindFile = 'v10m.labsea79' # 10-m meridional wind
25			atempFile = 'tair.labsea1979' # 2-m air temperature
26			aqhFile = 'qa.labsea1979' # 2-m specific humidity
27	dimitri	1.5	lwdownFile = 'flo.labsea1979' # downward longwave radiation
28			swdownFile = 'fsh.labsea1979' # downward shortwave radiation
29	dimitri	1.4	precipFile = 'prate.labsea1979' # precipitation
30	dimitri	1.11
31			The experiment uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
32			The test is a 1-cpu, 10-hour integration. Both the atmospheric
33			state and the open-water surface fluxes are provided by pkg/exf.
34	heimbach	1.2
35
36	dimitri	1.3	Using testscript to test sea-ice code
37			=====================================
38
39			Running the testscript experiment:
40	dimitri	1.10	cd MITgcm/verification
41			./testreport -t lab_sea
42	dimitri	1.4	The default experiment is Experiment 8, below.
43			It uses pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf.
44	dimitri	1.3
45			Note that fairly large differences in accuracy occur across different
46			platforms. For example, testscript comparisons between g77 (Linux)
47			and f77 (SGI) generated output gives:
48
49	dimitri	1.4	T S U V
50			C D M c m s m s m s m s
51			n p a R g m m e . m m e . m m e . m m e .
52			f n k u 2 i a a d i a a d i a a d i a a d
53			g d e n d n x n . n x n . n x n . n x n .
54	dimitri	1.3
55	dimitri	1.4	Y Y Y Y 8 10 9 11 10 9 11 13 10 9 8 8 9 8 9 8 9 FAIL lab_sea