Return-Path: Received: from [149.220.60.110] (account huebl@hzdr.de [149.220.60.110] verified) by hzdr.de (CommuniGate Pro SMTP 6.1.12) with ESMTPSA id 15447230 for picongpu-users@hzdr.de; Fri, 13 Jan 2017 16:17:03 +0100 Subject: Re: Current density To: picongpu-users@hzdr.de References: From: Axel Huebl Organization: HZDR Message-ID: <2fb02a83-0461-b47c-2d6c-17f9700bcbcb@hzdr.de> Date: Fri, 13 Jan 2017 16:17:02 +0100 MIME-Version: 1.0 In-Reply-To: Content-Type: text/plain; charset=utf-8 Content-Transfer-Encoding: 7bit Hi Danila, yep, unfortunately - sorry! Or you do a rough estimate from the particle data (position, positionOffset, momentum & weighting), but that is more work then just getting the solver-exact answer directly. Best, Axel On 13.01.2017 15:52, Khikhlukha Danila wrote: > Hi Axel, > thank you for your clear answer! I will try it -- I need to redo the simulation then. > > Best, > Danila. > ________________________________________ > From: picongpu-users@hzdr.de [picongpu-users@hzdr.de] on behalf of Axel Huebl [a.huebl@hzdr.de] > Sent: Friday, January 13, 2017 3:41 PM > To: picongpu-users@hzdr.de > Subject: Re: [PIConGPU-Users] Current density > > Dear Danila, > > thank you for the question! > PIConGPU in 0.2+ writes HDF5 and ADIOS output in the openPMD format. The > best way to get started is via: > > Introduction: http://openPMD.org > Full format: https://github.com/openPMD/openPMD-standard > > and then using the files directly. For data set interpretation and > processing you have several options: > a) investigate the structure just like a file system with tools such as: > - HDF5Viewer > - HDFCompass > - h5py > > b) Use frameworks such as: > - openPMD-viewer https://github.com/openPMD/openPMD-viewer > - postPIC https://github.com/skuschel/postpic > - yt-project (develop branch!) http://yt-project.org/ > > or more tools, see > https://github.com/openPMD/openPMD-projects > > for a full list. > > > The current density itself, our field "J", can be dumped directly. Just > add "FieldJ" in the "fileOutput.param" at the following position: > > > https://github.com/ComputationalRadiationPhysics/picongpu/blob/0.2.1/src/picongpu/include/simulation_defines/param/fileOutput.param#L95-L97 > > typedef MakeSeq< > FieldE, > FieldB, > FieldJ >> ::type NativeFileOutputFields; > > > Best, > Axel > > On 12.01.2017 11:03, Khikhlukha Danila wrote: >> Dear all, >> >> Currently I'm using PoG version 0.2.1. My simulation finished >> successfully and I'm trying to calculate a current density field in the >> given cross section of the computational grid. >> >> However hdf5 files have fields only as follows: >> >> * B >> * E >> * e_chargeDensity >> * e_energyDensity >> * e_particleCounter >> * e_particleMomentumComponent >> >> Is there a way how to get current density field? I guess one can try to >> use e_chargeDensity and e_particleMomentumComponent to calculate it, >> however it is unclear to me which units are used for the momentum field... >> >> Thank you in advance, >> Danila. >> >> >> > > -- > > Axel Huebl > Phone +49 351 260 3582 > https://www.hzdr.de/crp > Computational Radiation Physics > Laser Particle Acceleration Division > Helmholtz-Zentrum Dresden - Rossendorf e.V. > > Bautzner Landstrasse 400, 01328 Dresden > POB 510119, D-01314 Dresden > Vorstand: Prof. Dr.Dr.h.c. R. Sauerbrey > Prof. Dr.Dr.h.c. P. Joehnk > VR 1693 beim Amtsgericht Dresden > > ############################################################# > This message is sent to you because you are subscribed to > the mailing list . > To unsubscribe, E-mail to: > To switch to the DIGEST mode, E-mail to > To switch to the INDEX mode, E-mail to > Send administrative queries to > > > ############################################################# > This message is sent to you because you are subscribed to > the mailing list . > To unsubscribe, E-mail to: > To switch to the DIGEST mode, E-mail to > To switch to the INDEX mode, E-mail to > Send administrative queries to > -- Axel Huebl Phone +49 351 260 3582 https://www.hzdr.de/crp Computational Radiation Physics Laser Particle Acceleration Division Helmholtz-Zentrum Dresden - Rossendorf e.V. Bautzner Landstrasse 400, 01328 Dresden POB 510119, D-01314 Dresden Vorstand: Prof. Dr.Dr.h.c. R. Sauerbrey Prof. Dr.Dr.h.c. P. Joehnk VR 1693 beim Amtsgericht Dresden