#######################################################################################################################################
#                                                                                                                                     #
#         M O D U L E  :                                                                                                              #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xiirv/.IXEO.7.23.$U' :                                                                                          #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

$c                  #include  <stdio.h>
$c
$c                  extern    int       atoi();
$c                  extern    double    atof();
$c                  extern    char      *getenv();
$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                  #define   RAYONb    (atof(getenv("RayonB_Sphere")))
$c                  #define   RAYONs    (atof(getenv("RayonS_Sphere")))
$c
$c                  #define   dROUGE    (atoi(getenv("R_Sphere_D")))
$c                  #define   dVERTE    (atoi(getenv("V_Sphere_D")))
$c                  #define   dBLEUE    (atoi(getenv("B_Sphere_D")))
$c
$c                  #define   aROUGE    (atoi(getenv("R_Sphere_A")))
$c                  #define   aVERTE    (atoi(getenv("V_Sphere_A")))
$c                  #define   aBLEUE    (atoi(getenv("B_Sphere_A")))
$c
$c                  #define   N0        1
$c
$c                  extern    double    cos();
$c                  extern    double    sin();
$c
$c                  int       main()
$c                            {
$c                            int       n=N0;
$c                            int       N=((Ntheta+1)*(Nphi+1));
$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                                                double    X=RAYONs*cos(theta)*cos(phi),Y=RAYONs*sin(phi),Z=RAYONs*sin(theta)*cos(phi);
$c                                                double    dX=0,dY=0,dZ=0;
$c
$c                                                int       ROUGE=((aROUGE*(n-N0))+(dROUGE*(N-n)))/(N-N0);
$c                                                int       VERTE=((aVERTE*(n-N0))+(dVERTE*(N-n)))/(N-N0);
$c                                                int       BLEUE=((aBLEUE*(n-N0))+(dBLEUE*(N-n)))/(N-N0);
$c
$c                                                printf("\necho ' X=%f Y=%f Z=%f"
$c                                                      ,X,Y,Z
$c                                                       );
$c                                                printf(" dX=%f dY=%f dZ=%f"
$c                                                      ,dX,dY,dZ
$c                                                       );
$c                                                printf(" R=%f",RAYONb);
$c                                                printf(" ROUGE=%d VERTE=%d BLEUE=%d '"
$c                                                      ,ROUGE,VERTE,BLEUE
$c                                                       );
$c                                                printf(" >>! $xTV/lX2000");
$c                                                printf("\necho ' X=%f Y=%f Z=%f"
$c                                                      ,X,Y,Z
$c                                                       );
$c                                                printf(" dX=%f dY=%f dZ=%f"
$c                                                      ,dX,dY,dZ
$c                                                       );
$c                                                printf(" R=%f",RAYONb);
$c                                                printf(" ROUGE=%d VERTE=%d BLEUE=%d '"
$c                                                      ,ROUGE,VERTE,BLEUE
$c                                                       );
$c                                                printf(" >>! $xTV/lX2000");
$c
$c                                                n++;
$c                                                }
$c                                      }
$c
$c                            printf("\n");
$c                            }