TALK=T;RUN(1,1)
  PHOTON USE
  p;;;;;
  gr ou z 1 1
  gr ou y m z 1 1;gr ou z 2;gr ou x 1 z 1 1 ;gr ou x m z 1 1
  gr ou y m z 2 2;gr ou z 3;gr ou x 1 z 2 2 ;gr ou x m z 2 2
  gr ou y m z 3 3;gr ou z 4;gr ou x 1 z 3 3 ;gr ou x m z 3 3
  gr ou y m z 4 4;gr ou z 5;gr ou x 1 z 4 4 ;gr ou x m z 4 4
  vec z 1 sh
  msg
  msg        In slab of IZ=1 the flow is 45 degree
  msg        from right bottom corner.
  msg
  pause
  vec off;con off;red
  vec z 2 sh
  msg
  msg        In slab IZ=2 the flow is a solid body
  msg        counter clockwise rotation.
  msg
  pause
  vec off;con off;red
  con h1 z 3 fil;.001
  vec z 3 sh
  msg
  msg        In slab IZ=3 the flow is the superposition
  msg        of those in slabs 1 and 2.
  msg
  msg        The source dispersion plume is as shown.
  msg
  msg        Hit Enter to get the result of convection
  msg        fluxes transformation.
  pause
  vec off;con off;red
  con h1 z 4 fil;.001
  vec z 4 sh
  msg
  MSG
  msg        In slab IZ=4 the dispersion is exactly the same as
  msg        before in spite of background flow being in opposite
  msg        direction.
  ENDUSE
  DISPLAY
    The four problems are arranged in  four  slab  cartesian
    box. Each slab is devoted to the particular problem.

    The PLANTing  of  initial velocity fields,  manipulation
    with them and the alteration  of  convection  fluxes  by
    add-extra-velocity option are demonstrated.
  ENDDIS
    PLANTBEGIN
PATCH(INITZ2,INIVAL,1,NX,1,NY,2,2,1,1)
   VAL= -(YG2D-15.)
COVAL(INITZ2,U1,zero,GRND)
   VAL= XG2D-15.
COVAL(INITZ2,V1,zero,GRND)
     The velocity field in the second slab is initialised as
     solid body rotation.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PATCH(INITZ3,INIVAL,1,NX,1,NY,3,3,1,1)
   VAL= U1[,,1]+U1[,,2]+YG2D
COVAL(INITZ3,U1,zero,GRND)
   VAL= V1[,,1]+V1[,,2]-XG2D
COVAL(INITZ3,V1,zero,GRND)
     The velocity field for the third sub-domain,  IZ=3,  is
     the superposition of the  velocity  components  in  two
     previous  slabs  and  the  cell  node  coordinates.  It
     results in 45 degree flow of  10  m/s  from  north-west
     edge of the domain.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
NAMSAT=MOSG
U1AD=GRND
   VELAD=U1[,,2]+YG2D
  REGION(,NX-1) /IZ.EQ.4
V1AD=GRND
   VELAD=V1[,,2]-XG2D
  REGION(,,,NY-1)
  IF(IZ.EQ.4)
      Although, the velocity field at the fourth slab is the
      same  as for first one,  the add-extra-velocity option
      is  activated  as  above  pointers  tell.  The   extra
      velocities  added  to  the  main  velocity  components
      alters the convection fluxes to be equal  to  ones  at
      third  slab.  The  expected  distribution of convected
      property H1 should be identical to the one at IZ=3.

      Please note the differences in REGION  commands.  They
      are  attributed  to  the  staggered nature of velocity
      nodes and the usage of either switch or IF command  to
      limit the Z-direction extent of velocity alterations.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
      PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT( Alteration of convection fluxes        )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 621
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 20
 XULAST =20.
 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 =20.
 YFRAC(1)=0.05 ;YFRAC(2)=0.1
 YFRAC(3)=0.15 ;YFRAC(4)=0.2
 YFRAC(5)=0.25 ;YFRAC(6)=0.3
 YFRAC(7)=0.35 ;YFRAC(8)=0.4
 YFRAC(9)=0.45 ;YFRAC(10)=0.5
 YFRAC(11)=0.55 ;YFRAC(12)=0.6
 YFRAC(13)=0.65 ;YFRAC(14)=0.7
 YFRAC(15)=0.75 ;YFRAC(16)=0.8
 YFRAC(17)=0.85 ;YFRAC(18)=0.9
 YFRAC(19)=0.95 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 4
 ZWLAST =4.
 ZFRAC(1)=0.25 ;ZFRAC(2)=0.5
 ZFRAC(3)=0.75 ;ZFRAC(4)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(3)=U1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(14)=H1
 NAME(150)=HPOR
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(U1,Y,N,N,N,N,Y)
 SOLUTN(V1,Y,N,N,N,N,Y)
 SOLUTN(W1,Y,N,N,N,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(HPOR,Y,N,N,N,N,Y)
 EPOR = 0 ;HPOR = 150 ;NPOR = 0 ;VPOR = 0
 ************************************************************
  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(H1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 U1AD = GRND ;V1AD = GRND
 W1AD =0.
 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. ;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-04 ;ENUT =0.
 PRNDTL(H1)=1.
 PRT(H1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(U1)=-5. ;FIINIT(V1)=5.
 FIINIT(W1)=1.0E-10 ;FIINIT(H1)=1.0E-10
 FIINIT(HPOR)=0.
 
 PATCH(INITZ2 ,INIVAL, 1, 20, 1, 20, 2, 2, 1, 1)
 INIT(INITZ2 ,U1 ,0. , GRND )
 INIT(INITZ2 ,V1 ,0. , GRND )
 
 PATCH(INITZ3 ,INIVAL, 1, 20, 1, 20, 3, 3, 1, 1)
 INIT(INITZ3 ,U1 ,0. , GRND )
 INIT(INITZ3 ,V1 ,0. , GRND )
 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(FIXSOR ,CELL , 10, 10, 10, 10, 3, 4, 1, 1)
 COVAL(FIXSOR ,H1 , FIXVAL ,1. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 10 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(H1)=20
 ENDIT(H1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(U1,LINRLX,1.)
 RELAX(V1,LINRLX,1.)
 RELAX(W1,LINRLX,1.)
 RELAX(H1,FALSDT,1.0E+09)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 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(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(U1)=6.563 ;EX(V1)=6.563
 EX(W1)=1.0E-10 ;EX(H1)=0.02162
 EX(HPOR)=1.0E-11
 ************************************************************
  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(U1,Y,N,Y,N,N,N)
 OUTPUT(V1,Y,N,Y,N,N,N)
 OUTPUT(W1,Y,N,Y,N,N,N)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(HPOR,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 5 ;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
 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
STOP