for profil in seqProfs:
Stream_1.AddProfile(profil)
-Stream_1.SetDDZ( 0.200 )
+Stream_1.SetDDZ( 0.050 )
Stream_1.SetSpatialStep( 0.5 )
Stream_1.Update()
# --- basic properties control: edges length, surfaces
-controlGeomProps(geompy, riveGauche, 355.795218, 7765.650632)
-controlGeomProps(geompy, litMineur, 383.674511, 3410.00875 )
-controlGeomProps(geompy, riveDroite, 419.322218, 2549.630118)
+controlGeomProps(geompy, riveGauche, 355.239855, 7755.628676)
+controlGeomProps(geompy, litMineur, 383.395884, 3432.296956)
+controlGeomProps(geompy, riveDroite, 419.214176, 2537.363869)
#----------------------
# --- Meshing
controlMeshStats(plaine, 1680, 227, 3251)
controlSubMeshStats(litMineur_2, 2400)
-controlSubMeshStats(riveDroite_1, 271)
+controlSubMeshStats(riveDroite_1, 257)
controlSubMeshStats(riveGauche_1, 580)
#----------------------
# --- Z interpolation on the bathymety/altimetry on the mesh nodes
statz = interpolZ(nomCas, fichierMaillage, dicoGroupeRegion, zUndef)
-refstatz = {'riveDroite': (99.8278, 100.8), 'riveGauche': (100.0, 100.0), 'litMineur': (80.4541, 100.8)}
+#print statz
+refstatz = {'riveDroite': (100.0, 100.0, 100.0, 0.0, 100.0, 100.0),
+ 'riveGauche': (100.0, 100.0, 100.0, 0.0, 100.0, 100.0),
+ 'litMineur': (80.35, 100.0, 92.88, 4.82, 84.44, 100.0)}
controlStatZ(statz, refstatz)
# --- add a field on nodes of type double with z values, named "BOTTOM"