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