[CIG-SEISMO] 1000s of receivers SPECFEM Cartesian

[Andrea Colombi]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-seismo/attachments/20130605/f0fc177d/attachment.html>