TALK=T;RUN(1,1)
  PHOTON USE
  p
  parphi
  1 5000 1

  gr ou x 1
  con f2 x 1 fi;0.001
  upause 1
  con f5 x 1 fi;0.001
  upause 1
  con f8 x 1 fi;0.001
  upause 1
  con f11 x 1 fi;0.001
  upause 1
  con f14 x 1 fi;0.001
  upause 1
  con f17 x 1 fi;0.001
  upause 1
  con f20 x 1 fi;0.001
  upause 1
  con f23 x 1 fi;0.001
  upause 1
  con f26 x 1 fi;0.001
  upause 1
  con f29 x 1 fi;0.001
  upause 1
  ENDUSE
  DISPLAY
    This Q1 uses the multi-fluid model to simulate turbulence in
    the surface layer of a body of water subjected to a uniform
    shear stress. The process is time-dependent.

    A 2D population of fluids is used, the population-defining
    attributes being vertical velocity, v,  and horizontal
    velocity, w.

    The geometrical grid is 1D, the discretized direction being
    vertical, of 1 meter length.

    The shear is simulated by filling the bottom cell with the
    fluid having the attributes v=0, w=0, and the top cell with
    the fluid having the attributes v=0, w=1 meter/second.

    The micro-mixing rate is set equal to RATE, prescribed as
    1.0 hertz, times a constant: CONMIX.

    Mixing is simulated by setting the effective viscosity to 1.0
    meter ** 2 /seconds.
    The geometrical grid is as shown below on the left and the 2D
    population grid as shown on the right.

    |--------|
    |   T    |          v                                        NV
    |--------|        +0.5 |------|------|------|------|------| int
    |        |          ^  |      |      |      |      |      | erv
    |--------|          |  |------|------|------|------|------| als
    |        |         IV  |      |      |      |      |      |
    |--------|          |  |------|------|------|------|------|
    |        |          |  |   IFLUID=IV+(IW-1)*NV     |      |
    |--------|          |  |------|------|------|------|------|
    |        |         0.0 |  b   |      |      |      |  t   |
    |--------|             |------|------|------|------|------|
    |        |             |      |      |      |      |      |
    |--------| y           |------|------|------|------|------|
    |        |             |      |      |      |      |      |
    |--------| ^           |------|------|------|------|------|
    |        | |           |      |      |      |      |      |
    |--------| |       -0.5|------|------|------|------|------|
    |   B    | |                                IW
    |--------| |             0.0   ------------> w        1.0
     ---> z                        NW intervals
  ENDDIS
  sharbegin
  index 1
  power 2
  factor 0.00001
  gradfac 0
  sharend
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(1d ocean layer; 2d pop. ; v & w as PDAs )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =10.
 FSTEP = 1 ;LSTEP = 10
 TFRAC(1)=0.1 ;TFRAC(2)=0.2
 TFRAC(3)=0.3 ;TFRAC(4)=0.4
 TFRAC(5)=0.5 ;TFRAC(6)=0.6
 TFRAC(7)=0.7 ;TFRAC(8)=0.8
 TFRAC(9)=0.9 ;TFRAC(10)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 9
 YVLAST =1.
 YFRAC(1)=0.111111 ;YFRAC(2)=0.222222
 YFRAC(3)=0.333333 ;YFRAC(4)=0.444444
 YFRAC(5)=0.555556 ;YFRAC(6)=0.666667
 YFRAC(7)=0.777778 ;YFRAC(8)=0.888889
 YFRAC(9)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(7)=W1 ;NAME(19)=VIST
 NAME(20)=F1 ;NAME(21)=F2
 NAME(22)=F3 ;NAME(23)=F4
 NAME(24)=F5 ;NAME(25)=F6
 NAME(26)=F7 ;NAME(27)=F8
 NAME(28)=F9 ;NAME(29)=F10
 NAME(30)=F11 ;NAME(31)=F12
 NAME(32)=F13 ;NAME(33)=F14
 NAME(34)=F15 ;NAME(35)=F16
 NAME(36)=F17 ;NAME(37)=F18
 NAME(38)=F19 ;NAME(39)=F20
 NAME(40)=F21 ;NAME(41)=F22
 NAME(42)=F23 ;NAME(43)=F24
 NAME(44)=F25 ;NAME(45)=F26
 NAME(46)=F27 ;NAME(47)=F28
 NAME(48)=F29 ;NAME(49)=F30
 NAME(50)=F31 ;NAME(51)=F32
 NAME(52)=F33 ;NAME(53)=F34
 NAME(54)=F35 ;NAME(55)=F36
 NAME(56)=F37 ;NAME(57)=F38
 NAME(58)=F39 ;NAME(59)=F40
 NAME(60)=F41 ;NAME(61)=F42
 NAME(62)=F43 ;NAME(63)=F44
 NAME(64)=F45 ;NAME(65)=F46
 NAME(66)=F47 ;NAME(67)=F48
 NAME(68)=F49 ;NAME(69)=F50
 NAME(70)=F51 ;NAME(71)=F52
 NAME(72)=F53 ;NAME(73)=F54
 NAME(74)=F55 ;NAME(75)=F56
 NAME(76)=F57 ;NAME(77)=F58
 NAME(78)=F59 ;NAME(79)=F60
 NAME(80)=F61 ;NAME(81)=F62
 NAME(82)=F63 ;NAME(83)=F64
 NAME(84)=F65 ;NAME(85)=F66
 NAME(86)=F67 ;NAME(87)=F68
 NAME(88)=F69 ;NAME(89)=F70
 NAME(90)=F71 ;NAME(91)=F72
 NAME(92)=F73 ;NAME(93)=F74
 NAME(94)=F75 ;NAME(95)=F76
 NAME(96)=F77 ;NAME(97)=F78
 NAME(98)=F79 ;NAME(99)=F80
 NAME(100)=F81 ;NAME(101)=F82
 NAME(102)=F83 ;NAME(103)=F84
 NAME(104)=F85 ;NAME(105)=F86
 NAME(106)=F87 ;NAME(107)=F88
 NAME(108)=F89 ;NAME(109)=F90
 NAME(110)=F91 ;NAME(111)=F92
 NAME(112)=F93 ;NAME(113)=F94
 NAME(114)=F95 ;NAME(115)=F96
 NAME(116)=F97 ;NAME(117)=F98
 NAME(118)=F99 ;NAME(119)=F100
 NAME(120)=F101 ;NAME(121)=F102
 NAME(122)=F103 ;NAME(123)=F104
 NAME(124)=F105 ;NAME(125)=F106
 NAME(126)=F107 ;NAME(127)=F108
 NAME(128)=F109 ;NAME(129)=F110
 NAME(130)=F111 ;NAME(131)=F112
 NAME(132)=F113 ;NAME(133)=F114
 NAME(134)=F115 ;NAME(135)=F116
 NAME(136)=F117 ;NAME(137)=F118
 NAME(138)=F119 ;NAME(139)=F120
 NAME(140)=F121 ;NAME(141)=FRSM
 NAME(142)=NETM ;NAME(143)=PGRD
 NAME(144)=FLSM ;NAME(145)=RMOA
 NAME(146)=RMSA ;NAME(147)=RMSO
 NAME(148)=AVER ;NAME(149)=AVEF
 NAME(150)=RATE
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(VIST,Y,N,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(F11,Y,Y,N,N,N,Y)
 SOLUTN(F12,Y,Y,N,N,N,Y)
 SOLUTN(F13,Y,Y,N,N,N,Y)
 SOLUTN(F14,Y,Y,N,N,N,Y)
 SOLUTN(F15,Y,Y,N,N,N,Y)
 SOLUTN(F16,Y,Y,N,N,N,Y)
 SOLUTN(F17,Y,Y,N,N,N,Y)
 SOLUTN(F18,Y,Y,N,N,N,Y)
 SOLUTN(F19,Y,Y,N,N,N,Y)
 SOLUTN(F20,Y,Y,N,N,N,Y)
 SOLUTN(F21,Y,Y,N,N,N,Y)
 SOLUTN(F22,Y,Y,N,N,N,Y)
 SOLUTN(F23,Y,Y,N,N,N,Y)
 SOLUTN(F24,Y,Y,N,N,N,Y)
 SOLUTN(F25,Y,Y,N,N,N,Y)
 SOLUTN(F26,Y,Y,N,N,N,Y)
 SOLUTN(F27,Y,Y,N,N,N,Y)
 SOLUTN(F28,Y,Y,N,N,N,Y)
 SOLUTN(F29,Y,Y,N,N,N,Y)
 SOLUTN(F30,Y,Y,N,N,N,Y)
 SOLUTN(F31,Y,Y,N,N,N,Y)
 SOLUTN(F32,Y,Y,N,N,N,Y)
 SOLUTN(F33,Y,Y,N,N,N,Y)
 SOLUTN(F34,Y,Y,N,N,N,Y)
 SOLUTN(F35,Y,Y,N,N,N,Y)
 SOLUTN(F36,Y,Y,N,N,N,Y)
 SOLUTN(F37,Y,Y,N,N,N,Y)
 SOLUTN(F38,Y,Y,N,N,N,Y)
 SOLUTN(F39,Y,Y,N,N,N,Y)
 SOLUTN(F40,Y,Y,N,N,N,Y)
 SOLUTN(F41,Y,Y,N,N,N,Y)
 SOLUTN(F42,Y,Y,N,N,N,Y)
 SOLUTN(F43,Y,Y,N,N,N,Y)
 SOLUTN(F44,Y,Y,N,N,N,Y)
 SOLUTN(F45,Y,Y,N,N,N,Y)
 SOLUTN(F46,Y,Y,N,N,N,Y)
 SOLUTN(F47,Y,Y,N,N,N,Y)
 SOLUTN(F48,Y,Y,N,N,N,Y)
 SOLUTN(F49,Y,Y,N,N,N,Y)
 SOLUTN(F50,Y,Y,N,N,N,Y)
 SOLUTN(F51,Y,Y,N,N,N,Y)
 SOLUTN(F52,Y,Y,N,N,N,Y)
 SOLUTN(F53,Y,Y,N,N,N,Y)
 SOLUTN(F54,Y,Y,N,N,N,Y)
 SOLUTN(F55,Y,Y,N,N,N,Y)
 SOLUTN(F56,Y,Y,N,N,N,Y)
 SOLUTN(F57,Y,Y,N,N,N,Y)
 SOLUTN(F58,Y,Y,N,N,N,Y)
 SOLUTN(F59,Y,Y,N,N,N,Y)
 SOLUTN(F60,Y,Y,N,N,N,Y)
 SOLUTN(F61,Y,Y,N,N,N,Y)
 SOLUTN(F62,Y,Y,N,N,N,Y)
 SOLUTN(F63,Y,Y,N,N,N,Y)
 SOLUTN(F64,Y,Y,N,N,N,Y)
 SOLUTN(F65,Y,Y,N,N,N,Y)
 SOLUTN(F66,Y,Y,N,N,N,Y)
 SOLUTN(F67,Y,Y,N,N,N,Y)
 SOLUTN(F68,Y,Y,N,N,N,Y)
 SOLUTN(F69,Y,Y,N,N,N,Y)
 SOLUTN(F70,Y,Y,N,N,N,Y)
 SOLUTN(F71,Y,Y,N,N,N,Y)
 SOLUTN(F72,Y,Y,N,N,N,Y)
 SOLUTN(F73,Y,Y,N,N,N,Y)
 SOLUTN(F74,Y,Y,N,N,N,Y)
 SOLUTN(F75,Y,Y,N,N,N,Y)
 SOLUTN(F76,Y,Y,N,N,N,Y)
 SOLUTN(F77,Y,Y,N,N,N,Y)
 SOLUTN(F78,Y,Y,N,N,N,Y)
 SOLUTN(F79,Y,Y,N,N,N,Y)
 SOLUTN(F80,Y,Y,N,N,N,Y)
 SOLUTN(F81,Y,Y,N,N,N,Y)
 SOLUTN(F82,Y,Y,N,N,N,Y)
 SOLUTN(F83,Y,Y,N,N,N,Y)
 SOLUTN(F84,Y,Y,N,N,N,Y)
 SOLUTN(F85,Y,Y,N,N,N,Y)
 SOLUTN(F86,Y,Y,N,N,N,Y)
 SOLUTN(F87,Y,Y,N,N,N,Y)
 SOLUTN(F88,Y,Y,N,N,N,Y)
 SOLUTN(F89,Y,Y,N,N,N,Y)
 SOLUTN(F90,Y,Y,N,N,N,Y)
 SOLUTN(F91,Y,Y,N,N,N,Y)
 SOLUTN(F92,Y,Y,N,N,N,Y)
 SOLUTN(F93,Y,Y,N,N,N,Y)
 SOLUTN(F94,Y,Y,N,N,N,Y)
 SOLUTN(F95,Y,Y,N,N,N,Y)
 SOLUTN(F96,Y,Y,N,N,N,Y)
 SOLUTN(F97,Y,Y,N,N,N,Y)
 SOLUTN(F98,Y,Y,N,N,N,Y)
 SOLUTN(F99,Y,Y,N,N,N,Y)
 SOLUTN(F100,Y,Y,N,N,N,Y)
 SOLUTN(F101,Y,Y,N,N,N,Y)
 SOLUTN(F102,Y,Y,N,N,N,Y)
 SOLUTN(F103,Y,Y,N,N,N,Y)
 SOLUTN(F104,Y,Y,N,N,N,Y)
 SOLUTN(F105,Y,Y,N,N,N,Y)
 SOLUTN(F106,Y,Y,N,N,N,Y)
 SOLUTN(F107,Y,Y,N,N,N,Y)
 SOLUTN(F108,Y,Y,N,N,N,Y)
 SOLUTN(F109,Y,Y,N,N,N,Y)
 SOLUTN(F110,Y,Y,N,N,N,Y)
 SOLUTN(F111,Y,Y,N,N,N,Y)
 SOLUTN(F112,Y,Y,N,N,N,Y)
 SOLUTN(F113,Y,Y,N,N,N,Y)
 SOLUTN(F114,Y,Y,N,N,N,Y)
 SOLUTN(F115,Y,Y,N,N,N,Y)
 SOLUTN(F116,Y,Y,N,N,N,Y)
 SOLUTN(F117,Y,Y,N,N,N,Y)
 SOLUTN(F118,Y,Y,N,N,N,Y)
 SOLUTN(F119,Y,Y,N,N,N,Y)
 SOLUTN(F120,Y,Y,N,N,N,Y)
 SOLUTN(F121,Y,Y,N,N,N,Y)
 SOLUTN(FRSM,Y,N,N,N,N,Y)
 SOLUTN(NETM,Y,N,N,N,N,Y)
 SOLUTN(PGRD,Y,N,N,N,N,Y)
 SOLUTN(FLSM,Y,N,N,N,N,Y)
 SOLUTN(RMOA,Y,N,N,N,N,Y)
 SOLUTN(RMSA,Y,N,N,N,N,Y)
 SOLUTN(RMSO,Y,N,N,N,N,Y)
 SOLUTN(AVER,Y,N,N,N,N,Y)
 SOLUTN(AVEF,Y,N,N,N,N,Y)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 VIST = 19
 ************************************************************
  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(W1,N,N,Y,Y,Y,Y)
 TERMS(F1,N,N,Y,Y,Y,N)
 TERMS(F2,N,N,Y,Y,Y,N)
 TERMS(F3,N,N,Y,Y,Y,N)
 TERMS(F4,N,N,Y,Y,Y,N)
 TERMS(F5,N,N,Y,Y,Y,N)
 TERMS(F6,N,N,Y,Y,Y,N)
 TERMS(F7,N,N,Y,Y,Y,N)
 TERMS(F8,N,N,Y,Y,Y,N)
 TERMS(F9,N,N,Y,Y,Y,N)
 TERMS(F10,N,N,Y,Y,Y,N)
 TERMS(F11,N,N,Y,Y,Y,N)
 TERMS(F12,N,N,Y,Y,Y,N)
 TERMS(F13,N,N,Y,Y,Y,N)
 TERMS(F14,N,N,Y,Y,Y,N)
 TERMS(F15,N,N,Y,Y,Y,N)
 TERMS(F16,N,N,Y,Y,Y,N)
 TERMS(F17,N,N,Y,Y,Y,N)
 TERMS(F18,N,N,Y,Y,Y,N)
 TERMS(F19,N,N,Y,Y,Y,N)
 TERMS(F20,N,N,Y,Y,Y,N)
 TERMS(F21,N,N,Y,Y,Y,N)
 TERMS(F22,N,N,Y,Y,Y,N)
 TERMS(F23,N,N,Y,Y,Y,N)
 TERMS(F24,N,N,Y,Y,Y,N)
 TERMS(F25,N,N,Y,Y,Y,N)
 TERMS(F26,N,N,Y,Y,Y,N)
 TERMS(F27,N,N,Y,Y,Y,N)
 TERMS(F28,N,N,Y,Y,Y,N)
 TERMS(F29,N,N,Y,Y,Y,N)
 TERMS(F30,N,N,Y,Y,Y,N)
 TERMS(F31,N,N,Y,Y,Y,N)
 TERMS(F32,N,N,Y,Y,Y,N)
 TERMS(F33,N,N,Y,Y,Y,N)
 TERMS(F34,N,N,Y,Y,Y,N)
 TERMS(F35,N,N,Y,Y,Y,N)
 TERMS(F36,N,N,Y,Y,Y,N)
 TERMS(F37,N,N,Y,Y,Y,N)
 TERMS(F38,N,N,Y,Y,Y,N)
 TERMS(F39,N,N,Y,Y,Y,N)
 TERMS(F40,N,N,Y,Y,Y,N)
 TERMS(F41,N,N,Y,Y,Y,N)
 TERMS(F42,N,N,Y,Y,Y,N)
 TERMS(F43,N,N,Y,Y,Y,N)
 TERMS(F44,N,N,Y,Y,Y,N)
 TERMS(F45,N,N,Y,Y,Y,N)
 TERMS(F46,N,N,Y,Y,Y,N)
 TERMS(F47,N,N,Y,Y,Y,N)
 TERMS(F48,N,N,Y,Y,Y,N)
 TERMS(F49,N,N,Y,Y,Y,N)
 TERMS(F50,N,N,Y,Y,Y,N)
 TERMS(F51,N,N,Y,Y,Y,N)
 TERMS(F52,N,N,Y,Y,Y,N)
 TERMS(F53,N,N,Y,Y,Y,N)
 TERMS(F54,N,N,Y,Y,Y,N)
 TERMS(F55,N,N,Y,Y,Y,N)
 TERMS(F56,N,N,Y,Y,Y,N)
 TERMS(F57,N,N,Y,Y,Y,N)
 TERMS(F58,N,N,Y,Y,Y,N)
 TERMS(F59,N,N,Y,Y,Y,N)
 TERMS(F60,N,N,Y,Y,Y,N)
 TERMS(F61,N,N,Y,Y,Y,N)
 TERMS(F62,N,N,Y,Y,Y,N)
 TERMS(F63,N,N,Y,Y,Y,N)
 TERMS(F64,N,N,Y,Y,Y,N)
 TERMS(F65,N,N,Y,Y,Y,N)
 TERMS(F66,N,N,Y,Y,Y,N)
 TERMS(F67,N,N,Y,Y,Y,N)
 TERMS(F68,N,N,Y,Y,Y,N)
 TERMS(F69,N,N,Y,Y,Y,N)
 TERMS(F70,N,N,Y,Y,Y,N)
 TERMS(F71,N,N,Y,Y,Y,N)
 TERMS(F72,N,N,Y,Y,Y,N)
 TERMS(F73,N,N,Y,Y,Y,N)
 TERMS(F74,N,N,Y,Y,Y,N)
 TERMS(F75,N,N,Y,Y,Y,N)
 TERMS(F76,N,N,Y,Y,Y,N)
 TERMS(F77,N,N,Y,Y,Y,N)
 TERMS(F78,N,N,Y,Y,Y,N)
 TERMS(F79,N,N,Y,Y,Y,N)
 TERMS(F80,N,N,Y,Y,Y,N)
 TERMS(F81,N,N,Y,Y,Y,N)
 TERMS(F82,N,N,Y,Y,Y,N)
 TERMS(F83,N,N,Y,Y,Y,N)
 TERMS(F84,N,N,Y,Y,Y,N)
 TERMS(F85,N,N,Y,Y,Y,N)
 TERMS(F86,N,N,Y,Y,Y,N)
 TERMS(F87,N,N,Y,Y,Y,N)
 TERMS(F88,N,N,Y,Y,Y,N)
 TERMS(F89,N,N,Y,Y,Y,N)
 TERMS(F90,N,N,Y,Y,Y,N)
 TERMS(F91,N,N,Y,Y,Y,N)
 TERMS(F92,N,N,Y,Y,Y,N)
 TERMS(F93,N,N,Y,Y,Y,N)
 TERMS(F94,N,N,Y,Y,Y,N)
 TERMS(F95,N,N,Y,Y,Y,N)
 TERMS(F96,N,N,Y,Y,Y,N)
 TERMS(F97,N,N,Y,Y,Y,N)
 TERMS(F98,N,N,Y,Y,Y,N)
 TERMS(F99,N,N,Y,Y,Y,N)
 TERMS(F100,N,N,Y,Y,Y,N)
 TERMS(F101,N,N,Y,Y,Y,N)
 TERMS(F102,N,N,Y,Y,Y,N)
 TERMS(F103,N,N,Y,Y,Y,N)
 TERMS(F104,N,N,Y,Y,Y,N)
 TERMS(F105,N,N,Y,Y,Y,N)
 TERMS(F106,N,N,Y,Y,Y,N)
 TERMS(F107,N,N,Y,Y,Y,N)
 TERMS(F108,N,N,Y,Y,Y,N)
 TERMS(F109,N,N,Y,Y,Y,N)
 TERMS(F110,N,N,Y,Y,Y,N)
 TERMS(F111,N,N,Y,Y,Y,N)
 TERMS(F112,N,N,Y,Y,Y,N)
 TERMS(F113,N,N,Y,Y,Y,N)
 TERMS(F114,N,N,Y,Y,Y,N)
 TERMS(F115,N,N,Y,Y,Y,N)
 TERMS(F116,N,N,Y,Y,Y,N)
 TERMS(F117,N,N,Y,Y,Y,N)
 TERMS(F118,N,N,Y,Y,Y,N)
 TERMS(F119,N,N,Y,Y,Y,N)
 TERMS(F120,N,N,Y,Y,Y,N)
 TERMS(F121,N,N,Y,Y,Y,N)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1000. ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =1.
 PRNDTL(W1)=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(F11)=1.
 PRNDTL(F12)=1. ;PRNDTL(F13)=1.
 PRNDTL(F14)=1. ;PRNDTL(F15)=1.
 PRNDTL(F16)=1. ;PRNDTL(F17)=1.
 PRNDTL(F18)=1. ;PRNDTL(F19)=1.
 PRNDTL(F20)=1. ;PRNDTL(F21)=1.
 PRNDTL(F22)=1. ;PRNDTL(F23)=1.
 PRNDTL(F24)=1. ;PRNDTL(F25)=1.
 PRNDTL(F26)=1. ;PRNDTL(F27)=1.
 PRNDTL(F28)=1. ;PRNDTL(F29)=1.
 PRNDTL(F30)=1. ;PRNDTL(F31)=1.
 PRNDTL(F32)=1. ;PRNDTL(F33)=1.
 PRNDTL(F34)=1. ;PRNDTL(F35)=1.
 PRNDTL(F36)=1. ;PRNDTL(F37)=1.
 PRNDTL(F38)=1. ;PRNDTL(F39)=1.
 PRNDTL(F40)=1. ;PRNDTL(F41)=1.
 PRNDTL(F42)=1. ;PRNDTL(F43)=1.
 PRNDTL(F44)=1. ;PRNDTL(F45)=1.
 PRNDTL(F46)=1. ;PRNDTL(F47)=1.
 PRNDTL(F48)=1. ;PRNDTL(F49)=1.
 PRNDTL(F50)=1. ;PRNDTL(F51)=1.
 PRNDTL(F52)=1. ;PRNDTL(F53)=1.
 PRNDTL(F54)=1. ;PRNDTL(F55)=1.
 PRNDTL(F56)=1. ;PRNDTL(F57)=1.
 PRNDTL(F58)=1. ;PRNDTL(F59)=1.
 PRNDTL(F60)=1. ;PRNDTL(F61)=1.
 PRNDTL(F62)=1. ;PRNDTL(F63)=1.
 PRNDTL(F64)=1. ;PRNDTL(F65)=1.
 PRNDTL(F66)=1. ;PRNDTL(F67)=1.
 PRNDTL(F68)=1. ;PRNDTL(F69)=1.
 PRNDTL(F70)=1. ;PRNDTL(F71)=1.
 PRNDTL(F72)=1. ;PRNDTL(F73)=1.
 PRNDTL(F74)=1. ;PRNDTL(F75)=1.
 PRNDTL(F76)=1. ;PRNDTL(F77)=1.
 PRNDTL(F78)=1. ;PRNDTL(F79)=1.
 PRNDTL(F80)=1. ;PRNDTL(F81)=1.
 PRNDTL(F82)=1. ;PRNDTL(F83)=1.
 PRNDTL(F84)=1. ;PRNDTL(F85)=1.
 PRNDTL(F86)=1. ;PRNDTL(F87)=1.
 PRNDTL(F88)=1. ;PRNDTL(F89)=1.
 PRNDTL(F90)=1. ;PRNDTL(F91)=1.
 PRNDTL(F92)=1. ;PRNDTL(F93)=1.
 PRNDTL(F94)=1. ;PRNDTL(F95)=1.
 PRNDTL(F96)=1. ;PRNDTL(F97)=1.
 PRNDTL(F98)=1. ;PRNDTL(F99)=1.
 PRNDTL(F100)=1. ;PRNDTL(F101)=1.
 PRNDTL(F102)=1. ;PRNDTL(F103)=1.
 PRNDTL(F104)=1. ;PRNDTL(F105)=1.
 PRNDTL(F106)=1. ;PRNDTL(F107)=1.
 PRNDTL(F108)=1. ;PRNDTL(F109)=1.
 PRNDTL(F110)=1. ;PRNDTL(F111)=1.
 PRNDTL(F112)=1. ;PRNDTL(F113)=1.
 PRNDTL(F114)=1. ;PRNDTL(F115)=1.
 PRNDTL(F116)=1. ;PRNDTL(F117)=1.
 PRNDTL(F118)=1. ;PRNDTL(F119)=1.
 PRNDTL(F120)=1. ;PRNDTL(F121)=1.
 PRT(W1)=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(F11)=1.
 PRT(F12)=1. ;PRT(F13)=1.
 PRT(F14)=1. ;PRT(F15)=1.
 PRT(F16)=1. ;PRT(F17)=1.
 PRT(F18)=1. ;PRT(F19)=1.
 PRT(F20)=1. ;PRT(F21)=1.
 PRT(F22)=1. ;PRT(F23)=1.
 PRT(F24)=1. ;PRT(F25)=1.
 PRT(F26)=1. ;PRT(F27)=1.
 PRT(F28)=1. ;PRT(F29)=1.
 PRT(F30)=1. ;PRT(F31)=1.
 PRT(F32)=1. ;PRT(F33)=1.
 PRT(F34)=1. ;PRT(F35)=1.
 PRT(F36)=1. ;PRT(F37)=1.
 PRT(F38)=1. ;PRT(F39)=1.
 PRT(F40)=1. ;PRT(F41)=1.
 PRT(F42)=1. ;PRT(F43)=1.
 PRT(F44)=1. ;PRT(F45)=1.
 PRT(F46)=1. ;PRT(F47)=1.
 PRT(F48)=1. ;PRT(F49)=1.
 PRT(F50)=1. ;PRT(F51)=1.
 PRT(F52)=1. ;PRT(F53)=1.
 PRT(F54)=1. ;PRT(F55)=1.
 PRT(F56)=1. ;PRT(F57)=1.
 PRT(F58)=1. ;PRT(F59)=1.
 PRT(F60)=1. ;PRT(F61)=1.
 PRT(F62)=1. ;PRT(F63)=1.
 PRT(F64)=1. ;PRT(F65)=1.
 PRT(F66)=1. ;PRT(F67)=1.
 PRT(F68)=1. ;PRT(F69)=1.
 PRT(F70)=1. ;PRT(F71)=1.
 PRT(F72)=1. ;PRT(F73)=1.
 PRT(F74)=1. ;PRT(F75)=1.
 PRT(F76)=1. ;PRT(F77)=1.
 PRT(F78)=1. ;PRT(F79)=1.
 PRT(F80)=1. ;PRT(F81)=1.
 PRT(F82)=1. ;PRT(F83)=1.
 PRT(F84)=1. ;PRT(F85)=1.
 PRT(F86)=1. ;PRT(F87)=1.
 PRT(F88)=1. ;PRT(F89)=1.
 PRT(F90)=1. ;PRT(F91)=1.
 PRT(F92)=1. ;PRT(F93)=1.
 PRT(F94)=1. ;PRT(F95)=1.
 PRT(F96)=1. ;PRT(F97)=1.
 PRT(F98)=1. ;PRT(F99)=1.
 PRT(F100)=1. ;PRT(F101)=1.
 PRT(F102)=1. ;PRT(F103)=1.
 PRT(F104)=1. ;PRT(F105)=1.
 PRT(F106)=1. ;PRT(F107)=1.
 PRT(F108)=1. ;PRT(F109)=1.
 PRT(F110)=1. ;PRT(F111)=1.
 PRT(F112)=1. ;PRT(F113)=1.
 PRT(F114)=1. ;PRT(F115)=1.
 PRT(F116)=1. ;PRT(F117)=1.
 PRT(F118)=1. ;PRT(F119)=1.
 PRT(F120)=1. ;PRT(F121)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(W1)=1.0E-10 ;FIINIT(VIST)=1.
 FIINIT(F1)=0. ;FIINIT(F2)=0.
 FIINIT(F3)=0. ;FIINIT(F4)=0.
 FIINIT(F5)=0. ;FIINIT(F6)=0.
 FIINIT(F7)=0. ;FIINIT(F8)=0.
 FIINIT(F9)=0. ;FIINIT(F10)=0.
 FIINIT(F11)=0. ;FIINIT(F12)=0.
 FIINIT(F13)=0. ;FIINIT(F14)=0.
 FIINIT(F15)=0. ;FIINIT(F16)=0.
 FIINIT(F17)=0. ;FIINIT(F18)=0.
 FIINIT(F19)=0. ;FIINIT(F20)=0.
 FIINIT(F21)=0. ;FIINIT(F22)=0.
 FIINIT(F23)=0. ;FIINIT(F24)=0.
 FIINIT(F25)=0. ;FIINIT(F26)=0.
 FIINIT(F27)=0. ;FIINIT(F28)=0.
 FIINIT(F29)=0. ;FIINIT(F30)=0.
 FIINIT(F31)=0. ;FIINIT(F32)=0.
 FIINIT(F33)=0. ;FIINIT(F34)=0.
 FIINIT(F35)=0. ;FIINIT(F36)=0.
 FIINIT(F37)=0. ;FIINIT(F38)=0.
 FIINIT(F39)=0. ;FIINIT(F40)=0.
 FIINIT(F41)=0. ;FIINIT(F42)=0.
 FIINIT(F43)=0. ;FIINIT(F44)=0.
 FIINIT(F45)=0. ;FIINIT(F46)=0.
 FIINIT(F47)=0. ;FIINIT(F48)=0.
 FIINIT(F49)=0. ;FIINIT(F50)=0.
 FIINIT(F51)=0. ;FIINIT(F52)=0.
 FIINIT(F53)=0. ;FIINIT(F54)=0.
 FIINIT(F55)=0. ;FIINIT(F56)=0.
 FIINIT(F57)=0. ;FIINIT(F58)=0.
 FIINIT(F59)=0. ;FIINIT(F60)=0.
 FIINIT(F61)=0. ;FIINIT(F62)=0.
 FIINIT(F63)=0. ;FIINIT(F64)=0.
 FIINIT(F65)=0. ;FIINIT(F66)=0.
 FIINIT(F67)=0. ;FIINIT(F68)=0.
 FIINIT(F69)=0. ;FIINIT(F70)=0.
 FIINIT(F71)=0. ;FIINIT(F72)=0.
 FIINIT(F73)=0. ;FIINIT(F74)=0.
 FIINIT(F75)=0. ;FIINIT(F76)=0.
 FIINIT(F77)=0. ;FIINIT(F78)=0.
 FIINIT(F79)=0. ;FIINIT(F80)=0.
 FIINIT(F81)=0. ;FIINIT(F82)=0.
 FIINIT(F83)=0. ;FIINIT(F84)=0.
 FIINIT(F85)=0. ;FIINIT(F86)=0.
 FIINIT(F87)=0. ;FIINIT(F88)=0.
 FIINIT(F89)=0. ;FIINIT(F90)=0.
 FIINIT(F91)=0. ;FIINIT(F92)=0.
 FIINIT(F93)=0. ;FIINIT(F94)=0.
 FIINIT(F95)=0. ;FIINIT(F96)=0.
 FIINIT(F97)=0. ;FIINIT(F98)=0.
 FIINIT(F99)=0. ;FIINIT(F100)=0.
 FIINIT(F101)=0. ;FIINIT(F102)=0.
 FIINIT(F103)=0. ;FIINIT(F104)=0.
 FIINIT(F105)=0. ;FIINIT(F106)=0.
 FIINIT(F107)=0. ;FIINIT(F108)=0.
 FIINIT(F109)=0. ;FIINIT(F110)=0.
 FIINIT(F111)=0. ;FIINIT(F112)=0.
 FIINIT(F113)=0. ;FIINIT(F114)=0.
 FIINIT(F115)=0. ;FIINIT(F116)=0.
 FIINIT(F117)=0. ;FIINIT(F118)=0.
 FIINIT(F119)=0. ;FIINIT(F120)=0.
 FIINIT(F121)=0. ;FIINIT(FRSM)=1.0E-10
 FIINIT(NETM)=1.0E-10 ;FIINIT(PGRD)=1.0E-10
 FIINIT(FLSM)=1.0E-10 ;FIINIT(RMOA)=1.0E-10
 FIINIT(RMSA)=1.0E-10 ;FIINIT(RMSO)=1.0E-10
 FIINIT(AVER)=1.0E-10 ;FIINIT(AVEF)=1.0E-10
 FIINIT(RATE)=1.
 
 PATCH(TOPINIT ,INIVAL, 1, 1, 9, 9, 1, 1, 1, 1)
 INIT(TOPINIT ,F116, FIXVAL ,1. )
 
 PATCH(BTMINIT ,INIVAL, 1, 1, 1, 1, 1, 1, 1, 1)
 INIT(BTMINIT ,F6 , FIXVAL ,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(BOTTOM ,CELL , 1, 1, 1, 1, 1, 1, 1, 10)
 COVAL(BOTTOM ,W1 , FIXVAL ,0. )
 COVAL(BOTTOM ,F1 , FIXVAL ,0. )
 COVAL(BOTTOM ,F2 , FIXVAL ,0. )
 COVAL(BOTTOM ,F3 , FIXVAL ,0. )
 COVAL(BOTTOM ,F4 , FIXVAL ,0. )
 COVAL(BOTTOM ,F5 , FIXVAL ,0. )
 COVAL(BOTTOM ,F6 , FIXVAL ,1. )
 COVAL(BOTTOM ,F7 , FIXVAL ,0. )
 COVAL(BOTTOM ,F8 , FIXVAL ,0. )
 COVAL(BOTTOM ,F9 , FIXVAL ,0. )
 COVAL(BOTTOM ,F10 , FIXVAL ,0. )
 COVAL(BOTTOM ,F11 , FIXVAL ,0. )
 COVAL(BOTTOM ,F12 , FIXVAL ,0. )
 COVAL(BOTTOM ,F13 , FIXVAL ,0. )
 COVAL(BOTTOM ,F14 , FIXVAL ,0. )
 COVAL(BOTTOM ,F15 , FIXVAL ,0. )
 COVAL(BOTTOM ,F16 , FIXVAL ,0. )
 COVAL(BOTTOM ,F17 , FIXVAL ,0. )
 COVAL(BOTTOM ,F18 , FIXVAL ,0. )
 COVAL(BOTTOM ,F19 , FIXVAL ,0. )
 COVAL(BOTTOM ,F20 , FIXVAL ,0. )
 COVAL(BOTTOM ,F21 , FIXVAL ,0. )
 COVAL(BOTTOM ,F22 , FIXVAL ,0. )
 COVAL(BOTTOM ,F23 , FIXVAL ,0. )
 COVAL(BOTTOM ,F24 , FIXVAL ,0. )
 COVAL(BOTTOM ,F25 , FIXVAL ,0. )
 COVAL(BOTTOM ,F26 , FIXVAL ,0. )
 COVAL(BOTTOM ,F27 , FIXVAL ,0. )
 COVAL(BOTTOM ,F28 , FIXVAL ,0. )
 COVAL(BOTTOM ,F29 , FIXVAL ,0. )
 COVAL(BOTTOM ,F30 , FIXVAL ,0. )
 COVAL(BOTTOM ,F31 , FIXVAL ,0. )
 COVAL(BOTTOM ,F32 , FIXVAL ,0. )
 COVAL(BOTTOM ,F33 , FIXVAL ,0. )
 COVAL(BOTTOM ,F34 , FIXVAL ,0. )
 COVAL(BOTTOM ,F35 , FIXVAL ,0. )
 COVAL(BOTTOM ,F36 , FIXVAL ,0. )
 COVAL(BOTTOM ,F37 , FIXVAL ,0. )
 COVAL(BOTTOM ,F38 , FIXVAL ,0. )
 COVAL(BOTTOM ,F39 , FIXVAL ,0. )
 COVAL(BOTTOM ,F40 , FIXVAL ,0. )
 COVAL(BOTTOM ,F41 , FIXVAL ,0. )
 COVAL(BOTTOM ,F42 , FIXVAL ,0. )
 COVAL(BOTTOM ,F43 , FIXVAL ,0. )
 COVAL(BOTTOM ,F44 , FIXVAL ,0. )
 COVAL(BOTTOM ,F45 , FIXVAL ,0. )
 COVAL(BOTTOM ,F46 , FIXVAL ,0. )
 COVAL(BOTTOM ,F47 , FIXVAL ,0. )
 COVAL(BOTTOM ,F48 , FIXVAL ,0. )
 COVAL(BOTTOM ,F49 , FIXVAL ,0. )
 COVAL(BOTTOM ,F50 , FIXVAL ,0. )
 COVAL(BOTTOM ,F51 , FIXVAL ,0. )
 COVAL(BOTTOM ,F52 , FIXVAL ,0. )
 COVAL(BOTTOM ,F53 , FIXVAL ,0. )
 COVAL(BOTTOM ,F54 , FIXVAL ,0. )
 COVAL(BOTTOM ,F55 , FIXVAL ,0. )
 COVAL(BOTTOM ,F56 , FIXVAL ,0. )
 COVAL(BOTTOM ,F57 , FIXVAL ,0. )
 COVAL(BOTTOM ,F58 , FIXVAL ,0. )
 COVAL(BOTTOM ,F59 , FIXVAL ,0. )
 COVAL(BOTTOM ,F60 , FIXVAL ,0. )
 COVAL(BOTTOM ,F61 , FIXVAL ,0. )
 COVAL(BOTTOM ,F62 , FIXVAL ,0. )
 COVAL(BOTTOM ,F63 , FIXVAL ,0. )
 COVAL(BOTTOM ,F64 , FIXVAL ,0. )
 COVAL(BOTTOM ,F65 , FIXVAL ,0. )
 COVAL(BOTTOM ,F66 , FIXVAL ,0. )
 COVAL(BOTTOM ,F67 , FIXVAL ,0. )
 COVAL(BOTTOM ,F68 , FIXVAL ,0. )
 COVAL(BOTTOM ,F69 , FIXVAL ,0. )
 COVAL(BOTTOM ,F70 , FIXVAL ,0. )
 COVAL(BOTTOM ,F71 , FIXVAL ,0. )
 COVAL(BOTTOM ,F72 , FIXVAL ,0. )
 COVAL(BOTTOM ,F73 , FIXVAL ,0. )
 COVAL(BOTTOM ,F74 , FIXVAL ,0. )
 COVAL(BOTTOM ,F75 , FIXVAL ,0. )
 COVAL(BOTTOM ,F76 , FIXVAL ,0. )
 COVAL(BOTTOM ,F77 , FIXVAL ,0. )
 COVAL(BOTTOM ,F78 , FIXVAL ,0. )
 COVAL(BOTTOM ,F79 , FIXVAL ,0. )
 COVAL(BOTTOM ,F80 , FIXVAL ,0. )
 COVAL(BOTTOM ,F81 , FIXVAL ,0. )
 COVAL(BOTTOM ,F82 , FIXVAL ,0. )
 COVAL(BOTTOM ,F83 , FIXVAL ,0. )
 COVAL(BOTTOM ,F84 , FIXVAL ,0. )
 COVAL(BOTTOM ,F85 , FIXVAL ,0. )
 COVAL(BOTTOM ,F86 , FIXVAL ,0. )
 COVAL(BOTTOM ,F87 , FIXVAL ,0. )
 COVAL(BOTTOM ,F88 , FIXVAL ,0. )
 COVAL(BOTTOM ,F89 , FIXVAL ,0. )
 COVAL(BOTTOM ,F90 , FIXVAL ,0. )
 COVAL(BOTTOM ,F91 , FIXVAL ,0. )
 COVAL(BOTTOM ,F92 , FIXVAL ,0. )
 COVAL(BOTTOM ,F93 , FIXVAL ,0. )
 COVAL(BOTTOM ,F94 , FIXVAL ,0. )
 COVAL(BOTTOM ,F95 , FIXVAL ,0. )
 COVAL(BOTTOM ,F96 , FIXVAL ,0. )
 COVAL(BOTTOM ,F97 , FIXVAL ,0. )
 COVAL(BOTTOM ,F98 , FIXVAL ,0. )
 COVAL(BOTTOM ,F99 , FIXVAL ,0. )
 COVAL(BOTTOM ,F100, FIXVAL ,0. )
 COVAL(BOTTOM ,F101, FIXVAL ,0. )
 COVAL(BOTTOM ,F102, FIXVAL ,0. )
 COVAL(BOTTOM ,F103, FIXVAL ,0. )
 COVAL(BOTTOM ,F104, FIXVAL ,0. )
 COVAL(BOTTOM ,F105, FIXVAL ,0. )
 COVAL(BOTTOM ,F106, FIXVAL ,0. )
 COVAL(BOTTOM ,F107, FIXVAL ,0. )
 COVAL(BOTTOM ,F108, FIXVAL ,0. )
 COVAL(BOTTOM ,F109, FIXVAL ,0. )
 COVAL(BOTTOM ,F110, FIXVAL ,0. )
 COVAL(BOTTOM ,F111, FIXVAL ,0. )
 COVAL(BOTTOM ,F112, FIXVAL ,0. )
 COVAL(BOTTOM ,F113, FIXVAL ,0. )
 COVAL(BOTTOM ,F114, FIXVAL ,0. )
 COVAL(BOTTOM ,F115, FIXVAL ,0. )
 COVAL(BOTTOM ,F116, FIXVAL ,0. )
 COVAL(BOTTOM ,F117, FIXVAL ,0. )
 COVAL(BOTTOM ,F118, FIXVAL ,0. )
 COVAL(BOTTOM ,F119, FIXVAL ,0. )
 COVAL(BOTTOM ,F120, FIXVAL ,0. )
 COVAL(BOTTOM ,F121, FIXVAL ,0. )
 
 PATCH(TOP ,CELL , 1, 1, 9, 9, 1, 1, 1, 10)
 COVAL(TOP ,W1 , FIXVAL ,1. )
 COVAL(TOP ,F1 , FIXVAL ,0. )
 COVAL(TOP ,F2 , FIXVAL ,0. )
 COVAL(TOP ,F3 , FIXVAL ,0. )
 COVAL(TOP ,F4 , FIXVAL ,0. )
 COVAL(TOP ,F5 , FIXVAL ,0. )
 COVAL(TOP ,F6 , FIXVAL ,0. )
 COVAL(TOP ,F7 , FIXVAL ,0. )
 COVAL(TOP ,F8 , FIXVAL ,0. )
 COVAL(TOP ,F9 , FIXVAL ,0. )
 COVAL(TOP ,F10 , FIXVAL ,0. )
 COVAL(TOP ,F11 , FIXVAL ,0. )
 COVAL(TOP ,F12 , FIXVAL ,0. )
 COVAL(TOP ,F13 , FIXVAL ,0. )
 COVAL(TOP ,F14 , FIXVAL ,0. )
 COVAL(TOP ,F15 , FIXVAL ,0. )
 COVAL(TOP ,F16 , FIXVAL ,0. )
 COVAL(TOP ,F17 , FIXVAL ,0. )
 COVAL(TOP ,F18 , FIXVAL ,0. )
 COVAL(TOP ,F19 , FIXVAL ,0. )
 COVAL(TOP ,F20 , FIXVAL ,0. )
 COVAL(TOP ,F21 , FIXVAL ,0. )
 COVAL(TOP ,F22 , FIXVAL ,0. )
 COVAL(TOP ,F23 , FIXVAL ,0. )
 COVAL(TOP ,F24 , FIXVAL ,0. )
 COVAL(TOP ,F25 , FIXVAL ,0. )
 COVAL(TOP ,F26 , FIXVAL ,0. )
 COVAL(TOP ,F27 , FIXVAL ,0. )
 COVAL(TOP ,F28 , FIXVAL ,0. )
 COVAL(TOP ,F29 , FIXVAL ,0. )
 COVAL(TOP ,F30 , FIXVAL ,0. )
 COVAL(TOP ,F31 , FIXVAL ,0. )
 COVAL(TOP ,F32 , FIXVAL ,0. )
 COVAL(TOP ,F33 , FIXVAL ,0. )
 COVAL(TOP ,F34 , FIXVAL ,0. )
 COVAL(TOP ,F35 , FIXVAL ,0. )
 COVAL(TOP ,F36 , FIXVAL ,0. )
 COVAL(TOP ,F37 , FIXVAL ,0. )
 COVAL(TOP ,F38 , FIXVAL ,0. )
 COVAL(TOP ,F39 , FIXVAL ,0. )
 COVAL(TOP ,F40 , FIXVAL ,0. )
 COVAL(TOP ,F41 , FIXVAL ,0. )
 COVAL(TOP ,F42 , FIXVAL ,0. )
 COVAL(TOP ,F43 , FIXVAL ,0. )
 COVAL(TOP ,F44 , FIXVAL ,0. )
 COVAL(TOP ,F45 , FIXVAL ,0. )
 COVAL(TOP ,F46 , FIXVAL ,0. )
 COVAL(TOP ,F47 , FIXVAL ,0. )
 COVAL(TOP ,F48 , FIXVAL ,0. )
 COVAL(TOP ,F49 , FIXVAL ,0. )
 COVAL(TOP ,F50 , FIXVAL ,0. )
 COVAL(TOP ,F51 , FIXVAL ,0. )
 COVAL(TOP ,F52 , FIXVAL ,0. )
 COVAL(TOP ,F53 , FIXVAL ,0. )
 COVAL(TOP ,F54 , FIXVAL ,0. )
 COVAL(TOP ,F55 , FIXVAL ,0. )
 COVAL(TOP ,F56 , FIXVAL ,0. )
 COVAL(TOP ,F57 , FIXVAL ,0. )
 COVAL(TOP ,F58 , FIXVAL ,0. )
 COVAL(TOP ,F59 , FIXVAL ,0. )
 COVAL(TOP ,F60 , FIXVAL ,0. )
 COVAL(TOP ,F61 , FIXVAL ,0. )
 COVAL(TOP ,F62 , FIXVAL ,0. )
 COVAL(TOP ,F63 , FIXVAL ,0. )
 COVAL(TOP ,F64 , FIXVAL ,0. )
 COVAL(TOP ,F65 , FIXVAL ,0. )
 COVAL(TOP ,F66 , FIXVAL ,0. )
 COVAL(TOP ,F67 , FIXVAL ,0. )
 COVAL(TOP ,F68 , FIXVAL ,0. )
 COVAL(TOP ,F69 , FIXVAL ,0. )
 COVAL(TOP ,F70 , FIXVAL ,0. )
 COVAL(TOP ,F71 , FIXVAL ,0. )
 COVAL(TOP ,F72 , FIXVAL ,0. )
 COVAL(TOP ,F73 , FIXVAL ,0. )
 COVAL(TOP ,F74 , FIXVAL ,0. )
 COVAL(TOP ,F75 , FIXVAL ,0. )
 COVAL(TOP ,F76 , FIXVAL ,0. )
 COVAL(TOP ,F77 , FIXVAL ,0. )
 COVAL(TOP ,F78 , FIXVAL ,0. )
 COVAL(TOP ,F79 , FIXVAL ,0. )
 COVAL(TOP ,F80 , FIXVAL ,0. )
 COVAL(TOP ,F81 , FIXVAL ,0. )
 COVAL(TOP ,F82 , FIXVAL ,0. )
 COVAL(TOP ,F83 , FIXVAL ,0. )
 COVAL(TOP ,F84 , FIXVAL ,0. )
 COVAL(TOP ,F85 , FIXVAL ,0. )
 COVAL(TOP ,F86 , FIXVAL ,0. )
 COVAL(TOP ,F87 , FIXVAL ,0. )
 COVAL(TOP ,F88 , FIXVAL ,0. )
 COVAL(TOP ,F89 , FIXVAL ,0. )
 COVAL(TOP ,F90 , FIXVAL ,0. )
 COVAL(TOP ,F91 , FIXVAL ,0. )
 COVAL(TOP ,F92 , FIXVAL ,0. )
 COVAL(TOP ,F93 , FIXVAL ,0. )
 COVAL(TOP ,F94 , FIXVAL ,0. )
 COVAL(TOP ,F95 , FIXVAL ,0. )
 COVAL(TOP ,F96 , FIXVAL ,0. )
 COVAL(TOP ,F97 , FIXVAL ,0. )
 COVAL(TOP ,F98 , FIXVAL ,0. )
 COVAL(TOP ,F99 , FIXVAL ,0. )
 COVAL(TOP ,F100, FIXVAL ,0. )
 COVAL(TOP ,F101, FIXVAL ,0. )
 COVAL(TOP ,F102, FIXVAL ,0. )
 COVAL(TOP ,F103, FIXVAL ,0. )
 COVAL(TOP ,F104, FIXVAL ,0. )
 COVAL(TOP ,F105, FIXVAL ,0. )
 COVAL(TOP ,F106, FIXVAL ,0. )
 COVAL(TOP ,F107, FIXVAL ,0. )
 COVAL(TOP ,F108, FIXVAL ,0. )
 COVAL(TOP ,F109, FIXVAL ,0. )
 COVAL(TOP ,F110, FIXVAL ,0. )
 COVAL(TOP ,F111, FIXVAL ,0. )
 COVAL(TOP ,F112, FIXVAL ,0. )
 COVAL(TOP ,F113, FIXVAL ,0. )
 COVAL(TOP ,F114, FIXVAL ,0. )
 COVAL(TOP ,F115, FIXVAL ,0. )
 COVAL(TOP ,F116, FIXVAL ,1. )
 COVAL(TOP ,F117, FIXVAL ,0. )
 COVAL(TOP ,F118, FIXVAL ,0. )
 COVAL(TOP ,F119, FIXVAL ,0. )
 COVAL(TOP ,F120, FIXVAL ,0. )
 COVAL(TOP ,F121, FIXVAL ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 200 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = F
 RESREF(W1)=1.0E-08 ;RESREF(F1)=0.
 RESREF(F2)=0. ;RESREF(F3)=1.0E-08
 RESREF(F4)=1.0E-08 ;RESREF(F5)=1.0E-08
 RESREF(F6)=1.0E-08 ;RESREF(F7)=1.0E-08
 RESREF(F8)=1.0E-08 ;RESREF(F9)=1.0E-08
 RESREF(F10)=1.0E-08 ;RESREF(F11)=1.0E-08
 RESREF(F12)=1.0E-08 ;RESREF(F13)=1.0E-08
 RESREF(F14)=1.0E-08 ;RESREF(F15)=1.0E-08
 RESREF(F16)=1.0E-08 ;RESREF(F17)=1.0E-08
 RESREF(F18)=1.0E-08 ;RESREF(F19)=1.0E-08
 RESREF(F20)=1.0E-08 ;RESREF(F21)=1.0E-08
 RESREF(F22)=1.0E-08 ;RESREF(F23)=1.0E-08
 RESREF(F24)=1.0E-08 ;RESREF(F25)=1.0E-08
 RESREF(F26)=1.0E-08 ;RESREF(F27)=1.0E-08
 RESREF(F28)=1.0E-08 ;RESREF(F29)=1.0E-08
 RESREF(F30)=1.0E-08 ;RESREF(F31)=1.0E-08
 RESREF(F32)=1.0E-08 ;RESREF(F33)=1.0E-08
 RESREF(F34)=1.0E-08 ;RESREF(F35)=1.0E-08
 RESREF(F36)=1.0E-08 ;RESREF(F37)=1.0E-08
 RESREF(F38)=1.0E-08 ;RESREF(F39)=1.0E-08
 RESREF(F40)=1.0E-08 ;RESREF(F41)=1.0E-08
 RESREF(F42)=1.0E-08 ;RESREF(F43)=1.0E-08
 RESREF(F44)=1.0E-08 ;RESREF(F45)=1.0E-08
 RESREF(F46)=1.0E-08 ;RESREF(F47)=1.0E-08
 RESREF(F48)=1.0E-08 ;RESREF(F49)=1.0E-08
 RESREF(F50)=1.0E-08 ;RESREF(F51)=1.0E-08
 RESREF(F52)=1.0E-08 ;RESREF(F53)=1.0E-08
 RESREF(F54)=1.0E-08 ;RESREF(F55)=1.0E-08
 RESREF(F56)=1.0E-08 ;RESREF(F57)=1.0E-08
 RESREF(F58)=1.0E-08 ;RESREF(F59)=1.0E-08
 RESREF(F60)=1.0E-08 ;RESREF(F61)=1.0E-08
 RESREF(F62)=1.0E-08 ;RESREF(F63)=1.0E-08
 RESREF(F64)=1.0E-08 ;RESREF(F65)=1.0E-08
 RESREF(F66)=1.0E-08 ;RESREF(F67)=1.0E-08
 RESREF(F68)=1.0E-08 ;RESREF(F69)=1.0E-08
 RESREF(F70)=1.0E-08 ;RESREF(F71)=1.0E-08
 RESREF(F72)=1.0E-08 ;RESREF(F73)=1.0E-08
 RESREF(F74)=1.0E-08 ;RESREF(F75)=1.0E-08
 RESREF(F76)=1.0E-08 ;RESREF(F77)=1.0E-08
 RESREF(F78)=1.0E-08 ;RESREF(F79)=1.0E-08
 RESREF(F80)=1.0E-08 ;RESREF(F81)=1.0E-08
 RESREF(F82)=1.0E-08 ;RESREF(F83)=1.0E-08
 RESREF(F84)=1.0E-08 ;RESREF(F85)=1.0E-08
 RESREF(F86)=1.0E-08 ;RESREF(F87)=1.0E-08
 RESREF(F88)=1.0E-08 ;RESREF(F89)=1.0E-08
 RESREF(F90)=1.0E-08 ;RESREF(F91)=1.0E-08
 RESREF(F92)=1.0E-08 ;RESREF(F93)=1.0E-08
 RESREF(F94)=1.0E-08 ;RESREF(F95)=1.0E-08
 RESREF(F96)=1.0E-08 ;RESREF(F97)=1.0E-08
 RESREF(F98)=1.0E-08 ;RESREF(F99)=1.0E-08
 RESREF(F100)=1.0E-08 ;RESREF(F101)=1.0E-08
 RESREF(F102)=1.0E-08 ;RESREF(F103)=1.0E-08
 RESREF(F104)=1.0E-08 ;RESREF(F105)=1.0E-08
 RESREF(F106)=1.0E-08 ;RESREF(F107)=1.0E-08
 RESREF(F108)=1.0E-08 ;RESREF(F109)=1.0E-08
 RESREF(F110)=1.0E-08 ;RESREF(F111)=1.0E-08
 RESREF(F112)=1.0E-08 ;RESREF(F113)=1.0E-08
 RESREF(F114)=1.0E-08 ;RESREF(F115)=1.0E-08
 RESREF(F116)=1.0E-08 ;RESREF(F117)=1.0E-08
 RESREF(F118)=1.0E-08 ;RESREF(F119)=1.0E-08
 RESREF(F120)=1.0E-08 ;RESREF(F121)=1.0E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(W1)=10 ;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(F11)=20
 LITER(F12)=20 ;LITER(F13)=20
 LITER(F14)=20 ;LITER(F15)=20
 LITER(F16)=20 ;LITER(F17)=20
 LITER(F18)=20 ;LITER(F19)=20
 LITER(F20)=20 ;LITER(F21)=20
 LITER(F22)=20 ;LITER(F23)=20
 LITER(F24)=20 ;LITER(F25)=20
 LITER(F26)=20 ;LITER(F27)=20
 LITER(F28)=20 ;LITER(F29)=20
 LITER(F30)=20 ;LITER(F31)=20
 LITER(F32)=20 ;LITER(F33)=20
 LITER(F34)=20 ;LITER(F35)=20
 LITER(F36)=20 ;LITER(F37)=20
 LITER(F38)=20 ;LITER(F39)=20
 LITER(F40)=20 ;LITER(F41)=20
 LITER(F42)=20 ;LITER(F43)=20
 LITER(F44)=20 ;LITER(F45)=20
 LITER(F46)=20 ;LITER(F47)=20
 LITER(F48)=20 ;LITER(F49)=20
 LITER(F50)=20 ;LITER(F51)=20
 LITER(F52)=20 ;LITER(F53)=20
 LITER(F54)=20 ;LITER(F55)=20
 LITER(F56)=20 ;LITER(F57)=20
 LITER(F58)=20 ;LITER(F59)=20
 LITER(F60)=20 ;LITER(F61)=20
 LITER(F62)=20 ;LITER(F63)=20
 LITER(F64)=20 ;LITER(F65)=20
 LITER(F66)=20 ;LITER(F67)=20
 LITER(F68)=20 ;LITER(F69)=20
 LITER(F70)=20 ;LITER(F71)=20
 LITER(F72)=20 ;LITER(F73)=20
 LITER(F74)=20 ;LITER(F75)=20
 LITER(F76)=20 ;LITER(F77)=20
 LITER(F78)=20 ;LITER(F79)=20
 LITER(F80)=20 ;LITER(F81)=20
 LITER(F82)=20 ;LITER(F83)=20
 LITER(F84)=20 ;LITER(F85)=20
 LITER(F86)=20 ;LITER(F87)=20
 LITER(F88)=20 ;LITER(F89)=20
 LITER(F90)=20 ;LITER(F91)=20
 LITER(F92)=20 ;LITER(F93)=20
 LITER(F94)=20 ;LITER(F95)=20
 LITER(F96)=20 ;LITER(F97)=20
 LITER(F98)=20 ;LITER(F99)=20
 LITER(F100)=20 ;LITER(F101)=20
 LITER(F102)=20 ;LITER(F103)=20
 LITER(F104)=20 ;LITER(F105)=20
 LITER(F106)=20 ;LITER(F107)=20
 LITER(F108)=20 ;LITER(F109)=20
 LITER(F110)=20 ;LITER(F111)=20
 LITER(F112)=20 ;LITER(F113)=20
 LITER(F114)=20 ;LITER(F115)=20
 LITER(F116)=20 ;LITER(F117)=20
 LITER(F118)=20 ;LITER(F119)=20
 LITER(F120)=20 ;LITER(F121)=20
 ENDIT(W1)=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(F11)=1.0E-03
 ENDIT(F12)=1.0E-03 ;ENDIT(F13)=1.0E-03
 ENDIT(F14)=1.0E-03 ;ENDIT(F15)=1.0E-03
 ENDIT(F16)=1.0E-03 ;ENDIT(F17)=1.0E-03
 ENDIT(F18)=1.0E-03 ;ENDIT(F19)=1.0E-03
 ENDIT(F20)=1.0E-03 ;ENDIT(F21)=1.0E-03
 ENDIT(F22)=1.0E-03 ;ENDIT(F23)=1.0E-03
 ENDIT(F24)=1.0E-03 ;ENDIT(F25)=1.0E-03
 ENDIT(F26)=1.0E-03 ;ENDIT(F27)=1.0E-03
 ENDIT(F28)=1.0E-03 ;ENDIT(F29)=1.0E-03
 ENDIT(F30)=1.0E-03 ;ENDIT(F31)=1.0E-03
 ENDIT(F32)=1.0E-03 ;ENDIT(F33)=1.0E-03
 ENDIT(F34)=1.0E-03 ;ENDIT(F35)=1.0E-03
 ENDIT(F36)=1.0E-03 ;ENDIT(F37)=1.0E-03
 ENDIT(F38)=1.0E-03 ;ENDIT(F39)=1.0E-03
 ENDIT(F40)=1.0E-03 ;ENDIT(F41)=1.0E-03
 ENDIT(F42)=1.0E-03 ;ENDIT(F43)=1.0E-03
 ENDIT(F44)=1.0E-03 ;ENDIT(F45)=1.0E-03
 ENDIT(F46)=1.0E-03 ;ENDIT(F47)=1.0E-03
 ENDIT(F48)=1.0E-03 ;ENDIT(F49)=1.0E-03
 ENDIT(F50)=1.0E-03 ;ENDIT(F51)=1.0E-03
 ENDIT(F52)=1.0E-03 ;ENDIT(F53)=1.0E-03
 ENDIT(F54)=1.0E-03 ;ENDIT(F55)=1.0E-03
 ENDIT(F56)=1.0E-03 ;ENDIT(F57)=1.0E-03
 ENDIT(F58)=1.0E-03 ;ENDIT(F59)=1.0E-03
 ENDIT(F60)=1.0E-03 ;ENDIT(F61)=1.0E-03
 ENDIT(F62)=1.0E-03 ;ENDIT(F63)=1.0E-03
 ENDIT(F64)=1.0E-03 ;ENDIT(F65)=1.0E-03
 ENDIT(F66)=1.0E-03 ;ENDIT(F67)=1.0E-03
 ENDIT(F68)=1.0E-03 ;ENDIT(F69)=1.0E-03
 ENDIT(F70)=1.0E-03 ;ENDIT(F71)=1.0E-03
 ENDIT(F72)=1.0E-03 ;ENDIT(F73)=1.0E-03
 ENDIT(F74)=1.0E-03 ;ENDIT(F75)=1.0E-03
 ENDIT(F76)=1.0E-03 ;ENDIT(F77)=1.0E-03
 ENDIT(F78)=1.0E-03 ;ENDIT(F79)=1.0E-03
 ENDIT(F80)=1.0E-03 ;ENDIT(F81)=1.0E-03
 ENDIT(F82)=1.0E-03 ;ENDIT(F83)=1.0E-03
 ENDIT(F84)=1.0E-03 ;ENDIT(F85)=1.0E-03
 ENDIT(F86)=1.0E-03 ;ENDIT(F87)=1.0E-03
 ENDIT(F88)=1.0E-03 ;ENDIT(F89)=1.0E-03
 ENDIT(F90)=1.0E-03 ;ENDIT(F91)=1.0E-03
 ENDIT(F92)=1.0E-03 ;ENDIT(F93)=1.0E-03
 ENDIT(F94)=1.0E-03 ;ENDIT(F95)=1.0E-03
 ENDIT(F96)=1.0E-03 ;ENDIT(F97)=1.0E-03
 ENDIT(F98)=1.0E-03 ;ENDIT(F99)=1.0E-03
 ENDIT(F100)=1.0E-03 ;ENDIT(F101)=1.0E-03
 ENDIT(F102)=1.0E-03 ;ENDIT(F103)=1.0E-03
 ENDIT(F104)=1.0E-03 ;ENDIT(F105)=1.0E-03
 ENDIT(F106)=1.0E-03 ;ENDIT(F107)=1.0E-03
 ENDIT(F108)=1.0E-03 ;ENDIT(F109)=1.0E-03
 ENDIT(F110)=1.0E-03 ;ENDIT(F111)=1.0E-03
 ENDIT(F112)=1.0E-03 ;ENDIT(F113)=1.0E-03
 ENDIT(F114)=1.0E-03 ;ENDIT(F115)=1.0E-03
 ENDIT(F116)=1.0E-03 ;ENDIT(F117)=1.0E-03
 ENDIT(F118)=1.0E-03 ;ENDIT(F119)=1.0E-03
 ENDIT(F120)=1.0E-03 ;ENDIT(F121)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(W1,FALSDT,1.)
 RELAX(VIST,LINRLX,1.)
 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(F11,LINRLX,0.5)
 RELAX(F12,LINRLX,0.5)
 RELAX(F13,LINRLX,0.5)
 RELAX(F14,LINRLX,0.5)
 RELAX(F15,LINRLX,0.5)
 RELAX(F16,LINRLX,0.5)
 RELAX(F17,LINRLX,0.5)
 RELAX(F18,LINRLX,0.5)
 RELAX(F19,LINRLX,0.5)
 RELAX(F20,LINRLX,0.5)
 RELAX(F21,LINRLX,0.5)
 RELAX(F22,LINRLX,0.5)
 RELAX(F23,LINRLX,0.5)
 RELAX(F24,LINRLX,0.5)
 RELAX(F25,LINRLX,0.5)
 RELAX(F26,LINRLX,0.5)
 RELAX(F27,LINRLX,0.5)
 RELAX(F28,LINRLX,0.5)
 RELAX(F29,LINRLX,0.5)
 RELAX(F30,LINRLX,0.5)
 RELAX(F31,LINRLX,0.5)
 RELAX(F32,LINRLX,0.5)
 RELAX(F33,LINRLX,0.5)
 RELAX(F34,LINRLX,0.5)
 RELAX(F35,LINRLX,0.5)
 RELAX(F36,LINRLX,0.5)
 RELAX(F37,LINRLX,0.5)
 RELAX(F38,LINRLX,0.5)
 RELAX(F39,LINRLX,0.5)
 RELAX(F40,LINRLX,0.5)
 RELAX(F41,LINRLX,0.5)
 RELAX(F42,LINRLX,0.5)
 RELAX(F43,LINRLX,0.5)
 RELAX(F44,LINRLX,0.5)
 RELAX(F45,LINRLX,0.5)
 RELAX(F46,LINRLX,0.5)
 RELAX(F47,LINRLX,0.5)
 RELAX(F48,LINRLX,0.5)
 RELAX(F49,LINRLX,0.5)
 RELAX(F50,LINRLX,0.5)
 RELAX(F51,LINRLX,0.5)
 RELAX(F52,LINRLX,0.5)
 RELAX(F53,LINRLX,0.5)
 RELAX(F54,LINRLX,0.5)
 RELAX(F55,LINRLX,0.5)
 RELAX(F56,LINRLX,0.5)
 RELAX(F57,LINRLX,0.5)
 RELAX(F58,LINRLX,0.5)
 RELAX(F59,LINRLX,0.5)
 RELAX(F60,LINRLX,0.5)
 RELAX(F61,LINRLX,0.5)
 RELAX(F62,LINRLX,0.5)
 RELAX(F63,LINRLX,0.5)
 RELAX(F64,LINRLX,0.5)
 RELAX(F65,LINRLX,0.5)
 RELAX(F66,LINRLX,0.5)
 RELAX(F67,LINRLX,0.5)
 RELAX(F68,LINRLX,0.5)
 RELAX(F69,LINRLX,0.5)
 RELAX(F70,LINRLX,0.5)
 RELAX(F71,LINRLX,0.5)
 RELAX(F72,LINRLX,0.5)
 RELAX(F73,LINRLX,0.5)
 RELAX(F74,LINRLX,0.5)
 RELAX(F75,LINRLX,0.5)
 RELAX(F76,LINRLX,0.5)
 RELAX(F77,LINRLX,0.5)
 RELAX(F78,LINRLX,0.5)
 RELAX(F79,LINRLX,0.5)
 RELAX(F80,LINRLX,0.5)
 RELAX(F81,LINRLX,0.5)
 RELAX(F82,LINRLX,0.5)
 RELAX(F83,LINRLX,0.5)
 RELAX(F84,LINRLX,0.5)
 RELAX(F85,LINRLX,0.5)
 RELAX(F86,LINRLX,0.5)
 RELAX(F87,LINRLX,0.5)
 RELAX(F88,LINRLX,0.5)
 RELAX(F89,LINRLX,0.5)
 RELAX(F90,LINRLX,0.5)
 RELAX(F91,LINRLX,0.5)
 RELAX(F92,LINRLX,0.5)
 RELAX(F93,LINRLX,0.5)
 RELAX(F94,LINRLX,0.5)
 RELAX(F95,LINRLX,0.5)
 RELAX(F96,LINRLX,0.5)
 RELAX(F97,LINRLX,0.5)
 RELAX(F98,LINRLX,0.5)
 RELAX(F99,LINRLX,0.5)
 RELAX(F100,LINRLX,0.5)
 RELAX(F101,LINRLX,0.5)
 RELAX(F102,LINRLX,0.5)
 RELAX(F103,LINRLX,0.5)
 RELAX(F104,LINRLX,0.5)
 RELAX(F105,LINRLX,0.5)
 RELAX(F106,LINRLX,0.5)
 RELAX(F107,LINRLX,0.5)
 RELAX(F108,LINRLX,0.5)
 RELAX(F109,LINRLX,0.5)
 RELAX(F110,LINRLX,0.5)
 RELAX(F111,LINRLX,0.5)
 RELAX(F112,LINRLX,0.5)
 RELAX(F113,LINRLX,0.5)
 RELAX(F114,LINRLX,0.5)
 RELAX(F115,LINRLX,0.5)
 RELAX(F116,LINRLX,0.5)
 RELAX(F117,LINRLX,0.5)
 RELAX(F118,LINRLX,0.5)
 RELAX(F119,LINRLX,0.5)
 RELAX(F120,LINRLX,0.5)
 RELAX(F121,LINRLX,0.5)
 RELAX(FRSM,LINRLX,1.)
 RELAX(NETM,LINRLX,1.)
 RELAX(PGRD,LINRLX,1.)
 RELAX(FLSM,LINRLX,1.)
 RELAX(RMOA,LINRLX,1.)
 RELAX(RMSA,LINRLX,1.)
 RELAX(RMSO,LINRLX,1.)
 RELAX(AVER,LINRLX,1.)
 RELAX(AVEF,LINRLX,1.)
 RELAX(RATE,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(VIST)=1.0E+10 ;VARMIN(VIST)=-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(F11)=1. ;VARMIN(F11)=0.
 VARMAX(F12)=1. ;VARMIN(F12)=0.
 VARMAX(F13)=1. ;VARMIN(F13)=0.
 VARMAX(F14)=1. ;VARMIN(F14)=0.
 VARMAX(F15)=1. ;VARMIN(F15)=0.
 VARMAX(F16)=1. ;VARMIN(F16)=0.
 VARMAX(F17)=1. ;VARMIN(F17)=0.
 VARMAX(F18)=1. ;VARMIN(F18)=0.
 VARMAX(F19)=1. ;VARMIN(F19)=0.
 VARMAX(F20)=1. ;VARMIN(F20)=0.
 VARMAX(F21)=1. ;VARMIN(F21)=0.
 VARMAX(F22)=1. ;VARMIN(F22)=0.
 VARMAX(F23)=1. ;VARMIN(F23)=0.
 VARMAX(F24)=1. ;VARMIN(F24)=0.
 VARMAX(F25)=1. ;VARMIN(F25)=0.
 VARMAX(F26)=1. ;VARMIN(F26)=0.
 VARMAX(F27)=1. ;VARMIN(F27)=0.
 VARMAX(F28)=1. ;VARMIN(F28)=0.
 VARMAX(F29)=1. ;VARMIN(F29)=0.
 VARMAX(F30)=1. ;VARMIN(F30)=0.
 VARMAX(F31)=1. ;VARMIN(F31)=0.
 VARMAX(F32)=1. ;VARMIN(F32)=0.
 VARMAX(F33)=1. ;VARMIN(F33)=0.
 VARMAX(F34)=1. ;VARMIN(F34)=0.
 VARMAX(F35)=1. ;VARMIN(F35)=0.
 VARMAX(F36)=1. ;VARMIN(F36)=0.
 VARMAX(F37)=1. ;VARMIN(F37)=0.
 VARMAX(F38)=1. ;VARMIN(F38)=0.
 VARMAX(F39)=1. ;VARMIN(F39)=0.
 VARMAX(F40)=1. ;VARMIN(F40)=0.
 VARMAX(F41)=1. ;VARMIN(F41)=0.
 VARMAX(F42)=1. ;VARMIN(F42)=0.
 VARMAX(F43)=1. ;VARMIN(F43)=0.
 VARMAX(F44)=1. ;VARMIN(F44)=0.
 VARMAX(F45)=1. ;VARMIN(F45)=0.
 VARMAX(F46)=1. ;VARMIN(F46)=0.
 VARMAX(F47)=1. ;VARMIN(F47)=0.
 VARMAX(F48)=1. ;VARMIN(F48)=0.
 VARMAX(F49)=1. ;VARMIN(F49)=0.
 VARMAX(F50)=1. ;VARMIN(F50)=0.
 VARMAX(F51)=1. ;VARMIN(F51)=0.
 VARMAX(F52)=1. ;VARMIN(F52)=0.
 VARMAX(F53)=1. ;VARMIN(F53)=0.
 VARMAX(F54)=1. ;VARMIN(F54)=0.
 VARMAX(F55)=1. ;VARMIN(F55)=0.
 VARMAX(F56)=1. ;VARMIN(F56)=0.
 VARMAX(F57)=1. ;VARMIN(F57)=0.
 VARMAX(F58)=1. ;VARMIN(F58)=0.
 VARMAX(F59)=1. ;VARMIN(F59)=0.
 VARMAX(F60)=1. ;VARMIN(F60)=0.
 VARMAX(F61)=1. ;VARMIN(F61)=0.
 VARMAX(F62)=1. ;VARMIN(F62)=0.
 VARMAX(F63)=1. ;VARMIN(F63)=0.
 VARMAX(F64)=1. ;VARMIN(F64)=0.
 VARMAX(F65)=1. ;VARMIN(F65)=0.
 VARMAX(F66)=1. ;VARMIN(F66)=0.
 VARMAX(F67)=1. ;VARMIN(F67)=0.
 VARMAX(F68)=1. ;VARMIN(F68)=0.
 VARMAX(F69)=1. ;VARMIN(F69)=0.
 VARMAX(F70)=1. ;VARMIN(F70)=0.
 VARMAX(F71)=1. ;VARMIN(F71)=0.
 VARMAX(F72)=1. ;VARMIN(F72)=0.
 VARMAX(F73)=1. ;VARMIN(F73)=0.
 VARMAX(F74)=1. ;VARMIN(F74)=0.
 VARMAX(F75)=1. ;VARMIN(F75)=0.
 VARMAX(F76)=1. ;VARMIN(F76)=0.
 VARMAX(F77)=1. ;VARMIN(F77)=0.
 VARMAX(F78)=1. ;VARMIN(F78)=0.
 VARMAX(F79)=1. ;VARMIN(F79)=0.
 VARMAX(F80)=1. ;VARMIN(F80)=0.
 VARMAX(F81)=1. ;VARMIN(F81)=0.
 VARMAX(F82)=1. ;VARMIN(F82)=0.
 VARMAX(F83)=1. ;VARMIN(F83)=0.
 VARMAX(F84)=1. ;VARMIN(F84)=0.
 VARMAX(F85)=1. ;VARMIN(F85)=0.
 VARMAX(F86)=1. ;VARMIN(F86)=0.
 VARMAX(F87)=1. ;VARMIN(F87)=0.
 VARMAX(F88)=1. ;VARMIN(F88)=0.
 VARMAX(F89)=1. ;VARMIN(F89)=0.
 VARMAX(F90)=1. ;VARMIN(F90)=0.
 VARMAX(F91)=1. ;VARMIN(F91)=0.
 VARMAX(F92)=1. ;VARMIN(F92)=0.
 VARMAX(F93)=1. ;VARMIN(F93)=0.
 VARMAX(F94)=1. ;VARMIN(F94)=0.
 VARMAX(F95)=1. ;VARMIN(F95)=0.
 VARMAX(F96)=1. ;VARMIN(F96)=0.
 VARMAX(F97)=1. ;VARMIN(F97)=0.
 VARMAX(F98)=1. ;VARMIN(F98)=0.
 VARMAX(F99)=1. ;VARMIN(F99)=0.
 VARMAX(F100)=1. ;VARMIN(F100)=0.
 VARMAX(F101)=1. ;VARMIN(F101)=0.
 VARMAX(F102)=1. ;VARMIN(F102)=0.
 VARMAX(F103)=1. ;VARMIN(F103)=0.
 VARMAX(F104)=1. ;VARMIN(F104)=0.
 VARMAX(F105)=1. ;VARMIN(F105)=0.
 VARMAX(F106)=1. ;VARMIN(F106)=0.
 VARMAX(F107)=1. ;VARMIN(F107)=0.
 VARMAX(F108)=1. ;VARMIN(F108)=0.
 VARMAX(F109)=1. ;VARMIN(F109)=0.
 VARMAX(F110)=1. ;VARMIN(F110)=0.
 VARMAX(F111)=1. ;VARMIN(F111)=0.
 VARMAX(F112)=1. ;VARMIN(F112)=0.
 VARMAX(F113)=1. ;VARMIN(F113)=0.
 VARMAX(F114)=1. ;VARMIN(F114)=0.
 VARMAX(F115)=1. ;VARMIN(F115)=0.
 VARMAX(F116)=1. ;VARMIN(F116)=0.
 VARMAX(F117)=1. ;VARMIN(F117)=0.
 VARMAX(F118)=1. ;VARMIN(F118)=0.
 VARMAX(F119)=1. ;VARMIN(F119)=0.
 VARMAX(F120)=1. ;VARMIN(F120)=0.
 VARMAX(F121)=1. ;VARMIN(F121)=0.
 VARMAX(FRSM)=1.0E+10 ;VARMIN(FRSM)=-1.0E+10
 VARMAX(NETM)=1.0E+10 ;VARMIN(NETM)=-1.0E+10
 VARMAX(PGRD)=1.0E+10 ;VARMIN(PGRD)=-1.0E+10
 VARMAX(FLSM)=1.0E+10 ;VARMIN(FLSM)=-1.0E+10
 VARMAX(RMOA)=1.0E+10 ;VARMIN(RMOA)=-1.0E+10
 VARMAX(RMSA)=1.0E+10 ;VARMIN(RMSA)=-1.0E+10
 VARMAX(RMSO)=1.0E+10 ;VARMIN(RMSO)=-1.0E+10
 VARMAX(AVER)=1. ;VARMIN(AVER)=0.
 VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,MFM,MFMMOD,C,MFM)
 SPEDAT(SET,MFM,NFLUIDS,I,121)
 SPEDAT(SET,MFM,NFLF,I,11)
 SPEDAT(SET,MFM,NFLR,I,11)
 SPEDAT(SET,MFM,CONMIX,R,100.)
 SPEDAT(SET,REALREAD,2DPOP,C,-5&5&0&10)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = T ;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(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(VIST,N,N,N,N,N,N)
 OUTPUT(F1,Y,Y,Y,Y,Y,Y)
 OUTPUT(F2,Y,Y,N,N,Y,Y)
 OUTPUT(F3,Y,Y,N,N,Y,Y)
 OUTPUT(F4,Y,Y,N,N,Y,Y)
 OUTPUT(F5,Y,Y,N,N,Y,Y)
 OUTPUT(F6,Y,Y,N,N,Y,Y)
 OUTPUT(F7,Y,Y,N,N,Y,Y)
 OUTPUT(F8,Y,Y,N,N,Y,Y)
 OUTPUT(F9,Y,Y,N,N,Y,Y)
 OUTPUT(F10,Y,Y,N,N,Y,Y)
 OUTPUT(F11,Y,Y,N,N,Y,Y)
 OUTPUT(F12,Y,Y,N,N,Y,Y)
 OUTPUT(F13,Y,Y,N,N,Y,Y)
 OUTPUT(F14,Y,Y,N,N,Y,Y)
 OUTPUT(F15,Y,Y,N,N,Y,Y)
 OUTPUT(F16,Y,Y,N,N,Y,Y)
 OUTPUT(F17,Y,Y,N,N,Y,Y)
 OUTPUT(F18,Y,Y,N,N,Y,Y)
 OUTPUT(F19,Y,Y,N,N,Y,Y)
 OUTPUT(F20,Y,Y,N,N,Y,Y)
 OUTPUT(F21,Y,Y,N,N,Y,Y)
 OUTPUT(F22,Y,Y,N,N,Y,Y)
 OUTPUT(F23,Y,Y,N,N,Y,Y)
 OUTPUT(F24,Y,Y,N,N,Y,Y)
 OUTPUT(F25,Y,Y,N,N,Y,Y)
 OUTPUT(F26,Y,Y,N,N,Y,Y)
 OUTPUT(F27,Y,Y,N,N,Y,Y)
 OUTPUT(F28,Y,Y,N,N,Y,Y)
 OUTPUT(F29,Y,Y,N,N,Y,Y)
 OUTPUT(F30,Y,Y,N,N,Y,Y)
 OUTPUT(F31,Y,Y,N,N,Y,Y)
 OUTPUT(F32,Y,Y,N,N,Y,Y)
 OUTPUT(F33,Y,Y,N,N,Y,Y)
 OUTPUT(F34,Y,Y,N,N,Y,Y)
 OUTPUT(F35,Y,Y,N,N,Y,Y)
 OUTPUT(F36,Y,Y,N,N,Y,Y)
 OUTPUT(F37,Y,Y,N,N,Y,Y)
 OUTPUT(F38,Y,Y,N,N,Y,Y)
 OUTPUT(F39,Y,Y,N,N,Y,Y)
 OUTPUT(F40,Y,Y,N,N,Y,Y)
 OUTPUT(F41,Y,Y,N,N,Y,Y)
 OUTPUT(F42,Y,Y,N,N,Y,Y)
 OUTPUT(F43,Y,Y,N,N,Y,Y)
 OUTPUT(F44,Y,Y,N,N,Y,Y)
 OUTPUT(F45,Y,Y,N,N,Y,Y)
 OUTPUT(F46,Y,Y,N,N,Y,Y)
 OUTPUT(F47,Y,Y,N,N,Y,Y)
 OUTPUT(F48,Y,Y,N,N,Y,Y)
 OUTPUT(F49,Y,Y,N,N,Y,Y)
 OUTPUT(F50,Y,Y,N,N,Y,Y)
 OUTPUT(F51,Y,Y,N,N,Y,Y)
 OUTPUT(F52,Y,Y,N,N,Y,Y)
 OUTPUT(F53,Y,Y,N,N,Y,Y)
 OUTPUT(F54,Y,Y,N,N,Y,Y)
 OUTPUT(F55,Y,Y,N,N,Y,Y)
 OUTPUT(F56,Y,Y,N,N,Y,Y)
 OUTPUT(F57,Y,Y,N,N,Y,Y)
 OUTPUT(F58,Y,Y,N,N,Y,Y)
 OUTPUT(F59,Y,Y,N,N,Y,Y)
 OUTPUT(F60,Y,Y,N,N,Y,Y)
 OUTPUT(F61,Y,Y,N,N,Y,Y)
 OUTPUT(F62,Y,Y,N,N,Y,Y)
 OUTPUT(F63,Y,Y,N,N,Y,Y)
 OUTPUT(F64,Y,Y,N,N,Y,Y)
 OUTPUT(F65,Y,Y,N,N,Y,Y)
 OUTPUT(F66,Y,Y,N,N,Y,Y)
 OUTPUT(F67,Y,Y,N,N,Y,Y)
 OUTPUT(F68,Y,Y,N,N,Y,Y)
 OUTPUT(F69,Y,Y,N,N,Y,Y)
 OUTPUT(F70,Y,Y,N,N,Y,Y)
 OUTPUT(F71,Y,Y,N,N,Y,Y)
 OUTPUT(F72,Y,Y,N,N,Y,Y)
 OUTPUT(F73,Y,Y,N,N,Y,Y)
 OUTPUT(F74,Y,Y,N,N,Y,Y)
 OUTPUT(F75,Y,Y,N,N,Y,Y)
 OUTPUT(F76,Y,Y,N,N,Y,Y)
 OUTPUT(F77,Y,Y,N,N,Y,Y)
 OUTPUT(F78,Y,Y,N,N,Y,Y)
 OUTPUT(F79,Y,Y,N,N,Y,Y)
 OUTPUT(F80,Y,Y,N,N,Y,Y)
 OUTPUT(F81,Y,Y,N,N,Y,Y)
 OUTPUT(F82,Y,Y,N,N,Y,Y)
 OUTPUT(F83,Y,Y,N,N,Y,Y)
 OUTPUT(F84,Y,Y,N,N,Y,Y)
 OUTPUT(F85,Y,Y,N,N,Y,Y)
 OUTPUT(F86,Y,Y,N,N,Y,Y)
 OUTPUT(F87,Y,Y,N,N,Y,Y)
 OUTPUT(F88,Y,Y,N,N,Y,Y)
 OUTPUT(F89,Y,Y,N,N,Y,Y)
 OUTPUT(F90,Y,Y,N,N,Y,Y)
 OUTPUT(F91,Y,Y,N,N,Y,Y)
 OUTPUT(F92,Y,Y,N,N,Y,Y)
 OUTPUT(F93,Y,Y,N,N,Y,Y)
 OUTPUT(F94,Y,Y,N,N,Y,Y)
 OUTPUT(F95,Y,Y,N,N,Y,Y)
 OUTPUT(F96,Y,Y,N,N,Y,Y)
 OUTPUT(F97,Y,Y,N,N,Y,Y)
 OUTPUT(F98,Y,Y,N,N,Y,Y)
 OUTPUT(F99,Y,Y,N,N,Y,Y)
 OUTPUT(F100,Y,Y,N,N,Y,Y)
 OUTPUT(F101,Y,Y,N,N,Y,Y)
 OUTPUT(F102,Y,Y,N,N,Y,Y)
 OUTPUT(F103,Y,Y,N,N,Y,Y)
 OUTPUT(F104,Y,Y,N,N,Y,Y)
 OUTPUT(F105,Y,Y,N,N,Y,Y)
 OUTPUT(F106,Y,Y,N,N,Y,Y)
 OUTPUT(F107,Y,Y,N,N,Y,Y)
 OUTPUT(F108,Y,Y,N,N,Y,Y)
 OUTPUT(F109,Y,Y,N,N,Y,Y)
 OUTPUT(F110,Y,Y,N,N,Y,Y)
 OUTPUT(F111,Y,Y,N,N,Y,Y)
 OUTPUT(F112,Y,Y,N,N,Y,Y)
 OUTPUT(F113,Y,Y,N,N,Y,Y)
 OUTPUT(F114,Y,Y,N,N,Y,Y)
 OUTPUT(F115,Y,Y,N,N,Y,Y)
 OUTPUT(F116,Y,Y,N,N,Y,Y)
 OUTPUT(F117,Y,Y,N,N,Y,Y)
 OUTPUT(F118,Y,Y,N,N,Y,Y)
 OUTPUT(F119,Y,Y,N,N,Y,Y)
 OUTPUT(F120,Y,Y,N,N,Y,Y)
 OUTPUT(F121,Y,Y,N,N,Y,Y)
 OUTPUT(FRSM,Y,N,Y,N,N,N)
 OUTPUT(NETM,Y,N,Y,N,N,N)
 OUTPUT(PGRD,Y,N,Y,N,N,N)
 OUTPUT(FLSM,Y,N,Y,N,N,N)
 OUTPUT(RMOA,Y,N,Y,N,N,N)
 OUTPUT(RMSA,Y,N,Y,N,N,N)
 OUTPUT(RMSO,Y,N,Y,N,N,N)
 OUTPUT(AVER,Y,N,Y,N,N,N)
 OUTPUT(AVEF,Y,N,Y,N,N,N)
 OUTPUT(RATE,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 5 ;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
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 1 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.2
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(Y ,PROFIL, 1, 1, 1, 9, 1, 1, 1, 1)
 
 PATCH(PLOT ,PROFIL, 1, 1, 1, 9, 1, 1, 1, 10)
 PLOT(PLOT ,W1 ,0. ,0. )
 PLOT(PLOT ,AVER,0. ,0. )
 PLOT(PLOT ,AVEF,0. ,0. )
 
 PATCH(PLOT2 ,PROFIL, 1, 1, 1, 9, 1, 1, 1, 10)
 PLOT(PLOT2 ,F3 ,0. ,0. )
 PLOT(PLOT2 ,F8 ,0. ,0. )
 PLOT(PLOT2 ,F13 ,0. ,0. )
 PLOT(PLOT2 ,F18 ,0. ,0. )
 PLOT(PLOT2 ,F23 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
STOP