Commit 263708f3 authored by Benoit Urruty's avatar Benoit Urruty

figure reunion

parent 1f5ef171
......@@ -64,14 +64,11 @@
" )\n",
"\n",
"\n",
"# save figure\n",
"urlfigure = '/home/urrutyb/Documents/simu/'\n",
"figure = urlfigure\n",
"# load data 1\n",
"url = '/home/urrutyb/Documents/simu/bedmap_20yrs/'\n",
"url1 = '/home/urrutyb/Documents/simu/bedmachine_MODIFIED/'\n",
"files = glob.glob(url + 'run*.pvtu')\n",
"files_no = sorted(glob.glob(url + 'run*1.pvtu'))\n",
"url='/home/urrutyb/Documents/PhD_TiPACCS/inversion/inversion_Ronnie/step_3_inversion/OPTIM_sansDHDt_avecVisco_R1/mesh_24/'\n",
"files = sorted(glob.glob(url + 'OPTIM_sansDHDT_avecVisco_R1_*.pvtu'))\n",
"figure='/home/urrutyb/Documents/PhD_TiPACCS/inversion/inversion_Ronnie/step_3_inversion/OPTIM_sansDHDt_avecVisco_R1/figure/'\n",
"\n",
"\n",
"# figure=url + 'figure/'\n",
"files = sorted(files)\n",
......@@ -82,64 +79,64 @@
"# for i, block in enumerate(blocks):\n",
"# block[\"node_value\"] = np.full(block.n_points, i)\n",
"\n",
"files = np.sort(glob.glob(url + 'Scalars_run*_.dat'))\n",
"files_1 = np.sort(glob.glob(url + 'INITMIP_Scalar_OUTPUT_run*_.dat'))\n",
"# variable INITMIP_scalar\n",
"Time = []\n",
"Volume = []\n",
"Volume_Floatation = []\n",
"Volume_change = []\n",
"SMB = []\n",
"BMB = []\n",
"residual_flux = []\n",
"Ice_Discharge = []\n",
"Iceflux_GL = []\n",
"GIA = []\n",
"FIA = []\n",
"Free = []\n",
"# variable scalar\n",
"timebis = []\n",
"volumebis = []\n",
"volumegrounded = []\n",
"h = []\n",
"dhdt = []\n",
"smbbis = []\n",
"pdc_melt = []\n",
"h_residual = []\n",
"convective_flux = []\n",
"min_conflux = []\n",
"max_conflux = []\n",
"max_dhdt = []\n",
"min_dhdt = []\n",
"for i in range(0, len(files)):\n",
" var1, var2, var3, var4, var5, var6, var7, var8, var9, var10, var11, var12 = np.loadtxt(\n",
" files_1[i], unpack=True)\n",
" var13, var14, var15, var16, var17, var18, var19, var20, var21, var22, var23, var24, var25 = np.loadtxt(\n",
" files[i], unpack=True)\n",
" Time = np.concatenate((Time, var1))\n",
" Volume = np.concatenate((Volume, var2))\n",
" Volume_Floatation = np.concatenate((Volume_Floatation, var3))\n",
" Volume_change = np.concatenate((Volume_change, var4))\n",
" SMB = np.concatenate((SMB, var5))\n",
" BMB = np.concatenate((BMB, var6))\n",
" residual_flux = np.concatenate((residual_flux, var7))\n",
" Ice_Discharge = np.concatenate((Ice_Discharge, var8))\n",
" Iceflux_GL = np.concatenate((Iceflux_GL, var9))\n",
" GIA = np.concatenate((GIA, var10))\n",
" FIA = np.concatenate((FIA, var11))\n",
" Free = np.concatenate((Free, var12))\n",
" #\n",
" timebis = np.concatenate((timebis, var13))\n",
" volumebis = np.concatenate((volumebis, var14))\n",
" volumegrounded = np.concatenate((volumegrounded, var15))\n",
" h = np.concatenate((h, var16))\n",
" dhdt = np.concatenate((dhdt, var17))\n",
" h_residual = np.concatenate((h_residual, var20))\n",
" convective_flux = np.concatenate((convective_flux, var21))\n",
" min_conflux = np.concatenate((min_conflux, var22))\n",
" max_conflux = np.concatenate((max_conflux, var23))\n",
" max_dhdt = np.concatenate((max_dhdt, var24))\n",
" min_dhdt = np.concatenate((min_dhdt, var25))\n",
"# files = np.sort(glob.glob(url + 'Scalars_run*_.dat'))\n",
"# files_1 = np.sort(glob.glob(url + 'INITMIP_Scalar_OUTPUT_run*_.dat'))\n",
"# # variable INITMIP_scalar\n",
"# Time = []\n",
"# Volume = []\n",
"# Volume_Floatation = []\n",
"# Volume_change = []\n",
"# SMB = []\n",
"# BMB = []\n",
"# residual_flux = []\n",
"# Ice_Discharge = []\n",
"# Iceflux_GL = []\n",
"# GIA = []\n",
"# FIA = []\n",
"# Free = []\n",
"# # variable scalar\n",
"# timebis = []\n",
"# volumebis = []\n",
"# volumegrounded = []\n",
"# h = []\n",
"# dhdt = []\n",
"# smbbis = []\n",
"# pdc_melt = []\n",
"# h_residual = []\n",
"# convective_flux = []\n",
"# min_conflux = []\n",
"# max_conflux = []\n",
"# max_dhdt = []\n",
"# min_dhdt = []\n",
"# for i in range(0, len(files)):\n",
"# var1, var2, var3, var4, var5, var6, var7, var8, var9, var10, var11, var12 = np.loadtxt(\n",
"# files_1[i], unpack=True)\n",
"# var13, var14, var15, var16, var17, var18, var19, var20, var21, var22, var23, var24, var25 = np.loadtxt(\n",
"# files[i], unpack=True)\n",
"# Time = np.concatenate((Time, var1))\n",
"# Volume = np.concatenate((Volume, var2))\n",
"# Volume_Floatation = np.concatenate((Volume_Floatation, var3))\n",
"# Volume_change = np.concatenate((Volume_change, var4))\n",
"# SMB = np.concatenate((SMB, var5))\n",
"# BMB = np.concatenate((BMB, var6))\n",
"# residual_flux = np.concatenate((residual_flux, var7))\n",
"# Ice_Discharge = np.concatenate((Ice_Discharge, var8))\n",
"# Iceflux_GL = np.concatenate((Iceflux_GL, var9))\n",
"# GIA = np.concatenate((GIA, var10))\n",
"# FIA = np.concatenate((FIA, var11))\n",
"# Free = np.concatenate((Free, var12))\n",
"# #\n",
"# timebis = np.concatenate((timebis, var13))\n",
"# volumebis = np.concatenate((volumebis, var14))\n",
"# volumegrounded = np.concatenate((volumegrounded, var15))\n",
"# h = np.concatenate((h, var16))\n",
"# dhdt = np.concatenate((dhdt, var17))\n",
"# h_residual = np.concatenate((h_residual, var20))\n",
"# convective_flux = np.concatenate((convective_flux, var21))\n",
"# min_conflux = np.concatenate((min_conflux, var22))\n",
"# max_conflux = np.concatenate((max_conflux, var23))\n",
"# max_dhdt = np.concatenate((max_dhdt, var24))\n",
"# min_dhdt = np.concatenate((min_dhdt, var25))\n",
"\n",
"\n",
"# load data 2\n",
This diff is collapsed.
......@@ -50,7 +50,7 @@ $$
$$
Mean Viscosity Derivative :
$$
\partial\eta=2\times prefactor \times\muth
\partial\eta=2\times prefactor \times\mu
$$
Friction Parameter :
$$
......@@ -62,11 +62,6 @@ $$
$$
Grounded area coefficient : give the value 1 if it's grounded or 0 if it’s at the grounding line
Flux:
$$
$$
### SSA
......@@ -96,7 +91,7 @@ $$
Derivative of the cost function
and the derivative of the cost function.
$$
velocityb = \sum_1^{Nobs}{(U_{SSA} - U_{Obs})}
$$
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment