scalaire T_0 ??
scalaire PasDeTemps ??
scalaire NbPasDeTemps ??
scalaire ToleranceNlin 1e-6
scalaire ToleranceLin 1e-11
scalaire MaxItNlin 10
scalaire MaxItLin 5000
scalaire Divergence 1e+10

##  Solveur PETSc  ##

solveurlinpetsc      typesolveur "SolveurLin" direct
solveurlinpetsc      suiviconvergence "SolveurLin" defaut
solveurlinpetsc      configprecondsor  "SolveurLin" [vrai]

solveurinstnlinpetsc typesolveur "SolveurInst"
solveurinstnlinpetsc solveurlin  "SolveurInst" SolveurLin
solveurinstnlinpetsc parametresnlin "SolveurInst" [ToleranceNlin,MaxItNlin]
solveurinstnlinpetsc parametresinst "SolveurInst" [PasDeTemps, NbPasDeTemps, T_0, eulerimplicit]

solveurlinpetsc      typesolveur "SolveurMec" direct
solveurlinpetsc      tolerance   "SolveurMec" [ToleranceLin,ToleranceLin,Divergence,MaxItLin]

##  Definition des champs du probleme

champgeolin "ChampGeo"
scalaire "M_0" 7

scalquad "M" 7
v3dquad "U" [0, 0, 0]

##  Pour les CL  ##

scalaire  "zero"  0.
scalaire "h" -3.2e-4
scalaire "M_infini" 10.


##   Hydrique  ##

{ Densite
	scallin "dessus" 	9.40
	scallin "centre"	7.20
	scallin "dessous"	9.40
}
scalsurgeom "d_b" DESSUS=Densite::dessus DESSOUS=Densite::dessous MILIEU=Densite::centre defaut=Densite::centre

{ Diffusion
	to2symaconstant "dessus" [4.7e-10, 4.7e-10, 4.7e-10, 0, 0, 0]
	to2symaconstant "centre" [5.04e-10, 5.04e-10, 5.04e-10, 0, 0, 0]
	to2symaconstant "dessous" [4.7e-10, 4.7e-10, 4.7e-10, 0, 0, 0]
}
to2surgeom "K_M" DESSUS=Diffusion::dessus DESSOUS=Diffusion::dessous MILIEU=Diffusion::centre defaut=Diffusion::centre

##  Mecanique  ##

{ Coeff
	to4elasticite "dessus" [4.05e+9, 4.05e+9, 4.05e+9, 0.3, 0.3, 0.3, 1.557692308e+9, 1.557692308e+9, 1.557692308e+9]
	to4elasticite "centre" [2.5e+9, 2.5e+9, 2.5e+9, 0.3, 0.3, 0.3, 9.61538461e+8, 9.61538461e+8, 9.61538461e+8]
	to4elasticite "dessous" [4.05e+9, 4.05e+9, 4.05e+9, 0.3, 0.3, 0.3, 1.557692308e+9, 1.557692308e+9, 1.557692308e+9]
}

to4surgeom "Eijkl" DESSUS=Coeff::dessus DESSOUS=Coeff::dessous MILIEU=Coeff::centre defaut=Coeff::centre

{ Dilat
	to2symaconstant "dessus" [7.5e-4, 7.5e-4, 7.5e-4, 0 ,0 ,0]
	to2symaconstant "centre" [4.5e-4, 4.5e-4, 4.5e-4, 0 ,0 ,0]
	to2symaconstant "dessous" [7.5e-4, 7.5e-4, 7.5e-4, 0 ,0 ,0]
}
to2surgeom "Beta" DESSUS=Dilat::dessus DESSOUS=Dilat::dessous MILIEU=Dilat::centre defaut=Dilat::centre