TALK=T;RUN( 1, 1)
  Q1 created by TACT menu, Version 3.6, Date 30/06/05
CPVNAM=TACT
   ***********************************************************
   *
   *  GROUP 1.  Run identifiers and other preliminaries.
   *
TEXT(COOLING_TOWER_SIMULATION 3D
   *
   *  List of input geometrical data
REAL(RAINH);  RAINH=5.8350E+00
REAL(INLETH); INLETH=6.9325E+00
REAL(FILLH);  FILLH=2.1950E+00
REAL(ELIMH);  ELIMH=9.7500E+00
REAL(TOWERH); TOWERH=1.1540E+02
REAL(BTMFR);  BTMFR=3.8800E+01
REAL(TOPFR);  TOPFR=4.0280E+01
REAL(ELIMR);  ELIMR=3.9750E+01
REAL(EXTR);   EXTR=24
   *  List of initial data
   *  Ambient conditions and inlet water
REAL(PAMB);  PAMB=1.0050E+05
REAL(TAMB);  TAMB=24
REAL(TWAMB); TWAMB=20
REAL(TINI);  TINI=4.0700E+01
REAL(WFLR);  WFLR=1.6890E+00
   *  Heat-transfer characteristics
REAL(HSLDA); HSLDA=2.5000E-02
REAL(HSN);   HSN=-5.0000E-01
REAL(HPLDA); HPLDA=8.6300E-01
REAL(HPN);   HPN=-8.9300E-01
REAL(HRLDA); HRLDA=6.4600E-02
REAL(HRN);   HRN=-5.0000E-01
   *  Resistance characteristics
REAL(VSLDA); VSLDA=5.3000E-01
REAL(VSN);   VSN=1.3200E+00
REAL(VPLDA); VPLDA=4
REAL(VPN);   VPN=30
REAL(VRLDA); VRLDA=3.0000E-02
REAL(VRN);   VRN=1.0000E-01
REAL(VELIM); VELIM=3.4500E+00
REAL(VSTRT); VSTRT=1.7200E-01
   *  Wind velocity
REAL(WAMB); WAMB=2.5000E+00
   *  List of auxiliary integer and real variables
INTEGER(NYM,NZM,IHZ1,IHZ2,ILY1)
NYM=19; NZM=21
IHZ1=10; IHZ2=15; ILY1=1
REAL(VLM);   VLM=1.9847E+07
REAL(LEWNO); LEWNO=0.9074
   *
   ***********************************************************
   *
   *  GROUP 2.  Time-dependence and related parameters.
   *
   *  Steady-state simulation
STEADY=T
   *
   ***********************************************************
   *
   *  GROUP 3.  x-directon grid specification.
   *
   *  Grid is Body-fitted, X is angular direction
NX=12
XULAST=3.1400E+00
   ***********************************************************
   *
   *  GROUP 4.  y-directon grid specification.
   *
   *  Grid is Body-fitted, Y is radial direction
NY=26
YVLAST=2.0028E+02
   ***********************************************************
   *
   *  GROUP 5.  z-directon grid specification.
   *
   *  Grid is Body-fitted, Z is axial direction
NZ=42
ZWLAST=315
   ***********************************************************
   *
   *  GROUP 6.  Body-fitting and other grid distortions.
   *
   *  Non-orthogonal BFC option selected
BFC=T;NONORT=T
READCO(*../d_earth/d_spe/d_tact/inplib/tac3)
STORE(UCRT,VCRT,WCRT)
   *
   ***********************************************************
   *
   *  GROUP 7.  Variables (including porosities) named,
   *            stored & solved.
   *
NPOR=21;NAME(NPOR)=NPOR
SOLVE(V1,W1,H1,C1);SOLUTN(P1,Y,Y,Y,N,N,N)
STORE(DEN1,NPOR,VISL,TWTR,TAIR)
STORE(WFLR,GHST,GCST,BETA,AIRF,GQ)
STORE(ACST,WEVP,ALFA,GHFS,HTCF,MTCF)
   *  The simulation is 3D so solve for U1 as well
SOLVE(U1)
   *
   ***********************************************************
   *
   *  GROUP 8.  Terms (in differential equations) and devices.
   *
TERMS(H1,Y,Y,Y,P,P,P)
   *
   ***********************************************************
   *
   *  GROUP 9.  Properties of the medium (or media).
   *
RHO1=GRND;PRNDTL(H1)=0.72;ENUL=GRND
   *
   ***********************************************************
   *
   *  GROUP 10. Interphase-transfer processes and proerties.
   *
   *
   ***********************************************************
   *
   *  GROUP 11. Initialization of fields of variables,
   *            porosities, etc.
   *
INIADD=F;FIINIT(P1)=0.0;FIINIT(H1)=0.0
FIINIT(V1)=0.0;FIINIT(W1)=0.0
FIINIT(VISL)=1.E-5;FIINIT(NPOR)=1.
FIINIT(U1)=0.0
PATCH(W,INIVAL,1,NX,1,NYM,1,34+1,1,1)
COVAL(W,W1,0.0,0.1)
PATCH(V,INIVAL,1,NX/2,1,NYM,1,10,1,1)
COVAL(V,V1,0.0,-0.1)
  ** Tower wall
PATCH(TW,INIVAL,1,NX,NYM,NYM,15+1,34,1,1)
COVAL(TW,NPOR,0.0,0.0)
   *
   ***********************************************************
   *
   *  GROUP 12. Convection and diffusion adjustments
   *
   *
   ***********************************************************
   *
   *  GROUP 13. Boundary conditions and special sources
   *
  ** Inlet struts **
PATCH(INLET,NORTH,1,NX,NYM,NYM,1,15,1,1)
COVAL(INLET,V1,GRND6,0.0)
  ** North boundary - Exit**
PATCH(NORBOU,NORTH,NX/2+1,NX,NY,NY,1,NZ,1,1)
COVAL(NORBOU,P1,FIXP,0.0)
COVAL(NORBOU,W1,ONLYMS,SAME);COVAL(NORBOU,V1,ONLYMS,SAME)
COVAL(NORBOU,C1,ONLYMS,SAME);COVAL(NORBOU,H1,ONLYMS,SAME)
COVAL(NORBOU,U1,ONLYMS,SAME)
  ** Top boundary **
PATCH(TOP,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(TOP,P1,FIXP,0.0)
COVAL(TOP,H1,ONLYMS,SAME);COVAL(TOP,V1,ONLYMS,SAME)
COVAL(TOP,C1,ONLYMS,SAME);COVAL(TOP,W1,ONLYMS,SAME)
COVAL(TOP,U1,ONLYMS,SAME)
  ** Inside wall **
PATCH(WALL,NWALL,1,NX,NYM,NYM,15+1,34,1,1)
COVAL(WALL,W1,1.0,0.0)
COVAL(WALL,U1,1.0,0.0)
  ** Outside wall **
PATCH(OUTWL,SWALL,1,NX,NYM+1,NYM+1,15+1,34,1,1)
COVAL(OUTWL,W1,1.0,0.0)
COVAL(OUTWL,U1,1.0,0.0)
  ** Ground no-slip **
PATCH(GRDNSL,LWALL,1,NX,1,NY,1,1,1,1)
COVAL(GRDNSL,V1,1.0,0.0)
COVAL(GRDNSL,U1,1.0,0.0)
  ** Buoyancy **
PATCH(BUOY,PHASEM,1,NX,1,NY,1,NZ,1,1)
COVAL(BUOY,W1,FIXFLU,GRND2)
  ** Eliminator **
PATCH(ELIM,HIGH,1,NX,1,NYM,NZM,NZM,1,1)
COVAL(ELIM,W1,GRND7,0.0)
  ** Rain **
PATCH(RAIN,VOLUME,1,NX,1,NYM,1,10,1,1)
COVAL(RAIN,U1,GRND5,0.0)
COVAL(RAIN,V1,GRND5,0.0);COVAL(RAIN,W1,GRND5,0.0)
COVAL(RAIN,H1,GRND1,GRND6);COVAL(RAIN,C1,GRND2,GRND7)
   * Fill 1st row *
PATCH(FILL,VOLUME,1,NX,1,18,10+1,15,1,1)
COVAL(FILL,U1,FIXVAL,0.0)
COVAL(FILL,V1,FIXVAL,0.0);COVAL(FILL,W1,GRND3,0.0)
COVAL(FILL,H1,GRND4,GRND6);COVAL(FILL,C1,GRND8,GRND7)
   * Corner rain *
PATCH(RAINCORN,VOLUME,1,NX,19,NYM,10+1,15,1,1)
COVAL(RAINCORN,U1,GRND5,0.0)
COVAL(RAINCORN,V1,GRND5,0.0);COVAL(RAINCORN,W1,GRND5,0.0)
COVAL(RAINCORN,H1,GRND1,GRND6);COVAL(RAINCORN,C1,GRND2,GRND7)
   * Fill 2nd row *
PATCH(FILLSR,VOLUME,1,NX,1,NYM,15+1,NZM-1,1,1)
COVAL(FILLSR,U1,FIXVAL,0.0)
COVAL(FILLSR,V1,FIXVAL,0.0);COVAL(FILLSR,W1,GRND3,0.0)
COVAL(FILLSR,H1,GRND4,GRND6);COVAL(FILLSR,C1,GRND8,GRND7)
LG(4)=F
  ** Water spray **
PATCH(SPRAY,VOLUME,1,NX,1,NYM,NZM,NZM,1,1)
COVAL(SPRAY,W1,GRND,0.0)
COVAL(SPRAY,H1,GRND10,GRND6);COVAL(SPRAY,C1,GRND9,GRND7)
   * Mass source due to evaporation
PATCH(MASRC,VOLUME,1,NX,1,NYM,1,NZM-1,1,1)
COVAL(MASRC,P1,FIXFLU,GRND4)
COVAL(MASRC,U1,0.0,0.0)
COVAL(MASRC,V1,0.0,0.0);COVAL(MASRC,W1,0.0,0.0)
COVAL(MASRC,H1,0.0,0.0);COVAL(MASRC,C1,0.0,0.0)
   * North boundary-wind inlet *
PATCH(BFCWN,NORTH,1,NX/2,NY,NY,1,NZ,1,1)
COVAL(BFCWN,P1,FIXFLU,GRND1);COVAL(BFCWN,U1,ONLYMS,GRND1)
COVAL(BFCWN,W1,ONLYMS,0.0);COVAL(BFCWN,V1,ONLYMS,GRND1)
COVAL(BFCWN,C1,ONLYMS,GRND3);COVAL(BFCWN,H1,ONLYMS,0.0)
COVAL(BFCWN,UCRT,0.0,0.0);COVAL(BFCWN,VCRT,0.0,WAMB*-1.)
COVAL(BFCWN,WCRT,0.0,0.0)
   * Inlet struts U1 *
PATCH(INLU1,EAST,1,NX,NYM,NYM,1,15,1,1)
COVAL(INLU1,U1,GRND6,0.0)
   *
   ***********************************************************
   *
   *  GROUP 14. Downstream pressure (for free parabolic flow).
   *
   *
   ***********************************************************
   *
   *  GROUP 15. Termination criteria for sweeps and
   *            outer iterations.
   *
   *  Number of sweeps
LSWEEP=1000
   *
   ***********************************************************
   *
   *  GROUP 16. Termination criteria for inner iterations.
   *
   *  Estimate of the maximum velocity within domain
REAL(MAXV)
MAXV=2.5000E+00
RESREF(P1)=VLM*1.E-12
RESREF(U1)=VLM*1.15*1.E-12*MAXV
RESREF(V1)=VLM*1.15*1.E-12*MAXV
RESREF(W1)=VLM*1.15*1.E-12*MAXV
RESREF(H1)=VLM*1.15*1.E-12*5.E4
RESREF(C1)=VLM*1.15*1.E-12*0.02
LITER(P1)=50
   *
   ***********************************************************
   *
   *  GROUP 17. Under-relaxation and related devices.
   *
   *  Relaxation applied
RELAX(U1,FALSDT,1.E-1)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,0.1)
RELAX(H1,FALSDT,1.E3)
RELAX(C1,FALSDT,1.E2)
   *
   ***********************************************************
   *
   *  GROUP 18. Limits on variables values or increments
   *            to them.
   *
   *
   ***********************************************************
   *
   *  GROUP 19. Data communicated by SATELLITE to GROUND
   *
NAMGRD=TACT
IG(1)=NZM;IG(2)=NYM;IG(3)=IHZ1;IG(4)=IHZ2;IG(5)=ILY1
IG(6)=35; BUOYC=-9.81
RG(1)=WFLR;RG(2)=PAMB;RG(3)=TAMB;RG(4)=TWAMB;RG(5)=WAMB
RG(6)=HSLDA;RG(7)=HSN;RG(8)=VSN;RG(9)=VSLDA;RG(10)=HPLDA
RG(11)=HPN;RG(12)=VPLDA;RG(13)=VPN;RG(14)=HRLDA
RG(15)=HRN;RG(16)=VRLDA;RG(17)=VRN;RG(18)=TINI
RG(19)=VELIM;RG(20)=VSTRT
RG(21)=FILLH;RG(22)=RAINH;RG(23)=INLETH;RG(24)=ELIMH
RG(25)=TOWERH;RG(26)=BTMFR;RG(27)=TOPFR;RG(28)=ELIMR
RG(29)=EXTR;RG(30)=LEWNO;RG(40)=ZWLAST
   * Viscosity values
RG(35)=0.11;RG(36)=1.E-5
LG(3)=F
   *
   ***********************************************************
   *
   *  GROUP 20. Control of preliminary printout
   *
   *1-input-2-grid-input
LG(1)=F;LG(2)=F
ECHO=F
   *
   ***********************************************************
   *
   *  GROUP 21. Frequency and extent of field printout.
   *
LG(2)=F
OUTPUT(P1,Y,P,P,Y,Y,Y);OUTPUT(U1,Y,P,P,Y,Y,Y)
OUTPUT(V1,Y,P,P,Y,Y,Y);OUTPUT(W1,Y,P,P,Y,Y,Y)
OUTPUT(H1,Y,P,P,Y,Y,Y);OUTPUT(C1,Y,P,P,Y,Y,Y)
OUTPUT(TWTR,Y,P,P,Y,Y,Y);OUTPUT(TAIR,Y,P,P,Y,Y,Y)
OUTPUT(DEN1,Y,P,P,Y,Y,Y)
OUTPUT(VISL,N,N,N,N,N,N);OUTPUT(NPOR,N,N,N,N,N,N)
OUTPUT(WFLR,N,N,N,N,N,N);OUTPUT(BETA,N,N,N,N,N,N)
OUTPUT(GHST,N,N,N,N,N,N);OUTPUT(GCST,N,N,N,N,N,N)
OUTPUT(MTCF,N,N,N,N,N,N);OUTPUT(HTCF,N,N,N,N,N,N)
OUTPUT(AIRF,N,N,N,N,N,N);OUTPUT(GQ,N,N,N,N,N,N)
OUTPUT(WEVP,N,N,N,N,N,N);OUTPUT(ALFA,N,N,N,N,N,N)
OUTPUT(GHFS,N,N,N,N,N,N)
OUTPUT(UCRT,N,N,N,N,N,N);OUTPUT(VCRT,N,N,N,N,N,N)
OUTPUT(WCRT,N,N,N,N,N,N)
   *
   ***********************************************************
   *
   *  GROUP 22. Location of spot-value & frequency of
   *            residual printout.
   *
   *  Assign cell-indices of spot-point monitoring location
IXMON=3
IYMON=19
IZMON=25
   *  Set frequency of spot-values/residuals print-outs
NPRMNT=1
TSTSWP=-1
   *
   ***********************************************************
   *
   *  GROUP 23. Variable-by-variable field printout and plot
   *            and/or tabulation of spot-values and residuals.
   *
ITABL=1;NPLT=1;IPLTL=1000
NPRINT=LSWEEP
   *
   ***********************************************************
   *
   *  GROUP 24. Preperation for continuation runs.
   *
   *
   ***********************************************************
   *
   *  MENU ITEMS ( DO NOT ALTER SETTINGS )
   *
INTEGER(MNZ,YNTCN); MNZ=200
INTEGER(RAINCN,FILCN1,FILCN2,ZTWRCN,ZNTCN,YTWRCN,CRAINCN)
RAINCN=10;FILCN1=5;FILCN2=5;ZTWRCN=13
ZNTCN=4;YTWRCN=18;CRAINCN=1;YNTCN=4
ARRAY(SHELLR,REAL,MNZ)
SHELLR(1)=3.9750E+01;SHELLR(2)=3.7750E+01;SHELLR(3)=3.5450E+01
SHELLR(4)=3.3250E+01;SHELLR(5)=3.0750E+01;SHELLR(6)=2.8725E+01
SHELLR(7)=2.6500E+01;SHELLR(8)=2.4850E+01;SHELLR(9)=2.3400E+01
SHELLR(10)=2.2375E+01;SHELLR(11)=2.2000E+01;SHELLR(12)=2.2550E+01
SHELLR(13)=2.3225E+01;SHELLR(14)=2.4000E+01;
   *
   ***********************************************************
STOP