Open liaohaiyang1534 opened 1 year ago
hi Haiyang,
i think you just have to switch the X and Y positions in your STATIONS file.
the STATIONS file format comes originally with a "lat/lon" position in mind (as it all started with regional seismology), so the format typically is something like
#sta #net #LAT #LON #elevation #burial_depth
with a mesh created in Cartesian coordinates, this then translates to LATITUDES going in Y-direction (like "North" being at higher Y values) and LONGITUDES in X-direction (like "going East" is towards positive X values). well, and burial_depth along the Z-coordinate.
in your case, you don't use lat/lon (so likely the UTM conversion is suppressed with SUPPRESS_UTM_PROJECTION == .true.
?), and you specify directly the x/y/z position of the stations. therefore, the #LAT
becomes #Y
and #LON
becomes #X
, i.e., the STATIONS file has a format like
#sta #net #Y #X #elevation #burial_depth
hope that explains it. you can also see the explanation in the doc/wiki: https://github.com/SPECFEM/specfem3d/wiki/A_reference_frame#reference-frame-convention
also, in case the solver cannot place any stations onto the mesh (given they are too far away from the actual mesh location), it will stop running, assuming that it makes no sense in computing waves when there is nothing to record at the end - so, don't waste your compute cycles :)
Hi, daniel
This is helpful, thank you! Although it seems like a very simple question now, I indeed overlooked it due to inertia in my thinking. Thanks again for your patience and attention to detail in your explanation.
the UTM conversion should be suppressed with SUPPRESS_UTM_PROJECTION == .true.
Hi,
I encountered this error: "need at least one receiver"
I am not using the in-house mesher
xmeshfem3D
. I am using thecubit
mesher andCPML
conditions.The size of the mesh is as follows:
CPML layer thickness:
These are the parameters of my source:
The format of my stations is correct, with X direction ranging from 100 to 250.
However, when I change the X coordinate of the station to start from 50, there will be station output files again, but only stations with X = 50 to 80 have data, and the other stations do not.
This is very confusing, and I have searched through the issues and documentation PDF on the GitHub page but still have not found a solution.
I'm very appreciated to have any advice and suggestions for this issue. (If you need any more clarifications, please let me know)
Regards
Haiyang Liao haiyangliao@smail.nju.edu.cn