dbs 03.06.01 The link featuren in this does not work. I am therefore scrapping the whole feature, because In-Form can do the job much better 871

 
  PHOTON USE
  p;;;
 
  msg Temperature contours
  con tem1 z m fi;0.001;  con tem1 y m fi;0.001
  con tem1 x m fi;0.001
  msg
  msg Press e to END
  ENDUSE
    GROUP 1. Run title and other preliminaries
TEXT(3D Steady Heat Conduction In A Cube
TITLE
  DISPLAY
  This modification of Library case 100 exemplifies the use of
  the Linear Temperature Link.
 
  Linear Temperature link allows 'pipe-like' heat transfer by
  providing sub-grid-scale links within the computational domain. If
  there are thermal leads, thermal bridges or thin metal plates
  embedded in a poorly-conducting medium that would prevent use of
  coarse grid, this feature can directly link originating and
  destination cells.
  enddis
 
REAL(XLENGTH,YLENGTH,ZLENGTH)
XLENGTH=1.0;YLENGTH=1.0;ZLENGTH=1.0
NX=5;NY=5;NZ=5
  domain size and grid
NX=10;NY=10;NZ=10;xulast=xlength;yvlast=ylength;zwlast=zlength
#unigrid
 
    GROUP 7. Variables stored, solved & named
solve(tem1);SOLUTN(TEM1,Y,Y,Y,N,N,Y)
store(prps)
    GROUP 8. Terms (in differential equations) & devices
   **For pure conduction, cut out built-in source and convection
     terms
TERMS(TEM1,N,N,Y,N,Y,Y)
    GROUP 9. Properties of the medium (or media)

    GROUP 12. Local adjustments of convection and diffusion
    GROUP 13. Boundary conditions and special sources
   **Corner at IX=IY=IZ=1
PATCH(COLD,CELL,1,1,1,1,1,1,1,1)
   **Fix temperature to zero
COVAL(COLD,TEM1,1.E2,0.0)
   **Corner at IX=NX, IY=NY, IZ=NZ
PATCH(HOT,CELL,NX,NX,NY,NY,NZ,NZ,1,1)
   **Fix temperature to 1.0
COVAL(HOT,TEM1,1.E2,1.0)
L($021)
mesg(insert linear-link patches? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y) then
    linear link patches
 patch(-a,west,1,1,2,2,2,4,1,1)
 coval(-a,tem1,1.0e2,0.0)
 
 patch(-b,east,nx,nx,2,2,3,3,1,1)
 coval(-b,tem1,1.0e2,0.0)
 
 patch(-c,south,1,2,ny,ny,1,2,1,1)
 coval(-c,tem1,1.0,0.0)
 
 patch(-d,north,4,5,ny,ny,4,5,1,1)
 coval(-d,tem1,1.0,0.0)
 
 patch(-e,cell,3,4,3,4,nz,nz,1,1)
 coval(-e,tem1,1.e5,0.0)
 
 patch(-f,cell,3,4,ny,ny,3,4,1,1)
 coval(-f,tem1,1.0e5,0.0)
 
 patch(-g,cell,1,1,ny,ny,nz,nz,1,1)
 coval(-g,tem1,1.0e2,0.0)
 
 patch(-h,cell,1,1,1,1,nz,nz,1,1)
 coval(-h,tem1,1.e2,0.0)
 
 patch(heater,north,3,4,ny,ny,3,4,1,1)
 coval(heater,tem1,fixval,10.0)
 
mesg(linear-link patches are active
endif
L($021)
    GROUP 15. Termination of sweeps
OUTPUT(PRPS,N,N,N,N,N,N);LSWEEP=100
    GROUP 16. Termination of iterations
LITER(TEM1)=100; OVRRLX= 1.7
ISOLX=1;ISOLY=1;ISOLZ=1
    GROUP 17. Under-relaxation devices
RELAX(TEM1,LINRLX,0.5)
    GROUP 21. Print-out of variables
OUTPUT(TEM1,Y,Y,Y,Y,Y,Y); OUTPUT(PRPS,N,N,N,N,N,N)
    GROUP 22. Spot-value print-out
IXMON=NX/2+1;IYMON=NY/2+1;IZMON=NZ/2+1;UWATCH=T
    GROUP 23. Field print-out and plot control
NXPRIN=NX/5;NYPRIN=NY/5;NZPRIN=NZ/5;nplt=1;uwatch=t
    GROUP 24. Dumps for restarts