28 |
|
|
29 |
! Local: |
! Local: |
30 |
DOUBLE PRECISION Y |
DOUBLE PRECISION Y |
31 |
|
DOUBLE PRECISION h ! step of the uniform grid |
32 |
integer i, it |
integer i, it |
33 |
|
|
34 |
DOUBLE PRECISION xvrai(iim), Gvrai(iim) |
DOUBLE PRECISION xvrai(iim), Gvrai(iim) |
37 |
!------------------------------------------------------------------ |
!------------------------------------------------------------------ |
38 |
|
|
39 |
it = 0 ! initial guess |
it = 0 ! initial guess |
40 |
|
h = twopi_d / iim |
41 |
|
|
42 |
DO i = 1, iim |
DO i = 1, iim |
43 |
Y = - pi_d + (i + xuv - 0.75d0) * twopi_d / iim |
Y = - pi_d + (i + xuv - 0.75d0) * h |
44 |
! - pi <= y < pi |
! - pi <= y < pi |
45 |
abs_y = abs(y) |
abs_y = abs(y) |
46 |
|
|
63 |
END DO |
END DO |
64 |
|
|
65 |
xlon = xvrai + clon |
xlon = xvrai + clon |
66 |
xprim = twopi_d / (iim * Gvrai) |
xprim = h / Gvrai |
67 |
|
|
68 |
end subroutine invert_zoom_x |
end subroutine invert_zoom_x |
69 |
|
|