```
photon use
p

up y
msg the grid. press return for temperature contours
gr z 1;pause; gr off;red;gr ou z 1
msg temperature contours. press return for velocity vectors
con tem1 z 1 fi;0.0002;pause; con off;red;se re ve 2;vec z 1 sh
msg velocity vectors. press return for reduced-pressure contours
pause;vec off;red;con p1 z 1 fi;0.0002
enduse

************************************************************
Group 1. Run Title
TEXT(NON-ISOTHERMAL 2D LAMINAR MIXING LAYER
************************************************************
Group 2. Transience
DISPLAY
This case is a illustratiom of the Colocated-Covariant
Method (CCM). The geometry is with two inlets at differents
temperatures. Both sides walls are symetry planes. The
velocities difference and the temperatures entail a
recirculation zone behind the solid edge.

The solid conductivity is temperature dependent:
lenda=14.17+6.25E-3 T
(T in K)
The kinematic viscosity of air follows the Sutherland law.

ENDDIS
L(PAUSE
DISPLAY
Geometry:                 outlet
--------
P1=P2=1.013E5 Pa

|               |               ^ H2=0.1 m
|               |               |
|               |               |
|               |               |
|               |               |
|               |               |
|    |\         |               - H1=0.1 m
|    |//\       |               |
Symetry    |    |///|      |    Symetry    |
plane    |    |///|      |    plane      |
|^^^^|///|      |               |
| V1 |///|^^^^^^|               |
|    |///|  V2  |               |

air         air
T1=600K     T2=300K
V1=2m/s     V2=0.5m/s

ENDDIS
L(PAUSE
**********************|**************************************
Groups 3, 4, 5  Grid Information
INTEGER(NX1,NX2,NX3,NXTOT)
NX1=20;NX2=10;NX3=25;NXTOT=NX1+NX2+NX3
INTEGER(NY1,NY2,NY3,NYTOT)
NY1=25;NY2=20;NY3=30;NYTOT=NY1+NY2+NY3
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,NXTOT,NYTOT,1)
* Set overall domain extent:
*        xulast  yvlast  zwlast    name
XSI= 8.000E-03;YSI= 2.000E-01;ZSI= 1.000E-03;RSET(D,CHAM    )
************************************************************
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= 0.0000E+00;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P2  )
XPO= 0.0000E+00;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P3  )
XPO= 0.0000E+00;YPO= 2.0000E-01;ZPO= 0.0000E+00;GSET(P,P4  )
XPO= 1.5000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P5  )
XPO= 1.5000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P6  )
XPO= 1.5000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P7  )
XPO= 3.0000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P9  )
XPO= 3.0000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P10 )
XPO= 2.0000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P11 )
XPO= 8.0000E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P13 )
XPO= 8.0000E-03;YPO= 9.2000E-02;ZPO= 0.0000E+00;GSET(P,P14 )
XPO= 8.0000E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P15 )
XPO= 8.0000E-03;YPO= 2.0000E-01;ZPO= 0.0000E+00;GSET(P,P16 )
* Set lines/arcs
GSET(L,L1,P6,P7,NY2,1.0)
GSET(L,L2,P10,P11,NY2,1.0)
GSET(L,L3,P14,P15,NY2,1.0)
GSET(L,L4,P2,P3,NY2,1.0)
GSET(L,L5,P7,P11,NX2,1.0)
GSET(L,L6,P1,P5,NX1,-1.3)
GSET(L,L7,P5,P9,NX2,1.0)
GSET(L,L8,P9,P13,NX3,1.3)
GSET(L,L9,P13,P14,NY1,S1.5)
GSET(L,L10,P9,P10,NY1,S1.5)
GSET(L,L11,P5,P6,NY1,S1.5)
GSET(L,L12,P1,P2,NY1,S1.5)
GSET(L,L13,P15,P16,NY3,1.8)
GSET(L,L14,P3,P4,NY3,1.8)
GSET(L,L18,P3,P7,NX1,-1.3)
GSET(L,L19,P11,P15,NX3,1.3)
GSET(L,L20,P10,P14,NX3,1.3)
GSET(L,L21,P4,P16,NX1+NX2+NX3,1.0)
* Set frames
GSET(F,F1,P1,-,P5,P6,P7,-,P3,P2)
GSET(F,F2,P5,-,P9,P10,P11,-,P7,P6)
GSET(F,F3,P9,-,P13,-,P14,-,P10,-)
GSET(F,F4,P10,-,P14,-,P15,-,P11,-)
GSET(F,F5,P3,P7.P11,P15,-,P16,-,P4,-)
* Match a grid mesh
GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(M,F2,+I+J,NX1+1,1,1,TRANS)
GSET(M,F3,+I+J,NX1+NX2+1,1,1,TRANS)
GSET(M,F4,+I+J,NX1+NX2+1,NY1+1,1,TRANS)
GSET(M,F5,+I+J,1,NY1+NY2+1,1,TRANS)
* Copy/Transfer/Block grid planes
GSET(C,K2,F,K1,1,NXTOT,1,NYTOT,+,0,0,1.0000E-03,INC,1)
**********
NONORT  =    T
* X-cyclic boundaries switched
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME( 42) =ENUL ; NAME( 43) =KOND
NAME( 44) =BLOK ; NAME( 45) =TEM1
NAME( 46) =WCRT ; NAME( 47) =VCRT
NAME( 48) =UCRT ; NAME( 49) =DEN1
NAME( 50) =PRPS
* Solved variables list
SOLVE(P1  ,U1  ,V1  ,TEM1)
* Stored variables list
STORE(PRPS,DEN1,UCRT,VCRT,WCRT,BLOK,KOND,ENUL)
SOLUTN(P1  ,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)

CCM=T
IVARBK  =      -1 ;ISOLBK =       1
************************************************************
Group 8. Terms & Devices
NEWRH1  =    T
NEWENL  =    T
ISOLX   =       0 ;ISOLY  =       0 ;ISOLZ  =       0
************************************************************
Group 9. Properties
PRESS0  = 1.013E+05
ENUT   = 0.000E+00
*  List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=2
*  Name
*Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
*
152 8000 0 500 GRND4 0 0
*       constants for GRND option no 1
1.4170E+01 6.2500E-03
*
31 GRND5 GRND6 1004.0 0.0263 3.33E-3 GRND5
*       constants for GRND option no 1 (Ideal gas law)
0 3.4843E-03
*       constants for GRND option no 2 (Sutherland law)
1.4600E-06 110
*       constants for GRND option no 3
0 0 7.1429E-01
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(BLOK) =  1.000E+00 ;FIINIT(TEM1) =  4.500E+02
FIINIT(PRPS) =  3.100E+01

CONPOR(STEEL   , -1.00,CELL  ,-#2,-#2,-#1,-#2,-#1,-#1)
INIT(STEEL   ,BLOK, 0.000E+00, 2.000E+00)
INIT(STEEL   ,PRPS, 0.000E+00, 1.520E+02)

************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources

INLET (BFCIN1  ,SOUTH ,#1,#1,#1,#1,#1,#1,#1,#1)
VALUE (BFCIN1  ,P1  , GRND1    )
VALUE (BFCIN1  ,U1  , GRND1    )
VALUE (BFCIN1  ,V1  , GRND1    )
VALUE (BFCIN1  ,TEM1, 6.000E+02)
VALUE (BFCIN1  ,VCRT, 2.000E+00)
VALUE (BFCIN1  ,DEN1, 5.883E-01)

INLET (BFCIN2  ,SOUTH ,#3,#3,#1,#1,#1,#1,#1,#1)
VALUE (BFCIN2  ,P1  , GRND1    )
VALUE (BFCIN2  ,U1  , GRND1    )
VALUE (BFCIN2  ,V1  , GRND1    )
VALUE (BFCIN2  ,TEM1, 3.000E+02)
VALUE (BFCIN2  ,VCRT, 5.000E-01)
VALUE (BFCIN2  ,DEN1, 1.176E+00)

OUTLET(OUT     ,NORTH ,#1,#3,#3,#3,#1,#1,#1,#1)
VALUE (OUT     ,P1  , 0.000E+00)
VALUE (OUT     ,U1  , 0.000E+00)
VALUE (OUT     ,V1  , 0.000E+00)

BFCA    = 1.189E+00
EGWF    =    T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =     600
SELREF  =    T
RESFAC  = 1.000E-02
************************************************************
Group 16. Terminate Iterations
LITER (U1  ) =    2 ;LITER (V1  ) =    2
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 9.000E-01)
RELAX(U1  ,FALSDT, 5.062E-02)
RELAX(V1  ,FALSDT, 5.062E-02)
RELAX(TEM1,FALSDT, 1.062E+03)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
CSG10   =Q1
************************************************************
Group 22. Monitor Print-Out
IXMON   = NX1+NX2-1 ;IYMON  = NY1+NY2+2 ;IZMON  =       1
TSTSWP  =      -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN  =       1; NYPRIN  =       1
IPROF   =       3; ITABL   =       3
```