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 1-CPU and 2-CPU executables
=======================================================

Generating 1-CPU executable:
  cd MITgcm/verification/lab_sea
  mkdir build
  cd build
  ../../../tools/genmake2
  make clean
  make depend
  make
  mv mitgcmuv mitgcmuv_1

Generating 2-CPU executable:
  cd ../code
  \mv SIZE.h SIZE.h_1
  \mv SIZE.h_mpi SIZE.h
  ../../../tools/genmake2 -mpi
  make clean
  make depend
  make
  mv mitgcmuv mitgcmuv_2x1


Instructions for running Experiment 1
=====================================

This is a 1-cpu, 10-hour integration used to make sure that all
the files are available and that the sea-ice model compiles and
runs.  The reading of atmospheric forcing files and the computation
of open-water bulk fluxes are carried out within pkg/seaice.

To run Experiment 1:
  cd ../input
  ./mitgcmuv_1 >&! 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


Instructions for running Experiment 2
=====================================

This is a 2-cpu, 10-hour integration.

To run Experiment 2:
  cd ../../../verification/lab_sea/input
  mpirun -np 2 mitgcmuv_2x1
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 1-CPU and 2-CPU executables
61	=======================================================
62
63	Generating 1-CPU executable:
64	cd MITgcm/verification/lab_sea
65	mkdir build
66	cd build
67	../../../tools/genmake2
68	make clean
69	make depend
70	make
71	mv mitgcmuv mitgcmuv_1
72
73	Generating 2-CPU executable:
74	cd ../code
75	\mv SIZE.h SIZE.h_1
76	\mv SIZE.h_mpi SIZE.h
77	../../../tools/genmake2 -mpi
78	make clean
79	make depend
80	make
81	mv mitgcmuv mitgcmuv_2x1
82
83
84	Instructions for running Experiment 1
85	=====================================
86
87	This is a 1-cpu, 10-hour integration used to make sure that all
88	the files are available and that the sea-ice model compiles and
89	runs. The reading of atmospheric forcing files and the computation
90	of open-water bulk fluxes are carried out within pkg/seaice.
91
92	To run Experiment 1:
93	cd ../input
94	./mitgcmuv_1 >&! output.txt
95
96	Use matlab script lookat_exp1.m to compare the output
97	of exp1 with that from checkpoint51f sea-ice code:
98	cd ../../../verification/lab_sea/matlab
99	matlab
100	lookat_ice
101
102
103	Instructions for running Experiment 2
104	=====================================
105
106	This is a 2-cpu, 10-hour integration.
107
108	To run Experiment 2:
109	cd ../../../verification/lab_sea/input
110	mpirun -np 2 mitgcmuv_2x1