[CIG-SHORT] Reg: Time varying load

[Sylvain Barbot]

Hi Anil,

This is quite simple. You need to use multiple events. The load in
subsequent events corresponds to the change of load (the rate of
change of the load multiplied by the time interval). The Relax code
will also ask you for other fault and dyke motion in that time period.
You will have to add zero faults and dyke to subsequent time steps.

The tutorial example to look at is examples/tutorials/run1.sh

Imagine the time series of loads is organized in a file called order.dat

> cat order.dat
# time load_filename
    0.000000    20020729_20020807_km.dat
    0.027397    20020808_20020817_km.dat
    0.054794    20020818_20020827_km.dat
    0.082191    20020828_20020906_km.dat
    0.109589    20020907_20020916_km.dat
...

and each time is given a spatial distribution of loads, for example

> cat 20021007_20021016_km.dat
# n x1 x2 length width t3 T phi
1 -1513.66 -2220.8 100.973 110.83 -1.5582e-05 0 0
2 -1522.32 -2103.36 100.973 110.83 -0.000107212 0 0
3 -1530.49 -1986.67 100.973 110.83 -0.00019306 0 0
4 -1538.19 -1870.68 100.973 110.83 -0.000265286 0 0
5 -1545.43 -1755.36 100.973 110.83 -0.000319284 0 0
6 -1552.22 -1640.67 100.973 110.83 -0.000354564 0 0
7 -1558.59 -1526.57 100.973 110.83 -0.000373968 0 0
8 -1564.52 -1413.02 100.973 110.83 -0.000382494 0 0
...

Then you can calculate the time series of deformation with

> ./relax <<EOF
...
# number of coseismic events (when slip distribution is prescribed)
`wc -l order.dat`
`awk '{if (0!=$1){print $1};print 0;print 0;print 0;system("wc
$2);system("cat $2)}' order.dat`
EOF

Cheers,
Sylvain

On Fri, Aug 4, 2017 at 10:05 AM, Anil Earnest <anilearnest at gmail.com> wrote:
> Dear Sylvain,
>          Is it possible to introduce a time varying load (like hydrological)
> into RELAX runs? i am able to introduce loads using the "number of surface
> loads" option and the
>
> "n xs ys length width t3 T phi " parameters. But how to make this load time
> varying as my load data varies over time and space?
>
> Could you please help.
>
> regards
>
> - Anil