read_prt5 python

[info]

Dear all,

I know the read_prt5.m file, and I wonder if someone has developed a version suitable for Python.

Thank you in adavnce

[Salah Benkorichi]

Not up to my knowledge,

I will put it to my queue list to do --Although I'm a bit busy these days.

--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+unsubscribe@googlegroups.com.
To post to this group, send email to fds...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/fds-smv/4f755733-d7ff-4b40-825e-39b5b869f06e%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Topi Sikanen]

This is not exactly the same thing as read_prt.m  but this is a python script that I have used from time to time. Does not work with eight-byte part files.

[info]

Thank you very much!

regards

[WP]

I have extended Topi's python script to be a simple visualizer of evac prt5 file.  May I upload it here?  Or share it on Github?  

Now evac output data (prt5) can be alternatively visualized by pygame (SDL for python) in 2D style.  This could be supplementary to smokeview for visualizing evac output data, and it is useful to draw the agents' compete tracjactory as many users required.  Thank Salah Benkorichi for sending this thread to me.  

In addition, I modified xyz.shape = (7,nplim) for evac prt5 data.  In Topi's script it was initially given by xyz.shape = (3,nplim).  Maybe xyz.shape = (3,nplim) is for Lagrangian Particles.  We use xyz.shape = (7,nplim) for evac prt5 data.  As below is the data format specified in dump_evac() in evac.f90.  Comments are much welcome.  

          NPP = NPP + 1

          TA(NPP) = HR%ILABEL

          XP(NPP) = REAL(HR%X,FB)

          YP(NPP) = REAL(HR%Y,FB)

          ZP(NPP) = MIN( MAX(REAL(HR%Z,FB), EVAC_Z_MIN), EVAC_Z_MAX)

          AP(NPP,1) = 180.0_FB*REAL(HR%Angle/Pi,FB)

          AP(NPP,2) =   2.0_FB*REAL(HR%Radius,FB) 

          AP(NPP,3) =   2.0_FB*REAL(HR%r_torso,FB) 

          ! Height of a human scaled by radius, default male 1.80 m

          AP(NPP,4) =  1.80_FB*REAL(HR%Radius/0.27_EB,FB)

  

Maybe this small evac tool can be integrated in BlenderFDS or any other open-source fds tools in the future.  

Suggestions are much appreciated.  

[Gregor]

I have extended Topi's python script to be a simple visualizer of evac prt5 file.  May I upload it here?  Or share it on Github? 

Share it on GitHub. You can link your script on this wiki page.

Do you know this project (prt5parser)?

[WP]

Yes.  I know Chraibi's work.  In fact we discussed several of his work one year ago, not only prt5parsers, but also severl others.  I notice that his prt5parsers is not listed on the webpage.  Probably it is a small tool by Fortran.  Maybe it is better to first integrate different works together and fully test them with enough test cases.  Anyway, I am glad to see our works to be considered as third-party tools.  Thank you.