The OPAL Discussion Forum

[Dirk Vandeplassche <dirkvdpla AT skynet.be>]

Dear Daniel,

Sorry for the late reply. Your py-script has been extremely useful, thank you very much for it !

Best regards,

Dirk

On 16 Jun 2025, at 21:06, Daniel Winklehner <winklehn AT mit.edu> wrote:

Hi Dirk,

 

I have attached a python script that I wrote to replace the C++ program. It may not cover all pertinent cases, but it should give you a starting point. I zipped it because some email providers do not allow sending .py files.

 

Best,

  Daniel

 

 

 

From: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> On Behalf Of Zhang Hui
Sent: Monday, June 16, 2025 16:12
To: Dirk Vandeplassche <dirkvdpla AT skynet.be>; opal AT lists.psi.ch
Cc: opal-developers AT lists.psi.ch
Subject: RE: [Opal] Electric field map in BANDRF cyclotron

 

Dear Dirk,

 

I was transforming ascii file to *.h5part map four years ago. The format of my ascii file looks like the following

 

nx ny nz

nx ny nz

nx ny nz

x1 y1 z1 Ex1 Ey1 Ez1

x1 y1 z2 Ex2 Ey2 Ez2

x1 y1 z3 Ex3 Ey3 Ez3

…

 

The original c++ code doing the transformation must be from Andreas. But I am not sure about that.

 

I would think, that the head could be different. However, the numbers of grids in x, y, y-directions have to be defined in head. And this information has to be correctly transferred into .h5part file.

 

The .h5part map can be visualized by a program called VisIt. If the .h5part file cannot be opened, or it looks strange, then something must be wrong, in the ascii data, and/or during the transformation.

 

Best regards,

 

Hui

 

 

 

__________________________________________
Paul Scherrer Institut 
Hui Zhang
WBGA/C34
Forschungsstrasse 111
5232 Villigen PSI
Schweiz

Telefon: +41563105398
E-Mail: hui.zhang AT psi.ch

From: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> On Behalf Of Dirk Vandeplassche
Sent: Monday, June 16, 2025 1:05 PM
To: opal AT lists.psi.ch
Cc: opal-developers AT lists.psi.ch
Subject: [Opal] Electric field map in BANDRF cyclotron

 

Dear all,

 

It would really be very helpful for me to obtain a detailed description of the ASCII-version of the electric field map to be used in a BANDRF-type cyclotron. I am especially interested in the header of this file, hoping that it may solve the error in OPAL "Couldn't determine type of fieldmap in file “efield.h5part””. The conversion from the ASCII file to the binary file with the “ascii2h5block” tool seems to work well — anyway it does not give any error message and returns a correct number of data with a header like

 

0

101

101

61

0.0005 0.0005 0.0001

-0.025         -0.025         -0.003         502.1066965                 342.215468                   -1294.568941

…

 

Its output is

 

Combine efield2.txt and  to efield2.h5part

Frequency 7.2615e+07 [Hz]

H5block grid

H5gridPx 101

H5gridPy 101

H5gridPz 61

number Edata 622261

gridPx 101 stepEx 0.0005

gridPy 101 stepEy 0.0005

gridPz 61 stepEz 0.0001

number Bdata 0

Done bye ...

 

Any hint that may solve the problem is welcome !

 

Many thanks in advance,

Dirk Vandeplassche

 

<ascii2h5block.zip>