search

UniMainzGeo / LaMEM

LaMEM (Lithosphere and Mantle Evolution Model) - a 3D parallel code to simulate geodynamic and geomechanical processes
MIT License
50 stars 18 forks source link

Sedimentation not active when run LaMEM using Julia #20

Open MarcGuardia opened 3 months ago

MarcGuardia commented 3 months ago

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 _sedimentmodel = 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.

Untittled-1

Thank you

MArc

boriskaus commented 2 months ago

@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?