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/gendata.m )

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

 Test Exp. |  Momentum   |  Viscosity  | useFullLeith | Biharmonic | side-drag   |
           | formulation | formulation |              | vs harmonic| (no_slip BC)|
----------------------------------------------------------------------------------
 standard  | Vector-Inv. |  Vort-Div   |  FullLeith   |   viscC4   |    No       |
(dir=input)|             |             |              |            |             |
----------------------------------------------------------------------------------
 .A4FlxF   |  Flux-Form  |  FLux-Form  |  FullLeith   |   viscC4   |    Yes      |
----------------------------------------------------------------------------------
 .AhFlxF   |  Flux-Form  |  Flux-Form  |    No        |   viscC2   |    No       |
----------------------------------------------------------------------------------
 .AhVrDv   | Vector-Inv. |  Vort-Div   |  FullLeith   |   viscC2   |    Yes      |
----------------------------------------------------------------------------------
 .AhStTn   | Vector-Inv. | Strain-Tens |  FullLeith   |   viscC2   |    Yes      |
----------------------------------------------------------------------------------
 
Notes: 
1) Stain-Tension viscosity formulation is used when setting 
     useStrainTensionVisc=.TRUE.,
   and currently only available with Vector-Invariant momentum.
   Default is .False., to use Vorticity & Divergence formulation
   (pkg/mom_vecinv) or Flux-Form formulation (pkg/mom_fluxform)
2) test .A4FlxF (input.A4FlxF) 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	jmc	1.2	input files are real*8 (see matlab script input/gendata.m )
9	jmc	1.1
10			start from initial density field (given by initial Temp), no forcing
11
12	jmc	1.2	Test Exp. \| Momentum \| Viscosity \| useFullLeith \| Biharmonic \| side-drag \|
13			\| formulation \| formulation \| \| vs harmonic\| (no_slip BC)\|
14			----------------------------------------------------------------------------------
15			standard \| Vector-Inv. \| Vort-Div \| FullLeith \| viscC4 \| No \|
16			(dir=input)\| \| \| \| \| \|
17			----------------------------------------------------------------------------------
18			.A4FlxF \| Flux-Form \| FLux-Form \| FullLeith \| viscC4 \| Yes \|
19			----------------------------------------------------------------------------------
20			.AhFlxF \| Flux-Form \| Flux-Form \| No \| viscC2 \| No \|
21			----------------------------------------------------------------------------------
22			.AhVrDv \| Vector-Inv. \| Vort-Div \| FullLeith \| viscC2 \| Yes \|
23			----------------------------------------------------------------------------------
24			.AhStTn \| Vector-Inv. \| Strain-Tens \| FullLeith \| viscC2 \| Yes \|
25			----------------------------------------------------------------------------------
26	jmc	1.1
27			Notes:
28			1) Stain-Tension viscosity formulation is used when setting
29			useStrainTensionVisc=.TRUE.,
30	jmc	1.2	and currently only available with Vector-Invariant momentum.
31			Default is .False., to use Vorticity & Divergence formulation
32			(pkg/mom_vecinv) or Flux-Form formulation (pkg/mom_fluxform)
33			2) test .A4FlxF (input.A4FlxF) starts from a pickup-file,
34	jmc	1.1	all other test exp. start from iter0=0
35			3) test with Biharmonic visc. are using:
36			viscC4leith = viscC4leithD = 1.85,
37			viscC4smag = 1.,
38			viscA4GridMax = 1.,
39			4) test with Harmonic visc. are using:
40			viscC2leith = viscC2leithD=1.85,
41			viscC2smag = 1.,
42			viscAhGridMax = 1.,
43
44			-------------------------------------