TALK=T;RUN(1,1)
  DISPLAY

  This Kelvin-Helmholz simulation can be run in the following
  modes:
                                  'pure k-e'  or     'pure MFM'
                                  ----------         ----------
  1. Effective viscosity from  (a) k**2 / e   or (b) mixl * mnsq
  2. Micro-mixing rate   from  (a) e / k      or (b) mnsq / mixl
  3. Micro-mixing length from  (e) k**1.5 / e or (b) source =
                                                     const * mnsq

  Eight combinations are obviously possible. One or other of the
  'pure' cases can be set by way of a single boolean variable.

  This Q1 comprises that created by Dr Hornby for the k-epsilon
  model, with addition of MFM features.

  It requires the use of a special ground (see appendix D).

  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(K-EP:KH instability                     )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =1.5
 FSTEP = 1 ;LSTEP = 600
 TFRAC(1)=1.666667E-03 ;TFRAC(2)=3.333333E-03
 TFRAC(3)=5.0E-03 ;TFRAC(4)=6.666667E-03
 TFRAC(5)=8.333335E-03 ;TFRAC(6)=1.0E-02
 TFRAC(7)=0.011667 ;TFRAC(8)=0.013333
 TFRAC(9)=0.015 ;TFRAC(10)=0.016667
 TFRAC(11)=0.018333 ;TFRAC(12)=0.02
 TFRAC(13)=0.021667 ;TFRAC(14)=0.023333
 TFRAC(15)=0.025 ;TFRAC(16)=0.026667
 TFRAC(17)=0.028333 ;TFRAC(18)=0.03
 TFRAC(19)=0.031667 ;TFRAC(20)=0.033333
 TFRAC(21)=0.035 ;TFRAC(22)=0.036667
 TFRAC(23)=0.038333 ;TFRAC(24)=0.04
 TFRAC(25)=0.041667 ;TFRAC(26)=0.043333
 TFRAC(27)=0.045 ;TFRAC(28)=0.046667
 TFRAC(29)=0.048333 ;TFRAC(30)=0.05
 TFRAC(31)=0.051667 ;TFRAC(32)=0.053333
 TFRAC(33)=0.055 ;TFRAC(34)=0.056667
 TFRAC(35)=0.058333 ;TFRAC(36)=0.06
 TFRAC(37)=0.061667 ;TFRAC(38)=0.063333
 TFRAC(39)=0.065 ;TFRAC(40)=0.066667
 TFRAC(41)=0.068333 ;TFRAC(42)=0.07
 TFRAC(43)=0.071667 ;TFRAC(44)=0.073333
 TFRAC(45)=0.075 ;TFRAC(46)=0.076667
 TFRAC(47)=0.078333 ;TFRAC(48)=0.08
 TFRAC(49)=0.081667 ;TFRAC(50)=0.083333
 TFRAC(51)=0.085 ;TFRAC(52)=0.086667
 TFRAC(53)=0.088333 ;TFRAC(54)=0.09
 TFRAC(55)=0.091667 ;TFRAC(56)=0.093333
 TFRAC(57)=0.095 ;TFRAC(58)=0.096667
 TFRAC(59)=0.098333 ;TFRAC(60)=0.1
 TFRAC(61)=0.101667 ;TFRAC(62)=0.103333
 TFRAC(63)=0.105 ;TFRAC(64)=0.106667
 TFRAC(65)=0.108333 ;TFRAC(66)=0.11
 TFRAC(67)=0.111667 ;TFRAC(68)=0.113333
 TFRAC(69)=0.115 ;TFRAC(70)=0.116667
 TFRAC(71)=0.118333 ;TFRAC(72)=0.12
 TFRAC(73)=0.121667 ;TFRAC(74)=0.123333
 TFRAC(75)=0.125 ;TFRAC(76)=0.126667
 TFRAC(77)=0.128333 ;TFRAC(78)=0.13
 TFRAC(79)=0.131667 ;TFRAC(80)=0.133333
 TFRAC(81)=0.135 ;TFRAC(82)=0.136667
 TFRAC(83)=0.138333 ;TFRAC(84)=0.14
 TFRAC(85)=0.141667 ;TFRAC(86)=0.143333
 TFRAC(87)=0.145 ;TFRAC(88)=0.146667
 TFRAC(89)=0.148333 ;TFRAC(90)=0.15
 TFRAC(91)=0.151667 ;TFRAC(92)=0.153333
 TFRAC(93)=0.155 ;TFRAC(94)=0.156667
 TFRAC(95)=0.158333 ;TFRAC(96)=0.16
 TFRAC(97)=0.161667 ;TFRAC(98)=0.163333
 TFRAC(99)=0.165 ;TFRAC(100)=0.166667
 TFRAC(101)=0.168333 ;TFRAC(102)=0.17
 TFRAC(103)=0.171667 ;TFRAC(104)=0.173333
 TFRAC(105)=0.175 ;TFRAC(106)=0.176667
 TFRAC(107)=0.178333 ;TFRAC(108)=0.18
 TFRAC(109)=0.181667 ;TFRAC(110)=0.183333
 TFRAC(111)=0.185 ;TFRAC(112)=0.186667
 TFRAC(113)=0.188333 ;TFRAC(114)=0.19
 TFRAC(115)=0.191667 ;TFRAC(116)=0.193333
 TFRAC(117)=0.195 ;TFRAC(118)=0.196667
 TFRAC(119)=0.198333 ;TFRAC(120)=0.2
 TFRAC(121)=0.201667 ;TFRAC(122)=0.203333
 TFRAC(123)=0.205 ;TFRAC(124)=0.206667
 TFRAC(125)=0.208333 ;TFRAC(126)=0.21
 TFRAC(127)=0.211667 ;TFRAC(128)=0.213333
 TFRAC(129)=0.215 ;TFRAC(130)=0.216667
 TFRAC(131)=0.218333 ;TFRAC(132)=0.22
 TFRAC(133)=0.221667 ;TFRAC(134)=0.223333
 TFRAC(135)=0.225 ;TFRAC(136)=0.226667
 TFRAC(137)=0.228333 ;TFRAC(138)=0.23
 TFRAC(139)=0.231667 ;TFRAC(140)=0.233333
 TFRAC(141)=0.235 ;TFRAC(142)=0.236667
 TFRAC(143)=0.238333 ;TFRAC(144)=0.24
 TFRAC(145)=0.241667 ;TFRAC(146)=0.243333
 TFRAC(147)=0.245 ;TFRAC(148)=0.246667
 TFRAC(149)=0.248333 ;TFRAC(150)=0.25
 TFRAC(151)=0.251667 ;TFRAC(152)=0.253333
 TFRAC(153)=0.255 ;TFRAC(154)=0.256667
 TFRAC(155)=0.258333 ;TFRAC(156)=0.26
 TFRAC(157)=0.261667 ;TFRAC(158)=0.263333
 TFRAC(159)=0.265 ;TFRAC(160)=0.266667
 TFRAC(161)=0.268333 ;TFRAC(162)=0.27
 TFRAC(163)=0.271667 ;TFRAC(164)=0.273333
 TFRAC(165)=0.275 ;TFRAC(166)=0.276667
 TFRAC(167)=0.278333 ;TFRAC(168)=0.28
 TFRAC(169)=0.281667 ;TFRAC(170)=0.283333
 TFRAC(171)=0.285 ;TFRAC(172)=0.286667
 TFRAC(173)=0.288333 ;TFRAC(174)=0.29
 TFRAC(175)=0.291667 ;TFRAC(176)=0.293333
 TFRAC(177)=0.295 ;TFRAC(178)=0.296667
 TFRAC(179)=0.298333 ;TFRAC(180)=0.3
 TFRAC(181)=0.301667 ;TFRAC(182)=0.303333
 TFRAC(183)=0.305 ;TFRAC(184)=0.306667
 TFRAC(185)=0.308333 ;TFRAC(186)=0.31
 TFRAC(187)=0.311667 ;TFRAC(188)=0.313333
 TFRAC(189)=0.315 ;TFRAC(190)=0.316667
 TFRAC(191)=0.318333 ;TFRAC(192)=0.32
 TFRAC(193)=0.321667 ;TFRAC(194)=0.323333
 TFRAC(195)=0.325 ;TFRAC(196)=0.326667
 TFRAC(197)=0.328333 ;TFRAC(198)=0.33
 TFRAC(199)=0.331667 ;TFRAC(200)=0.333333
 TFRAC(201)=0.335 ;TFRAC(202)=0.336667
 TFRAC(203)=0.338333 ;TFRAC(204)=0.34
 TFRAC(205)=0.341667 ;TFRAC(206)=0.343333
 TFRAC(207)=0.345 ;TFRAC(208)=0.346667
 TFRAC(209)=0.348333 ;TFRAC(210)=0.35
 TFRAC(211)=0.351667 ;TFRAC(212)=0.353333
 TFRAC(213)=0.355 ;TFRAC(214)=0.356667
 TFRAC(215)=0.358333 ;TFRAC(216)=0.36
 TFRAC(217)=0.361667 ;TFRAC(218)=0.363333
 TFRAC(219)=0.365 ;TFRAC(220)=0.366667
 TFRAC(221)=0.368333 ;TFRAC(222)=0.37
 TFRAC(223)=0.371667 ;TFRAC(224)=0.373333
 TFRAC(225)=0.375 ;TFRAC(226)=0.376667
 TFRAC(227)=0.378333 ;TFRAC(228)=0.38
 TFRAC(229)=0.381667 ;TFRAC(230)=0.383333
 TFRAC(231)=0.385 ;TFRAC(232)=0.386667
 TFRAC(233)=0.388333 ;TFRAC(234)=0.39
 TFRAC(235)=0.391667 ;TFRAC(236)=0.393333
 TFRAC(237)=0.395 ;TFRAC(238)=0.396667
 TFRAC(239)=0.398333 ;TFRAC(240)=0.4
 TFRAC(241)=0.401667 ;TFRAC(242)=0.403333
 TFRAC(243)=0.405 ;TFRAC(244)=0.406667
 TFRAC(245)=0.408333 ;TFRAC(246)=0.41
 TFRAC(247)=0.411667 ;TFRAC(248)=0.413333
 TFRAC(249)=0.415 ;TFRAC(250)=0.416667
 TFRAC(251)=0.418333 ;TFRAC(252)=0.42
 TFRAC(253)=0.421667 ;TFRAC(254)=0.423333
 TFRAC(255)=0.425 ;TFRAC(256)=0.426667
 TFRAC(257)=0.428333 ;TFRAC(258)=0.43
 TFRAC(259)=0.431667 ;TFRAC(260)=0.433333
 TFRAC(261)=0.435 ;TFRAC(262)=0.436667
 TFRAC(263)=0.438333 ;TFRAC(264)=0.44
 TFRAC(265)=0.441667 ;TFRAC(266)=0.443333
 TFRAC(267)=0.445 ;TFRAC(268)=0.446667
 TFRAC(269)=0.448333 ;TFRAC(270)=0.45
 TFRAC(271)=0.451667 ;TFRAC(272)=0.453333
 TFRAC(273)=0.455 ;TFRAC(274)=0.456667
 TFRAC(275)=0.458333 ;TFRAC(276)=0.46
 TFRAC(277)=0.461667 ;TFRAC(278)=0.463333
 TFRAC(279)=0.465 ;TFRAC(280)=0.466667
 TFRAC(281)=0.468333 ;TFRAC(282)=0.47
 TFRAC(283)=0.471667 ;TFRAC(284)=0.473333
 TFRAC(285)=0.475 ;TFRAC(286)=0.476667
 TFRAC(287)=0.478333 ;TFRAC(288)=0.48
 TFRAC(289)=0.481667 ;TFRAC(290)=0.483333
 TFRAC(291)=0.485 ;TFRAC(292)=0.486667
 TFRAC(293)=0.488333 ;TFRAC(294)=0.49
 TFRAC(295)=0.491667 ;TFRAC(296)=0.493333
 TFRAC(297)=0.495 ;TFRAC(298)=0.496667
 TFRAC(299)=0.498333 ;TFRAC(300)=0.5
 TFRAC(301)=0.501667 ;TFRAC(302)=0.503333
 TFRAC(303)=0.505 ;TFRAC(304)=0.506667
 TFRAC(305)=0.508333 ;TFRAC(306)=0.51
 TFRAC(307)=0.511667 ;TFRAC(308)=0.513333
 TFRAC(309)=0.515 ;TFRAC(310)=0.516667
 TFRAC(311)=0.518333 ;TFRAC(312)=0.52
 TFRAC(313)=0.521667 ;TFRAC(314)=0.523333
 TFRAC(315)=0.525 ;TFRAC(316)=0.526667
 TFRAC(317)=0.528333 ;TFRAC(318)=0.53
 TFRAC(319)=0.531667 ;TFRAC(320)=0.533333
 TFRAC(321)=0.535 ;TFRAC(322)=0.536667
 TFRAC(323)=0.538333 ;TFRAC(324)=0.54
 TFRAC(325)=0.541667 ;TFRAC(326)=0.543333
 TFRAC(327)=0.545 ;TFRAC(328)=0.546667
 TFRAC(329)=0.548333 ;TFRAC(330)=0.55
 TFRAC(331)=0.551667 ;TFRAC(332)=0.553333
 TFRAC(333)=0.555 ;TFRAC(334)=0.556667
 TFRAC(335)=0.558333 ;TFRAC(336)=0.56
 TFRAC(337)=0.561667 ;TFRAC(338)=0.563333
 TFRAC(339)=0.565 ;TFRAC(340)=0.566667
 TFRAC(341)=0.568333 ;TFRAC(342)=0.57
 TFRAC(343)=0.571667 ;TFRAC(344)=0.573333
 TFRAC(345)=0.575 ;TFRAC(346)=0.576667
 TFRAC(347)=0.578333 ;TFRAC(348)=0.58
 TFRAC(349)=0.581667 ;TFRAC(350)=0.583333
 TFRAC(351)=0.585 ;TFRAC(352)=0.586667
 TFRAC(353)=0.588333 ;TFRAC(354)=0.59
 TFRAC(355)=0.591667 ;TFRAC(356)=0.593333
 TFRAC(357)=0.595 ;TFRAC(358)=0.596667
 TFRAC(359)=0.598333 ;TFRAC(360)=0.6
 TFRAC(361)=0.601667 ;TFRAC(362)=0.603333
 TFRAC(363)=0.605 ;TFRAC(364)=0.606667
 TFRAC(365)=0.608333 ;TFRAC(366)=0.61
 TFRAC(367)=0.611667 ;TFRAC(368)=0.613333
 TFRAC(369)=0.615 ;TFRAC(370)=0.616667
 TFRAC(371)=0.618333 ;TFRAC(372)=0.62
 TFRAC(373)=0.621667 ;TFRAC(374)=0.623333
 TFRAC(375)=0.625 ;TFRAC(376)=0.626667
 TFRAC(377)=0.628333 ;TFRAC(378)=0.63
 TFRAC(379)=0.631667 ;TFRAC(380)=0.633333
 TFRAC(381)=0.635 ;TFRAC(382)=0.636667
 TFRAC(383)=0.638333 ;TFRAC(384)=0.64
 TFRAC(385)=0.641667 ;TFRAC(386)=0.643333
 TFRAC(387)=0.645 ;TFRAC(388)=0.646667
 TFRAC(389)=0.648333 ;TFRAC(390)=0.65
 TFRAC(391)=0.651667 ;TFRAC(392)=0.653333
 TFRAC(393)=0.655 ;TFRAC(394)=0.656667
 TFRAC(395)=0.658333 ;TFRAC(396)=0.66
 TFRAC(397)=0.661667 ;TFRAC(398)=0.663333
 TFRAC(399)=0.665 ;TFRAC(400)=0.666667
 TFRAC(401)=0.668333 ;TFRAC(402)=0.67
 TFRAC(403)=0.671667 ;TFRAC(404)=0.673333
 TFRAC(405)=0.675 ;TFRAC(406)=0.676667
 TFRAC(407)=0.678333 ;TFRAC(408)=0.68
 TFRAC(409)=0.681667 ;TFRAC(410)=0.683333
 TFRAC(411)=0.685 ;TFRAC(412)=0.686667
 TFRAC(413)=0.688333 ;TFRAC(414)=0.69
 TFRAC(415)=0.691667 ;TFRAC(416)=0.693333
 TFRAC(417)=0.695 ;TFRAC(418)=0.696667
 TFRAC(419)=0.698333 ;TFRAC(420)=0.7
 TFRAC(421)=0.701667 ;TFRAC(422)=0.703333
 TFRAC(423)=0.705 ;TFRAC(424)=0.706667
 TFRAC(425)=0.708333 ;TFRAC(426)=0.71
 TFRAC(427)=0.711667 ;TFRAC(428)=0.713333
 TFRAC(429)=0.715 ;TFRAC(430)=0.716667
 TFRAC(431)=0.718333 ;TFRAC(432)=0.72
 TFRAC(433)=0.721667 ;TFRAC(434)=0.723333
 TFRAC(435)=0.725 ;TFRAC(436)=0.726667
 TFRAC(437)=0.728333 ;TFRAC(438)=0.73
 TFRAC(439)=0.731667 ;TFRAC(440)=0.733333
 TFRAC(441)=0.735 ;TFRAC(442)=0.736667
 TFRAC(443)=0.738333 ;TFRAC(444)=0.74
 TFRAC(445)=0.741667 ;TFRAC(446)=0.743333
 TFRAC(447)=0.745 ;TFRAC(448)=0.746667
 TFRAC(449)=0.748333 ;TFRAC(450)=0.75
 TFRAC(451)=0.751667 ;TFRAC(452)=0.753333
 TFRAC(453)=0.755 ;TFRAC(454)=0.756667
 TFRAC(455)=0.758333 ;TFRAC(456)=0.76
 TFRAC(457)=0.761667 ;TFRAC(458)=0.763333
 TFRAC(459)=0.765 ;TFRAC(460)=0.766667
 TFRAC(461)=0.768333 ;TFRAC(462)=0.77
 TFRAC(463)=0.771667 ;TFRAC(464)=0.773333
 TFRAC(465)=0.775 ;TFRAC(466)=0.776667
 TFRAC(467)=0.778333 ;TFRAC(468)=0.78
 TFRAC(469)=0.781667 ;TFRAC(470)=0.783333
 TFRAC(471)=0.785 ;TFRAC(472)=0.786667
 TFRAC(473)=0.788333 ;TFRAC(474)=0.79
 TFRAC(475)=0.791667 ;TFRAC(476)=0.793333
 TFRAC(477)=0.795 ;TFRAC(478)=0.796667
 TFRAC(479)=0.798333 ;TFRAC(480)=0.8
 TFRAC(481)=0.801667 ;TFRAC(482)=0.803333
 TFRAC(483)=0.805 ;TFRAC(484)=0.806667
 TFRAC(485)=0.808333 ;TFRAC(486)=0.81
 TFRAC(487)=0.811667 ;TFRAC(488)=0.813333
 TFRAC(489)=0.815 ;TFRAC(490)=0.816667
 TFRAC(491)=0.818333 ;TFRAC(492)=0.82
 TFRAC(493)=0.821667 ;TFRAC(494)=0.823333
 TFRAC(495)=0.825 ;TFRAC(496)=0.826667
 TFRAC(497)=0.828333 ;TFRAC(498)=0.83
 TFRAC(499)=0.831667 ;TFRAC(500)=0.833333
 TFRAC(501)=0.835 ;TFRAC(502)=0.836667
 TFRAC(503)=0.838333 ;TFRAC(504)=0.84
 TFRAC(505)=0.841667 ;TFRAC(506)=0.843333
 TFRAC(507)=0.845 ;TFRAC(508)=0.846667
 TFRAC(509)=0.848333 ;TFRAC(510)=0.85
 TFRAC(511)=0.851667 ;TFRAC(512)=0.853333
 TFRAC(513)=0.855 ;TFRAC(514)=0.856667
 TFRAC(515)=0.858333 ;TFRAC(516)=0.86
 TFRAC(517)=0.861667 ;TFRAC(518)=0.863333
 TFRAC(519)=0.865 ;TFRAC(520)=0.866667
 TFRAC(521)=0.868333 ;TFRAC(522)=0.87
 TFRAC(523)=0.871667 ;TFRAC(524)=0.873333
 TFRAC(525)=0.875 ;TFRAC(526)=0.876667
 TFRAC(527)=0.878333 ;TFRAC(528)=0.88
 TFRAC(529)=0.881667 ;TFRAC(530)=0.883333
 TFRAC(531)=0.885 ;TFRAC(532)=0.886667
 TFRAC(533)=0.888333 ;TFRAC(534)=0.89
 TFRAC(535)=0.891667 ;TFRAC(536)=0.893333
 TFRAC(537)=0.895 ;TFRAC(538)=0.896667
 TFRAC(539)=0.898333 ;TFRAC(540)=0.9
 TFRAC(541)=0.901667 ;TFRAC(542)=0.903333
 TFRAC(543)=0.905 ;TFRAC(544)=0.906667
 TFRAC(545)=0.908333 ;TFRAC(546)=0.91
 TFRAC(547)=0.911667 ;TFRAC(548)=0.913333
 TFRAC(549)=0.915 ;TFRAC(550)=0.916667
 TFRAC(551)=0.918333 ;TFRAC(552)=0.92
 TFRAC(553)=0.921667 ;TFRAC(554)=0.923333
 TFRAC(555)=0.925 ;TFRAC(556)=0.926667
 TFRAC(557)=0.928333 ;TFRAC(558)=0.93
 TFRAC(559)=0.931667 ;TFRAC(560)=0.933333
 TFRAC(561)=0.935 ;TFRAC(562)=0.936667
 TFRAC(563)=0.938333 ;TFRAC(564)=0.94
 TFRAC(565)=0.941667 ;TFRAC(566)=0.943333
 TFRAC(567)=0.945 ;TFRAC(568)=0.946667
 TFRAC(569)=0.948333 ;TFRAC(570)=0.95
 TFRAC(571)=0.951667 ;TFRAC(572)=0.953333
 TFRAC(573)=0.955 ;TFRAC(574)=0.956667
 TFRAC(575)=0.958333 ;TFRAC(576)=0.96
 TFRAC(577)=0.961667 ;TFRAC(578)=0.963333
 TFRAC(579)=0.965 ;TFRAC(580)=0.966667
 TFRAC(581)=0.968333 ;TFRAC(582)=0.97
 TFRAC(583)=0.971667 ;TFRAC(584)=0.973333
 TFRAC(585)=0.975 ;TFRAC(586)=0.976667
 TFRAC(587)=0.978333 ;TFRAC(588)=0.98
 TFRAC(589)=0.981667 ;TFRAC(590)=0.983333
 TFRAC(591)=0.985 ;TFRAC(592)=0.986667
 TFRAC(593)=0.988333 ;TFRAC(594)=0.99
 TFRAC(595)=0.991667 ;TFRAC(596)=0.993333
 TFRAC(597)=0.995 ;TFRAC(598)=0.996667
 TFRAC(599)=0.998333 ;TFRAC(600)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 20
 XULAST =0.16
 XFRAC(1)=0.05 ;XFRAC(2)=0.1
 XFRAC(3)=0.15 ;XFRAC(4)=0.2
 XFRAC(5)=0.25 ;XFRAC(6)=0.3
 XFRAC(7)=0.35 ;XFRAC(8)=0.4
 XFRAC(9)=0.45 ;XFRAC(10)=0.5
 XFRAC(11)=0.55 ;XFRAC(12)=0.6
 XFRAC(13)=0.65 ;XFRAC(14)=0.7
 XFRAC(15)=0.75 ;XFRAC(16)=0.8
 XFRAC(17)=0.85 ;XFRAC(18)=0.9
 XFRAC(19)=0.95 ;XFRAC(20)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 30
 YVLAST =0.1
 YFRAC(1)=0.039728 ;YFRAC(2)=0.078893
 YFRAC(3)=0.117459 ;YFRAC(4)=0.155387
 YFRAC(5)=0.192631 ;YFRAC(6)=0.229136
 YFRAC(7)=0.264837 ;YFRAC(8)=0.299655
 YFRAC(9)=0.333489 ;YFRAC(10)=0.36621
 YFRAC(11)=0.397639 ;YFRAC(12)=0.427522
 YFRAC(13)=0.455445 ;YFRAC(14)=0.480606
 YFRAC(15)=0.5 ;YFRAC(16)=0.519394
 YFRAC(17)=0.544555 ;YFRAC(18)=0.572478
 YFRAC(19)=0.60236 ;YFRAC(20)=0.63379
 YFRAC(21)=0.666511 ;YFRAC(22)=0.700345
 YFRAC(23)=0.735163 ;YFRAC(24)=0.770864
 YFRAC(25)=0.807369 ;YFRAC(26)=0.844613
 YFRAC(27)=0.882541 ;YFRAC(28)=0.921107
 YFRAC(29)=0.960272 ;YFRAC(30)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
    * X-cyclic boundaries switched
 XCYIZ( 1, T )
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(14)=H1
 NAME(16)=SALT ;NAME(17)=C2
 NAME(132)=F1 ;NAME(133)=F2
 NAME(134)=F3 ;NAME(135)=F4
 NAME(136)=F5 ;NAME(137)=F6
 NAME(138)=F7 ;NAME(139)=F8
 NAME(140)=F9 ;NAME(141)=F10
 NAME(142)=AVEF ;NAME(143)=MNSQ
 NAME(144)=RATE ;NAME(145)=LEN1
 NAME(146)=MIXL ;NAME(147)=ENUT
 NAME(148)=EPKE ;NAME(149)=DEN1
 NAME(150)=TMP1
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(SALT,Y,Y,N,N,N,Y)
 SOLUTN(C2,Y,Y,N,N,N,Y)
 SOLUTN(F1,Y,Y,N,N,N,Y)
 SOLUTN(F2,Y,Y,N,N,N,Y)
 SOLUTN(F3,Y,Y,N,N,N,Y)
 SOLUTN(F4,Y,Y,N,N,N,Y)
 SOLUTN(F5,Y,Y,N,N,N,Y)
 SOLUTN(F6,Y,Y,N,N,N,Y)
 SOLUTN(F7,Y,Y,N,N,N,Y)
 SOLUTN(F8,Y,Y,N,N,N,Y)
 SOLUTN(F9,Y,Y,N,N,N,Y)
 SOLUTN(F10,Y,Y,N,N,N,Y)
 SOLUTN(AVEF,Y,N,N,N,N,Y)
 SOLUTN(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(MIXL,Y,Y,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 DEN1 = 149
 VIST = 147
 TEMP1 = 150
 LEN1 = 145
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(H1,N,Y,Y,Y,Y,Y)
 TERMS(SALT,N,Y,Y,Y,Y,Y)
 TERMS(C2,N,Y,Y,Y,N,Y)
 TERMS(F1,N,Y,Y,Y,Y,Y)
 TERMS(F2,N,Y,Y,Y,N,Y)
 TERMS(F3,N,Y,Y,Y,Y,Y)
 TERMS(F4,N,Y,Y,Y,N,Y)
 TERMS(F5,N,Y,Y,Y,Y,Y)
 TERMS(F6,N,Y,Y,Y,N,Y)
 TERMS(F7,N,Y,Y,Y,Y,Y)
 TERMS(F8,N,Y,Y,Y,N,Y)
 TERMS(F9,N,Y,Y,Y,Y,Y)
 TERMS(F10,N,Y,Y,Y,N,Y)
 TERMS(MIXL,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENT = T
 HUNIT =1.
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND1 ;TMP1 = GRND2 ;EL1 = GRND4
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =1086.97583 ;DRH1DP =0.
 RHO1A =1086.97583 ;RHO1B =-7.674451E-05 ;RHO1C =0.818824
 TMP1A =11. ;TMP1B =2.392344E-04 ;TMP1C =0.
 TMP2A =0. ;TMP2B =0. ;TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL =1.0E-06 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 0
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(H1)=7.572598 ;PRNDTL(SALT)=1000.
 PRNDTL(C2)=1. ;PRNDTL(F1)=1.
 PRNDTL(F2)=1. ;PRNDTL(F3)=1.
 PRNDTL(F4)=1. ;PRNDTL(F5)=1.
 PRNDTL(F6)=1. ;PRNDTL(F7)=1.
 PRNDTL(F8)=1. ;PRNDTL(F9)=1.
 PRNDTL(F10)=1. ;PRNDTL(MIXL)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(H1)=0.7 ;PRT(SALT)=0.7
 PRT(C2)=1. ;PRT(F1)=1.
 PRT(F2)=1. ;PRT(F3)=1.
 PRT(F4)=1. ;PRT(F5)=1.
 PRT(F6)=1. ;PRT(F7)=1.
 PRT(F8)=1. ;PRT(F9)=1.
 PRT(F10)=1. ;PRT(MIXL)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(V1)=1.0E-10 ;FIINIT(KE)=0.
 FIINIT(EP)=0. ;FIINIT(H1)=1.0E-10
 FIINIT(SALT)=1.0E-10 ;FIINIT(C2)=1.0E-10
 FIINIT(F1)=1.0E-10 ;FIINIT(F2)=1.0E-10
 FIINIT(F3)=1.0E-10 ;FIINIT(F4)=1.0E-10
 FIINIT(F5)=1.0E-10 ;FIINIT(F6)=1.0E-10
 FIINIT(F7)=1.0E-10 ;FIINIT(F8)=1.0E-10
 FIINIT(F9)=1.0E-10 ;FIINIT(F10)=1.0E-10
 FIINIT(AVEF)=1.0E-10 ;FIINIT(MNSQ)=1.0E-10
 FIINIT(RATE)=1.0E-10 ;FIINIT(LEN1)=1.0E-06
 FIINIT(MIXL)=1.0E-06 ;FIINIT(ENUT)=1.0E-10
 FIINIT(EPKE)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
 FIINIT(TMP1)=1.0E-10
 
 PATCH(SOL0 ,INIVAL, 1, 20, 1, 30, 1, 1, 1, 1)
 INIT(SOL0 ,U1 ,0. , GRND )
 INIT(SOL0 ,V1 ,0. , GRND )
 INIT(SOL0 ,H1 ,0. , GRND )
 INIT(SOL0 ,SALT,0. , GRND )
 
 PATCH(LOWER ,INIVAL, 1, 20, 1, 15, 1, 1, 1, 1)
 INIT(LOWER ,F1 ,0. ,1. )
 
 PATCH(UPPER ,INIVAL, 1, 20, 16, 30, 1, 1, 1, 1)
 INIT(UPPER ,F10 ,0. ,1. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 600)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(PREF ,CELL , 1, 1, 1, 1, 1, 1, 1, 600)
 COVAL(PREF ,P1 ,1. ,0. )
 
 PATCH(TRACE ,CELL , 1, 20, 15, 16, 1, 1, 1, 600)
 COVAL(TRACE ,C2 , FIXVAL ,1. )
 
 PATCH(BUOY ,PHASEM, 1, 20, 1, 30, 1, 1, 1, 600)
 COVAL(BUOY ,U1 , FIXFLU , GRND2 )
 COVAL(BUOY ,V1 , FIXFLU , GRND2 )
 
 PATCH(KEBUOY ,PHASEM, 1, 20, 1, 30, 1, 1, 1, 600)
 COVAL(KEBUOY ,KE , GRND4 , GRND4 )
 COVAL(KEBUOY ,EP , GRND4 , GRND4 )
 XCYCLE = T
 EGWF = T
 WALLCO = GRND2
 BUOYA =-0.342368 ; BUOYB =-9.804025 ; BUOYC =0.
 BUOYD =1086.97583
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 3 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(KE)=20
 LITER(EP)=20 ;LITER(H1)=20
 LITER(SALT)=20 ;LITER(C2)=20
 LITER(F1)=20 ;LITER(F2)=20
 LITER(F3)=20 ;LITER(F4)=20
 LITER(F5)=20 ;LITER(F6)=20
 LITER(F7)=20 ;LITER(F8)=20
 LITER(F9)=20 ;LITER(F10)=20
 LITER(MIXL)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
 ENDIT(SALT)=1.0E-03 ;ENDIT(C2)=1.0E-03
 ENDIT(F1)=1.0E-03 ;ENDIT(F2)=1.0E-03
 ENDIT(F3)=1.0E-03 ;ENDIT(F4)=1.0E-03
 ENDIT(F5)=1.0E-03 ;ENDIT(F6)=1.0E-03
 ENDIT(F7)=1.0E-03 ;ENDIT(F8)=1.0E-03
 ENDIT(F9)=1.0E-03 ;ENDIT(F10)=1.0E-03
 ENDIT(MIXL)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.7)
 RELAX(U1,FALSDT,5.333334E-03)
 RELAX(V1,FALSDT,5.333334E-03)
 RELAX(KE,LINRLX,0.7)
 RELAX(EP,LINRLX,0.7)
 RELAX(H1,FALSDT,5.333334E-03)
 RELAX(SALT,FALSDT,5.333334E-03)
 RELAX(C2,FALSDT,1.0E+09)
 RELAX(F1,LINRLX,0.5)
 RELAX(F2,LINRLX,0.5)
 RELAX(F3,LINRLX,0.5)
 RELAX(F4,LINRLX,0.5)
 RELAX(F5,LINRLX,0.5)
 RELAX(F6,LINRLX,0.5)
 RELAX(F7,LINRLX,0.5)
 RELAX(F8,LINRLX,0.5)
 RELAX(F9,LINRLX,0.5)
 RELAX(F10,LINRLX,0.5)
 RELAX(AVEF,LINRLX,1.)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(MIXL,FALSDT,1.0E+09)
 RELAX(ENUT,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(DEN1,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 KELIN = 3
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(SALT)=1.0E+10 ;VARMIN(SALT)=-1.0E+10
 VARMAX(C2)=1.0E+10 ;VARMIN(C2)=-1.0E+10
 VARMAX(F1)=1. ;VARMIN(F1)=0.
 VARMAX(F2)=1. ;VARMIN(F2)=0.
 VARMAX(F3)=1. ;VARMIN(F3)=0.
 VARMAX(F4)=1. ;VARMIN(F4)=0.
 VARMAX(F5)=1. ;VARMIN(F5)=0.
 VARMAX(F6)=1. ;VARMIN(F6)=0.
 VARMAX(F7)=1. ;VARMIN(F7)=0.
 VARMAX(F8)=1. ;VARMIN(F8)=0.
 VARMAX(F9)=1. ;VARMIN(F9)=0.
 VARMAX(F10)=1. ;VARMIN(F10)=0.
 VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
 VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IBUOYB = 14 ;IBUOYC = 16
 ISG62 = 1
 SPEDAT(SET,KEBUOY,GCEB,R,0.2)
 SPEDAT(SET,STORED,MIXL,C,=LEN1!ZSLFIN)
 SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN)
 SPEDAT(SET,MFM,MFMMOD,C,MNSQ)
 SPEDAT(SET,MFM,NFLUIDS,I,10)
 SPEDAT(SET,MFM,NFLR,I,1)
 SPEDAT(SET,MFM,NFLF,I,10)
 SPEDAT(SET,MFM,CONMIX,R,1.)
 SPEDAT(SET,MFM,VISCON,R,20.)
 SPEDAT(SET,MFM,POPMIN,R,-0.3)
 SPEDAT(SET,MFM,POPMAX,R,0.)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 IG( 1) = 2
 IG( 2) = 2
 IG( 3) = 2
 RG( 1) =0.3
 RG( 2) =1.2
 RG( 3) =1086.97583
 RG( 4) =4180.
 RG( 5) =4.598E+04
 RG( 6) =-2.951234E-04
 RG( 7) =5.016E+04
 RG( 8) =4.18E+04
 RG( 9) =9.81
 RG( 10) =7.53305E-04
 RG( 11) =3.
 RG( 12) =218.
 RG( 14) =110.5
 RG( 15) =5.0E-03
 RG( 16) =7.999999E-03
 RG( 17) =5.0E-03
 RG( 18) =8.333333E-03
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(SALT,Y,Y,Y,Y,Y,N)
 OUTPUT(C2,Y,N,Y,Y,Y,Y)
 OUTPUT(F1,Y,N,Y,Y,Y,Y)
 OUTPUT(F2,Y,N,Y,Y,Y,Y)
 OUTPUT(F3,Y,N,Y,Y,Y,Y)
 OUTPUT(F4,Y,N,Y,Y,Y,Y)
 OUTPUT(F5,Y,N,Y,Y,Y,Y)
 OUTPUT(F6,Y,N,Y,Y,Y,Y)
 OUTPUT(F7,Y,N,Y,Y,Y,Y)
 OUTPUT(F8,Y,N,Y,Y,Y,Y)
 OUTPUT(F9,Y,N,Y,Y,Y,Y)
 OUTPUT(F10,Y,N,Y,Y,Y,Y)
 OUTPUT(AVEF,Y,N,Y,N,N,N)
 OUTPUT(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(RATE,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(MIXL,Y,N,Y,Y,Y,Y)
 OUTPUT(ENUT,Y,Y,Y,Y,Y,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 15 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(PROF1 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF1 ,F1 ,0. ,0. )
 
 PATCH(PROF2 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF2 ,F2 ,0. ,0. )
 
 PATCH(PROF3 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF3 ,F3 ,0. ,0. )
 
 PATCH(PROF4 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF4 ,F4 ,0. ,0. )
 
 PATCH(PROF5 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF5 ,F5 ,0. ,0. )
 
 PATCH(PROF6 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF6 ,F6 ,0. ,0. )
 
 PATCH(PROF7 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF7 ,F7 ,0. ,0. )
 
 PATCH(PROF8 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF8 ,F8 ,0. ,0. )
 
 PATCH(PROF9 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF9 ,F9 ,0. ,0. )
 
 PATCH(PROF10 ,PROFIL, 10, 10, 1, 30, 1, 1, 1, 1)
 PLOT(PROF10 ,F10 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 50 ;IDISPB = 0 ;IDISPC = 0
STOP