Setting up for extract velocity field from an outflow plane(failed)

[potatoslc]

Hi everyone, I am trying to use branch 502 to extract velocity field for my simulation. I looked through the changes made in the codes and made modification on my input field but still failed to have this feature to work. Here is my setup. I am trying to extract a velocity plane at x=0.1(total simulation length is 0.2) and perpendicular to x-axis. Any comments and suggestions are welcomed. Thanks for the help.
peleLM.write_outflow_plane_dir = 1
peleLM.write_outflow_plane_loc = 0.1
peleLM.write_outflow_plane_int = 0.05
peleLM.write_outflow_plane_folder= ./turb_slice_0
peleLM.write_outflow_plane_file_prefix = turb_
peleLM.write_outflow_plane_verbose= 1

[baperry2]

Can you be more specific about what is not working for you? Note that peleLM.write_outflow_plane_loc should be an int corresponding the coarse grid index of the plane you want to save.

[potatoslc]

Hi, I am trying to extract the velocity field of my simulation flow in the rectangular channel (x=0.16 * y=0.02 *z=0.02 meter). The velocity field velocity plane is expected to be extracted is at x=0.08 meter. The new setup like this:
peleLM.write_outflow_plane_dir = 0
peleLM.write_outflow_plane_loc = 128 (I set 256 cells in the x direction)
peleLM.write_outflow_plane_int = 1
Also, I setup check points at PeleLMeX_Evolvo.cpp/PeleLMeX_Advance.cpp by using amrex::Print(). If function go through these points, I should able to see the message in the terminal.
I try to run a test for 5 time steps in total and the get messages:
Unused ParmParse Variables:
[TOP]::amr.chk_int(nvals = 1) :: [500]
[TOP]::amr.plt_per(nvals = 1) :: [0.0001]
[TOP]::peleLM.write_outflow_plane_dir(nvals = 1) :: [0]
[TOP]::peleLM.write_outflow_plane_int(nvals = 1) :: [1]
[TOP]::peleLM.write_outflow_plane_loc(nvals = 1) :: [1]
[TOP]::peleLM.write_outflow_plane_verbose(nvals = 1) :: [1]
Also I did not see any messages be printed out in the terminal, so I assume the function WriteOutflowPlane(m_nstep, m_cur_time); in the PeleLMeX_Evolve.cpp is not triggered for unknown reasons.

[baperry2]

This is weird, those ParmParse variables should always be queried, so they shouldn't be showing up in the unused variables list. Verify that the git hash of PeleLMeX that you are running matches 6ac6dd4. This should be printed at the top of the terminal output when you run the code, should look something like this (may or may not show the v23.11 stuff, the important thing is that the hash matches).

 ================= Build infos =================
 PeleLMeX    git hash: v23.11-130-g6ac6dd4-dirty

Might also be a good idea to cleanup and rebuild just to be sure, make TPLrealclean && make realclean && make TPL && make -j

[potatoslc]

I have PeleLMeX git hash: v23.11-134-g4f842c6-dirty .That's from the Build infos. I think I used the code directly download from turbslice branch. I want to make sure I am using the right branch.

[baperry2]

The turbSlice branch from @drummerdoc's fork is the correct one to be using. But the git hash indicates you are using the development branch. Do git status to see what branch you are on, then probably you need to do git checkout turbSlice to get the correct branch

[potatoslc]

That works after I switch to the turbSlice branch. Thanks for the help.