| 1 | c
|
|---|
| 2 | c Main storage of simulation variables
|
|---|
| 3 | c
|
|---|
| 4 | integer lvt1,lvt2,lbt1,lbt2,lorder2
|
|---|
| 5 | parameter (lvt1 = lx1*ly1*lz1*lelv)
|
|---|
| 6 | parameter (lvt2 = lx2*ly2*lz2*lelv)
|
|---|
| 7 | parameter (lbt1 = lbx1*lby1*lbz1*lbelv)
|
|---|
| 8 | parameter (lbt2 = lbx2*lby2*lbz2*lbelv)
|
|---|
| 9 |
|
|---|
| 10 | parameter (lorder2 = max(1,lorder-2) )
|
|---|
| 11 | c
|
|---|
| 12 | c Solution and data
|
|---|
| 13 | c
|
|---|
| 14 | real bq(lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 15 | common /bqcb/ bq
|
|---|
| 16 |
|
|---|
| 17 | c Can be used for post-processing runs (SIZE .gt. 10+3*LDIMT flds)
|
|---|
| 18 | real vxlag (lx1,ly1,lz1,lelv,2)
|
|---|
| 19 | $ ,vylag (lx1,ly1,lz1,lelv,2)
|
|---|
| 20 | $ ,vzlag (lx1,ly1,lz1,lelv,2)
|
|---|
| 21 | $ ,tlag (lx1,ly1,lz1,lelt,lorder-1,ldimt)
|
|---|
| 22 | $ ,vgradt1(lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 23 | $ ,vgradt2(lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 24 | $ ,abx1 (lx1,ly1,lz1,lelv)
|
|---|
| 25 | $ ,aby1 (lx1,ly1,lz1,lelv)
|
|---|
| 26 | $ ,abz1 (lx1,ly1,lz1,lelv)
|
|---|
| 27 | $ ,abx2 (lx1,ly1,lz1,lelv)
|
|---|
| 28 | $ ,aby2 (lx1,ly1,lz1,lelv)
|
|---|
| 29 | $ ,abz2 (lx1,ly1,lz1,lelv)
|
|---|
| 30 | $ ,vdiff_e(lx1,ly1,lz1,lelt)
|
|---|
| 31 |
|
|---|
| 32 | c Solution data
|
|---|
| 33 | real vx (lx1,ly1,lz1,lelv)
|
|---|
| 34 | $ ,vy (lx1,ly1,lz1,lelv)
|
|---|
| 35 | $ ,vz (lx1,ly1,lz1,lelv)
|
|---|
| 36 | $ ,vx_e (lx1,ly1,lz1,lelv)
|
|---|
| 37 | $ ,vy_e (lx1,ly1,lz1,lelv)
|
|---|
| 38 | $ ,vz_e (lx1,ly1,lz1,lelv)
|
|---|
| 39 | $ ,t (lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 40 | $ ,vtrans (lx1,ly1,lz1,lelt,ldimt1)
|
|---|
| 41 | $ ,vdiff (lx1,ly1,lz1,lelt,ldimt1)
|
|---|
| 42 | $ ,bfx (lx1,ly1,lz1,lelv)
|
|---|
| 43 | $ ,bfy (lx1,ly1,lz1,lelv)
|
|---|
| 44 | $ ,bfz (lx1,ly1,lz1,lelv)
|
|---|
| 45 | $ ,cflf (lx1,ly1,lz1,lelv)
|
|---|
| 46 | $ ,bmnv (lx1*ly1*lz1*lelv*ldim,lorder+1) ! binv*mask
|
|---|
| 47 | $ ,bmass (lx1*ly1*lz1*lelv*ldim,lorder+1) ! bmass
|
|---|
| 48 | $ ,bdivw (lx1*ly1*lz1*lelv*ldim,lorder+1) ! bdivw*mask
|
|---|
| 49 | $ ,c_vx (lxd*lyd*lzd*lelv*ldim,lorder+1) ! characteristics
|
|---|
| 50 | $ ,fw (2*ldim,lelt) ! face weights for DG
|
|---|
| 51 |
|
|---|
| 52 | common /vptsol/ vxlag, vylag, vzlag, tlag, vgradt1, vgradt2,
|
|---|
| 53 | $ abx1, aby1, abz1, abx2, aby2, abz2, vdiff_e,
|
|---|
| 54 | $ vx, vy, vz, t, vtrans, vdiff, bfx, bfy, bfz, cflf, c_vx,fw,
|
|---|
| 55 | $ bmnv, bmass, bdivw,
|
|---|
| 56 | $ vx_e,vy_e,vz_e
|
|---|
| 57 |
|
|---|
| 58 | c Solution data for magnetic field
|
|---|
| 59 | real bx (lbx1,lby1,lbz1,lbelv)
|
|---|
| 60 | $ ,by (lbx1,lby1,lbz1,lbelv)
|
|---|
| 61 | $ ,bz (lbx1,lby1,lbz1,lbelv)
|
|---|
| 62 | $ ,pm (lbx2,lby2,lbz2,lbelv)
|
|---|
| 63 | $ ,bmx (lbx1,lby1,lbz1,lbelv) ! magnetic field rhs
|
|---|
| 64 | $ ,bmy (lbx1,lby1,lbz1,lbelv)
|
|---|
| 65 | $ ,bmz (lbx1,lby1,lbz1,lbelv)
|
|---|
| 66 | $ ,bbx1 (lbx1,lby1,lbz1,lbelv) ! extrapolation terms for
|
|---|
| 67 | $ ,bby1 (lbx1,lby1,lbz1,lbelv) ! magnetic field rhs
|
|---|
| 68 | $ ,bbz1 (lbx1,lby1,lbz1,lbelv)
|
|---|
| 69 | $ ,bbx2 (lbx1,lby1,lbz1,lbelv)
|
|---|
| 70 | $ ,bby2 (lbx1,lby1,lbz1,lbelv)
|
|---|
| 71 | $ ,bbz2 (lbx1,lby1,lbz1,lbelv)
|
|---|
| 72 | $ ,bxlag (lbx1*lby1*lbz1*lbelv,lorder-1)
|
|---|
| 73 | $ ,bylag (lbx1*lby1*lbz1*lbelv,lorder-1)
|
|---|
| 74 | $ ,bzlag (lbx1*lby1*lbz1*lbelv,lorder-1)
|
|---|
| 75 | $ ,pmlag (lbx2*lby2*lbz2*lbelv,lorder2)
|
|---|
| 76 |
|
|---|
| 77 | common /vptsolm/
|
|---|
| 78 | $ bx, by, bz, pm, bmx, bmy, bmz,
|
|---|
| 79 | $ bbx1, bby1, bbz1, bbx2, bby2, bbz2, bxlag, bylag, bzlag,
|
|---|
| 80 | $ pmlag
|
|---|
| 81 |
|
|---|
| 82 | real nu_star
|
|---|
| 83 | common /expvis/ nu_star
|
|---|
| 84 |
|
|---|
| 85 | real pr(lx2,ly2,lz2,lelv), prlag(lx2,ly2,lz2,lelv,lorder2)
|
|---|
| 86 | common /cbm2/ pr, prlag
|
|---|
| 87 |
|
|---|
| 88 | real qtl(lx2,ly2,lz2,lelt), usrdiv(lx2,ly2,lz2,lelt)
|
|---|
| 89 | common /diverg/ qtl, usrdiv
|
|---|
| 90 |
|
|---|
| 91 | real p0th, dp0thdt, gamma0, p0thn, p0thlag(2)
|
|---|
| 92 | common /p0therm/ p0th, dp0thdt, gamma0, p0thn, p0thlag
|
|---|
| 93 |
|
|---|
| 94 | real v1mask (lx1,ly1,lz1,lelv)
|
|---|
| 95 | $ ,v2mask (lx1,ly1,lz1,lelv)
|
|---|
| 96 | $ ,v3mask (lx1,ly1,lz1,lelv)
|
|---|
| 97 | $ ,pmask (lx1,ly1,lz1,lelv)
|
|---|
| 98 | $ ,tmask (lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 99 | $ ,omask (lx1,ly1,lz1,lelt)
|
|---|
| 100 | $ ,vmult (lx1,ly1,lz1,lelv)
|
|---|
| 101 | $ ,tmult (lx1,ly1,lz1,lelt,ldimt)
|
|---|
| 102 | $ ,b1mask (lbx1,lby1,lbz1,lbelv) ! masks for mag. field
|
|---|
| 103 | $ ,b2mask (lbx1,lby1,lbz1,lbelv)
|
|---|
| 104 | $ ,b3mask (lbx1,lby1,lbz1,lbelv)
|
|---|
| 105 | $ ,bpmask (lbx1,lby1,lbz1,lbelv) ! magnetic pressure
|
|---|
| 106 | common /vptmsk/ v1mask,v2mask,v3mask,pmask,tmask,omask,vmult,
|
|---|
| 107 | $ tmult,b1mask,b2mask,b3mask,bpmask
|
|---|
| 108 | c
|
|---|
| 109 | c Solution and data for perturbation fields
|
|---|
| 110 | c
|
|---|
| 111 | real vxp (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 112 | $ ,vyp (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 113 | $ ,vzp (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 114 | $ ,prp (lpx2*lpy2*lpz2*lpelv,lpert)
|
|---|
| 115 | $ ,tp (lpx1*lpy1*lpz1*lpelt,ldimt,lpert)
|
|---|
| 116 | $ ,bqp (lpx1*lpy1*lpz1*lpelt,ldimt,lpert)
|
|---|
| 117 | $ ,bfxp (lpx1*lpy1*lpz1*lpelv,lpert) ! perturbation field rhs
|
|---|
| 118 | $ ,bfyp (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 119 | $ ,bfzp (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 120 | $ ,vxlagp (lpx1*lpy1*lpz1*lpelv,lorder-1,lpert)
|
|---|
| 121 | $ ,vylagp (lpx1*lpy1*lpz1*lpelv,lorder-1,lpert)
|
|---|
| 122 | $ ,vzlagp (lpx1*lpy1*lpz1*lpelv,lorder-1,lpert)
|
|---|
| 123 | $ ,prlagp (lpx2*lpy2*lpz2*lpelv,lorder2,lpert)
|
|---|
| 124 | $ ,tlagp (lpx1*lpy1*lpz1*lpelt,ldimt,lorder-1,lpert)
|
|---|
| 125 | $ ,exx1p (lpx1*lpy1*lpz1*lpelv,lpert) ! extrapolation terms for
|
|---|
| 126 | $ ,exy1p (lpx1*lpy1*lpz1*lpelv,lpert) ! perturbation field rhs
|
|---|
| 127 | $ ,exz1p (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 128 | $ ,exx2p (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 129 | $ ,exy2p (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 130 | $ ,exz2p (lpx1*lpy1*lpz1*lpelv,lpert)
|
|---|
| 131 | $ ,vgradt1p(lpx1*lpy1*lpz1*lpelt,ldimt,lpert)
|
|---|
| 132 | $ ,vgradt2p(lpx1*lpy1*lpz1*lpelt,ldimt,lpert)
|
|---|
| 133 | common /pvptsl/ vxp, vyp, vzp, prp, tp, bqp, bfxp, bfyp, bfzp,
|
|---|
| 134 | $ vxlagp, vylagp, vzlagp, prlagp, tlagp,
|
|---|
| 135 | $ exx1p, exy1p, exz1p, exx2p, exy2p, exz2p,
|
|---|
| 136 | $ vgradt1p, vgradt2p
|
|---|
| 137 |
|
|---|
| 138 | integer jp
|
|---|
| 139 | common /ppointr/ jp
|
|---|
