Metric computation for cfd mmg3D and mshmet

Hy all,

I’m trying to to generate an anisotropic mesh in order to capture a shock around an airfoil in a CFD simulation made with SU2.
The result that I’m trying to get is something similar to the following picture (my case is 3D):
mesh1 sol1

I’m using the density flow field in order to compute the metric, using mshmet. I haven’t found much documentation about this software, I hope I’m using it correctly.

the command that I use is:
mshmet -temp.mesh -eps 0.0001 -hmin 0.001 -hmax 10
The dentisy flow field is stored in the file .sol with the same name of the mesh. (temp.sol)
Mshmet saves a file containing the metric. It has 6 columns for each row (vertex).

Then, the new mesh is generate by mmg3d with the followiwng command:
mmg3d_O3 temp.mesh -sol -hgrad 1.5 -hmin 0.001 -hmax 10

Iterating this procedure, I got this result:

It captures only the density gradient near the stagnation point, but doest not capture the shock.
I’ve already tried varying all the parameters, including the flow field for the computation (I also tryied Mach number), but my result is not close to the expected one.

Can anyone tell me what I’m doing wrong?
I was also wondering the difference between hmin and hmax of mshmet and mmg3d.