Issue regarding the reading of Output files when the rmax0 is set for a GENESIS simulation (wrong size of the grid) about ocelot

[sserkez]

Dear Hector,

Thank you very much for helping us to make ocelot great! (again) =)
At the moment the development of the Ocelot adaptor is focused in iagapov/ocelot/dev_gen branch.
The code discussed looks like this:

    if out('dgrid') == 0:
        rbeam = sqrt(out('rxbeam')**2 + out('rybeam')**2)
        ray = sqrt(out('zrayl') * out('xlamds') / np.pi * (1 + (out('zwaist') / out('zrayl')))**2)
        out.leng = out('rmax0') * (rbeam + ray)
    else:
        out.leng = 2 * out('dgrid')
    out.ncar = int(out('ncar'))  # number of mesh points
    
    #universal solution?
    out.leng=out('meshsize')*(out.ncar-1)

We have also found this bug only recently (I always use dgrid!=0) and decided to go for the most universal solution possible, that is to use the explicit mesh size value, provided in genesis *.out:

flags for output parameter
1 1 1 1 1 0 1 1 1 1 1 0 0 1 0 0 0 0 0
342 entries per record
762 history records
0.1831E-08 wavelength
0.1831E-07 seperation of output slices
101 number of gridpoints
0.1000E-04 meshsize
8192 number of particles
0 particle: records in z
0 particle: records in t

so that out.leng=out('meshsize')*(out.ncar-1)

Thanks again for letting us know, I'll fix the Rayleigh length calculation for the future!

Svitozar