[WIP] Fastscape c++ integration with adapter

source/mesh_deformation/fastscapecc.cc

@@ -485,7 +485,7 @@ namespace aspect
           fastscapelib::spl_eroder<FlowGraphType> spl_eroder
           =
             // fastscapelib::make_spl_eroder(flow_graph, kff, n, m, 1e-5)
-            fastscapelib::make_spl_eroder(flow_graph, kff, n, m, 1e-5);
+            fastscapelib::make_spl_eroder(flow_graph, 1e-6, 0.4, 1, 1e-5);
           std::cout << "KFF : " << kff<< std::endl;
 
 
@@ -528,113 +528,10 @@ namespace aspect
 // the other plugins where we don't output information in *every* time step, but only
 // when requested.
 
-    // Create output directory
-    // const std::string output_directory = (this->get_output_directory() + "/fastscapeCC");
-    // if (!std::filesystem::exists(output_directory))
-    //   std::filesystem::create_directories(output_directory);
-
-    // // Vector to track output files
-    // std::vector<std::string> vtu_filenames;
-
-    // const auto write_output = [&](const unsigned int step_index,
-    //                               const xt::xarray<double> &elevation,
-    //                               const xt::xarray<double> &uplift_rate,
-    //                               const xt::xarray<double> &erosion,
-    //                               const xt::xarray<double> &drainage_area)
-    // {
-    //   dealii::Vector<double> elevation_output(n_grid_nodes);
-    //   dealii::Vector<double> uplift_rate_output(n_grid_nodes);
-    //   dealii::Vector<double> erosion_output(n_grid_nodes);
-    //   dealii::Vector<double> drainage_area_output(n_grid_nodes);
-
-    //   for (unsigned int j = 0; j < n_grid_nodes; ++j)
-    //     {
-    //       elevation_output[j]      = elevation[j];
-    //       uplift_rate_output[j]    = uplift_rate[j];
-    //       erosion_output[j]        = erosion[j];
-    //       drainage_area_output[j]  = drainage_area[j];
-    //     }
-
-    //   dealii::DataOut<dim-1, dim> data_out;
-    //   data_out.attach_dof_handler(surface_mesh_dof_handler);
-    //   data_out.add_data_vector(elevation_output, "Elevation");
-    //   data_out.add_data_vector(uplift_rate_output, "UpliftRate");
-    //   data_out.add_data_vector(erosion_output, "Erosion");
-    //   data_out.add_data_vector(drainage_area_output, "DrainageArea");
-    //   data_out.build_patches();
-
-    //   const std::string filename = output_directory + "/fastscape_surface_iteration_" +
-    //                               Utilities::int_to_string(step_index, 4) + ".vtu";
-    //   std::ofstream output(filename);
-    //   data_out.write_vtu(output);
-    //   this->get_pcout() << "➤ Wrote VTK file: " << filename << std::endl;
-
-    //   vtu_filenames.push_back(filename);
-    // };
-
-    // // Output initial state
-    // if (output_internal_fastscape_steps)
-    //   write_output(0, elevation, uplift_rate, xt::zeros<double>({n_grid_nodes}), xt::zeros<double>({n_grid_nodes}));
-
-    // // === FASTSCAPE LOOP ===
-    // for (unsigned int i = 1; i <= fastscape_iterations; ++i)
-    //   {
-    //     std::cout << "\nFastScape iteration " << i << "/" << fastscape_iterations << std::endl;
-
-    //         xt::xarray<double> uplifted_elevation = elevation + fastscape_timestep_in_years * uplift_rate;
-
-    //         std::cout << "[DEBUG] uplifted_elevation max = " << xt::amax(uplifted_elevation)()
-    //                   << ", min = " << xt::amin(uplifted_elevation)() << std::endl;
-
-    //         flow_graph.update_routes(uplifted_elevation);
-    //         flow_graph.accumulate(drainage_area, 1.0);
-
-    //         std::cout << "[DEBUG] drainage_area max = " << xt::amax(drainage_area)()
-    //                   << ", min = " << xt::amin(drainage_area)() << std::endl;
-
-    //         auto spl_erosion = spl_eroder.erode(uplifted_elevation, drainage_area, fastscape_timestep_in_years);
-
-    //         std::cout << "[DEBUG] spl_erosion max = " << xt::amax(spl_erosion)()
-    //                   << ", min = " << xt::amin(spl_erosion)() << std::endl;
-
-    //         //Working for raster grid only
-    //         // auto diff_erosion = diffusion_eroder.erode(uplifted_elevation - spl_erosion, fastscape_timestep_in_years);
-
-    //         sediment_flux = flow_graph.accumulate(spl_erosion);
-    //         elevation = uplifted_elevation - spl_erosion ; //- diff_erosion for raster grid
-
-    //         std::cout << "[DEBUG] updated elevation max = " << xt::amax(elevation)()
-    //                   << ", min = " << xt::amin(elevation)() << std::endl;
-
-    //     if (output_internal_fastscape_steps)
-    //       write_output(i, elevation, uplift_rate, spl_erosion, drainage_area);
-    //   }
-
-    // // Final output (even if internal output was disabled)
-    // if (!output_internal_fastscape_steps)
-    //   write_output(fastscape_iterations, elevation, uplift_rate, xt::zeros<double>({n_grid_nodes}), drainage_area);
-
-    // // === WRITE PVD FILE ===
-    // {
-    //   const std::string pvd_filename = output_directory + "/fastscape_iterations.pvd";
-    //   std::ofstream pvd_file(pvd_filename);
-
-    //   pvd_file << "<?xml version=\"1.0\"?>\n"
-    //           << "<VTKFile type=\"Collection\" version=\"0.1\" byte_order=\"LittleEndian\">\n"
-    //           << "  <Collection>\n";
-
-    //   for (unsigned int i = 0; i < vtu_filenames.size(); ++i)
-    //     {
-    //       pvd_file << "    <DataSet timestep=\"" << i << "\" group=\"\" part=\"0\"\n"
-    //               << "             file=\""
-    //               << "fastscape_surface_iteration_" << Utilities::int_to_string(i, 4) << ".vtu\"/>\n";
-    //     }
-
-    //   pvd_file << "  </Collection>\n"
-    //           << "</VTKFile>\n";
-
-    //   this->get_pcout() << "➤ Wrote PVD file: " << pvd_filename << std::endl;
-    // }
+      // Read parameter
+
+
+
 
 
           // Create output directory