#######################################################################################################################################
# #
# M O U V E M E N T B R O W N I E N E T S A F R O N T I E R E #
# E N C O O R D O N N E E S S P H E R I Q U E S : #
# #
# #
# Author of '$xiirv/.BROW.31.1.$U' : #
# #
# Jean-Francois Colonna (LACTAMME, 20110320105910). #
# #
#######################################################################################################################################
$Z setParam _____NPoints 50000
$Z setParam _____Loi "$K_VIDE"
$Z # Argument introduit le 20110408104112... #
$Z setParam _____BorneInfR 0
$Z setParam _____BorneSupR +0.01
$Z setParam _____BorneInfP 0
$Z setParam _____BorneSupP $dpi
$Z setParam _____BorneInfT 0
$Z setParam _____BorneSupT $pi
$Z SETParam _____Couleur $NEXIST
$Z setParam _____Rayon 0.004
$Z setParam _____Zoom 1.15
$Z setParam _____Seuil 9
$Z $DELETE $xTV/POINTS.?.* >& $nul
$Z set PremierPoint=1
$Z set PointOrigine=0
$Z $xci/valeurs_alea$X \
$Z premiere=$PremierPoint \
$Z derniere=$_____NPoints \
$Z inf=$_____BorneInfR sup=$_____BorneSupR \
$Z $_____Loi \
$Z graine=$_____Graine_R \
$Z > $xTV/POINTS.1.RHO
$Z $xci/valeurs_alea$X \
$Z premiere=$PremierPoint \
$Z derniere=$_____NPoints \
$Z inf=$_____BorneInfP sup=$_____BorneSupP \
$Z graine=$_____Graine_P \
$Z > $xTV/POINTS.1.PHI
$Z $xci/valeurs_alea$X \
$Z premiere=$PremierPoint \
$Z derniere=$_____NPoints \
$Z inf=$_____BorneInfT sup=$_____BorneSupT \
$Z graine=$_____Graine_T \
$Z > $xTV/POINTS.1.THETA
$Z $xrv/RPT_XYZ.01$X \
$Z ne=0 \
$Z LISTE_RHO=$xTV/POINTS.1.RHO \
$Z LISTE_PHI=$xTV/POINTS.1.PHI \
$Z LISTE_THETA=$xTV/POINTS.1.THETA \
$Z pX=1 pY=0 pZ=0 \
$Z > $xTV/POINTS.1$COORD_X
$Z $xrv/RPT_XYZ.01$X \
$Z ne=0 \
$Z LISTE_RHO=$xTV/POINTS.1.RHO \
$Z LISTE_PHI=$xTV/POINTS.1.PHI \
$Z LISTE_THETA=$xTV/POINTS.1.THETA \
$Z pX=0 pY=1 pZ=0 \
$Z > $xTV/POINTS.1$COORD_Y
$Z $xrv/RPT_XYZ.01$X \
$Z ne=0 \
$Z LISTE_RHO=$xTV/POINTS.1.RHO \
$Z LISTE_PHI=$xTV/POINTS.1.PHI \
$Z LISTE_THETA=$xTV/POINTS.1.THETA \
$Z pX=0 pY=0 pZ=1 \
$Z > $xTV/POINTS.1$COORD_Z
$Z # On notera le 20110320112407 que si : #
$Z # #
$Z # THETA = pi/2 #
$Z # #
$Z # pour tout {RHO,PHI}, alors on se retrouve dans le cas bidimensionnel avec : #
$Z # #
$Z # Z = 0 #
$Z # #
$Z $xrv/accumule.01$X \
$Z ne=0 \
$Z fichier=$xTV/POINTS.1$COORD_X \
$Z translater=VRAI \
$Z > $xTV/POINTS.2$COORD_X
$Z $xrv/accumule.01$X \
$Z ne=0 \
$Z fichier=$xTV/POINTS.1$COORD_Y \
$Z translater=VRAI \
$Z > $xTV/POINTS.2$COORD_Y
$Z $xrv/accumule.01$X \
$Z ne=0 \
$Z fichier=$xTV/POINTS.1$COORD_Z \
$Z translater=VRAI \
$Z > $xTV/POINTS.2$COORD_Z