Commit f16437f2 authored by Benoit Urruty's avatar Benoit Urruty

new_mesh

parent 263708f3
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
......@@ -105,7 +105,7 @@
"plt.figure()\n",
"plt.scatter(x,y,c=np.sqrt(blocks[0].point_arrays['uobs'][:,0]**2+blocks[0].point_arrays['uobs'][:,1]**2))\n",
"for n in range(0,len(flowlines)):\n",
" x_fl,y_fl=np.loadtxt(flowlines[n],skiprows=1,delimiter=',',usecols=(32,33),unpack=True)\n",
" x_fl,y_fl=np.loadtxt(flowlines[n],skiprows=1,delimiter=',',usecols=(22,23),unpack=True)\n",
" plt.scatter(x_fl,y_fl)\n",
" eta=np.zeros([len(blocks),len(x_fl)])\n",
" beta=np.zeros([len(blocks),len(x_fl)])\n",
......@@ -256,9 +256,9 @@
" ax1.set_ylabel('SSA velocity (m/an)') \n",
" \n",
" ax2.semilogy(globals()[xx]['dist'],globals()[xx]['eta'][n,],color=colors[n])\n",
" ax2.set_ylabel('eta ($Pa.s$)') \n",
" ax2.set_ylabel('eta ($Pa.a.m$)') \n",
" ax3.semilogy(globals()[xx]['dist'],abs(globals()[xx]['beta'][n,]),color=colors[n])\n",
" ax3.set_ylabel('beta ($Pa. m^{-1/3}.s^{1/3}$)')\n",
" ax3.set_ylabel('beta')\n",
"\n",
" \n",
" ax5.semilogy(globals()[xx]['dist'],abs(globals()[xx]['djdbeta'][n,]),color=colors[n])\n",
This diff is collapsed.
......@@ -98,7 +98,7 @@
"t_inf = time.time()\n",
"\n",
"\n",
"meth='nearest'\n",
"meth='linear'\n",
"for n in range(0,len(flowlines)):\n",
" x_fl,y_fl=np.loadtxt(flowlines[n],skiprows=1,delimiter=',',usecols=(22,23),unpack=True)\n",
" groundedmask=np.zeros([len(blocks),len(x_fl)])\n",
{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"toc": true
},
"source": [
"<h1>Table of Contents<span class=\"tocSkip\"></span></h1>\n",
"<div class=\"toc\"><ul class=\"toc-item\"></ul></div>"
]
},
{
"cell_type": "code",
"execution_count": 10,
......@@ -933,6 +943,19 @@
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.6"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": true,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
......
......@@ -58,7 +58,7 @@ $$
$$
Friction parameter derivative :
$$
\frac{\partial }{\partial\alpha}=\frac{\partial J}{\partial\beta}\times 10 ^ {\alpha}\times ln 10
\frac{\partial J}{\partial\alpha}=\frac{\partial J}{\partial\beta}\times 10 ^ {\alpha}\times ln 10
$$
Grounded area coefficient : give the value 1 if it's grounded or 0 if it’s at the grounding line
......
......@@ -2792,3 +2792,27 @@ ne pas tout faire tourner jusque qq itérations
Entered on [2020-05-26 mar. 14:58]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*regularisation][regularisation]]
** 2020-06 juin
*** 2020-06-04 jeudi
**** reunion
***** mesh là ou le mask est 0 on force output positif
mask = -1 ou 1 output negatif
avec boucle if
distancesolver1: fait le calcul
FEM src module
voit si epaisseur nul
***** sif inversion
faire interpolation sur un noeud exactement
***** linear solver
grep ssa
***** presentation
données immportante pour l'initialisation
Entered on [2020-06-04 jeu. 15:44]
[[file:~/Documents/PhD_TiPACCS/TiPACCS_git/journal/journal.org::*reunion][reunion]]
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