#######################################################################################################################################
# #
# V I S U A L I S A T I O N D E S V E C T E U R S V I T E S S E S : #
# #
# #
# Nota : #
# #
# On notera que ci-apres dans #
# l'utilisation de '$xrv/particule.10$X' #
# il faut l'option : #
# #
# coordonnees_dans_0_1=VRAI #
# #
# car, en effet, ce ne sont pas des coordonnees #
# {X,Y,Z} issues de '$xrr/N_corps.11$X' que l'on #
# traite (et qui seraient deja dans [0,1], mais #
# des vitesses {VX,VY,VZ} qui n'ont donc pas subies #
# cette normalisation dans '$xrr/N_corps.11$X' #
# #
# #
# Author of '$xivP/disk.000000009/face.2/.NCOR.6.22.$U' : #
# #
# Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss). #
# #
#######################################################################################################################################
$Z SETParam _____Pas 1
$Z # Introduit le 20200517095212... #
$Z setParam _____Jonction VRAI
$Z # Introduit le 20031210145214 dans l'idee de visualiser la trajectoire des vitesses (dans #
$Z # l'espace {VX,VY,VZ}) et non pas l'hodographe des vitesses... #
$Z setParam _____Points 64
$Z setParam _____ParamParticule "$K_VIDE"
$Z # Introduit le 20200515132345... #
$c #define N0 Premiere
$c #define N Derniere
$c #define PAS _____Pas
$c
$c int main()
$c {
$c int n;
$c
$c for (n=N0 ; n<=N ; n=n+PAS)
$c {
$c printf("$DELETE $xTV/COORDONNEES\n");
$c printf("$DELETE $xTV/COORDONNEES$COORD_X\n");
$c printf("$DELETE $xTV/COORDONNEES$COORD_Y\n");
$c printf("$DELETE $xTV/COORDONNEES$COORD_Z\n");
$c printf("$DELETE $xTV/COULEURS$ROUGE\n");
$c printf("$DELETE $xTV/COULEURS$VERTE\n");
$c printf("$DELETE $xTV/COULEURS$BLEUE\n");
$c
$c printf("$CA $xTV/LISTE_POINTS |
$c $GRE 'periode=%d ' |
$c $GRE -v 'CENTRE_DE_GRAVITE='
$c > $xTV/COORDONNEES\n"
$c ,n,'"','"'
$c );
$c
$c printf("set Particules=`$WC $xTV/COORDONNEES`\n");
$c printf("set Particules=$Particules[$WC1]\n");
$c
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* VX=//%c
$c -e %cs/ .*//%c
$c > $xTV/COORDONNEES$COORD_X\n"
$c ,'"','"','"','"'
$c );
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* VY=//%c
$c -e %cs/ .*//%c
$c > $xTV/COORDONNEES$COORD_Y\n"
$c ,'"','"','"','"'
$c );
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* VZ=//%c
$c -e %cs/ .*//%c
$c > $xTV/COORDONNEES$COORD_Z\n"
$c ,'"','"','"','"'
$c );
$c
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* ROUGE=//%c
$c -e %cs/ .*//%c
$c > $xTV/COULEURS$ROUGE\n"
$c ,'"','"','"','"'
$c );
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* VERTE=//%c
$c -e %cs/ .*//%c
$c > $xTV/COULEURS$VERTE\n"
$c ,'"','"','"','"'
$c );
$c printf("$CA $xTV/COORDONNEES |
$c $SE -e %cs/^.* BLEUE=//%c
$c -e %cs/ .*//%c
$c > $xTV/COULEURS$BLEUE\n"
$c ,'"','"','"','"'
$c );
$c
$c printf("$xrv/particule.10$X
$c np=1
$c iterations=$Particules
$c isoles=VRAI joindre_central=$_____Jonction
$c jonction_a_un_corps=VRAI corps_de_jonction=$_____Reference
$c Apoints=$_____Points
$c LISTE_X=$xTV/COORDONNEES$COORD_X
$c LISTE_Y=$xTV/COORDONNEES$COORD_Y
$c LISTE_Z=$xTV/COORDONNEES$COORD_Z
$c LISTE_ROUGE=$xTV/COULEURS$ROUGE
$c LISTE_VERTE=$xTV/COULEURS$VERTE
$c LISTE_BLEUE=$xTV/COULEURS$BLEUE
$c Ar=8.0e-14
$c LISTE_RAYON=$xTV/RAYONS
$c Lz=100
$c editer_extrema_hors=FAUX
$c zoom_automatique=FAUX ZOOM=1.0
$c coordonnees_dans_0_1=VRAI
$c AXYZ=2.0e-5 BXYZ=0.5
$c TRANSLATION_OX=0.0 TRANSLATION_OY=0.0
$c fond=FAUX
$c centrer=FAUX reference=$_____Reference
$c attenuation_au_bord=0.25
$c $_____ParamParticule
$c R=$xTV/PART.1.
$c $formatI\n"
$c );
$c /* L'argument "editer_extrema_hors=FAUX" a ete introduit le 20200515112729... */
$c
$c printf("$xci/neutre$X
$c A=$xTV/PART.1.%04d$ROUGE
$c R=$_____imagesR2.%04d$ROUGE
$c $formatI\n"
$c ,N0,n
$c );
$c printf("$xci/neutre$X
$c A=$xTV/PART.1.%04d$VERTE
$c R=$_____imagesR2.%04d$VERTE
$c $formatI\n"
$c ,N0,n
$c );
$c printf("$xci/neutre$X
$c A=$xTV/PART.1.%04d$BLEUE
$c R=$_____imagesR2.%04d$BLEUE
$c $formatI\n"
$c ,N0,n
$c );
$c }
$c }