49 |
c |
c |
50 |
|
|
51 |
REAL dzoom |
REAL dzoom |
52 |
REAL*8 ylat(jjp1), yprim(jjp1) |
DOUBLE PRECISION ylat(jjp1), yprim(jjp1) |
53 |
REAL*8 yuv |
DOUBLE PRECISION yuv |
54 |
REAL*8 yt(0:nmax2) |
DOUBLE PRECISION yt(0:nmax2) |
55 |
REAL*8 fhyp(0:nmax2),beta,Ytprim(0:nmax2),fxm(0:nmax2) |
DOUBLE PRECISION fhyp(0:nmax2),beta,Ytprim(0:nmax2),fxm(0:nmax2) |
56 |
SAVE Ytprim, yt,Yf |
SAVE Ytprim, yt,Yf |
57 |
REAL*8 Yf(0:nmax2),yypr(0:nmax2) |
DOUBLE PRECISION Yf(0:nmax2),yypr(0:nmax2) |
58 |
REAL*8 yvrai(jjp1), yprimm(jjp1),ylatt(jjp1) |
DOUBLE PRECISION yvrai(jjp1), yprimm(jjp1),ylatt(jjp1) |
59 |
REAL*8 pi,depi,pis2,epsilon,y0,pisjm |
DOUBLE PRECISION pi,depi,pis2,epsilon,y0,pisjm |
60 |
REAL*8 yo1,yi,ylon2,ymoy,Yprimin,champmin,champmax |
DOUBLE PRECISION yo1,yi,ylon2,ymoy,Yprimin,champmin,champmax |
61 |
REAL*8 yfi,Yf1,ffdy |
DOUBLE PRECISION yfi,Yf1,ffdy |
62 |
REAL*8 ypn,deply,y00 |
DOUBLE PRECISION ypn,deply,y00 |
63 |
SAVE y00, deply |
SAVE y00, deply |
64 |
|
|
65 |
INTEGER i,j,it,ik,iter,jlat |
INTEGER i,j,it,ik,iter,jlat |
66 |
INTEGER jpn,jjpn |
INTEGER jpn,jjpn |
67 |
SAVE jpn |
SAVE jpn |
68 |
REAL*8 a0,a1,a2,a3,yi2,heavyy0,heavyy0m |
DOUBLE PRECISION a0,a1,a2,a3,yi2,heavyy0,heavyy0m |
69 |
REAL*8 fa(0:nmax2),fb(0:nmax2) |
DOUBLE PRECISION fa(0:nmax2),fb(0:nmax2) |
70 |
REAL y0min,y0max |
REAL y0min,y0max |
71 |
|
|
72 |
REAL*8 heavyside |
DOUBLE PRECISION heavyside |
73 |
|
|
74 |
pi = 2. * ASIN(1.) |
pi = 2. * ASIN(1.) |
75 |
depi = 2. * pi |
depi = 2. * pi |