+++ /dev/null
-
-c---------------------------------------------------------------------
-c---------------------------------------------------------------------
-c--- applu.incl
-c---------------------------------------------------------------------
-c---------------------------------------------------------------------
-
-c---------------------------------------------------------------------
-c npbparams.h defines parameters that depend on the class and
-c number of nodes
-c---------------------------------------------------------------------
-
- include 'npbparams.h'
-
-c---------------------------------------------------------------------
-c parameters which can be overridden in runtime config file
-c (in addition to size of problem - isiz01,02,03 give the maximum size)
-c ipr = 1 to print out verbose information
-c omega = 2.0 is correct for all classes
-c tolrsd is tolerance levels for steady state residuals
-c---------------------------------------------------------------------
- integer ipr_default
- parameter (ipr_default = 1)
- double precision omega_default
- parameter (omega_default = 1.2d0)
- double precision tolrsd1_def, tolrsd2_def, tolrsd3_def,
- > tolrsd4_def, tolrsd5_def
- parameter (tolrsd1_def=1.0e-08,
- > tolrsd2_def=1.0e-08, tolrsd3_def=1.0e-08,
- > tolrsd4_def=1.0e-08, tolrsd5_def=1.0e-08)
-
- double precision c1, c2, c3, c4, c5
- parameter( c1 = 1.40d+00, c2 = 0.40d+00,
- > c3 = 1.00d-01, c4 = 1.00d+00,
- > c5 = 1.40d+00 )
-
-c---------------------------------------------------------------------
-c grid
-c---------------------------------------------------------------------
- integer nx, ny, nz
- integer nx0, ny0, nz0
- integer ipt, ist, iend
- integer jpt, jst, jend
- integer ii1, ii2
- integer ji1, ji2
- integer ki1, ki2
- double precision dxi, deta, dzeta
- double precision tx1, tx2, tx3
- double precision ty1, ty2, ty3
- double precision tz1, tz2, tz3
-
- common/cgcon/ dxi, deta, dzeta,
- > tx1, tx2, tx3,
- > ty1, ty2, ty3,
- > tz1, tz2, tz3,
- > nx, ny, nz,
- > nx0, ny0, nz0,
- > ipt, ist, iend,
- > jpt, jst, jend,
- > ii1, ii2,
- > ji1, ji2,
- > ki1, ki2
-
-c---------------------------------------------------------------------
-c dissipation
-c---------------------------------------------------------------------
- double precision dx1, dx2, dx3, dx4, dx5
- double precision dy1, dy2, dy3, dy4, dy5
- double precision dz1, dz2, dz3, dz4, dz5
- double precision dssp
-
- common/disp/ dx1,dx2,dx3,dx4,dx5,
- > dy1,dy2,dy3,dy4,dy5,
- > dz1,dz2,dz3,dz4,dz5,
- > dssp
-
-c---------------------------------------------------------------------
-c field variables and residuals
-c---------------------------------------------------------------------
- double precision u(5,-1:isiz1+2,-1:isiz2+2,isiz3),
- > rsd(5,-1:isiz1+2,-1:isiz2+2,isiz3),
- > frct(5,-1:isiz1+2,-1:isiz2+2,isiz3),
- > flux(5,0:isiz1+1,0:isiz2+1,isiz3)
-
- common/cvar/ u,
- > rsd,
- > frct,
- > flux
-
-
-c---------------------------------------------------------------------
-c output control parameters
-c---------------------------------------------------------------------
- integer ipr, inorm
-
- common/cprcon/ ipr, inorm
-
-c---------------------------------------------------------------------
-c newton-raphson iteration control parameters
-c---------------------------------------------------------------------
- integer itmax, invert
- double precision dt, omega, tolrsd(5),
- > rsdnm(5), errnm(5), frc, ttotal
-
- common/ctscon/ dt, omega, tolrsd,
- > rsdnm, errnm, frc, ttotal,
- > itmax, invert
-
- double precision a(5,5,isiz1,isiz2),
- > b(5,5,isiz1,isiz2),
- > c(5,5,isiz1,isiz2),
- > d(5,5,isiz1,isiz2)
-
- common/cjac/ a, b, c, d
-
-c---------------------------------------------------------------------
-c coefficients of the exact solution
-c---------------------------------------------------------------------
- double precision ce(5,13)
-
- common/cexact/ ce
-
-c---------------------------------------------------------------------
-c multi-processor common blocks
-c---------------------------------------------------------------------
- integer id, ndim, num, xdim, ydim, row, col
- common/dim/ id,ndim,num,xdim,ydim,row,col
-
- integer north,south,east,west
- common/neigh/ north,south,east, west
-
- integer from_s,from_n,from_e,from_w
- parameter (from_s=1,from_n=2,from_e=3,from_w=4)
-
- integer npmax
- parameter (npmax=isiz01+isiz02)
-
- logical icommn(npmax+1),icomms(npmax+1),
- > icomme(npmax+1),icommw(npmax+1)
- double precision buf(5,2*isiz2*isiz3),
- > buf1(5,2*isiz2*isiz3)
-
- common/comm/ buf, buf1,
- > icommn,icomms,
- > icomme,icommw
-
- double precision maxtime
- common/timer/maxtime
-
-
-c---------------------------------------------------------------------
-c end of include file
-c---------------------------------------------------------------------