Open lpk-py opened 7 years ago
Is this community still active? I noticed the last post was in Nov. 2015.
I think fully fluid planets are currently not supported, because we did not yet implement the free surface (which is implicit for the elastic case). I guess we also assume the source is in the elastic part of the planet, hence the error you are getting. Both should not be terribly hard to add though.
@martinvandriel Thanks for the clarification. Look forward to the new features in next versions. If you are planning the support to the fluid Earth, I suggest to add a controlling parameter in inparam_mesh
to make the builtin models to be acoustic automatically.
BTW, is there a limit for the parallel subprocesses, related to the dominant period of meshes? When I tried to divide the model into slices, there could be an error in the OUTPUT
of MESHER
define search sflobal index....
partition sflobal index....
ERORR: having more then 8 neighbours (+myself)
check mesh decomposition in the solid)
iproct = 0
or sometimes
define search sflobal index....
ERROR: too many neighbours in solid
I replaced with zeros all the vs values in
<model>_axisem.bm
output byAxiSEM
in/MESHER/Diags/
, to make an acoustic Earth model. Then assign the file as external model ininparam_mesh
.However, the
SOLVER
exited due to some errors related to source location. Below is theOUTPUT_<RUNNAME>
file on my PC:On the cluster, the
SOLVER
also raise errors due to source (OUTPUT
tells that multiple processes have the source).