Can you try explicitly loading it as an AMReXDataset, as in the snippet in my last email?
Andrew On Wed, Mar 27, 2024, 8:40 PM <[email protected]> wrote: > 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] >
_______________________________________________ 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]
