I received some great assistance last time I posted in the forum, so I would like to try again.
I am working on setting up SPECFEM3D_Globe to run with different models, so I am currently working on the /get_timestep_and_layers.f90 file in the /src/shared directory. I can see this is where DT is being calculated based on the model being used and parameters related to the mesh dimensions (NEX_XI and NEX_ETA). However, it looks like there are only some values that are acceptable. When I attempt to run it with say a DT = 1 sec., I receive an mpi error due to a floating point exception. If DT is decreased to less the < 0.16 sec, this problem goes away but I am left with an unreasonably high run time on my machine. Are more processors needed to allow for higher DT values or, if anyone else is familiar with this issue, could this be related to something else?
As always, any advice is greatly appreciated.