#######################################################################################################################################
#                                                                                                                                     #
#         D I S P O S I T I O N   " H A R M O N I E U S E "   D E   P O I N T S   S U R   U N E   S P H E R E                         #
#         A V E C   V I S U A L I S A T I O N   D E   L A   S P H E R E                                                               #
#         E T   V I S U A L I S A T I O N   D E   L A   D Y N A M I Q U E   D U   P R O C E S S U S  :                                #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirs/.PSPH.31.2.1.$U' :                                                                                        #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20030401114526).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$c                  #include  <stdio.h>
$c
$c                  #define   RAYON     _____Rayon
$c
$c                  #define   PI        3.141592653589793
$c
$c                  #define   PARITE    4
$c
$c                  #define   Ntheta    (16*PARITE)
$c                  #define   theta0    0.0
$c                  #define   thetaN    (2*PI)
$c
$c                  #define   Nphi      (8*PARITE)
$c                  #define   phi0      (-PI/2)
$c                  #define   phiN      (+PI/2)
$c
$c                  extern    double    cos();
$c                  extern    double    sin();
$c
$c                  int       main()
$c                            {
$c                            double    theta,phi;
$c
$c                            for       (theta=theta0 ;
$c                                       theta<=thetaN ;
$c                                       theta=theta+((thetaN-theta0)/Ntheta)
$c                                       )
$c                                      {
$c                                      for       (phi=phi0 ;
$c                                                 phi<=phiN ;
$c                                                 phi=phi+((phiN-phi0)/Nphi)
$c                                                 )
$c                                                {
$c                                                int       n;
$c                                                for       (n=1 ; n<=2 ; n++)
$c                                      /* En effet, il faut "faire deux tours" a cause de l'option "isoles=FAUX" ci-apres dans      */
$c                                      /* '$xrv/particule.10$X'...                                                                  */
$c                                                          {
$c                                                          printf("echo %f >>! $xTV/COORDONNEES_02$COORD_X\n"
$c                                                                ,RAYON*cos(theta)*cos(phi)
$c                                                                 );
$c                                                          printf("echo %f >>! $xTV/COORDONNEES_02$COORD_Y\n"
$c                                                                ,RAYON*sin(phi)
$c                                                                 );
$c                                                          printf("echo %f >>! $xTV/COORDONNEES_02$COORD_Z\n"
$c                                                                ,RAYON*sin(theta)*cos(phi)
$c                                                                 );
$c
$c                                                          printf("echo $_____RayonS >>! $xTV/RAYONS_02\n");
$c
$c                                                          printf("echo $_____CouleurSR >>! $xTV/COULEURS_02$ROUGE\n");
$c                                                          printf("echo $_____CouleurSV >>! $xTV/COULEURS_02$VERTE\n");
$c                                                          printf("echo $_____CouleurSB >>! $xTV/COULEURS_02$BLEUE\n");
$c                                                          }
$c                                                }
$c                                      }
$c                            }