geodynamics/axisem

velocity seismograms output broken?

martinvandriel opened this issue · 1 comments

Reported by Wangxin:

But when we use use post processing to generate the velocity output. We can change the parameter "SEISTYPE disp" to "SEISTYPE velo" in param_post_processing. When we change the "SEISTYPE velo", there is something wrong. For example: "frottl: severe (24): end-of-file during read, unit 60,...file.../Data/ABU_FNet_velo.dat". I check the source code post_processing.F90. In post_processing.F90, Line 338, the code will load seismograms from directories. If we change the "SEISTYPE velo", it will read data "/Data/*_velo.dat". Actually, after use command "./submit.csh IASP91_exp20s", it will only generate the *_disp.dat in ./IASP91_exp20s/Data.

there is none