How to extract the uplift rate from the ASPECT-FastScape coupled simulation results

ASPECT
1

Hello,everyone

I am currently working with the ASPECT-FastScape coupling for simulations, and I’ve encountered an issue. I would like to visualize the uplift rate distribution from the coupled simulation results. However, the .VTK output files from FastScape only contain the topography variable. Could you please advise on how I can extract or calculate the uplift rate and create a corresponding visualization?

Thank you very much for your assistance!

Yuan

2

Hi @YUAN,

Apologies for the slow reply to your message and thank you for posting to the forum.

Could you please advise on how I can extract or calculate the uplift rate and create a corresponding visualization?

Would the velocity distribution from the surface of ASPECT provide the information you are looking for?

Cheers,

John

3

Hi Yuan,

If you want to visualize it directly on the VTK files, in the latest version on line 547 in mesh_deformation/fastscape.cc you’ll see we send the bedrock_river_incision_rate into the execute_fastscape command.

      // Find timestep size, run FastScape, and make visualizations.
      execute_fastscape(elevation,
                        bedrock_river_incision_rate_array,  // HHHHH VTk field
                        velocity_x,
                        velocity_y,
                        velocity_z,
                        bedrock_transport_coefficient_array,
                        fastscape_timestep_in_years,
                        fastscape_iterations);

This is actually only used for the visualization, as when you call fastscape’s VTK function, you can send in one array to be visualized on the “HHHHH” field within e.g., paraview, you should be able to change “bedrock_river_incision_rate” to “velocity_z”,

      // Find timestep size, run FastScape, and make visualizations.
      execute_fastscape(elevation,
                        velocity_z,   // HHHHH VTK field
                        velocity_x,
                        velocity_y,
                        velocity_z,
                        bedrock_transport_coefficient_array,
                        fastscape_timestep_in_years,
                        fastscape_iterations);

and then I believe that should work.

Best,

Derek.

4

Sorry for my late reply. Thank you very much for your help.

5

Thank you so much for your reply. I have modified fastscape.cc, and it’s now working.