'fig11_JOSAA_2004' !The file names resulting from this computation will contain this chain
0.3     !         longnorm: length of the normals, computed in prof_norm.XXX. Useless!
400     !Number of sources N refering to section 5 of G. Tayeb and S. Enoch, "Combined fictitious-sources�scattering-matrix method," J. Opt. Soc. Am. A 21, 1417-1423 (2004)
13.944  !a: half width (along x) of the jacket
2.6996  !b: half height (along y) of the jacket 
1.      !coef4: coef4>1 put the sources farther from the profile; coef4<1 put them closer. Begin with 1 if you do not know, then experiment...
200     !nfpoly: number of Fourier coefficients used to describe the profile, called n in Eq. (10) of G. Tayeb and S. Enoch, "Combined fictitious-sources�scattering-matrix method," J. Opt. Soc. Am. A 21, 1417-1423 (2004)
1.,0.1,1 !re_eps1,inc,nbre:      real part of the relative permittivity of the external medium (initial, step, number of computations)
0.,0.1,1 !im_eps1,inc,nbre: imaginary part of the relative permittivity of the external medium (initial, step, number of computations)
(3.4641,0.) !cnumater2 : complex optical index of the jacket [time dependency in exp(-iwt)]
2.02,0.,1 !wlvi,incwlv,nwlv: wavelength in vacuum (initial, step, number of computations)
0       !npol: POLARISATION, 0=E// to the rods, 1=H// to the rods
2       !nchampinc: [0=plane wave] [1=wire source outside the jacket] [2=Gaussian beam]
30.     !tetad : incidence angle (degres), taken from Oy (only used for plane wave or Gaussian beam)
(25.,0.) !zfilsou: position (x,y) of the wire source (only used if nchampinc=1)
5.      !dpaq: Gaussian beam width (only used for Gaussian beam)
3.,0.   !x0paq,y0paq: waist position of the Gaussian beam
60      !nbop: number of plane waves used to build the Gaussian beam
.false. !ldl: curvilinear abscissa for scaler products, see below
.false. !lsvd: (true:SVD) (false:Eq normales de Gauss)
180     !nntot: Fourier coefficients of the diffracted field: from -nntot to + nntot
720     !npp: nbre of pts for g(teta). Take >= 2*(nntot+1)
3     !nntig: truncation of the diffracted field (et eventuellement interieur) for each rod (from -nntig to +nntig)
1.d-4 !relconv: recommanded value: 1.d-4. Display warnings if the bn convergence for the worse rod is less more bad than relconv
0.    !centx: center of map abscissa (for field maps)
0.    !centy: center of map ordinate (for field maps)
30.   !lgrilx: horizontal side of the map 
30.   !lgrily:   vertical side of the map    
0.3   !quad: grid spacing (for field maps) (if=0., the map is not computed)
40    !Number of rods along x axis, on the first horizontal line of rods starting from the bottom
9     !Number of rows of rods, along the y axis
-2.355589  !y coordinate of the centers of the first line of rods
0.68  !Distance between centers of rods
(1.,0.) !Optical index of the rods
0.294 !Radius of the rods
.false. !lacfich: intermediate file results accumulation
1       !nresu: file results (0=ser&gtet only)(1=all)

---------------------------------------------------------------------------
Indices optiques:
Pour cnumater2 et les indices lus dans nomfipos:
  Si Im(cnumater2)=-1. , l'indice varie en fonction de wlv, et le mat�riau est choisi est fonction de Re(cnumater2) selon la correspondance suivante (� compl�ter dans subroutine calcindice au fur et � mesure des besoins):
    cnumater2=(1.,-1.) => Au
    cnumater2=(2.,-1.) => Ag
  Sinon, on lit directement l'indice, qui est alors ind�pendant de wlv
---------------------------------------------------------------------------
ATTENTION: ce programme a �t� con�u en pensant aux cas:
 nprofil = 3 ou 4, et devrait aussi fonctionner avec 1 et 2
 opts = 4
Pour les autres cas, il y a probablement des ajustements � faire
---------------------------------------------------------------------------
Differents profils pour la gaine:
nprofil=1: Ellipse, on peut choisir la taille (a et b) et toutes options opts ci-dessous.
        2: Haricot, profil d�fini dans subroutine haricot, opts= 4 ou 20 seulement.
        3: Polygone � angles arrondis. Profil d�fini dans subroutine polygone, opts= 4 ou 20 seulement. 
           Le cas actuellement programm� est un rectangle centr� � l'origine. Les sommets sont calcul�s � partir de a et b.
        4: Polygone quelconque � angles arrondis. Les sommets sont lus dans nomfipoly. opts= 4 ou 20 seulement (� v�rifier).
---------------------------------------------------------------------------
Differentes options de sources:
opts=0: comme MVDE; e1 et e2 sont alors les coeff d'homothetie (e1<1 et e2>1)
     1: selon teta; e1 et e2 sont alors des distances
     2: vers centre ellipse; e1 et e2 sont alors des distances
     3: sur normale; e1 et e2 sont alors des distances
     4: fonction du rayon de courbure (nd non utilise; nd=2*nb). 
        ATTENTION: visualiser la position des sources, car elles peuvent se placer du mauvais c�t� du profil.
    10: nappes
    20: pour un point du profil sur 2, on porte sur la normale 
         coef4 * la distance entre 2 points cons�cutifs du profil
        ATTENTION: visualiser la position des sources, car elles peuvent se placer du mauvais c�t� du profil.
---------------------------------------------------------------------------
nprosc et ldl:
nprosc : le produit scalaire de la colonne (u,v) est calcule en integrant sur P:
          |u|**2 avec ponderation cp1  (u=champ)
          |v|**2 avec ponderation cp3  (v=d�riv�e normale)
          si nprosc=0 , cp3=1
          si nprosc=5 , cp3 obtenu a partir de moyennes sur tous les vect. de base, et depend de la polarisation
          NOTA: nprosc n'est pas utilise pour les tiges IC
          NOTA: nprosc n'est pas utilise si la resolution se fait au moyen de SVD
                qui utilise son propre produit scalaire (equivalent a nprosc=0)
ldl    : T = on pondere les produits scalaires en tenant compte de la distance entre les pts sur P
         F = pas de ponderation 
---------------------------------------------------------------------------
