#######################################################################################################################################
# #
# D I S C R I M I N A T I O N E N T R E L E S P A R T I C U L E S A B S O R B A B L E S #
# E T C E L L E S Q U I N E L E S O N T P A S #
# D A N S U N M I L I E U D E P R O P A G A T I O N C O R R E S P O N D A N T #
# A L A C A R T E D E D E N S I T E T R I D I M E N S I O N N E L L E #
# D E S B R O N C H I O L E S T E R M I N A L E S D E S P O U M O N S #
# A V E C A J U S T E M E N T A U T O M A T I Q U E : #
# #
# #
# Author of '$xiird/.ACIN.c.16.$U' : #
# #
# Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss). #
# #
#######################################################################################################################################
$Z setParam _____AbsorbableMin 1
$Z setParam _____AbsorbableMax 100
$Z # ATTENTION : il n'est pas possible d'utiliser ici '5' (et '1' pour '$_____Absorbable_s') #
$Z # car, en effet, les nombres entiers generes par 'GenR' sont alors mal equi-repartis, ce que #
$Z # l'on compense donc en "elargissant" l'intervalle de [1,5] a [1,100]... #
$Z # #
$Z # On notera le 20160506125321 que la valeur "100" a aussi ete choisie car, en effet, elle #
$Z # represente en fait "100%"... #
$Z GenR $xTV/ABSORBABLE $_____AbsorbableMin $_____AbsorbableMax g=5353 entiers=VRAI epsilon=0.5 compatibilite_20030926=VRAI
$Z # L'option "compatibilite_20030926=VRAI" a ete introduite le 20041001102308... #
$Z # #
$Z # On genere un fichier d'absorbabilite dans {$_____AbsorbableMin,$_____AbsorbableMax}. #
$Z setParam _____Absorbable_i $_____AbsorbableMin
$Z setParam _____Absorbable_s 20
$Z # ATTENTION : il n'est pas possible d'utiliser ici '1' (et '5' pour '$_____AbsorbableMax') #
$Z # car, en effet, les nombres entiers generes par 'GenR' sont alors mal equi-repartis, ce que #
$Z # l'on compense donc en "elargissant" l'intervalle de [1,5] a [1,100]... #
$Z # #
$Z # On notera le 20160506125321 que "20" represente en fait "20%" en rappelant que #
$Z # '$_____AbsorbableMax' vaut "100", soit "100%". Ainsi, il y aura donc 20% de particules #
$Z # absorbables, ce qui correspond a la proportion d'Oxygene dans l'air... #
$Z setParam _____Marqueur_Abs 1
$Z setParam _____Marqueur_NAb 2
$Z $DELETE $xTV/SELECTEUR_Abs
$Z $xrv/passe_ban.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z minimum=$_____Absorbable_i maximum=$_____Absorbable_s \
$Z hors_bande=VRAI \
$Z marquer=VRAI \
$Z valeur_inferieure=$_____Marqueur_NAb \
$Z valeur_intermediaire=$_____Marqueur_Abs \
$Z valeur_superieure=$_____Marqueur_NAb \
$Z fichier=$xTV/ABSORBABLE \
$Z formater=VRAI entiers=VRAI signe="$K_VIDE" \
$Z > $xTV/SELECTEUR_Abs
$Z # Les particules dans {$_____Absorbable_i,$_____Absorbable_s} sont "Abs" ("ABSorbables"), #
$Z # avec "i" pour "Inferieur" et "s" pour "Superieur". Les autres sont dites "NAb" (soit #
$Z # "Non ABsorbables"). #
$Z SETParam _____Absorbable_Racine $NEXIST
$Z if ($_____Absorbable_Racine == $EXIST) then
$Z # Cas ou seules des particules de type 'Abs' seront emises depuis la racine 'mmm' de l'arbre. #
$Z $DELETE $xTV/SELECTEUR_mmm
$Z
$Z if ($_____Dimension == 3) then
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$NEXIST/" \
$Z > $xTV/SELECTEUR_mmm
$Z else
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$NEXIST/" \
$Z > $xTV/SELECTEUR_mmm
$Z endif
$Z
$Z $DELETE $xTV/SELECTEUR_Abs_1
$Z
$Z if ($_____Dimension == 3) then
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$_____Marqueur_Abs/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$_____Marqueur_NAb/" \
$Z > $xTV/SELECTEUR_Abs_1
$Z # Ainsi les particules '$_____Marqueur_mmm' seront obligatoirement '$_____Marqueur_Abs'... #
$Z else
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$_____Marqueur_Abs/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$_____Marqueur_Abs/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$_____Marqueur_NAb/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$_____Marqueur_NAb/" \
$Z > $xTV/SELECTEUR_Abs_1
$Z # Ainsi les particules '$_____Marqueur_mmm' et '$_____Marqueur_mmM' seront obligatoirement #
$Z # '$_____Marqueur_Abs'... #
$Z endif
$Z
$Z $DELETE $xTV/SELECTEUR_Abs_2
$Z
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$xTV/SELECTEUR_Abs_1 \
$Z fichierF=$xTV/SELECTEUR_Abs \
$Z selection=$xTV/SELECTEUR_mmm \
$Z formater=VRAI entiers=VRAI signe="$K_VIDE" \
$Z > $xTV/SELECTEUR_Abs_2
$Z # Ainsi les particules '$_____Marqueur_mmm' seront obligatoirement '$_____Marqueur_Abs', #
$Z # en plus de celles qui etaient deja marquees '$_____Marqueur_Abs'... #
$Z
$Z $DELETE $xTV/SELECTEUR_mmm
$Z
$Z $DELETE $xTV/SELECTEUR_Abs
$Z
$Z $xrv/neutre$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier=$xTV/SELECTEUR_Abs_2 \
$Z formater=VRAI entiers=VRAI signe="$K_VIDE" \
$Z > $xTV/SELECTEUR_Abs
$Z
$Z $DELETE $xTV/SELECTEUR_Abs_2
$Z else
$Z endif
$Z SETParam _____Vitesse_Racine $NEXIST
$Z if ($_____Vitesse_Racine == $EXIST) then
$Z # Cas ou les particules emises depuis la racine devront avoir une direction de vitesse #
$Z # alignee sur l'axe {mmm,MMM}. #
$Z $DELETE $xTV/SELECTEUR_mmm
$Z
$Z if ($_____Dimension == 3) then
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$NEXIST/" \
$Z > $xTV/SELECTEUR_mmm
$Z else
$Z $CA $xTV/SELECTEUR_Zones \
$Z | $SE -e "s/^$_____Marqueur_ext"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmm"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_Mmm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMm"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mmM"'$'"/$EXIST/" \
$Z -e "s/^$_____Marqueur_MmM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_mMM"'$'"/$NEXIST/" \
$Z -e "s/^$_____Marqueur_MMM"'$'"/$NEXIST/" \
$Z > $xTV/SELECTEUR_mmm
$Z endif
$Z
$Z
$Z $DELETE $xTV/VITESSES.rho
$Z $DELETE $xTV/VITESSES.theta
$Z $DELETE $xTV/VITESSES.phi
$Z
$Z $xrv/XYZ_RPT.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_X=$xTV/VITESSES$COORD_X \
$Z LISTE_Y=$xTV/VITESSES$COORD_Y \
$Z LISTE_Z=$xTV/VITESSES$COORD_Z \
$Z pRHO=1 pTHETA=0 pPHI=0 \
$Z > $xTV/VITESSES.rho
$Z $xrv/XYZ_RPT.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_X=$xTV/VITESSES$COORD_X \
$Z LISTE_Y=$xTV/VITESSES$COORD_Y \
$Z LISTE_Z=$xTV/VITESSES$COORD_Z \
$Z pRHO=0 pTHETA=1 pPHI=0 \
$Z > $xTV/VITESSES.theta
$Z $xrv/XYZ_RPT.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_X=$xTV/VITESSES$COORD_X \
$Z LISTE_Y=$xTV/VITESSES$COORD_Y \
$Z LISTE_Z=$xTV/VITESSES$COORD_Z \
$Z pRHO=0 pTHETA=0 pPHI=1 \
$Z > $xTV/VITESSES.phi
$Z
$Z set Delta_mmm_MMM_X=`$xcg/SOUS.01$X a=$MaximumX b=$MinimumX`
$Z set Delta_mmm_MMM_Y=`$xcg/SOUS.01$X a=$MaximumY b=$MinimumY`
$Z
$Z if ($_____Dimension == 3) then
$Z set Delta_mmm_MMM_Z=`$xcg/SOUS.01$X a=$MaximumZ b=$MinimumZ`
$Z else
$Z set Delta_mmm_MMM_Z=0
$Z endif
$Z
$Z set Params="$K_VIDE"
$Z set Params="$Params"" X=$Delta_mmm_MMM_X eX=FAUX"
$Z set Params="$Params"" Y=$Delta_mmm_MMM_Y eY=FAUX"
$Z set Params="$Params"" Z=$Delta_mmm_MMM_Z eZ=FAUX"
$Z set Rho__mmm_MMM=`$xcg/XYZ_RPT$X $Params eRHO=VRAI eTHETA=FAUX ePHI=FAUX | $SE -e "s/^.*=//"`
$Z set Theta_mmm_MMM=`$xcg/XYZ_RPT$X $Params eRHO=FAUX eTHETA=VRAI ePHI=FAUX | $SE -e "s/^.*=//"`
$Z set Phi__mmm_MMM=`$xcg/XYZ_RPT$X $Params eRHO=FAUX eTHETA=FAUX ePHI=VRAI | $SE -e "s/^.*=//"`
$Z
$Z setParam _____FMVitesse 2.0
$Z # Facteur d'augmentation du module de la vitesse des particules naissant a la racine... #
$Z # ATTENTION : le 20020322121300, j'ai essaye la valeur 8 ; celle-ci a declenche le probleme #
$Z # 'v $xrk/rdn_walk.52$I PRINT_ERREUR..le.Gradient.est.nul'. Il est important donc de ne pas #
$Z # trop augmenter le module de la vitesse des particules... #
$Z
$Z $DELETE $xTV/VITESSES_1.rho
$Z
$Z set Vitesse_Racine=`$xcg/MUL2.01$X a=$_____FMVitesse b=$_____MVitesse`
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$Vitesse_Racine \
$Z fichierF=$xTV/VITESSES.rho \
$Z selection=$xTV/SELECTEUR_mmm \
$Z > $xTV/VITESSES_1.rho
$Z
$Z $DELETE $xTV/VITESSES.rho
$Z
$Z $xrv/neutre$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier=$xTV/VITESSES_1.rho \
$Z > $xTV/VITESSES.rho
$Z
$Z $DELETE $xTV/VITESSES_1.rho
$Z
$Z $DELETE $xTV/VITESSES_1.theta
$Z
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$Theta_mmm_MMM \
$Z fichierF=$xTV/VITESSES.theta \
$Z selection=$xTV/SELECTEUR_mmm \
$Z > $xTV/VITESSES_1.theta
$Z
$Z $DELETE $xTV/VITESSES.theta
$Z
$Z setParam _____RdnThetaMin -0.05
$Z
$Z setParam _____RdnThetaMax +0.05
$Z
$Z GenR $xTV/RANDOM.theta $_____RdnThetaMin $_____RdnThetaMax g=3535
$Z
$Z $DELETE $xTV/RANDOM_1.theta
$Z
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$xTV/RANDOM.theta \
$Z fichierF=0 \
$Z selection=$xTV/SELECTEUR_mmm \
$Z > $xTV/RANDOM_1.theta
$Z
$Z $DELETE $xTV/VITESSES.theta
$Z
$Z $xrv/ADD2.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier1=$xTV/VITESSES_1.theta \
$Z fichier2=$xTV/RANDOM_1.theta \
$Z > $xTV/VITESSES.theta
$Z
$Z $DELETE $xTV/VITESSES_1.theta
$Z $DELETE $xTV/RANDOM_1.theta
$Z
$Z $DELETE $xTV/VITESSES_1.phi
$Z
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$Phi__mmm_MMM \
$Z fichierF=$xTV/VITESSES.phi \
$Z selection=$xTV/SELECTEUR_mmm \
$Z > $xTV/VITESSES_1.phi
$Z
$Z $DELETE $xTV/VITESSES.phi
$Z
$Z setParam _____RdnPhiMin -0.05
$Z
$Z setParam _____RdnPhiMax +0.05
$Z
$Z GenR $xTV/RANDOM.phi $_____RdnPhiMin $_____RdnPhiMax g=3939
$Z
$Z $DELETE $xTV/RANDOM_1.phi
$Z
$Z $xrv/COND.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichierV=$xTV/RANDOM.phi \
$Z fichierF=0 \
$Z selection=$xTV/SELECTEUR_mmm \
$Z > $xTV/RANDOM_1.phi
$Z
$Z $DELETE $xTV/VITESSES.phi
$Z
$Z $xrv/ADD2.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier1=$xTV/VITESSES_1.phi \
$Z fichier2=$xTV/RANDOM_1.phi \
$Z > $xTV/VITESSES.phi
$Z
$Z $DELETE $xTV/VITESSES_1.phi
$Z $DELETE $xTV/RANDOM_1.phi
$Z
$Z $DELETE $xTV/VITESSES$COORD_X
$Z $DELETE $xTV/VITESSES$COORD_Y
$Z $DELETE $xTV/VITESSES$COORD_Z
$Z
$Z $xrv/RPT_XYZ.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_RHO=$xTV/VITESSES.rho \
$Z LISTE_THETA=$xTV/VITESSES.theta \
$Z LISTE_PHI=$xTV/VITESSES.phi \
$Z pX=1 pY=0 pZ=0 \
$Z > $xTV/VITESSES$COORD_X
$Z $xrv/RPT_XYZ.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_RHO=$xTV/VITESSES.rho \
$Z LISTE_THETA=$xTV/VITESSES.theta \
$Z LISTE_PHI=$xTV/VITESSES.phi \
$Z pX=0 pY=1 pZ=0 \
$Z > $xTV/VITESSES$COORD_Y
$Z $xrv/RPT_XYZ.01$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z LISTE_RHO=$xTV/VITESSES.rho \
$Z LISTE_THETA=$xTV/VITESSES.theta \
$Z LISTE_PHI=$xTV/VITESSES.phi \
$Z pX=0 pY=0 pZ=1 \
$Z > $xTV/VITESSES$COORD_Z
$Z
$Z $DELETE $xTV/VITESSES.rho
$Z $DELETE $xTV/VITESSES.theta
$Z $DELETE $xTV/VITESSES.phi
$Z
$Z $DELETE $xTV/SELECTEUR_mmm
$Z else
$Z endif
$Z setParam _____Immobilis_Abs +1
$Z setParam _____Rayon_Abs 1.0
$Z setParam _____CouleurR_Abs $GRIS_8
$Z setParam _____CouleurV_Abs $GRIS_7
$Z setParam _____CouleurB_Abs $GRIS_0
$Z setParam _____Immobilis_NAb -1000000000
$Z # Le 20020327183038, je suis passe de '10000000000' a '1000000000' a cause du test #
$Z # 'v $xrk/rdn_walk.52$K N_EST_PAS_ENTIER.ACCES_IMMOBILISABLES.corps' et ce afin donc de #
$Z # ne pas avoir de valeurs entieres superieures a 'INFINI = +2147483647' en multipliant #
$Z # le contenu de '$xTV/IMMOBILISABLE' par '$_____Immobilis_NAb', en particulier en presence #
$Z # de probabilites egale a 1.000 comme cela s'est vu dans 'v _____xivPdf_14_1/026189_026700'. #
$Z setParam _____Rayon_NAb 1.0
$Z setParam _____CouleurR_NAb $GRIS_4
$Z setParam _____CouleurV_NAb $GRIS_4
$Z setParam _____CouleurB_NAb $GRIS_4
$Z $DELETE $xTV/IMMOBILISABLE_1
$Z $CA $xTV/SELECTEUR_Abs \
$Z | $SE -e "s/^$_____Marqueur_Abs"'$'"/$_____Immobilis_Abs/" \
$Z -e "s/^$_____Marqueur_NAb"'$'"/$_____Immobilis_NAb/" \
$Z > $xTV/IMMOBILISABLE_1
$Z $DELETE $xTV/IMMOBILISABLE_2
$Z $xrv/MUL2.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier1=$xTV/IMMOBILISABLE_1 \
$Z fichier2=$xTV/IMMOBILISABLE \
$Z > $xTV/IMMOBILISABLE_2
$Z $DELETE $xTV/IMMOBILISABLE_1
$Z $DELETE $xTV/IMMOBILISABLE
$Z $xrv/neutre$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier=$xTV/IMMOBILISABLE_2 \
$Z > $xTV/IMMOBILISABLE
$Z $DELETE $xTV/IMMOBILISABLE_2
$Z $DELETE $xTV/RAYONS_1
$Z $CA $xTV/SELECTEUR_Abs \
$Z | $SE -e "s/^$_____Marqueur_Abs"'$'"/$_____Rayon_Abs/" \
$Z -e "s/^$_____Marqueur_NAb"'$'"/$_____Rayon_NAb/" \
$Z > $xTV/RAYONS_1
$Z $DELETE $xTV/RAYONS_2
$Z $xrv/MUL2.11$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier1=$xTV/RAYONS_1 \
$Z fichier2=$xTV/RAYONS \
$Z > $xTV/RAYONS_2
$Z $DELETE $xTV/RAYONS_1
$Z $DELETE $xTV/RAYONS
$Z $xrv/neutre$X \
$Z CAL1=FAUX \
$Z ne=$PartN \
$Z fichier=$xTV/RAYONS_2 \
$Z > $xTV/RAYONS
$Z $DELETE $xTV/RAYONS_2
$Z $DELETE $xTV/COULEURS$ROUGE
$Z $DELETE $xTV/COULEURS$VERTE
$Z $DELETE $xTV/COULEURS$BLEUE
$Z $CA $xTV/SELECTEUR_Abs \
$Z | $SE -e "s/^$_____Marqueur_Abs"'$'"/$_____CouleurR_Abs/" \
$Z -e "s/^$_____Marqueur_NAb"'$'"/$_____CouleurR_NAb/" \
$Z > $xTV/COULEURS$ROUGE
$Z $CA $xTV/SELECTEUR_Abs \
$Z | $SE -e "s/^$_____Marqueur_Abs"'$'"/$_____CouleurV_Abs/" \
$Z -e "s/^$_____Marqueur_NAb"'$'"/$_____CouleurV_NAb/" \
$Z > $xTV/COULEURS$VERTE
$Z $CA $xTV/SELECTEUR_Abs \
$Z | $SE -e "s/^$_____Marqueur_Abs"'$'"/$_____CouleurB_Abs/" \
$Z -e "s/^$_____Marqueur_NAb"'$'"/$_____CouleurB_NAb/" \
$Z > $xTV/COULEURS$BLEUE
$Z saut
$Z echo "nombre de particules de type 'Abs'=`$GRE '^$_____Marqueur_Abs'"'$'" $xTV/SELECTEUR_Abs | $WCl`"
$Z echo "nombre de particules de type 'NAb'=`$GRE '^$_____Marqueur_NAb'"'$'" $xTV/SELECTEUR_Abs | $WCl`"
$Z saut
$Z $DELETE $xTV/SELECTEUR_Abs