Dear Andrew,
Thanks a lot!
I have tried it, but it does not read particle field(such as
particle_position_x does not appear),
and only read in the fliud field (p,rho,u,v,vof,w):
yt : [INFO ] 2024-03-28 11:33:30,934 Parameters: current_time
= 5.896837762381966e-07
yt : [INFO ] 2024-03-28 11:33:30,934 Parameters: domain_dimensions
= [128 128 16]
yt : [INFO ] 2024-03-28 11:33:30,934 Parameters: domain_left_edge
= [0. 0. 0.]
yt : [INFO ] 2024-03-28 11:33:30,934 Parameters: domain_right_edge
= [0.4 0.4 0.05]
[('boxlib', 'p'), ('boxlib', 'rho'), ('boxlib', 'u'), ('boxlib', 'v'),
('boxlib', 'vof'), ('boxlib', 'w')]
Below are the Header file in plt***/particles, the AMReX output codes and the
yt postprocess codes, respectively.
**************************************************************************
Here is the Header in the plt00001/particles
Version_Two_Dot_One_double
3
7
lev
u
v
w
ax
ay
az
.......
*****************************************************************************************
*****************************************************************************************
Here is the AMReX output codes:
void
AmrCoreEuler::DataWritePlot(amrex::Real time, int step)
{
const amrex::Vector<std::string> prmvarname = { "rho","u","v","w","p","vof"
};
const amrex::Vector<std::string> ptcvarname = { "lev", "u", "v", "w", "ax",
"ay", "az"};
const std::string& pfname = amrex::Concatenate("./Output/plt", step);
const int ref_ratio = 2;
const amrex::Vector<IntVect> outrf = Vector<IntVect>(finestLevel() + 1,
IntVect{ ref_ratio });
WriteMultiLevelPlotfile(pfname, finestLevel()+1,
amrex::GetVecOfConstPtrs(PrmVar_new),
prmvarname, geom, time, Vector<int>(step, 0),outrf );
if (do_tracers) {
PTCPC->Checkpoint(pfname, "particles", ptcvarname);
}
}
***********************************************************************************************************
*********************************************************
Here is yt postprocess code:
import yt
ds = yt.load(root_dir+'/Plt00001')
print(ds.field_list)
**********************************************************
best wishes,
Yifeng He
_______________________________________________
yt-users mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/yt-users.python.org/
Member address: [email protected]
