#######################################################################################################################################
#                                                                                                                                     #
#         F L O T   D E   P A R T I C U L E S   L E   L O N G                                                                         #
#         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  :                                              #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiird/.ACIN.A.11.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$Z                  setParam            _____Palette                  $xiP/trou_noir.84

$Z                  SET       Premiere=1
$Z                  SET       Derniere=240
$Z                                      # '$Derniere' doit etre un multiple de 'N' defini ci-apres avec '3*(DIMENSION-1)=3*(9-1)=24'. #

$Z                  Pal

$Z                  setenv    Acinus_Ajuster      FAUX
$Z                  setenv    Acinus_Gros         1

$c                  #define   DIMENSION 9
$c                  #define   GRAINE    11
$c
$c                  #define   N0        1
$c                  #define   N         (3*(DIMENSION-1))
$c
$c                  #define   L0        0.0
$c                  #define   L         1.0
$c
$c                  int       main()
$c                            {
$c                            int       n=N0;
$c                            int       niveau;
$c                            int       points;
$c                            int       Npoints=(Derniere-Premiere+1)/(N-N0+1)+1;
$c                                      /* Il faut un point de plus (+1) que ce l'on veut visualiser car, en effet, les extremites   */
$c                                      /* des segments ne doivent etre tracees qu'une fois bien que presque toutes appartiennent a  */
$c                                      /* deux segments pusique :                                                                   */
$c                                      /*                                                                                           */
$c                                      /*                  extremite(segment(n-1)) = origine(segment(n))                            */
$c                                      /*                                                                                           */
$c                                      /* Mais malgre tout, il faut ce '+1' afin que la definition du 'lambda' ci-apres soit        */
$c                                      /* compatible avec celle de 'v $xrv/particule.31$I DEFV.Float.INIT.lambda'.                  */
$c
$c                            printf("setenv    Acinus_APoints      %d\n"
$c                                  ,Npoints
$c                                   );
$c
$c                            for       (niveau=N0 ; niveau<=N ; niveau++)
$c                                      {
$c                                      for       (points=1 ; points<Npoints ; points++)
$c                                                {
$c                                                double    lambda=((double)(points-1))/((double)(Npoints-1));
$c                                                double    Dlambda=1.0/((double)2*(Npoints-1));
$c
$c                                                printf("setenv    Acinus_LambdaMin    %f\n"
$c                                                      ,lambda-Dlambda
$c                                                       );
$c                                                printf("setenv    Acinus_LambdaMax    %f\n"
$c                                                      ,lambda+Dlambda
$c                                                       );
$c
$c                                                printf("$xu/Felici.M/GenImages.01$Z
$c                                                                    $_____images.%04d
$c                                                                    %d %d %d
$c                                                                    0 0 0
$c                                                                    $_____Palette
$c                                                                    %d
$c                                                                    %d %d\n"
$c                                                      ,n
$c                                                      ,DIMENSION,DIMENSION,DIMENSION
$c                                                      ,GRAINE
$c                                                      ,niveau,niveau
$c                                                       );
$c
$c                                                n++;
$c                                                }
$c                                      }
$c                            }