[CIG-SHORT] Problem setting initial conditions

[Oliver Boyd]

Hi Brad,
Thanks for your help. I had initially tried using db_initial_strain as
well as several of other things. I did as you suggested below, but do not
see the initial strain field being set after running step2Kin. Thanks for
reminding me of the use of --component.help. I was trying to figure out
what variables I needed to set in db_initial_strain. Based on the manual I
thought it was strain-xx Š, but I wasn't sure. I thought I might be able
to get to it with help-properties, but I still haven't figured that out.
In any case, I'm still doing something wrong.

Thanks again for your help,
Oliver

On 2/2/14 11:07 AM, "Brad Aagaard" <baagaard at usgs.gov> wrote:

>Oliver,
>
>There are separate spatial databases for initial stress and initial
>strain.
>
>[pylithapp.timedependent.materials.upper_crust]
>db_initial_stress = spatialdata.spatialdb.SimpleDB
>db_initial_stress.label = Initial stress in upper crust
>db_initial_stress.iohandler.filename =
>spatialdb/afterslip_stress_uppercrust.spatialdb
>db_initial_stress.query_type = nearest
>
># ADD db_initial_strain
>db_initial_strain = spatialdata.spatialdb.SimpleDB
>db_initial_strain.label = Initial strain in upper crust
>db_initial_strain.iohandler.filename =
>spatialdb/afterslip_strain_uppercrust.spatialdb
>db_initial_strain.query_type = nearest
>
>Remember you can use pylith --COMPONENT.help and
>--COMPONENT.help-properties and --COMPONENT.help-components to
>interrogate what properties and components are available. In this case
>COMPONENT would be
>timedependent.materials.upper_crust.
>
>pylith step1Kin.cfg --timedependent.materials.upper_crust.help-components
>facilities of 'elasticisotropic3d':
>     db_initial_state=<component name>: Database for initial state
>variables.
>         current value: 'nullcomponent', from {default}
>         configurable as: nullcomponent, db_initial_state
>     db_initial_strain=<component name>: Database for initial strain.
>         current value: 'nullcomponent', from {default}
>         configurable as: nullcomponent, db_initial_strain
>     db_initial_stress=<component name>: Database for initial stress.
>         current value: 'nullcomponent', from {default}
>         configurable as: nullcomponent, db_initial_stress
>     db_properties=<component name>: Database for physical property
>parameters.
>         current value: 'simpledb', from {default}
>         configurable as: simpledb, db_properties
>     output=<component name>: Output manager for elastic material
>information.
>         current value: 'outputmatelastic', from {default}
>         configurable as: outputmatelastic, output
>     perf_logger=<component name>: Performance and memory logging.
>         current value: 'memory_logger', from {default}
>         configurable as: memory_logger, perf_logger
>     quadrature=<component name>: Quadrature object for numerical
>integration.
>         current value: 'quadrature', from {default}
>         configurable as: quadrature
>
>Brad
>
>On 1/31/14 11:26 AM, Oliver Boyd wrote:
>> Hello,
>> I'm trying to generate a set of initial conditions in one problem and
>> then use them in another. The first problem calculates the stress
>> changes due to an earthquake. Currently, the second problem only looks
>> at the viscoelastic response (which of course can be done in the first
>> problem). Once I get the viscoelastic response in the second problem to
>> equal the first (if I let the first progress through time), I will
>> change the fault interface model, but this is not yet implemented.
>>
>> I'm able to transfer over the initial stress through db_initial_stress
>> but I don't know how to transfer the initial strain. I've tried putting
>> the initial strain in the db_initial_stress.iohandler.filename file with
>> the components given by p71 of the manual, but they don't appear to have
>> an effect when I run the simulation. Attached are the files needed for
>> the simulation.
>>
>> Thanks for your assistance,
>> Oliver
>>
>> Flow
>> pylith step1Kin.cfg
>> python afterslip_initial_parms.py
>> pylith step2Kin.cfg
>>
>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>
>
>_______________________________________________
>CIG-SHORT mailing list
>CIG-SHORT at geodynamics.org
>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short