search

geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
223 stars 235 forks source link

Adiabatic boundary get_data_component() function assert gives an error #3970

Open alarshi opened 3 years ago

alarshi commented 3 years ago

I am working on a model in which adiabatic conditions are computed based on the material model and an input ascii profile (reference_profile.cc). The material model, (equiliblirum_grain_size.cc), calls the get_data_component() function in AdiabaticBoundary class defined in utilities.cc file. However, to initialize the adiabatic conditions, the material model equiliblirum_grain_size.cc is called, which throws an error in the assert of get_data_component(): this->get_time() - first_data_file_model_time >= 0.0 (ln 2450 in utilities.cc). Below is the error message I get: image

anne-glerum commented 3 years ago

It could be that this->get_time() is not set yet, see lines 834-837, where I had to add a specific case to circumvent this: https://github.com/geodynamics/aspect/blob/f62dedc6d3feee5680d10f756bc37f6dec18294a/source/structured_data.cc#L834