/[lmdze]/trunk/Sources/dyn3d/fyhyp.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/fyhyp.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 148 by guez, Tue May 26 17:46:03 2015 UTC revision 149 by guez, Thu Jun 18 12:23:44 2015 UTC
# Line 15  contains Line 15  contains
15    
16      ! Il vaut mieux avoir : grossismy * dzoom < pi / 2      ! Il vaut mieux avoir : grossismy * dzoom < pi / 2
17    
18      use coefpoly_m, only: coefpoly      use coefpoly_m, only: coefpoly, a0, a1, a2, a3
19      USE dimens_m, only: jjm      USE dimens_m, only: jjm
20      use dynetat0_m, only: clat, grossismy, dzoomy, tauy      use dynetat0_m, only: clat, grossismy, dzoomy, tauy
21      use heavyside_m, only: heavyside      use heavyside_m, only: heavyside
# Line 47  contains Line 47  contains
47      INTEGER i, j, it, ik, iter, jlat      INTEGER i, j, it, ik, iter, jlat
48      INTEGER jpn, jjpn      INTEGER jpn, jjpn
49      SAVE jpn      SAVE jpn
50      DOUBLE PRECISION a0, a1, a2, a3, yi2, heavyy0, heavyy0m      DOUBLE PRECISION yi2, heavyy0, heavyy0m
51      DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)      DOUBLE PRECISION fa(0:nmax2), fb(0:nmax2)
52      REAL y0min, y0max      REAL y0min, y0max
53    
# Line 187  contains Line 187  contains
187            ! et Y'(yi)            ! et Y'(yi)
188    
189            CALL coefpoly(yf(it), yf(it + 1), ytprim(it), ytprim(it + 1), &            CALL coefpoly(yf(it), yf(it + 1), ytprim(it), ytprim(it + 1), &
190                 yt(it), yt(it + 1), a0, a1, a2, a3)                 yt(it), yt(it + 1))
191    
192            yf1 = yf(it)            yf1 = yf(it)
193            yprimin = a1 + 2.*a2*yi + 3.*a3*yi*yi            yprimin = a1 + 2.*a2*yi + 3.*a3*yi*yi
Legend:
Removed from v.148  
changed lines
  Added in v.149
  ViewVC Help
Powered by ViewVC 1.1.21