[aspect-devel] The x, y coordinate in depth_average.cc?

[Ian Rose]

I agree that the formatting of the .gnuplot output is confusing to read. It
is specifically meant for reading with gnuplot, and is not particularly
easy to parse by hand. The negative time values have to do with
constructing a depth-time surface plot (I think?), but you shouldn't take
them too seriously.

The .txt output is kind of meant to remedy this parsing difficulty, so if
you are looking for an easier format to deal with for constructing more
general plots, I would recommend that.

On Mon, Dec 21, 2015 at 10:48 PM, Nan Zhang <Nan.Zhang at colorado.edu> wrote:

> Hi Ian,
>
> May I ask you what are the first two columns as x,y in aspect-1.3? As Timo
> said, the x is depth, and y is the time, why the second column shows some
> negative value?
>
> Which array in depth_average.cc carry on the depth and time, do you know?
>
> On Tue, Dec 22, 2015 at 2:41 PM, Ian Rose <ian.rose at berkeley.edu> wrote:
>
>> I think you would need the development version. Looking at the commit
>> logs, this output option is from October.
>>
>> On Mon, Dec 21, 2015 at 10:32 PM, Nan Zhang <Nan.Zhang at colorado.edu>
>> wrote:
>>
>>> Hi Ian,
>>>
>>> I am using aspect-1.3. So, which latest version I should download?
>>>
>>> Nan
>>>
>>> On Tue, Dec 22, 2015 at 1:34 PM, Ian Rose <ian.rose at berkeley.edu> wrote:
>>>
>>>> Hi Nan,
>>>>
>>>> What version of ASPECT are you using? The txt output format for the
>>>> depth average postprocessor is fairly new, so you may need to update it.
>>>>
>>>> Best,
>>>> Ian
>>>>
>>>> On Mon, Dec 21, 2015 at 7:21 PM, Nan Zhang <Nan.Zhang at colorado.edu>
>>>> wrote:
>>>>
>>>>> Hi Timo,
>>>>>
>>>>> I did a test in shell_simple_2d.prm with:
>>>>> *********************************
>>>>> subsection Depth average
>>>>>     set Time between graphical output = 1e6
>>>>>     set Output format                 = txt
>>>>>   end
>>>>> ...
>>>>> *********************************
>>>>>
>>>>> While, I got error like:
>>>>> **************************************
>>>>> ...
>>>>> Line <81> of file <input string>:
>>>>>     The entry value
>>>>>         txt
>>>>>     for the entry named
>>>>>         Output format
>>>>>     does not match the given pattern
>>>>>         [Selection
>>>>> none|dx|ucd|gnuplot|povray|eps|gmv|tecplot|tecplot_binary|vtk|vtu|hdf5|svg|deal.II
>>>>> intermediate ]
>>>>>
>>>>> ----------------------------------------------------
>>>>> Exception on processing:
>>>>>
>>>>> --------------------------------------------------------
>>>>> An error occurred in line <353> of file
>>>>> </group/director1026/znan/sw_packages/aspect_trilinos/aspect-1.3/source/main.cc>
>>>>> in function
>>>>> ...
>>>>>
>>>>> **************************************
>>>>>
>>>>> I checked the user manual. The format does not include txt.
>>>>>
>>>>>
>>>>>
>>>>> I want to go back to the x,y coordinate columns. The first couple
>>>>> lines in the depth_average.gnuplot:
>>>>> ***********************************
>>>>> ...
>>>>> #
>>>>> # <x> <y> <temperature> <adiabatic_temperature> <velocity_magnitude>
>>>>> <sinking_velocity> <Vs> <Vp> <viscosity>
>>>>> 0 -159126 1422.18 0 3.09299e-20 3.1557e-09 -1 -1 1e+22
>>>>> 285500 -159126 1422.18 0 3.09299e-20 3.1557e-09 -1 -1 1e+22
>>>>>
>>>>> 0 0 1422.18 0 3.09299e-20 3.1557e-09 -1 -1 1e+22
>>>>> 285500 0 1422.18 0 3.09299e-20 3.1557e-09 -1 -1 1e+22
>>>>> ...
>>>>> ***********************************
>>>>> If the x,y is the comment inside the .gnuplot file, what are the
>>>>> numbers "0 -159126 ...".
>>>>>
>>>>> My purpose was just to get a depth_average file for fast plotting.
>>>>>
>>>>>
>>>>> On Mon, Dec 21, 2015 at 11:12 PM, Timo Heister <heister at clemson.edu>
>>>>> wrote:
>>>>>
>>>>>> Nan,
>>>>>>
>>>>>> it is written as x,y in the comment inside the .gnuplot file, but this
>>>>>> is because we use the deal.II functionality to output the data to
>>>>>> different file formats. The first column is the depth while the second
>>>>>> column is the time.
>>>>>>
>>>>>> > Who could tell me where is the x,y array, in entries.push_back
>>>>>> (data_point);?
>>>>>>
>>>>>> entries[i] is the depth average at timestep i (or rather the i-th time
>>>>>> where we do output) and entries[i].values[j][k] is the value of field
>>>>>> j at depth index k.
>>>>>>
>>>>>> Are you trying to parse the data? If yes you might want to look into
>>>>>> "Output format = txt".
>>>>>>
>>>>>> On Mon, Dec 21, 2015 at 9:56 AM, Nan Zhang <Nan.Zhang at colorado.edu>
>>>>>> wrote:
>>>>>> > Hi all,
>>>>>> >
>>>>>> > I am trying to change the output format a little bit for
>>>>>> > depth_average.gnuplot file. I see the first two columns in
>>>>>> > depth_average.gnuplot are x,y. But, when I dug into
>>>>>> depth_average.cc, I
>>>>>> > could not find the x,y coordinate array.
>>>>>> >
>>>>>> > Who could tell me where is the x,y array, in entries.push_back
>>>>>> > (data_point);?
>>>>>> >
>>>>>> > Bests,
>>>>>> > Nan
>>>>>> >
>>>>>> >
>>>>>> >
>>>>>> > _______________________________________________
>>>>>> > Aspect-devel mailing list
>>>>>> > Aspect-devel at geodynamics.org
>>>>>> > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>>>>>>
>>>>>>
>>>>>>
>>>>>> --
>>>>>> Timo Heister
>>>>>> http://www.math.clemson.edu/~heister/
>>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> Aspect-devel mailing list
>>>>> Aspect-devel at geodynamics.org
>>>>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>>>>>
>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20151221/bed0b305/attachment-0001.html>