UniMainzGeo
Sedimentation not active when run LaMEM using Julia
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
@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?
