{"id":22697832,"url":"https://github.com/odinn-sciml/iceflow_sandbox","last_synced_at":"2025-04-13T04:50:32.762Z","repository":{"id":60964396,"uuid":"540474767","full_name":"ODINN-SciML/iceflow_sandbox","owner":"ODINN-SciML","description":"Ice flow model experiments for OGGM-ODINN","archived":false,"fork":false,"pushed_at":"2023-09-27T10:00:33.000Z","size":19228,"stargazers_count":5,"open_issues_count":5,"forks_count":2,"subscribers_count":3,"default_branch":"main","last_synced_at":"2025-03-26T21:46:59.819Z","etag":null,"topics":["cryosphere","differential-eq","glaciers","julia"],"latest_commit_sha":null,"homepage":null,"language":"Jupyter Notebook","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"mit","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/ODINN-SciML.png","metadata":{"files":{"readme":"README.md","changelog":null,"contributing":null,"funding":null,"license":"LICENSE","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null}},"created_at":"2022-09-23T14:18:10.000Z","updated_at":"2024-09-23T11:40:08.000Z","dependencies_parsed_at":"2023-09-27T15:56:37.662Z","dependency_job_id":null,"html_url":"https://github.com/ODINN-SciML/iceflow_sandbox","commit_stats":{"total_commits":109,"total_committers":4,"mean_commits":27.25,"dds":0.1834862385321101,"last_synced_commit":"38ad149711b3aa5a3ed33d92a622cca0960d38cc"},"previous_names":[],"tags_count":0,"template":false,"template_full_name":null,"repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ODINN-SciML%2Ficeflow_sandbox","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ODINN-SciML%2Ficeflow_sandbox/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ODINN-SciML%2Ficeflow_sandbox/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ODINN-SciML%2Ficeflow_sandbox/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/ODINN-SciML","download_url":"https://codeload.github.com/ODINN-SciML/iceflow_sandbox/tar.gz/refs/heads/main","host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":248665756,"owners_count":21142123,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2022-07-04T15:15:14.044Z","host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":["cryosphere","differential-eq","glaciers","julia"],"created_at":"2024-12-10T05:15:52.404Z","updated_at":"2025-04-13T04:50:32.740Z","avatar_url":"https://github.com/ODINN-SciML.png","language":"Jupyter Notebook","funding_links":[],"categories":[],"sub_categories":[],"readme":"# ODINN.jl iceflow sandbox\n\nIce flow model experiments for OGGM-ODINN. Currently under development by Lucille Gimenes.\n\nCurrent developments:\n\nAccelerating ice flow simulations in OGGM with Julia\n- Coupling a 1D Shallow Ice Approximation model in Julia to OGGM.\n- Coupling a 2D Shallow Ice Approximation from ODINN.jl to OGGM.\n\nAt the moment, the options available are :\n- Run a glacier simulation with any type of climate with the Julia SIA model (simply by chosing `cfg.PARAMS[\"evolution_model\"]=\"SIA1D\"` in OGGM) \n- Run a glacier spin up (fixed geometry or dynamic) \n\n## Overview \n\n### Ice flow model 1D \n\nThe new OGGM class (available in the forked repository of OGGM on this user profile : https://github.com/lucillegimenes/oggm ) `IceflowJuliaModel` and its child `SIA1D` calls Julia code (availble in oggm.core \u003e `SIA1D.jl`\u0026 `SIA1D_utils.jl`) and a Julia solver from `DifferentialEquations.jl` to solve the SIA equation on glaciers. \n\n`IceFlowJuliaModel` is a copy of the `FlowlineModel` class, with the methods `run_until` and `run_until_and_store` adapted to be used with the functions `glacier_evolution` and `glacier_evolution_store` from `SIA_1D.jl`. `SIA_1D` has no method `step` compared to `FluxBasedModel` or `SemiImplicitModel`. The new class is now functionnal only with elevation bands (only one flowline) and doesn't have calving implemented. \n\n\n\u003ccenter\u003e\u003cimg src=\"https://github.com/lucillegimenes/iceflow_sandbox/blob/main/plots/OGGM_plots/Siachen_profile.png\" width=\"850\"\u003e\u003c/center\u003e\n\n\u003e Example with Siachen glacier : profile at the end of a GCM climate data run (GCM MRI-ESM2-0 with ssp126). In green is the initial surface of the glacier, and in blue is the solution given by the SemiImplicitModel for the year 2300. In dashed orange is the surface given by the SIA_1D model.\n\n\u003ccenter\u003e\u003cimg src=\"https://github.com/lucillegimenes/iceflow_sandbox/blob/main/plots/OGGM_plots/Siachen_evolution.png\" width=\"900\"\u003e\u003c/center\u003e\n\n\u003e Evolution of volume, area and length of Siachen glacier, under the SSP126 of GCM MRI-ESM2-0\n\n### Benchmarking \n\nTo find which solver was the most efficient, a benchmark was conducted on 12 different glaciers (see notebook Benchmark_from_OGGM)\n\n\u003ccenter\u003e\u003cimg src=\"https://github.com/lucillegimenes/iceflow_sandbox/blob/main/plots/OGGM_benchmarks/random_300y_20r.png\" width=\"900\"\u003e\u003c/center\u003e\n\n\u003e Benchmark conducted with saving annual variables, as well as model geometry and flowline diagnostic variables.\n\nIt appears that the RDPK3Sp35() solver (\"5-stage, third order low-storage scheme with embedded error estimator, optimized for compressible fluid mechanics\")\nis the most efficient solver. \n\n### Assessment of the new model SIA1D with climate experiment \n\n\u003ccenter\u003e\u003cimg src=\"https://github.com/lucillegimenes/iceflow_sandbox/blob/main/plots/OGGM_plots/Experiments/sinus/khumbu_exp_sinus_spin_up_20r_length.png\" width=\"800\"\u003e\u003c/center\u003e\n\n\u003e Evolution of volume, area and length of glacier RGI60-15.03733 (also called Khumbu) under experimental sinusoidal climate represented in bottom subplot.\n\n## In this repository, you will find : \n\n- Notebooks illustrating the results from the iceflow model `SIA1D` : `OGGM_SIA1D_test_(Python).ipynb` is an exemple with simulated glacier Mer de glace, and `Benchmark_from_OGGM_(Python).ipynb` aims to benchmark overall the new iceflow model. Climate experiments are done with `Report_figures.ipynb`. These notebooks will work with `SIA1D` only if the used OGGM version is the one available at https://github.com/lucillegimenes/oggm.\n- Scripts : most of the Julia files available in this directory was used during the early developpements of the iceflow model, but is not useful anymore. In scripts \u003e AD is available a version of `SIA1D` that is compatible with Julia automatic differentiation. \n- Plots obtained from the notebooks\n\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fodinn-sciml%2Ficeflow_sandbox","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fodinn-sciml%2Ficeflow_sandbox","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fodinn-sciml%2Ficeflow_sandbox/lists"}