# K-H instability

param name kh
param nonhydrostatic true
param totalcount 6001
param saveinterval 40
param basetimestep 0.0005
param zdiffusionimplicit false
param iternonhydro 100
domain 33 4 67 0.0015 0.01 0.0009375 0.0 0.0 0.0
domain 33 4 67 0.0015 0.01 0.0009375 0.0435 0.0 0.0
domain 33 4 67 0.0015 0.01 0.0009375 0.087 0.0 0.0
domain 33 4 67 0.0015 0.01 0.0009375 0.1305 0.0 0.0

field 1 0 1 2 3
bcdomain top rigidslip
bcdomain down rigidslip
bcdomain east rigidslip
bcdomain west rigidslip
bcdomain north rigidslip
bcdomain south rigidslip

phys u uface true  true   10000000000.0 0.0
phys v vface false false  1.0 0.0
phys w wface true  true   10000000000.0 0.0
phys temp center false false 7.09 20.0
phys salinity center true true 10000000000.0 0.0
advect uq uq uw

connect west-east 0 1 0 0 false
connect west-east 1 2 0 0 false
connect west-east 2 3 0 0 false
connect west-east 3 0 0 0 false

depinit u velocity.inp 0
depinit salinity salinity.inp 0 

output xyz