[CIG-SHORT] Velocity value for nodes on free surface

[Brad Aagaard]

Christian,

I follow the general idea of what you are describing but not all of the 
details. Imposing a slip rate is different than imposing a velocity on 
one side of the fault. When imposing a slip rate, the geometry, BC, and 
elastic properties will control how that slip rate is partitioned into 
the relative velocities on each side of the fault. If you question your 
results, you can try simpler geometries, BC, etc to build up your 
intuition and confidence that everything is setup the way you want it.

Regards,
Brad



On 11/19/2013 02:52 AM, Christian Baillard wrote:
> Hi,
>
> I'm in the 2d subduction case, on one the west side I impose 0 displacement for Mantle+Crust whether on the east side I impose 3.5 cm/yr for the crust (mantle fixed), see below.
>
> # The creep rate and slip time are uniform, so use UniformDB for simplicity.
> [pylithapp.timedependent.interfaces.fault_slabbot.eq_srcs.rupture.slip_function]
> slip_rate = spatialdata.spatialdb.UniformDB
> slip_rate.label = Slip rate
> slip_rate.values = [left-lateral-slip, fault-opening]
> slip_rate.data = [-3.5*cm/year, 0.0*cm/year]
>
> When I check the velocities for the western bound of the crust I have velocities around 3 cm/year (the western bound is 250 km away from trench, I impose some locking patches for the subduction interface). Is that coherent to have this value of 3cm/yr instead of 3.5 cm/yr?
>
> Best regards,
>
> Christian
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>