| |
| 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 isiz1,isiz2,isiz3 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 ist, iend |
| integer 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, |
| > ist, iend, |
| > 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 to improve cache performance, second two dimensions padded by 1 |
| c for even number sizes only. |
| c Note: corresponding array (called "v") in routines blts, buts, |
| c and l2norm are similarly padded |
| c--------------------------------------------------------------------- |
| double precision u(5,isiz1/2*2+1, |
| > isiz2/2*2+1, |
| > isiz3), |
| > rsd(5,isiz1/2*2+1, |
| > isiz2/2*2+1, |
| > isiz3), |
| > frct(5,isiz1/2*2+1, |
| > isiz2/2*2+1, |
| > isiz3), |
| > flux(5,isiz1), |
| > qs(isiz1/2*2+1,isiz2/2*2+1,isiz3), |
| > rho_i(isiz1/2*2+1,isiz2/2*2+1,isiz3) |
| |
| common/cvar/ u, rsd, frct, flux, |
| > qs, rho_i |
| |
| |
| 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/2*2+1,isiz2), |
| > b(5,5,isiz1/2*2+1,isiz2), |
| > c(5,5,isiz1/2*2+1,isiz2), |
| > d(5,5,isiz1/2*2+1,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 timers |
| c--------------------------------------------------------------------- |
| integer t_rhsx,t_rhsy,t_rhsz,t_rhs,t_jacld,t_blts, |
| > t_jacu,t_buts,t_add,t_l2norm,t_last,t_total |
| parameter (t_total = 1) |
| parameter (t_rhsx = 2) |
| parameter (t_rhsy = 3) |
| parameter (t_rhsz = 4) |
| parameter (t_rhs = 5) |
| parameter (t_jacld = 6) |
| parameter (t_blts = 7) |
| parameter (t_jacu = 8) |
| parameter (t_buts = 9) |
| parameter (t_add = 10) |
| parameter (t_l2norm = 11) |
| parameter (t_last = 11) |
| logical timeron |
| double precision maxtime |
| |
| common/timer/maxtime,timeron |
| |
| |
| c--------------------------------------------------------------------- |
| c end of include file |
| c--------------------------------------------------------------------- |