'fig6_PRL_2006' !The file names resulting from this computation will contain this chain
0.3     !longnorm: longueur des normales (0.d0 -> automatique pour ellipses)
300     !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 (1/2 largeur gaine, si ellipse ou polygone rectangulaire)
2.65  !b (1/2 hauteur gaine, si ellipse ou polygone rectangulaire) 
1.      !coef4 (si opts=4 ou 20, coef4>1 eloigne les sources de P, coef4<1 les rapproche)
200     !nfpoly: coefficients de Fourier pour lissage polygone
5.7,0.1,1 !re_eps1,inc,nbre calculs: partie r�elle permittivit� relative milieu exterieur
0.,0.1,1 !im_eps1,inc,nbre calculs: partie imaginaire permittivit� relative milieu exterieur
(3.4641,0.) !cnumater2 : indice gaine
2.02,0.,1 !wlvi,incwlv,nwlv: longueur d'onde ds le vide. ATTENTION: doit etre en microns si on fait varier l'indice en fct de wlv
0       !npol: POLARISATION, 0=E//,1=H//
1       !nchampinc: [0=onde plane] [1=fil source � l'ext�rieur de la gaine] [2=paquet gaussien]
30.     !tetad : angle incidence (degres), � partir de Oy, comme pour un r�seau
(0.,5.3) !zfilsou: affixe du fil source si nchampinc=1
5.      !dpaq: largeur du paquet
3.,0.   !x0paq,y0paq: centre du paquet (position du waist)
60      !nbop: nbre d'ondes planes dans le paquet
.false. !ldl: abscisse curviligne pour produits scalaires, voir plus bas
.false. !lsvd: (true:SVD) (false:Eq normales de Gauss)
180     !nntot: coeff de Fourier du champ diffract� de -nntot � + nntot
720     !npp: nbre de pts pour g(teta). Prendre >= 2*(nntot+1)
6     !nntig: troncature champ diffracte (et eventuellement interieur) pour chaque tige (de -nntig a +nntig)
1.d-4 !relconv: valeur conseill�e: 1.d-4. Affiche warnings si la convergence des bn pour la pire des tiges est inf�rieure � relconv
0.    !centx: abscisse du centre de la grille (cartes de champ)
0.    !centy: ordonnee du centre de la grille
30.   !lgrilx: cote horizontal de la grille 
30.   !lgrily: cote vertical de la grille    
0.2   !quad: cote d'une maille de la grille (si=0., ne fait pas la carte)
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:accumulation des resultats intermediaires
1       !nresu:Fichiers de resultats (0=ser&gtet seult)(1=tous)

---------------------------------------------------------------------------
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 
---------------------------------------------------------------------------
