[aspect-devel] How to use Compositional initial conditions/Ascii data model

[Rene Gassmoeller]

Hello Liu,
most often this error indicates that something is wrong with the format
of the input file, e.g. the number of given points does not match the
number of points in the header line, the number of columns is wrong, or
one value in one column is missing. Please have a look and check in
detail what you changed compared to the original file. In case you do
not find anything, can you try to send a simple example model and file
so I can reproduce the problem? If I can reproduce the problem I can try
to improve the error checks inside the code and provide more useful output.

The third column indicates the value of the first compositional field at
the particular coordinate. In the example case it is between 0 and 1,
because it is a parametrization for how much material of a given
composition is in one place (and for 2 compositions this can be done
with one field between 0 and 1). If you have more different fields you
need to add additional columns to the file (one per field). One might
think adding additional compositions as larger values for one field
would work, but in practice it does not, because if you mix field 0 and
2, you would end up with composition 1, which is clearly not correct.
That is why one usually uses (n-1) compositional fields to represent n
compositions.

Best,
Rene

On 05/18/2016 04:36 AM, 刘梦雪 wrote:
> Hi all, 
> When I used the plugin of 'Compositional initial conditions/Ascii data
> model',an error occurred like this,
> ......
> An error occurred in line <634> of file
> </home/lmx/candi/aspect/source/utilities.cc> in function
>     void aspect::Utilities::AsciiDataLookup<dim>::load_file(const
> string&, ompi_communicator_t* const&) [with int dim = 2; std::string =
> std::basic_string<char>; MPI_Comm = ompi_communicator_t*]
> The violated condition was: 
>     i == (components + dim) * data_table.n_elements()
> The name and call sequence of the exception was:
>     ExcMessage (std::string("Number of read in points does not match
> number of expected points. File corrupted ?"))
> Additional Information: 
> Number of read in points does not match number of expected points.
> File corrupted?
> --------------------------------------------------------
>
> Aborting!
> ----------------------------------------------------
> --------------------------------------------------------------------------
> mpirun noticed that the job aborted, but has no info as to the process
> that caused that situation.
> .........
>
> What I did is first take the cookbooks/compositional-passive.prm as an
> example,then changed the geometry to 2m*2m and subsection
>  Compositional initial conditions.After that I copied the box_2d.txt
> from files 'data/compositional-initial-conditions/ascii-data/test',and
> changed the x,y coordinates but still the same number of points.Is
> this error related to the number of points as its too limited ?And
> there is also another thing I'm confu sed that what's the meaning of
> the numbers in the third column in this data file,such as 0 or 1.Does
> it mean that there are two compositional fields?That is,is the number
> of these different numbers the same as the number of compositional
> fields?Hope anyone can solve this ,thanks in advance,have a nice day!
>
> Liu
> 2016/5/18
>
>
>
>
>
> _______________________________________________
> 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/20160518/f6c3d991/attachment-0001.html>