geopsy.org

Hi Marc, yes I got the procedure - let me explain by example in the target I specify the CWD, let's say CURRDIR, and some callable script, let's name it 'forward.sh' for convenience. In forward.sh I call some binary with options, e.g. myforward -c config.cfg -i input.inp -o misfit usually input.inp ...

Hi Marc, thanks for checking - I'll probably have to see the example in order to understand how things work now - using a shell script in the target I expect a file 'parameters' to be created in the current working directory as indicated in the target. The script will read the file parameters callin...

Dear Marc and all, I just stumbled over the following problem using DinverExt plugin: The sampler creates hexadress like named subdirectories for each thread to output the 'parameters' file . Using the target I usually call then a bash or python script to read the parameters file into some code for ...

Hi Marc, ok - all fine - I'll go back to an earlier version - I should still have 2.0.3 or so somewhere... an ls in /usr/local gives me: /usr/local/sesarray /usr/local/sesarray-2.0.0-snapshot-200701005 /usr/local/sesarray-2.0.0-snapshot-20070801 /usr/local/sesarray-2.0.0-snapshot-20070810 /usr/local...

geopsy -db /home/mao/DATAPROC/LUENEBURG/lg20070904.gpy -tool geopsyarray -slot 0 -group ost-2007090400-z -param ost-00-batch-fk-28012009.log here I copy&paste the param file: ### Parameters ### WINDOW MIN LENGTH (s) = 80 WINDOW MAX LENGTH (s) = 80 WINDOW LENGTH TYPE (at least/exactly/freq. dep.)...

Hi Marc, I just noticed that something apparently changed in the last version for the batch processing - I try to run the same (let's better say a similar) script for batch processing as I did last November using version 2.0.3 (I think) - now I get: ----WARNING--- Loading signal ...---- Unknown form...

Hi Sylvette, I think you need to pipe the output of your command to gpell and then the ellipticites are computed for the models which produced a misfit less than xxx - I assume that the -e option can only work, if you had done an inversion with ellipticity target - if not, then the ellipticity is pr...

Hi Marc,

thanks - yes finally I noticed that there was an additional column with
either 0 or 1 - thanks for the reply - I couldn't remember that I put it here on
the forum.

tbye
Matthias

for whom it concerns: if scanning the minimum misfit models from several reports you need to avoid some additional header lines - I used: gpdcreport run*report -pm |\ awk '$1!="#" {nd=$2;print $(3+nd) " " $1}' |\ sort -n -k 1 | head -200 so the only difference is to suppress thos...

Hi Marc, I tried today the Filter button on the bottom of max2curve tool. It behaves as expected at first glance - it removes the histogram outside the specified k-bounds. Later I learned that this filtering is not permanent - it seems not possible to apply gird statistics on the filtered map... the...

Hi Giuseppe, the resampling messages was probably related to the autocorrelation curves. If you have computed them independently (e.g. from different data sets) with different frequency sampling (e.g. one data from 1 to 5 Hz with 50 samples log, another from 1.4 to 10. Hz with 35 samples in log freq...

Dear leedward, have you tried a small damping on your cross spectral matrix before inverting. I think you can try easily a simple diagonal loading and you will get more stable results then - it will come at the expense of reduced resolution capabilities. The reponse will look then in between the und...

you can simply pipe the output into figue for plotting if you like -
or once you have created a new report containing the best nnn models
you can also use gpdcreport and pipe the output to figue for plotting

bye
Matthias

One question which might be of more borad intrest to others, too. When using dinver with one's own forward problem, then it is possible to use the plugins as you describe in the documentation. What, if we would just like to use the importance sampling starting from an already computed set of models ...