How to obtain the gridpoint acceleration in 3DEC

3DEC
1

Hi,
How can I obtain the gridpoint acceleration in 3DEC, similar to what is shown in the FLAC below?
Thanks in advance!

image
image1680×1284 119 KB

2

I want to obtain the peak (maximum) acceleration of each gridpoint in the model during the earthquake process and visualize it using gp.extra.

3

It seems that the FISH function for block.gp.accel is missing. We will add it ASAP. (what version are you using?).

In the short term, you can get it the old fashioned way:

fish def get_acc
  $kxacc=41  ; fortran offset for x acceleration.   
  loop foreach gp block.gp.list
    gpi = block.gp.index(gp)
    ax = memory.fortran.float(gpi+$kxacc)
    ay = memory.fortran.float(gpi+$kxacc + 1)
    az = memory.fortran.float(gpi+$kxacc + 2)
    
    ; this is a vector.  Use math.mag if you want just the magnitude
    block.gp.extra(gp) = vector(ax,ay,az)
  end_loop
end
[get_acc]
5

Hi, Jhazzard
I’m using 3DEC 9.0.177.
Thank you very much for your help! I’ll try the FISH code you provided first.
Have a good day!
JLUJYP