MB-FGE Test: Heat conduction in a 5-pointed star.
  **************************************************************
  DISPLAY
   Heat conduction in a 5-pointed star (Chris J.Weire 16/3/94).
 
   This is a  simple example of  a heat conduction  problem,
   using 5  linked blocks.   The "star"  is heated  along  a
   vertical line at its base, and cooled along its edges.
 
   The links used  here are not  "natural", in the  sense of
   being NORTH-SOUTH,  EAST-WEST or  HIGH-LOW. This  type of
   link is currently only supported for scalar problems.
 
   Q1-file contains Photon Use information.
   ---------------------------------------------------------
  ENDDIS
L(PAUSE
  **************************************************************
  PHOTON USE
  p ; ; ; ; ;
 
  mgr 1 k 1 col 2
  mgr 2 k 1 col 3
  mgr 3 k 1 col 4
  mgr 4 k 1 col 5
  mgr 5 k 1 col 6
  msg This is the grid, showing 5 blocks.  Press Return
  pause
  cl
  mcon 1 temp k 1 fi
  0.001
  mcon 2 temp k 1 fi
  0.001
  mcon 3 temp k 1 fi
  0.001
  mcon 4 temp k 1 fi
  0.001
  mcon 5 temp k 1 fi
  0.001
  msg This is the Temperature field.
  ENDUSE
    GROUP 1. Run title and other preliminaries
TEXT(MB-FGE: Heat Conduction in a 5-pointed star
TITLE
INTEGER(NS);
REAL(PI,SPHT,COND,ANG,XP,YP)
PI= 3.14159
  ** NS is grid size
NS= 10
    GROUP 6. Body-fitted coordinates or grid distortion
BFC= T; GSET(D,NS,NS,1); GSET(P,PC,0.,0.,0.)
ANG=2.*PI/10; XP=SIN(ANG);   YP=COS(ANG);    GSET(P,P1,XP,YP,0.0)
ANG=0.*PI/10; XP=2*SIN(ANG); YP=2.*COS(ANG); GSET(P,P2,XP,YP,0.0)
ANG=-2*PI/10; XP=SIN(ANG);   YP=COS(ANG);    GSET(P,P3,XP,YP,0.0)
GSET(L,L1,PC,P1,NS); GSET(L,L2,P1,P2,NS)
GSET(L,L3,P2,P3,NS); GSET(L,L4,P3,PC,NS)
GSET(F,F1,PC,-,P1,-,P2,-,P3,-); GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NS,1,NS,+,0.0,0.0,0.1)
DUMPC(MBGR1)
GSET(D,NS,NS,1)
ANG= 2.*PI/10+2*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P1,XP,YP,0.0)
ANG= 0.*PI/10+2*PI/5.; XP= 2.*SIN(ANG); YP= 2.*COS(ANG)
GSET(P,P2,XP,YP,0.0)
ANG=-2.*PI/10+2*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P3,XP,YP,0.0)
GSET(L,L1,PC,P1,NS); GSET(L,L2,P1,P2,NS)
GSET(L,L3,P2,P3,NS); GSET(L,L4,P3,PC,NS)
GSET(F,F1,PC,-,P1,-,P2,-,P3,-); GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NS,1,NS,+,0.,0.,0.1)
DUMPC(MBGR2)
GSET(D,NS,NS,1)
ANG= 2.*PI/10+4*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P1,XP,YP,0.0)
ANG= 0.*PI/10+4*PI/5.; XP= 2.*SIN(ANG); YP= 2.*COS(ANG)
GSET(P,P2,XP,YP,0.0)
ANG=-2.*PI/10+4*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P3,XP,YP,0.0)
GSET(L,L1,PC,P1,NS); GSET(L,L2,P1,P2,NS)
GSET(L,L3,P2,P3,NS); GSET(L,L4,P3,PC,NS)
GSET(F,F1,PC,-,P1,-,P2,-,P3,-); GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NS,1,NS,+,0.,0.,0.1)
DUMPC(MBGR3)
GSET(D,NS,NS,1)
ANG= 2.*PI/10+6*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P1,XP,YP,0.0)
ANG= 0.*PI/10+6*PI/5.; XP= 2.*SIN(ANG); YP= 2.*COS(ANG)
GSET(P,P2,XP,YP,0.0)
ANG=-2.*PI/10+6*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P3,XP,YP,0.0)
GSET(L,L1,PC,P1,NS); GSET(L,L2,P1,P2,NS)
GSET(L,L3,P2,P3,NS); GSET(L,L4,P3,PC,NS)
GSET(F,F1,PC,-,P1,-,P2,-,P3,-); GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NS,1,NS,+,0.,0.,0.1)
DUMPC(MBGR4)
GSET(D,NS,NS,1)
ANG= 2.*PI/10+8*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P1,XP,YP,0.0)
ANG= 0.*PI/10+8*PI/5.; XP= 2.*SIN(ANG); YP= 2.*COS(ANG)
GSET(P,P2,XP,YP,0.0)
ANG=-2.*PI/10+8*PI/5.; XP= SIN(ANG);    YP= COS(ANG)
GSET(P,P3,XP,YP,0.0)
GSET(L,L1,PC,P1,NS); GSET(L,L2,P1,P2,NS)
GSET(L,L3,P2,P3,NS); GSET(L,L4,P3,PC,NS)
GSET(F,F1,PC,-,P1,-,P2,-,P3,-); GSET(M,F1,+I+J,1,1,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NS,1,NS,+,0.,0.,0.1)
DUMPC(MBGR5)
NUMBLK= 5;  READCO(MBGR+)
GVIEW(Z);  VIEW
  ** Links not natural, so generated manually
MPATCH(1,MBL1.2,SOUTH,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(2,MBL2.1,WEST, 1, 1,1,NS,1,NZ,1,LSTEP)
MPATCH(2,MBL2.3,SOUTH,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(3,MBL3.2,WEST, 1, 1,1,NS,1,NZ,1,LSTEP)
MPATCH(3,MBL3.4,SOUTH,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(4,MBL4.3,WEST, 1, 1,1,NS,1,NZ,1,LSTEP)
MPATCH(4,MBL4.5,SOUTH,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(5,MBL5.4,WEST, 1, 1,1,NS,1,NZ,1,LSTEP)
MPATCH(5,MBL5.1,SOUTH,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(1,MBL1.5,WEST, 1, 1,1,NS,1,NZ,1,LSTEP)
    GROUP 7. Variables stored, solved & named
SOLVE(C1); NAME(C1)= TEMP; STORE(VPOR)
    GROUP 8. Terms in differential equations & devices
TERMS(TEMP,N,N,Y,N,Y,Y)
    GROUP 9. Properties of the medium (or media)
SPHT= 500.0; COND= 400.0; ENUL= 1.0; RHO1= 1.0E4;
PRNDTL(TEMP)= SPHT*ENUL*RHO1/COND
    GROUP 11. Initialization of variable or porosity fields
FIINIT(TEMP)= 0.0
    GROUP 13. Boundary conditions and special sources
MPATCH(3,HEAT3,CELL,1,NS,1, 1,1,NZ,1,LSTEP)
MPATCH(4,HEAT4,CELL,1, 1,1,NS,1,NZ,1,LSTEP)
COVAL(HEAT3,TEMP,FIXVAL,100.); COVAL(HEAT4,TEMP,FIXVAL,100.)
MPATCH(1,SINK1,EWALL, NS,NS,  1,NS, 1,1, 1,LSTEP)
MPATCH(1,SINK2,NWALL,  1,NS, NS,NS, 1,1, 1,LSTEP)
MPATCH(2,SINK3,EWALL, NS,NS,  1,NS, 1,1, 1,LSTEP)
MPATCH(2,SINK4,NWALL,  1,NS, NS,NS, 1,1, 1,LSTEP)
MPATCH(3,SINK5,EWALL, NS,NS,  1,NS, 1,1, 1,LSTEP)
MPATCH(3,SINK6,NWALL,  1,NS, NS,NS, 1,1, 1,LSTEP)
MPATCH(4,SINK7,EWALL, NS,NS,  1,NS, 1,1, 1,LSTEP)
MPATCH(4,SINK8,NWALL,  1,NS, NS,NS, 1,1, 1,LSTEP)
MPATCH(5,SINK9,EWALL, NS,NS,  1,NS, 1,1, 1,LSTEP)
MPATCH(5,SINK10,NWALL,  1,NS, NS,NS, 1,1, 1,LSTEP)
DO II=1,10
+ COVAL(SINK:II:,TEMP,FIXVAL,0.)
ENDDO
    GROUP 15. Termination of sweeps
LSWEEP= 10; TSTSWP= -1;  LITER(TEMP)=-50
    GROUP 19. Data communicated by satellite to GROUND
    * LSG4 = T, activates nonorthogonality treatment.
CSG3= LCRU; LSG4= T
    GROUP 22. Spot-value print-out
IXMON = 3; IYMON = 3; IZMON = 1