[CIG-SHORT] Save pylith results only to H5 files

[Francisco Delgado]

Hi Brad, thanks for your help. Right now my code to save the h5 file is
[pylithapp.problem.formulation]
output = [subdomain]
output.subdomain = pylith.meshio.OutputSolnSubset

[pylithapp.problem.formulation.output.subdomain]
writer = pylith.meshio.DataWriterHDF5
label = face_pos #surface
vertex_data_fields = [displacement]
writer.filename = surface.h5

but it still save the vtk files. Is there something else I'm missing??

Thanks

On Tue, Oct 25, 2016 at 12:39 PM, Brad Aagaard <baagaard at usgs.gov> wrote:

> Francisco,
>
> I don't know what you mean by "the following do not turn off the vtk
> writer for the rest of the faces". Both sets of parameters yout listed
> refer to the same component: pylithapp.problem.formulation.output.subdomain.
> As a result, the second set will override the first set of parameters.
>
> The error message you are getting suggests there is something wrong with
> line 109 in the file visco_shell.cfg. The error message seems to indicate
> that line 109 is missing a leading square bracket [. The parameters you
> listed do have the square bracket but I don't know which line corresponds
> to line 109.
>
> If you are trying to turn off output for the entire domain while keeping
> output over the ground surface, then use:
>
> [pylithapp.problem.formulation]
> output = [subdomain]
> output.subdomain = pylith.meshio.OutputSolnSubset
>
> [pylithapp.problem.formulation.output.subdomain]
> writer = pylith.meshio.DataWriterHDF5
> label = face_zpos
> vertex_data_fields = [displacement]
> writer.filename = surface.h5
>
>
> This sets the array of output managers for the solution to an array with
> one item: subdomain. We then set the output manager called 'subdomain' to
> the output manager corresponding to a subset of the solution
> (pylith.meshio.OutputSolnSubset). As a result, there will not be output
> over the domain.
>
> Regards,
> Brad
>
>
>
>
> On 10/25/2016 09:15 AM, Francisco Delgado wrote:
>
>> Dear Pylith developers
>>
>> Some months ago I asked about how can I turn off the VTK writer and use
>> only the H5 one. The following lines work to save the ground
>> displacements in a h5 file
>>
>> # Surface
>> [pylithapp.problem.formulation.output.subdomain]
>> writer = pylith.meshio.DataWriterHDF5
>> label = face_zpos
>> vertex_data_fields = [displacement]
>> writer.filename = surface.h5
>>
>> but the following do not turn off the vtk writer for the rest of the faces
>>
>> [pylithapp.problem.formulation.output.subdomain]
>> writer = pylith.meshio.DataWriterVTK
>> label = face_res
>> vertex_info_fields = []
>> vertex_data_fields = []
>> cell_info_fields = []
>> cell_data_fields = []
>>
>> as they output the following error
>>
>> * >> {command line}::*
>>
>> * >> pyre.inventory(error)*
>>
>> * -- File contains parsing errors: visco_shell.cfg*
>>
>> *[line 109]: 'pylithapp.problem.formulation.output.subdomain]\n'*
>>
>> *pylithapp: configuration error(s)*
>>
>>
>>
>> How can I solve this issue??
>>
>> Thanks
>>
>> --
>> --------------------------------------------------------
>> Francisco Delgado
>> PhD student in Geological Sciences
>> Cornell University
>>
>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>
>>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short




-- 
--------------------------------------------------------
Francisco Delgado
PhD student in Geological Sciences
Cornell University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20161025/240521da/attachment.html>