| 1 | c
|
|---|
| 2 | c Points (z) and weights (w) on velocity, pressure
|
|---|
| 3 | c
|
|---|
| 4 | c zgl -- velocity points on Gauss-Lobatto points i = 1,...nx
|
|---|
| 5 | c zgp -- pressure points on Gauss points i = 1,...nxp (nxp = nx-2)
|
|---|
| 6 | c
|
|---|
| 7 |
|
|---|
| 8 | c integer lxm ! defined in HSMG
|
|---|
| 9 | c parameter (lxm = lx1)
|
|---|
| 10 | integer lxq
|
|---|
| 11 | parameter (lxq = lx2)
|
|---|
| 12 | c
|
|---|
| 13 | real zgl(lx1), wgl(lx1), zgp(lx1), wgp(lxq)
|
|---|
| 14 | common /wz1/ zgl , wgl , zgp , wgp
|
|---|
| 15 | c
|
|---|
| 16 | c Tensor- (outer-) product of 1D weights (for volumetric integration)
|
|---|
| 17 | c
|
|---|
| 18 | real wgl1(lx1*lx1), wgl2(lxq*lxq), wgli(lx1*lx1)
|
|---|
| 19 | common /wz2/ wgl1 , wgl2 , wgli
|
|---|
| 20 | c
|
|---|
| 21 | c
|
|---|
| 22 | c Frequently used derivative matrices:
|
|---|
| 23 | c
|
|---|
| 24 | c D1, D1t --- differentiate on mesh 1 (velocity mesh)
|
|---|
| 25 | c D2, D2t --- differentiate on mesh 2 (pressure mesh)
|
|---|
| 26 | c
|
|---|
| 27 | c DXd,DXdt --- differentiate from velocity mesh ONTO dealiased mesh
|
|---|
| 28 | c (currently the same as D1 and D1t...)
|
|---|
| 29 | c
|
|---|
| 30 | c
|
|---|
| 31 | real d1 (lx1*lx1) , d1t (lx1*lx1)
|
|---|
| 32 | $ , d2 (lx1*lx1) , b2p (lx1*lx1)
|
|---|
| 33 | $ , B1iA1 (lx1*lx1) , B1iA1t (lx1*lx1)
|
|---|
| 34 | $ , da (lx1*lx1) , dat (lx1*lx1)
|
|---|
| 35 | $ , iggl (lx1*lxq) , igglt (lx1*lxq)
|
|---|
| 36 | $ , dglg (lx1*lxq) , dglgt (lx1*lxq)
|
|---|
| 37 | $ , wglg (lx1*lxq) , wglgt (lx1*lxq)
|
|---|
| 38 | common /deriv/ d1,d1t,d2,b2p,B1iA1,B1iA1t
|
|---|
| 39 | $ ,da,dat,iggl,igglt,dglg,dglgt,wglg,wglgt
|
|---|
