examples/smpi/NAS/LU/applu.incl

   1
   2 c---------------------------------------------------------------------
   3 c---------------------------------------------------------------------
   4 c---  applu.incl
   5 c---------------------------------------------------------------------
   6 c---------------------------------------------------------------------
   7
   8 c---------------------------------------------------------------------
   9 c   npbparams.h defines parameters that depend on the class and
  10 c   number of nodes
  11 c---------------------------------------------------------------------
  12
  13       include 'npbparams.h'
  14
  15 c---------------------------------------------------------------------
  16 c   parameters which can be overridden in runtime config file
  17 c   (in addition to size of problem - isiz01,02,03 give the maximum size)
  18 c   ipr = 1 to print out verbose information
  19 c   omega = 2.0 is correct for all classes
  20 c   tolrsd is tolerance levels for steady state residuals
  21 c---------------------------------------------------------------------
  22       integer ipr_default
  23       parameter (ipr_default = 1)
  24       double precision omega_default
  25       parameter (omega_default = 1.2d0)
  26       double precision tolrsd1_def, tolrsd2_def, tolrsd3_def,
  27      >                 tolrsd4_def, tolrsd5_def
  28       parameter (tolrsd1_def=1.0e-08,
  29      >          tolrsd2_def=1.0e-08, tolrsd3_def=1.0e-08,
  30      >          tolrsd4_def=1.0e-08, tolrsd5_def=1.0e-08)
  31
  32       double precision c1, c2, c3, c4, c5
  33       parameter( c1 = 1.40d+00, c2 = 0.40d+00,
  34      >           c3 = 1.00d-01, c4 = 1.00d+00,
  35      >           c5 = 1.40d+00 )
  36
  37 c---------------------------------------------------------------------
  38 c   grid
  39 c---------------------------------------------------------------------
  40       integer nx, ny, nz
  41       integer nx0, ny0, nz0
  42       integer ipt, ist, iend
  43       integer jpt, jst, jend
  44       integer ii1, ii2
  45       integer ji1, ji2
  46       integer ki1, ki2
  47       double precision  dxi, deta, dzeta
  48       double precision  tx1, tx2, tx3
  49       double precision  ty1, ty2, ty3
  50       double precision  tz1, tz2, tz3
  51
  52       common/cgcon/ dxi, deta, dzeta,
  53      >              tx1, tx2, tx3,
  54      >              ty1, ty2, ty3,
  55      >              tz1, tz2, tz3,
  56      >              nx, ny, nz,
  57      >              nx0, ny0, nz0,
  58      >              ipt, ist, iend,
  59      >              jpt, jst, jend,
  60      >              ii1, ii2,
  61      >              ji1, ji2,
  62      >              ki1, ki2
  63
  64 c---------------------------------------------------------------------
  65 c   dissipation
  66 c---------------------------------------------------------------------
  67       double precision dx1, dx2, dx3, dx4, dx5
  68       double precision dy1, dy2, dy3, dy4, dy5
  69       double precision dz1, dz2, dz3, dz4, dz5
  70       double precision dssp
  71
  72       common/disp/ dx1,dx2,dx3,dx4,dx5,
  73      >             dy1,dy2,dy3,dy4,dy5,
  74      >             dz1,dz2,dz3,dz4,dz5,
  75      >             dssp
  76
  77 c---------------------------------------------------------------------
  78 c   field variables and residuals
  79 c---------------------------------------------------------------------
  80       double precision u(5,-1:isiz1+2,-1:isiz2+2,isiz3),
  81      >       rsd(5,-1:isiz1+2,-1:isiz2+2,isiz3),
  82      >       frct(5,-1:isiz1+2,-1:isiz2+2,isiz3),
  83      >       flux(5,0:isiz1+1,0:isiz2+1,isiz3)
  84
  85       common/cvar/ u,
  86      >             rsd,
  87      >             frct,
  88      >             flux
  89
  90
  91 c---------------------------------------------------------------------
  92 c   output control parameters
  93 c---------------------------------------------------------------------
  94       integer ipr, inorm
  95
  96       common/cprcon/ ipr, inorm
  97
  98 c---------------------------------------------------------------------
  99 c   newton-raphson iteration control parameters
 100 c---------------------------------------------------------------------
 101       integer itmax, invert
 102       double precision  dt, omega, tolrsd(5),
 103      >        rsdnm(5), errnm(5), frc, ttotal
 104
 105       common/ctscon/ dt, omega, tolrsd,
 106      >               rsdnm, errnm, frc, ttotal,
 107      >               itmax, invert
 108
 109       double precision a(5,5,isiz1,isiz2),
 110      >       b(5,5,isiz1,isiz2),
 111      >       c(5,5,isiz1,isiz2),
 112      >       d(5,5,isiz1,isiz2)
 113
 114       common/cjac/ a, b, c, d
 115
 116 c---------------------------------------------------------------------
 117 c   coefficients of the exact solution
 118 c---------------------------------------------------------------------
 119       double precision ce(5,13)
 120
 121       common/cexact/ ce
 122
 123 c---------------------------------------------------------------------
 124 c   multi-processor common blocks
 125 c---------------------------------------------------------------------
 126       integer id, ndim, num, xdim, ydim, row, col
 127       common/dim/ id,ndim,num,xdim,ydim,row,col
 128
 129       integer north,south,east,west
 130       common/neigh/ north,south,east, west
 131
 132       integer from_s,from_n,from_e,from_w
 133       parameter (from_s=1,from_n=2,from_e=3,from_w=4)
 134
 135       integer npmax
 136       parameter (npmax=isiz01+isiz02)
 137
 138       logical icommn(npmax+1),icomms(npmax+1),
 139      >        icomme(npmax+1),icommw(npmax+1)
 140       double precision  buf(5,2*isiz2*isiz3),
 141      >                  buf1(5,2*isiz2*isiz3)
 142
 143       common/comm/ buf, buf1,
 144      >             icommn,icomms,
 145      >             icomme,icommw
 146
 147       double precision maxtime
 148       common/timer/maxtime
 149
 150
 151 c---------------------------------------------------------------------
 152 c   end of include file
 153 c---------------------------------------------------------------------