[CIG-SEISMO] 1000s of receivers SPECFEM Cartesian

[schiemenz at geophysik.uni-muenchen.de]

I can recommend the Seismic Unix option (SU format) in the code.  I used
this option to efficiently handle ~50,000 receivers.

Regards
Alan

> Hello,
>
> when running simulation with SPECFEM cartesian I would like to have in
> output many receivers, lets say ~10000 for relatively
> long simulation (1000000 timestep), in parallel simulation that may
> feature
> 100-1000 cpu. However the program stales because of unclear memory reason
> (not a direct error of SPECFEM) related to the cluster.
>
> My strategy is at the moment to dump the fields at each corner point of
> the
> mesh through netcdf library that can be easily ported on other machines.
> For a single simulation everything works smooth and by adjusting the
> sampling rate the total size of the file remains under control.
> However, in this approach, if the mesh changes, the position of the point
> also changes, complicating eventually the comparison of the results.
>
> Having tested few configuration, I know that I/O of ASCII seismograms is a
> kind of heavy task that slows down the time integration but I wonder if a
> different strategy to do that is available. Not sure if generic binary
> output would help if the files must be read from other potentially
> different machines.
>
> Thanks for the help
> Andrea
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo