Sedimentation not active when run LaMEM using Julia

[MarcGuardia]

Hello,

I am trying to run simulations on salt tectonics with erosion-sedimentation active using a julia script. The problem is that when I start the simulation, the sedimentation was not happening.

I checked the julia script and the sedimentation is supposed to be active:

sediment_model = 1,
sed_num_layers = 9,
sed_time_delims = [0, 3, 6, 9, 12, 15, 18, 21, 24],
sed_rates = [0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004, 0.004],
sed_phases = [5, 6, 5, 6, 5, 6, 5, 6, 5],
sed_levels = [3, 3, 3, 3, 3, 3, 3, 3, 3, 3]),

But the .dat file that is generated/writted from the .jl file which uses LaMEM to run the simulation does not include the line sediment_model = 1. Then, LaMEM is not detecting active the sedimentation. I attached captures from the .jl original script and the generated .dat file used by LaMEM. It happens in both Linux PC and Windows PC.

It is not a big issue because the simulations can be run directly from the NAME.dat file, but I thought it was good to inform about it.

Thank you

MArc

[boriskaus]

@MarcGuardia thanks for pointing this out! I think that this is really an issue that concerns LaMEM.jl and not the LaMEM main directory (this one).
In any case, I looked at it and was unable to directly reproduce it. Can you perhaps come up with a minimum working example (MWE) that shows the problem and list that here or on the LaMEM.jl repo?