TALK=T;RUN(1,1)
  DISPLAY
  The purpose of creation of this case is a demonstration
  of transfer objects feature.

  This example explores the unsteady distribution of pollution
  on a reasonably large region of ground.

  The pollution source is set in the west-low corner and
  acts only at the first time step. Further the pollution
  cloud will be distributed by wind.

  Use of transfer objects permits the calculation of unidirectional
  flow on small parts one after another. The solved area in this
  example is divided into four parts. The solution in each part is
  made in separate runs.

               ---------------------
               !         !         !
               !         !         !
               ! 3rd run ! 4th run !
               !         !         !
       W1      !         !         !
       -->     ----------+----------
               !         !         !
               !         !         !
               ! 1st run ! 2nd run !
               !         !         !
        X !    !         !         !
          !    ---------------------
          !              ^
          !----- Z      /!\
                         !
                         !U1

  The HTR$
O1.. and ETRO1.. transfer objects on the high
  and east boundaries are formed at the end of each time step of the
  first run by means of two '(export' In-Form statements. They store
  the values of outlet mass flux and concentration on these
  boundaries in order to transfer them to the second and third runs.

  The second run reads the information at the low boundary from the
  HTRO1..$
 objects at the start of each time step by means of '(import
  In-Form statements and at the end of each time step dumps it at
  the east boundary in the ETRO2 object.

  The third run reads from the ETRO1 object and writes from the
  HTRO3 object.

  The fourth run reads the information from HTRO3 and ETRO2 
  import transfer objects at the low and west boundaries.

  The wind profile at the inlet boundaries is set by means of
  In-For$
m statements as a logarithmic velocity profile.

  The ground relief (HIG variable) is calculated by an
  I$
n-Form formula.

  The MARK variable defined by In-Form is used for the
  image of the ground relief in Photon.

  The ground roughness is simulated by a change in air density
  with height in an atmospheric layer. The density of the air is ca$
lculated
  from a b$
arometric formula by means of In-Form.

  The fifth and last run simulates the flow in the whole region
  without partitioning. It will be useful to compare with the previ$
ous runs.

  The Q1 contains PHOTON USE commands
  ENDDIS
  PHOTON USE
  p
  110
  1 5 1;;

  VI -1 1 1
  gr ou y 1
  SURF MARK Y .99
  SURF MARK X .99
  SURF MARK Z .99
  msg First run: ground geometry
  pause
  p
  110;;;

  vi y
  gr ou y 1
  con rho1 y m fi;.001
  msg density contours
  pause
  con cl;red
  ve y 1 sh
  msg velocity vectors
  pause
  con cl;red
  set con scale range on
  con conc y 1 fi;0 0.01;.001
  set con scale range off
  msg concentration contours
  pause
  p
  210
  1 5 1;;

  VI -1 1 1
  gr ou y 1
  SURF MARK Y .99
  SURF MARK X .99
  SURF MARK Z .99
  msg Second run: ground geometry
  pause
  p
  210;;;

  vi y
  gr ou y 1
  con rho1 y m fi;.001
  msg density contours
  pause
  con cl;red
  ve y 1 sh
  msg velocity vectors
  pause
  con cl;red
  set con scale range on
  con conc y 1 fi;0 0.01;.001
  set con scale range off
  msg concentration contours
  pause
  p
  310
  1 5 1;;

  VI -1 1 1
  gr ou y 1
  SURF MARK Y .99
  SURF MARK X .99
  SURF MARK Z .99
  msg Third run: ground geometry
  pause
  p
  310;;;

  vi y
  gr ou y 1
  con rho1 y m fi;.001
  msg density contours
  pause
  con cl;red
  ve y 1 sh
  msg velocity vectors
  pause
  con cl;red
  set con scale range on
  con conc y 1 fi;0 0.01;.001
  set con scale range off
  msg concentration contours
  pause
  p
  410
  1 5 1;;

  VI -1 1 1
  gr ou y 1
  SURF MARK Y .99
  SURF MARK X .99
  SURF MARK Z .99
  msg Fourth run: ground geometry
  pause
  p
  410;;;

  vi y
  gr ou y 1
  con rho1 y m fi;.001
  msg density contours
  pause
  con cl;red
  ve y 1 sh
  msg velocity vectors
  pause
  con cl;red
  set con scale range on
  con conc y 1 fi;0 0.01;.001
  set con scale range off
  msg concentration contours
  pause
  p
  510
  1 5 1;;

  VI -1 1 1
  gr ou y 1
  SURF MARK Y .99
  SURF MARK X .99
  SURF MARK Z .99
  msg First run: ground geometry
  pause
  p
  510;;;

  vi y
  gr ou y 1
  con rho1 y m fi;.001
  msg density contours
  pause
  con cl;red
  ve y 1 sh
  msg velocity vectors
  pause
  con cl;red
  set con scale range on
  con conc y 1 fi;0 0.01;.001
  set con scale range off
  msg concentration contours
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(ATMOSPHERIC B.L., first run             )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 859
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =5000.
 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 = 20
 XULAST =5000.
 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 = 20
 YVLAST =1000.
 YFRAC(1)=2.5E-03 ;YFRAC(2)=0.01
 YFRAC(3)=0.0225 ;YFRAC(4)=0.04
 YFRAC(5)=0.0625 ;YFRAC(6)=0.09
 YFRAC(7)=0.1225 ;YFRAC(8)=0.16
 YFRAC(9)=0.2025 ;YFRAC(10)=0.25
 YFRAC(11)=0.3025 ;YFRAC(12)=0.36
 YFRAC(13)=0.4225 ;YFRAC(14)=0.49
 YFRAC(15)=0.5625 ;YFRAC(16)=0.64
 YFRAC(17)=0.7225 ;YFRAC(18)=0.81
 YFRAC(19)=0.9025 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =5000.
 ZFRAC(1)=0.05 ;ZFRAC(2)=0.1
 ZFRAC(3)=0.15 ;ZFRAC(4)=0.2
 ZFRAC(5)=0.25 ;ZFRAC(6)=0.3
 ZFRAC(7)=0.35 ;ZFRAC(8)=0.4
 ZFRAC(9)=0.45 ;ZFRAC(10)=0.5
 ZFRAC(11)=0.55 ;ZFRAC(12)=0.6
 ZFRAC(13)=0.65 ;ZFRAC(14)=0.7
 ZFRAC(15)=0.75 ;ZFRAC(16)=0.8
 ZFRAC(17)=0.85 ;ZFRAC(18)=0.9
 ZFRAC(19)=0.95 ;ZFRAC(20)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(12)=KE ;NAME(13)=EP
 NAME(144)=RHO1 ;NAME(145)=MARK
 NAME(146)=HIG ;NAME(147)=EL1
 NAME(148)=ENUT ;NAME(149)=EPKE
 NAME(150)=CONC
    * 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(W1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 SOLUTN(MARK,Y,N,N,N,N,Y)
 SOLUTN(HIG,Y,N,N,N,N,Y)
 SOLUTN(EL1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,Y)
 SOLUTN(CONC,Y,Y,N,N,N,Y)
 DEN1 = 144
 VIST = 148
 LEN1 = 147
 ************************************************************
  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(W1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(CONC,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENT = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1.189 ;TMP1 =0. ;EL1 = GRND4
 TSURR =0. ;TEMP0 =273. ;PRESS0 =1.0E+05
 DVO1DT =3.41E-03 ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL =1.544E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 0
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(CONC)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(CONC)=1.
 CP1 =1005. ;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(W1)=1.0E-10
 FIINIT(KE)=1.0E-10 ;FIINIT(EP)=1.0E-10
 FIINIT(RHO1)=1.0E-10 ;FIINIT(MARK)=1.0E-10
 FIINIT(HIG)=1.0E-10 ;FIINIT(EL1)=1.0E-10
 FIINIT(ENUT)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
 FIINIT(CONC)=0.
 
 PATCH(INIT ,INIVAL, 1, 20, 1, 20, 1, 20, 1, 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, 10)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(HOUT ,HIGH , 1, 20, 1, 20, 20, 20, 1, 10)
 COVAL(HOUT ,P1 ,1. ,0. )
 
 PATCH(EOUT ,EAST , 20, 20, 1, 20, 1, 20, 1, 10)
 COVAL(EOUT ,P1 ,1. ,0. )
 
 PATCH(NOUT ,NORTH , 1, 20, 20, 20, 1, 20, 1, 10)
 COVAL(NOUT ,P1 ,1. ,0. )
 
 PATCH(SOURCE ,CELL , 1, 7, 1, 1, 1, 4, 1, 1)
 COVAL(SOURCE ,CONC, FIXFLU ,1.0E+05 )
 
 PATCH(GROUND ,SWALL , 1, 20, 1, 1, 1, 20, 1, 10)
 COVAL(GROUND ,U1 , GRND5 ,0. )
 COVAL(GROUND ,W1 , GRND5 ,0. )
 COVAL(GROUND ,KE , GRND5 , GRND5 )
 COVAL(GROUND ,EP , GRND5 , GRND5 )
 
 PATCH(LINLET ,LOW , 1, 20, 1, 20, 1, 1, 1, 10)
 COVAL(LINLET ,P1 ,In-Form:source - see Grp 19)
 COVAL(LINLET ,U1 ,In-Form:source - see Grp 19)
 COVAL(LINLET ,W1 ,In-Form:source - see Grp 19)
 COVAL(LINLET ,KE ,In-Form:source - see Grp 19)
 COVAL(LINLET ,EP ,In-Form:source - see Grp 19)
 
 PATCH(WINLET ,WEST , 1, 1, 1, 20, 1, 20, 1, 10)
 COVAL(WINLET ,P1 ,In-Form:source - see Grp 19)
 COVAL(WINLET ,U1 ,In-Form:source - see Grp 19)
 COVAL(WINLET ,W1 ,In-Form:source - see Grp 19)
 COVAL(WINLET ,KE ,In-Form:source - see Grp 19)
 COVAL(WINLET ,EP ,In-Form:source - see Grp 19)
 
 PATCH(HPAT1 ,HIGH , 1, 20, 1, 20, 20, 20, 1, 1)
 
 PATCH(EPAT1 ,EAST , 20, 20, 1, 20, 1, 20, 1, 1)
 
 PATCH(HPAT2 ,HIGH , 1, 20, 1, 20, 20, 20, 2, 2)
 
 PATCH(EPAT2 ,EAST , 20, 20, 1, 20, 1, 20, 2, 2)
 
 PATCH(HPAT3 ,HIGH , 1, 20, 1, 20, 20, 20, 3, 3)
 
 PATCH(EPAT3 ,EAST , 20, 20, 1, 20, 1, 20, 3, 3)
 
 PATCH(HPAT4 ,HIGH , 1, 20, 1, 20, 20, 20, 4, 4)
 
 PATCH(EPAT4 ,EAST , 20, 20, 1, 20, 1, 20, 4, 4)
 
 PATCH(HPAT5 ,HIGH , 1, 20, 1, 20, 20, 20, 5, 5)
 
 PATCH(EPAT5 ,EAST , 20, 20, 1, 20, 1, 20, 5, 5)
 
 PATCH(HPAT6 ,HIGH , 1, 20, 1, 20, 20, 20, 6, 6)
 
 PATCH(EPAT6 ,EAST , 20, 20, 1, 20, 1, 20, 6, 6)
 
 PATCH(HPAT7 ,HIGH , 1, 20, 1, 20, 20, 20, 7, 7)
 
 PATCH(EPAT7 ,EAST , 20, 20, 1, 20, 1, 20, 7, 7)
 
 PATCH(HPAT8 ,HIGH , 1, 20, 1, 20, 20, 20, 8, 8)
 
 PATCH(EPAT8 ,EAST , 20, 20, 1, 20, 1, 20, 8, 8)
 
 PATCH(HPAT9 ,HIGH , 1, 20, 1, 20, 20, 20, 9, 9)
 
 PATCH(EPAT9 ,EAST , 20, 20, 1, 20, 1, 20, 9, 9)
 
 PATCH(HPAT10 ,HIGH , 1, 20, 1, 20, 20, 20, 10, 10)
 
 PATCH(EPAT10 ,EAST , 20, 20, 1, 20, 1, 20, 10, 10)
 XCYCLE = F
 WALLA =0.02 ;WALLB =0.
 EGWF = T
 WALLCO = GRND5
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 50 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=20 ;LITER(EP)=20
 LITER(CONC)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(CONC)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(RHO1,LINRLX,1.)
 RELAX(MARK,LINRLX,1.)
 RELAX(HIG,LINRLX,1.)
 RELAX(EL1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(CONC,FALSDT,1.0E+09)
 KELIN = 0
 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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10
 VARMAX(HIG)=1.0E+10 ;VARMIN(HIG)=-1.0E+10
 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(CONC)=1.0E+10 ;VARMIN(CONC)=0.
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 CONWIZ = T
 ISG52 = 2
 ISG62 = 1
 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0)
 SPEDAT(SET,SOURCE,R1!LINLET,C,=0.733131*YG^0.21)
 SPEDAT(SET,SOURCE,U1!LINLET,C,=0.616595*YG^0.21!ONLYMS)
 SPEDAT(SET,SOURCE,W1!LINLET,C,=0.616595*YG^0.21!ONLYMS)
 SPEDAT(SET,SOURCE,KE!LINLET,C,=0.029928!ONLYMS)
 SPEDAT(SET,SOURCE,EP!LINLET,C,=2.074951E-03/(YG-0.)!ONLYMS)
 SPEDAT(SET,SOURCE,R1!WINLET,C,=0.733131*YG^0.21)
 SPEDAT(SET,SOURCE,U1!WINLET,C,=0.616595*YG^0.21!ONLYMS)
 SPEDAT(SET,SOURCE,W1!WINLET,C,=0.616595*YG^0.21!ONLYMS)
 SPEDAT(SET,SOURCE,KE!WINLET,C,=0.029928!ONLYMS)
 SPEDAT(SET,SOURCE,EP!WINLET,C,=2.074951E-03/(YG-0.)!ONLYMS)
 SPEDAT(SET,INITIAL,U1!INIT,C,=0.616595*YG^0.21)
 SPEDAT(SET,INITIAL,W1!INIT,C,=0.616595*YG^0.21)
 SPEDAT(SET,INITIAL,KE!INIT,C,=0.029928)
 SPEDAT(SET,INITIAL,EP!INIT,C,=2.074951E-03/(YG-0.))
 SPEDAT(SET,STORED,HIG!GROUND,C,=100.*(1-(0.+XG)/1.0E+04)*(1-(0.+Z$)
 SPEDAT(SET,STORED,HIG!GROUND,C,G)/1.0E+04)*(2.+SIN(6*(0.+XG)/1.0E$)
 SPEDAT(SET,STORED,HIG!GROUND,C,+04)+SIN(24*(0.+ZG)/1.0E+04))!TSTS$)
 SPEDAT(SET,STORED,HIG!GROUND,C,TR)
 SPEDAT(SET,STORED,MARK,C,=1!IF(YG.GT.HIG[&1])!TSTSTR)
 SPEDAT(SET,PROPERTY,RHO1,C,=1.189*EXP((YG+HIG[&1])/(-8000)))
 SPEDAT(SET,EXPORT,HPAT1,C,HTRO11)
 SPEDAT(SET,EXPORT,EPAT1,C,ETRO11)
 SPEDAT(SET,EXPORT,HPAT2,C,HTRO12)
 SPEDAT(SET,EXPORT,EPAT2,C,ETRO12)
 SPEDAT(SET,EXPORT,HPAT3,C,HTRO13)
 SPEDAT(SET,EXPORT,EPAT3,C,ETRO13)
 SPEDAT(SET,EXPORT,HPAT4,C,HTRO14)
 SPEDAT(SET,EXPORT,EPAT4,C,ETRO14)
 SPEDAT(SET,EXPORT,HPAT5,C,HTRO15)
 SPEDAT(SET,EXPORT,EPAT5,C,ETRO15)
 SPEDAT(SET,EXPORT,HPAT6,C,HTRO16)
 SPEDAT(SET,EXPORT,EPAT6,C,ETRO16)
 SPEDAT(SET,EXPORT,HPAT7,C,HTRO17)
 SPEDAT(SET,EXPORT,EPAT7,C,ETRO17)
 SPEDAT(SET,EXPORT,HPAT8,C,HTRO18)
 SPEDAT(SET,EXPORT,EPAT8,C,ETRO18)
 SPEDAT(SET,EXPORT,HPAT9,C,HTRO19)
 SPEDAT(SET,EXPORT,EPAT9,C,ETRO19)
 SPEDAT(SET,EXPORT,HPAT10,C,HTRO110)
 SPEDAT(SET,EXPORT,EPAT10,C,ETRO110)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.6653 ;EX(U1)=1.872
 EX(V1)=0.02486 ;EX(W1)=1.874
 EX(KE)=0.06406 ;EX(EP)=2.039E-04
 EX(RHO1)=1.12 ;EX(MARK)=0.6183
 EX(HIG)=7.655 ;EX(EL1)=125.
 EX(ENUT)=16.139999 ;EX(EPKE)=1.0E-10
 EX(CONC)=0.0108
 ************************************************************
  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(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 OUTPUT(MARK,Y,N,Y,N,N,N)
 OUTPUT(HIG,Y,N,Y,N,N,N)
 OUTPUT(EL1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(CONC,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 10 ;IZMON = 18
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 1 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 1 ;IDISPB = 1 ;IDISPC = 10
 CSG1    ='1'
STOP