Skip to content

Commit

Permalink
fixed sliding law in UQ flow
Browse files Browse the repository at this point in the history
  • Loading branch information
@aaschwanden
aaschwanden committed Nov 16, 2024
1 parent 77de183 commit faa54fe
Show file tree
Hide file tree
Showing 3 changed files with 24 additions and 20 deletions.
4 changes: 4 additions & 0 deletions uq/draw_ragis_samples.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -254,6 +254,10 @@
f"seasonal_calving_id_{int(id)}_1900_2025.nc"
for id in distributions[key].ppf(unif_sample[:, i])
]
elif key == "sliding_law":
dist_sample[:, i] = [
slidinglaw[id] for id in distributions[key].ppf(unif_sample[:, i])
]
elif key == "climate_file":
dist_sample[:, i] = [
climate[id] for id in distributions[key].ppf(unif_sample[:, i])
Expand Down
20 changes: 10 additions & 10 deletions uq/ensemble_gris_flow_lhs_10.csv
View file
Original file line number Diff line number Diff line change
@@ -1,11 +1,11 @@
id,basal_resistance.pseudo_plastic.q,basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden,basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max,basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min,basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min,basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max,sliding_law,stress_balance.sia.enhancement_factor,stress_balance.ssa.Glen_exponent,stress_balance.sia.Glen_exponent,calving.thickness_calving.threshold,calving.vonmises_calving.sigma_max,climate,climate_file,fractures,frontal_melt,ocean.th.gamma_T,hydrology,ocean.models,ocean_file,prescribed_retreat_file
0,0.523519292786216,0.031051901544783107,40.87145962473322,27.071334361188452,-558.8970433443315,1084.1675033975566,1.0,1.1805038984112035,3.1688408974469233,2.1387144889667837,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
1,0.39745400312385293,0.020218124230385135,58.87474175449745,30.138456538815678,-369.5632227845549,441.1377432878015,1.0,2.550645700419452,3.3632054947248182,3.994071123051578,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
2,0.5613539282359795,0.01028402211785213,53.97767491352718,32.15889023427513,-232.8452450726968,1208.3049048371252,0.0,2.9076050640379902,2.8055531177776802,3.398909921502095,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
3,0.7154886157522484,0.03614567165800017,44.436358547346565,11.12389718734676,-624.295232783998,869.5033318304487,0.0,3.1719911201972977,3.478412592710214,1.7891118488926492,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
4,0.6840939199729917,0.0290386671464724,46.25402972357654,14.414484258015072,-99.20435949132911,203.7034528178146,1.0,1.5171180368988817,3.209881121740347,3.527105241385735,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
5,0.9465658611268803,0.02783668174877839,56.32151884728766,24.019967331896623,-745.5653552754864,1042.4758605798747,0.0,3.9211130632614166,3.2916267464653317,1.3215368508858574,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
6,0.28291552014140076,0.022611741336849442,43.35598983838478,20.91837731575358,-188.36025459311304,492.0619367666505,0.0,2.150343033310853,3.0578876397103696,2.6278708703252667,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
7,0.4176089698087678,0.018685138297598523,55.02291365566566,6.613922407853845,-874.4477461801431,706.9273457798915,0.0,3.4078501074537484,2.8735578177624435,2.943780845882936,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
8,0.8887750228190243,0.013749901788512292,50.936621026049686,18.457933276387855,-488.08556807474,1367.7770430885207,1.0,1.8577354424749233,2.9088777021339984,2.355717288744213,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
9,0.795557452516821,0.03791948685808393,48.1664297614296,8.081034736636923,-975.2891917601597,137.4987943956447,1.0,2.3788076562881555,3.0488768043506598,1.0854815799026452,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
0,0.8174403266730802,0.015648408796911595,49.07043643693543,19.05850773343824,-458.5382474386829,659.3563290169802,regularized_coulomb,2.030844112219657,3.452721276321805,1.2118169414322837,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
1,0.8790914804605443,0.02977261716270392,55.93859998561837,33.381582703342836,-200.23042604144655,894.3619648624309,regularized_coulomb,2.586063281008711,3.4120004862959634,3.529703870721692,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
2,0.4801093022964081,0.02186321213093271,40.204449018810884,16.527418818008698,-646.6608554003333,420.00608316786384,pseudo_plastic,2.2073328875327376,3.058713999974761,3.3291975318640783,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
3,0.67756840872587,0.03896752054166095,45.20801086263051,13.800942724366497,-376.5216889902948,503.67461745251893,regularized_coulomb,3.0144983202929017,3.1915664994667985,3.0395215888473124,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
4,0.2651648848643861,0.03546783809019627,50.626873502015656,29.72336757132002,-510.052289787341,1214.7702320061228,pseudo_plastic,1.5626607176526452,3.2739637013104517,3.9548593049116674,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
5,0.42740874165042203,0.010601002747490592,42.40757086832328,7.16191151257691,-703.6260211884612,96.66535011070879,regularized_coulomb,1.0949229330255552,3.3021844030362026,2.3987711386484896,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
6,0.7659387836029681,0.02438227032348007,53.29082976224862,25.74905630769057,-869.1195419723383,219.54396428998086,pseudo_plastic,3.97602820072622,2.7769387397603786,2.1755087573748906,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
7,0.6182860477575036,0.031294808303140007,57.44262762099119,9.649225266845345,-911.5836341105037,1381.2276562767634,pseudo_plastic,1.7600314655494391,2.9370278306529607,1.5398041482121874,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
8,0.32689312320790365,0.02763063516902181,47.276171213345606,27.489279484868316,-121.02998772162539,1142.1074305512614,regularized_coulomb,3.499159086427073,2.9829747255061214,2.5667624334701955,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
9,0.9731021433928932,0.01677902859236864,59.11407946143228,21.900675376764735,-82.03815063337902,920.29966208376,pseudo_plastic,3.371016215069299,2.8534537733239147,1.8525705703791964,50,750000,given_smb,HIRHAM5-monthly-ERA5_1975_2021.nc,false,off,0.0001,diffuse,const,,pism_g450m_frontretreat_calfin_2007.nc
20 changes: 10 additions & 10 deletions uq/gris_flow_lhs_10.csv
View file
Original file line number Diff line number Diff line change
@@ -1,11 +1,11 @@
id,basal_resistance.pseudo_plastic.q,basal_yield_stress.mohr_coulomb.till_effective_fraction_overburden,basal_yield_stress.mohr_coulomb.topg_to_phi.phi_max,basal_yield_stress.mohr_coulomb.topg_to_phi.phi_min,basal_yield_stress.mohr_coulomb.topg_to_phi.topg_min,basal_yield_stress.mohr_coulomb.topg_to_phi.topg_max,sliding_law,stress_balance.sia.enhancement_factor,stress_balance.ssa.Glen_exponent,stress_balance.sia.Glen_exponent
0,0.523519292786216,0.031051901544783107,40.87145962473322,27.071334361188452,-558.8970433443315,1084.1675033975566,1.0,1.1805038984112035,3.1688408974469233,2.1387144889667837
1,0.39745400312385293,0.020218124230385135,58.87474175449745,30.138456538815678,-369.5632227845549,441.1377432878015,1.0,2.550645700419452,3.3632054947248182,3.994071123051578
2,0.5613539282359795,0.01028402211785213,53.97767491352718,32.15889023427513,-232.8452450726968,1208.3049048371252,0.0,2.9076050640379902,2.8055531177776802,3.398909921502095
3,0.7154886157522484,0.03614567165800017,44.436358547346565,11.12389718734676,-624.295232783998,869.5033318304487,0.0,3.1719911201972977,3.478412592710214,1.7891118488926492
4,0.6840939199729917,0.0290386671464724,46.25402972357654,14.414484258015072,-99.20435949132911,203.7034528178146,1.0,1.5171180368988817,3.209881121740347,3.527105241385735
5,0.9465658611268803,0.02783668174877839,56.32151884728766,24.019967331896623,-745.5653552754864,1042.4758605798747,0.0,3.9211130632614166,3.2916267464653317,1.3215368508858574
6,0.28291552014140076,0.022611741336849442,43.35598983838478,20.91837731575358,-188.36025459311304,492.0619367666505,0.0,2.150343033310853,3.0578876397103696,2.6278708703252667
7,0.4176089698087678,0.018685138297598523,55.02291365566566,6.613922407853845,-874.4477461801431,706.9273457798915,0.0,3.4078501074537484,2.8735578177624435,2.943780845882936
8,0.8887750228190243,0.013749901788512292,50.936621026049686,18.457933276387855,-488.08556807474,1367.7770430885207,1.0,1.8577354424749233,2.9088777021339984,2.355717288744213
9,0.795557452516821,0.03791948685808393,48.1664297614296,8.081034736636923,-975.2891917601597,137.4987943956447,1.0,2.3788076562881555,3.0488768043506598,1.0854815799026452
0,0.8174403266730802,0.015648408796911595,49.07043643693543,19.05850773343824,-458.5382474386829,659.3563290169802,regularized_coulomb,2.030844112219657,3.452721276321805,1.2118169414322837
1,0.8790914804605443,0.02977261716270392,55.93859998561837,33.381582703342836,-200.23042604144655,894.3619648624309,regularized_coulomb,2.586063281008711,3.4120004862959634,3.529703870721692
2,0.4801093022964081,0.02186321213093271,40.204449018810884,16.527418818008698,-646.6608554003333,420.00608316786384,pseudo_plastic,2.2073328875327376,3.058713999974761,3.3291975318640783
3,0.67756840872587,0.03896752054166095,45.20801086263051,13.800942724366497,-376.5216889902948,503.67461745251893,regularized_coulomb,3.0144983202929017,3.1915664994667985,3.0395215888473124
4,0.2651648848643861,0.03546783809019627,50.626873502015656,29.72336757132002,-510.052289787341,1214.7702320061228,pseudo_plastic,1.5626607176526452,3.2739637013104517,3.9548593049116674
5,0.42740874165042203,0.010601002747490592,42.40757086832328,7.16191151257691,-703.6260211884612,96.66535011070879,regularized_coulomb,1.0949229330255552,3.3021844030362026,2.3987711386484896
6,0.7659387836029681,0.02438227032348007,53.29082976224862,25.74905630769057,-869.1195419723383,219.54396428998086,pseudo_plastic,3.97602820072622,2.7769387397603786,2.1755087573748906
7,0.6182860477575036,0.031294808303140007,57.44262762099119,9.649225266845345,-911.5836341105037,1381.2276562767634,pseudo_plastic,1.7600314655494391,2.9370278306529607,1.5398041482121874
8,0.32689312320790365,0.02763063516902181,47.276171213345606,27.489279484868316,-121.02998772162539,1142.1074305512614,regularized_coulomb,3.499159086427073,2.9829747255061214,2.5667624334701955
9,0.9731021433928932,0.01677902859236864,59.11407946143228,21.900675376764735,-82.03815063337902,920.29966208376,pseudo_plastic,3.371016215069299,2.8534537733239147,1.8525705703791964

0 comments on commit faa54fe

Please sign in to comment.