PHOTON USE
   p ; ; ; ; ;
 
 
 
 
 
   msg Computational Grid
   mgr 1 z 1 x 2 6 y 2 6 col 4
   mgr 2 z 1 x 1 5 y 2 6 col 6
   mgr 3 z 1 x 2 6 y 2 6 col 9
   mgr 4 z 1 x 2 6 y 2 6 col 12
   mgr 5 z 1 x 2 6 y 1 5 col 14
   mgr 1 ou z m x 2 6 y 2 6
   mgr 2 ou z m x 1 5 y 2 6
   mgr 3 ou z m x 2 6 y 2 6
   mgr 4 ou z m x 2 6 y 2 6
   mgr 5 ou z m x 2 6 y 1 5
   mgr 2 ou x 1 y 2 6
   mgr 3 ou x m y 2 6
   mgr 4 ou y m x 2 6
   mgr 5 ou y 1 x 2 6
   msg Press return for Velocity Vectors
   pause
   mvec 1 x 4 sh
   mvec 4 x 4 sh
   mvec 5 x 4 sh
   red
   msg Vectors on Centre Plane
   msg Press return for Vectors at entry plane
   pause
   vec cl;vi z
   mvec 1 z 1 sh
   mvec 2 z 1 sh
   mvec 3 z 1 sh
   mvec 4 z 1 sh
   mvec 5 z 1 sh
   red
   msg Vectors at Entry Plane
   msg Press return for Vectors near exit plane
   pause
   vec cl;vi y;up -z
   mvec 1 z 13 sh
   mvec 2 z 13 sh
   mvec 3 z 13 sh
   mvec 4 z 13 sh
   mvec 5 z 13 sh
   red
   msg Vectors near Exit Plane
   msg Press return for Pressure Contours
   pause
   vec cl;vi 1 1 1;up y
   mcon 1 p1 x 4 fi
   .001
   mcon 4 p1 x 4 fi
   .001
   mcon 5 p1 x 4 fi
   .001
   msg Pressure contours
   msg Press return for LTLS contours near exit plane
   pause
   vec cl;con cl; vi y;up -z
   mcon 1 ltls z 13 fi
    .001
   mcon 2 ltls z 13 fi
    .001
   mcon 3 ltls z 13 fi
    .001
   mcon 4 ltls z 13 fi
    .001
   mcon 5 ltls z 13 fi
    .001
   red
   msg Press E then return to end
  ENDUSE
 
 ************************************************************
   Q1 created by GENERAL menu, Version 2.2, Date 17/09/96
 CPVNAM=GENERAL
 ************************************************************
  Group 1. Run Title
 TEXT(Pipe Flow with Multi-Block Grid and LVEL:585
 TITLE
 ************************************************************
  Group 6. Body-Fitted coordinates
 BFC=T
    * Set points
 XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P1  )
 XPO= 5.0000E-01;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P2  )
 XPO= 5.0000E-01;YPO= 5.0000E-01;ZPO= 0.0000E+00;GSET(P,P3  )
 XPO= 0.0000E+00;YPO= 5.0000E-01;ZPO= 0.0000E+00;GSET(P,P4  )
 XPO= 9.5710E-01;YPO= 9.5710E-01;ZPO= 0.0000E+00;GSET(P,P5  )
 XPO= 9.5710E-01;YPO=-4.5710E-01;ZPO= 0.0000E+00;GSET(P,P6  )
 XPO=-4.5710E-01;YPO=-4.5710E-01;ZPO= 0.0000E+00;GSET(P,P7  )
 XPO=-4.5710E-01;YPO= 9.5710E-01;ZPO= 0.0000E+00;GSET(P,P8  )
 XPO= 2.5000E-01;YPO= 1.2500E+00;ZPO= 0.0000E+00;GSET(P,P9  )
 XPO= 1.2500E+00;YPO= 2.5000E-01;ZPO= 0.0000E+00;GSET(P,P10 )
 XPO= 2.5000E-01;YPO=-7.5000E-01;ZPO= 0.0000E+00;GSET(P,P11 )
 XPO=-7.5000E-01;YPO= 2.5000E-01;ZPO= 0.0000E+00;GSET(P,P12 )
    * Set lines/arcs
 GSET(L,L1,P1,P2,5,1.0)
 GSET(L,L2,P2,P3,5,1.0)
 GSET(L,L3,P3,P4,5,1.0)
 GSET(L,L4,P4,P1,5,1.0)
 GSET(L,L5,P4,P8,5,1.0)
 GSET(L,L6,P3,P5,5,1.0)
 GSET(L,L7,P2,P6,5,1.0)
 GSET(L,L8,P1,P7,5,1.0)
 GSET(L,A9,P8,P5,5,1.0,ARC,P9)
 GSET(L,A10,P5,P6,5,1.0,ARC,P10)
 GSET(L,A11,P6,P7,5,1.0,ARC,P11)
 GSET(L,A12,P7,P8,5,1.0,ARC,P12)
    * Set frames
 GSET(F,F1,P1,-,P2,-,P3,-,P4,-)
 GSET(F,F2,P7,-,P1,-,P4,-,P8,-)
 GSET(F,F3,P2,-,P6,-,P5,-,P3,-)
 GSET(F,F4,P4,-,P3,-,P5,-,P8,-)
 GSET(F,F5,P7,-,P6,-,P2,-,P1,-)
    * Match a grid mesh
 GSET(D,7,7,15)
 GSET(M,F1,+I+J,2,2,1,TRANS)
 GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
 GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
 GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
 GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
 GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
 GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
 GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
 DUMPC(GRID1)
 
 GSET(D,6,7,15)
 GSET(M,F2,+I+J,1,2,1,TRANS)
 GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
 GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
 GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
 GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
 GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
 GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
 DUMPC(GRID2)
 
 GSET(D,6,7,15)
 GSET(M,F3,+I+J,2,2,1,TRANS)
 GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
 GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
 GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
 GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
 GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
 GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
 DUMPC(GRID3)
 
 GSET(D,7,6,15)
 GSET(M,F4,+I+J,2,2,1,TRANS)
 GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
 GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
 GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
 GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
 GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
 GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
 DUMPC(GRID4)
 
 GSET(D,7,6,15)
 GSET(M,F5,+I+J,2,1,1,TRANS)
 GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
 GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
 GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
 GSET(C,I1,F,I2,1,NY-1,1,NZ,+,0,0,0)
 GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
 GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
 DUMPC(GRID5)
    **********
 
 NUMBLK=5
   ### Make default grid BIG, so no cells are left-handed
 XSI=5; YSI=5; ZSI=5; RSET(D,CHAM)
 READCO(GRID+)
 MESG(Show the grid? (N/y)
 READVDU(ANS,CHAR,N)
 IF(:ANS:.EQ.Y) THEN
+  READCO(VIEWONLY)
+  VIEW
 ENDIF
 
 NONORT  =    T
    * X-cyclic boundaries switched
 
 MPATCH(1,MBL1-2,WEST ,   2,   2,2,NY-1,1,NZ,1,1)
 MPATCH(2,MBL2-1,EAST ,NX-1,NX-1,2,NY-1,1,NZ,1,1)
 
 MPATCH(1,MBL1-3,EAST ,NX-1,NX-1,2,NY-1,1,NZ,1,1)
 MPATCH(3,MBL3-1,WEST ,   2,   2,2,NY-1,1,NZ,1,1)
 
 MPATCH(1,MBL1-4,NORTH,2,NX-1,NY-1,NY-1,1,NZ,1,1)
 MPATCH(4,MBL4-1,SOUTH,2,NX-1,   2,   2,1,NZ,1,1)
 
 MPATCH(1,MBL1-5,SOUTH,2,NX-1,   2,   2,1,NZ,1,1)
 MPATCH(5,MBL5-1,NORTH,2,NX-1,NY-1,NY-1,1,NZ,1,1)
  
   
 MPATCH(2,MBL2-4,NORTH,1,NX-1,NY-1,NY-1,1,NZ,1,1)
 MPATCH(4,MBL4-2,WEST ,2,   2,   2,NY  ,1,NZ,1,1)
 SPEDAT(SET,GCV,MBL4-2,C,WNL)
 
 MPATCH(2,MBL2-5,SOUTH,1,NX-1,2,   2,1,NZ,1,1)
 MPATCH(5,MBL5-2,WEST ,2,   2,1,NY-1,1,NZ,1,1)
 SPEDAT(SET,GCV,MBL5-2,C,ESL)
 
 MPATCH(3,MBL3-4,NORTH,   2,NX  ,NY-1,NY-1,1,NZ,1,1)
 MPATCH(4,MBL4-3,EAST ,NX-1,NX-1,   2,NY  ,1,NZ,1,1)
 SPEDAT(SET,GCV,MBL4-3,C,ESL)
 
 MPATCH(3,MBL3-5,SOUTH,   2,  NX,2,   2,1,NZ,1,1)
 MPATCH(5,MBL5-3,EAST ,NX-1,NX-1,1,NY-1,1,NZ,1,1)
 SPEDAT(SET,GCV,MBL5-3,C,WNL)
 
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  )
    * Stored variables list
 STORE(PRPS,DEN1,UCRT,VCRT,WCRT,WDIS,WGAP,ENUT,EL1)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 TURMOD(LVEL)
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T ; NEWENL  =    T
 ************************************************************
  Group 9. Properties
 RHO1=1.189
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS) =  0.000E+00
 INIADD  =    T
 ************************************************************
  Group 13. Boundary & Special Sources
 REAL(WIN,FAC)
 WIN=5; FAC=1
 MPATCH(1,IN1,LOW,2,NX-1,2,NY-1,1,1,1,1)
 COVAL(IN1,P1,FIXFLU,RHO1*WIN)
 COVAL(IN1,W1,ONLYMS,WIN)
 
 MPATCH(2,IN2,LOW,1,NX-1,2,NY-1,1,1,1,1)
 COVAL(IN2,P1,FIXFLU,RHO1*WIN)
 COVAL(IN2,W1,ONLYMS,WIN)
 COVAL(IN2,V1,ONLYMS,WIN/FAC)
 
 MPATCH(3,IN3,LOW,2,NX,2,NY-1,1,1,1,1)
 COVAL(IN3,P1,FIXFLU,RHO1*WIN)
 COVAL(IN3,W1,ONLYMS,WIN)
 COVAL(IN3,V1,ONLYMS,-WIN/FAC)
 
 MPATCH(4,IN4,LOW,2,NX-1,2,NY,1,1,1,1)
 COVAL(IN4,P1,FIXFLU,RHO1*WIN)
 COVAL(IN4,W1,ONLYMS,WIN)
 COVAL(IN4,U1,ONLYMS,WIN/FAC)
 
 MPATCH(5,IN5,LOW,2,NX-1,1,NY-1,1,1,1,1)
 COVAL(IN5,P1,FIXFLU,RHO1*WIN)
 COVAL(IN5,W1,ONLYMS,WIN)
 COVAL(IN5,U1,ONLYMS,-WIN/FAC)
 
 MPATCH(1,OUT1,HIGH,2,NX-1,2,NY-1,NZ,NZ,1,1)
 COVAL(OUT1,P1,1E3,0)
 
 MPATCH(2,OUT2,HIGH,1,NX-1,2,NY-1,NZ,NZ,1,1)
 COVAL(OUT2,P1,1E3,0)
 
 MPATCH(3,OUT3,HIGH,2,NX  ,2,NY-1,NZ,NZ,1,1)
 COVAL(OUT3,P1,1E3,0)
 
 MPATCH(4,OUT4,HIGH,2,NX-1,2,NY  ,NZ,NZ,1,1)
 COVAL(OUT4,P1,1E3,0)
 
 MPATCH(5,OUT5,HIGH,2,NX-1,1,NY-1,NZ,NZ,1,1)
 COVAL(OUT5,P1,1E3,0)
 
 MPATCH(2,WALL1,WWALL,1,1,2,NY-1,1,NZ,1,1)
 COVAL(WALL1,W1,GRND2,0)
 
 MPATCH(3,WALL2,EWALL,NX,NX,2,NY-1,1,NZ,1,1)
 COVAL(WALL2,W1,GRND2,0)
 
 MPATCH(4,WALL3,NWALL,2,NX-1,NY,NY,1,NZ,1,1)
 COVAL(WALL3,W1,GRND2,0)
 
 MPATCH(5,WALL4,SWALL,2,NX-1,1,1,1,NZ,1,1)
 COVAL(WALL4,W1,GRND2,0)
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       100; SELREF  =    T; RESFAC  = 1.000E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(U1  ,FALSDT, 3/WIN)
 RELAX(V1  ,FALSDT, 3/WIN)
 RELAX(W1  ,FALSDT, 3/WIN)
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
GCV=T
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      18 ;IYMON  =       3 ;IZMON  =       7
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 24. Dumps For Restarts
 ************************************************************