Contents 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.

More pkg/seaice test experiments, configured for low and
high-resolution global cube-sphere domains are described
in MITgcm_contrib/high_res_cube/README_ice.


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

Running the testscript experiment:
  cd MITgcm/verification
  ./testreport -t lab_sea

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


Instructions for generating and running a 1-CPU experiment
==========================================================

Generating 1-CPU executable:
  cd MITgcm/verification/lab_sea
  mkdir build
  cd build
  cp ../code/*.h ../code/packages.conf .
  ../../../tools/genmake2
  make depend
  make

To run Experiment 1:
  cd ../input
  ../build/mitgcmuv > output.txt

Use matlab script lookat_exp1.m to compare the output
of exp1 with that from checkpoint51f sea-ice code:
  cd ../../../verification/lab_sea/matlab
  matlab
  lookat_ice
1	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	A brief description of the sea-ice model is in "seaice.ps".
7
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	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	lwdownFile = 'flo.labsea1979' # downward longwave radiation
28	swdownFile = 'fsh.labsea1979' # downward shortwave radiation
29	precipFile = 'prate.labsea1979' # precipitation
30
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
35	More pkg/seaice test experiments, configured for low and
36	high-resolution global cube-sphere domains are described
37	in MITgcm_contrib/high_res_cube/README_ice.
38
39
40	Using testscript to test sea-ice code
41	=====================================
42
43	Running the testscript experiment:
44	cd MITgcm/verification
45	./testreport -t lab_sea
46
47	Note that fairly large differences in accuracy occur across different
48	platforms. For example, testscript comparisons between g77 (Linux)
49	and f77 (SGI) generated output gives:
50
51	T S U V
52	C D M c m s m s m s m s
53	n p a R g m m e . m m e . m m e . m m e .
54	f n k u 2 i a a d i a a d i a a d i a a d
55	g d e n d n x n . n x n . n x n . n x n .
56
57	Y Y Y Y 8 10 9 11 10 9 11 13 10 9 8 8 9 8 9 8 9 FAIL lab_sea
58
59
60	Instructions for generating and running a 1-CPU experiment
61	==========================================================
62
63	Generating 1-CPU executable:
64	cd MITgcm/verification/lab_sea
65	mkdir build
66	cd build
67	cp ../code/*.h ../code/packages.conf .
68	../../../tools/genmake2
69	make depend
70	make
71
72	To run Experiment 1:
73	cd ../input
74	../build/mitgcmuv > output.txt
75
76	Use matlab script lookat_exp1.m to compare the output
77	of exp1 with that from checkpoint51f sea-ice code:
78	cd ../../../verification/lab_sea/matlab
79	matlab
80	lookat_ice