THESE - Université Ferhat Abbas de Sétif
THESE - Université Ferhat Abbas de Sétif
THESE - Université Ferhat Abbas de Sétif
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
FERM=-CDV*RAU*DEXP(-<br />
XL/2.0d0)*LAGO2(XL,LL)*DBSJ0(ANG)*CVAR<br />
RETURN<br />
end<br />
*calcul <strong>de</strong> la fonction a integrer en sin: avec signe moins<br />
FUNCTION FEIM(XX)<br />
IMPLICIT DOUBLE PRECISION (a-h,o-z)<br />
c REAL*8 DBSJ0,ER,LAGO2,FER<br />
DOUBLE PRECISION LAGO2<br />
EXTERNAL DBSJ0,LAGO2,RJB<br />
COMMON/PAM1/WL,Z,alambd,ft,DRL,LL<br />
COMMON/PAM3/R,RSZ,W0<br />
COMMON/LIST/MAXSUB<br />
COMMON/PHASE/GOUY<br />
RAU=XX<br />
c XL est la variable R normalise pour le poly <strong>de</strong> Laguerre<br />
XL=2.0d0*RAU*RAU/WL/WL<br />
pi=4.d0*datan(1.d0(<br />
C R est la variable radiale dans le plan z<br />
c RAU est la variable radiale dans le plan <strong>de</strong> l"EOD<br />
c calcul avec lentille<br />
AGRC=pi*RAU*RAU*(DRL+(1.0d0/Z)-(1.0d0/ft))/alambd<br />
c calcul sans lentille<br />
c AGRC=pi*RAU*RAU*(DRL+(1.0d0/Z))/alambd<br />
c SVAR=DSIN(GOUY*DFLOAT(2*LL-1)+PHI+AGRC)<br />
SVAR=DSIN(AGRC)<br />
!!!!IMPORTANT : l'introduction du déphasage <strong>de</strong> Gouy perturbe<br />
!!le calcul. En particulier on ne retrouve pas qu'une gaussienne qui<br />
cc traverse une lentille est une gaussienne transformée mais une<br />
cc gaussienne tout <strong>de</strong> même. De plus la phase <strong>de</strong> Gouy ne dépendant pas<br />
cc <strong>de</strong> la coordonnée radiale n'apporte rien au calcul <strong>de</strong> diffraction.<br />
CDV=(W0/WL)*2.0d0*pi/alambd/Z<br />
ANG=2.d0*pi*R*RAU/alambd/Z<br />
FEIM=-CDV*RAU*DEXP(-<br />
XL/2.0d0)*LAGO2(XL,LL)*DBSJ0(ANG)*SVAR<br />
RETURN<br />
end<br />
98