Annotation of /MITgcm/verification/MLAdjust/README

Simple set-up to test flow-dependent horizontal viscosity implementation.

 Domain size is 50 x 26 x 40 grid-cells, 
 with uniform resolution dx=dy= 1.km , dz = 5.m
 
 Zonally re-entrant, flat bottom channel (closed by Northern Wall @ j=26)

 input files are real*8 (see matlab script input/genit.m )

 start from initial density field (given by initial Temp), no forcing

 Test Exp. |  Momentum   | useFullLeith | Strain-Tension | Biharmonic | side-drag   |
           | formulation |              |  formulation   | vs harmonic| (no_slip BC)|
-------------------------------------------------------------------------------------
 standard  | Vector-Inv. | FullLeith    |       No       |   viscC4   |    No       |
(dir=input)|             |              |                |            |             |
-------------------------------------------------------------------------------------
 .0.smag   | Vector-Inv. | FullLeith    |       No       |   viscC2   |    No       |
-------------------------------------------------------------------------------------
 .0.leith  | Vector-Inv. | FullLeith    |       Yes      |   viscC2   |    No       |
-------------------------------------------------------------------------------------
 .0.leithD | Vector-Inv. |    No        |       No       |   viscC2   |    No       |
-------------------------------------------------------------------------------------
 .1.leith  | Vector-Inv. | FullLeith    |       No       |   viscC4   |    No       |
-------------------------------------------------------------------------------------
 
Notes: 
1) Stain-Tension viscosity formulation is used when setting 
     useStrainTensionVisc=.TRUE.,
   default is .False., to use vorticity & divergence formulation
2) test 1.leith (input.1.leith) starts from a pickup-file, 
   all other test exp. start from iter0=0
3) test with Biharmonic visc. are using:
    viscC4leith = viscC4leithD = 1.85, 
    viscC4smag  = 1.,
    viscA4GridMax = 1.,
4) test with Harmonic visc. are using:
    viscC2leith = viscC2leithD=1.85, 
    viscC2smag  = 1.,
    viscAhGridMax = 1.,

-------------------------------------
1	jmc	1.1	Simple set-up to test flow-dependent horizontal viscosity implementation.
2
3			Domain size is 50 x 26 x 40 grid-cells,
4			with uniform resolution dx=dy= 1.km , dz = 5.m
5
6			Zonally re-entrant, flat bottom channel (closed by Northern Wall @ j=26)
7
8			input files are real*8 (see matlab script input/genit.m )
9
10			start from initial density field (given by initial Temp), no forcing
11
12			Test Exp. \| Momentum \| useFullLeith \| Strain-Tension \| Biharmonic \| side-drag \|
13			\| formulation \| \| formulation \| vs harmonic\| (no_slip BC)\|
14			-------------------------------------------------------------------------------------
15			standard \| Vector-Inv. \| FullLeith \| No \| viscC4 \| No \|
16			(dir=input)\| \| \| \| \| \|
17			-------------------------------------------------------------------------------------
18			.0.smag \| Vector-Inv. \| FullLeith \| No \| viscC2 \| No \|
19			-------------------------------------------------------------------------------------
20			.0.leith \| Vector-Inv. \| FullLeith \| Yes \| viscC2 \| No \|
21			-------------------------------------------------------------------------------------
22			.0.leithD \| Vector-Inv. \| No \| No \| viscC2 \| No \|
23			-------------------------------------------------------------------------------------
24			.1.leith \| Vector-Inv. \| FullLeith \| No \| viscC4 \| No \|
25			-------------------------------------------------------------------------------------
26
27			Notes:
28			1) Stain-Tension viscosity formulation is used when setting
29			useStrainTensionVisc=.TRUE.,
30			default is .False., to use vorticity & divergence formulation
31			2) test 1.leith (input.1.leith) starts from a pickup-file,
32			all other test exp. start from iter0=0
33			3) test with Biharmonic visc. are using:
34			viscC4leith = viscC4leithD = 1.85,
35			viscC4smag = 1.,
36			viscA4GridMax = 1.,
37			4) test with Harmonic visc. are using:
38			viscC2leith = viscC2leithD=1.85,
39			viscC2smag = 1.,
40			viscAhGridMax = 1.,
41
42			-------------------------------------