| ... | @@ -969,9 +969,9 @@ where (*Â*<sub>*x*</sub>, *Â*<sub>*y*</sub>, *Â*<sub>*z*</sub>) denote | ... | @@ -969,9 +969,9 @@ where (*Â*<sub>*x*</sub>, *Â*<sub>*y*</sub>, *Â*<sub>*z*</sub>) denote | 
|  | the path integrals |  | the path integrals | 
|  |  |  |  | 
|  | over corresponding cell edges. The difference operator |  | over corresponding cell edges. The difference operator | 
|  | *Δ*<sub>`i`*y*</sub>, for instance, is given by |  | *Δ*<sub>`iy`</sub>, for instance, is given by | 
|  |  |  |  | 
|  | *Δ*<sub>`iy`</sub>*Â*<sub>*x*</sub> = *Â*<sub>*x*</sub>(`ix`,`iy`+1,` `iz`) − *Â*<sub>*x*</sub>(`ix`,`iy`,`iz`) |  | *Δ*<sub>`iy`</sub>*Â*<sub>*x*</sub> = *Â*<sub>*x*</sub>(`ix`,`iy`+1,`iz`) − *Â*<sub>*x*</sub>(`ix`,`iy`,`iz`) | 
|  |  |  |  | 
|  | and similarly for others. If the *Â*’s are *unambiguously* computed (or |  | and similarly for others. If the *Â*’s are *unambiguously* computed (or | 
|  | approximated) on coinciding cell edges in the grid hierarchy, the |  | approximated) on coinciding cell edges in the grid hierarchy, the | 
| ... | @@ -1011,15 +1011,15 @@ domain of length *L* (given by `_C.up[0]-_C.lo[0]`): | ... | @@ -1011,15 +1011,15 @@ domain of length *L* (given by `_C.up[0]-_C.lo[0]`): | 
|  |  |  |  | 
|  | 𝜚 = 𝜚<sub>0</sub> |  | 𝜚 = 𝜚<sub>0</sub> | 
|  |  |  |  | 
|  | **m** = 𝜚<sub>0</sub>(−sin(2*π*  *y*/*L*),sin(2*π*  *x*/*L*),0) |  | **m** = 𝜚<sub>0</sub>(−sin(2*π* *y*/*L*),sin(2*π* *x*/*L*),0) | 
|  |  |  |  | 
|  | **B** = *B*<sub>0</sub>(−sin(2*π*  *y*/*L*),sin(4*π*  *x*/*L*),0) |  | **B** = *B*<sub>0</sub>(−sin(2*π* *y*/*L*),sin(4*π* *x*/*L*),0) | 
|  |  |  |  | 
|  | *e* = *p*<sub>0</sub>/(*γ* − 1) + **m**<sup>2</sup>/(2𝜚) + **B**<sup>2</sup>/(2*μ*) |  | *e* = *p*<sub>0</sub>/(*γ*−1) + **m**<sup>2</sup>/(2𝜚) + **B**<sup>2</sup>/(2*μ*) | 
|  |  |  |  | 
|  | where 𝜚<sub>0</sub> = 2.77778, *p*<sub>0</sub> = 5/3 and |  | where 𝜚<sub>0</sub>=2.77778, *p*<sub>0</sub>=5/3 and | 
|  | *B*<sub>0</sub> = 1. The adiabatic index *γ* = 5/3 and the magnetic |  | *B*<sub>0</sub>=1. The adiabatic index *γ*=5/3 and the magnetic | 
|  | permeability *μ* = 1. |  | permeability *μ*=1. | 
|  |  |  |  | 
|  | void configUser(GRD **gm) |  | void configUser(GRD **gm) | 
|  | /** |  | /** | 
| ... | @@ -1092,10 +1092,7 @@ expression *a*<sup>2</sup> + *b*<sup>2</sup> + *c*<sup>2</sup>. | ... | @@ -1092,10 +1092,7 @@ expression *a*<sup>2</sup> + *b*<sup>2</sup> + *c*<sup>2</sup>. | 
|  |  |  |  | 
|  | IC for the shock-cloud interaction problem simulated in a Cartesian box |  | IC for the shock-cloud interaction problem simulated in a Cartesian box | 
|  | given by (*x*, *y*, *z*) ∈ \[ − 1/2, 1/2\]<sup>3</sup>: |  | given by (*x*, *y*, *z*) ∈ \[ − 1/2, 1/2\]<sup>3</sup>: | 
|  | $$(\\varrho,p,v\_x,v\_y,v\_z,B\_x,B\_y,B\_z)=\\left\\{\\begin{array}{ll} |  |  | 
|  | (3.86859,167.345,0,0,0,0,2.1826182,-2.1826182) & x<0.1\\\\ |  |  | 
|  | (1,1,-11.2536,0,0,0,0.56418958,0.56418958) & x\\ge 0.1 |  |  | 
|  | \\end{array}\\right.$$ |  |  | 
|  | At **x** = (0.3, 0, 0) a spherical clump with radius 0.15 and density of |  | At **x** = (0.3, 0, 0) a spherical clump with radius 0.15 and density of | 
|  | 10 is embedded and co-moving with its surrounding flow under the |  | 10 is embedded and co-moving with its surrounding flow under the | 
|  | assumption of pressure equilibrium. The adiabatic index *γ* = 5/3 and |  | assumption of pressure equilibrium. The adiabatic index *γ* = 5/3 and | 
| ... | @@ -1201,7 +1198,7 @@ magnetic permeability *μ* = 1. | ... | @@ -1201,7 +1198,7 @@ magnetic permeability *μ* = 1. | 
|  | return; |  | return; | 
|  | } |  | } | 
|  |  |  |  | 
|  | ### Defining boundary conditions |  | ## Defining boundary conditions | 
|  |  |  |  | 
|  | A superblock `g` includes up to 3 ghost cells per side per coordinate |  | A superblock `g` includes up to 3 ghost cells per side per coordinate | 
|  | direction beyond its active region. The number of ghost cells depends on |  | direction beyond its active region. The number of ghost cells depends on | 
| ... |  | ... |  |