- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_interpolation.f90
r5119 r5158 34 34 jl = 0 35 35 ju = n + 1 36 dowhile (ju - jl > 1)36 DO while (ju - jl > 1) 37 37 jm = (ju + jl) / 2 ! Compute a midpoint, 38 38 IF (ascnd .eqv. (x >= xx(jm))) THEN … … 120 120 end if ! Done hunting, value bracketed. 121 121 122 do! Hunt is done, so begin the final bisection phase:122 DO ! Hunt is done, so begin the final bisection phase: 123 123 IF (jhi - jlo <= 1) THEN 124 124 IF (x == xx(n)) jlo = n - 1
Note: See TracChangeset
for help on using the changeset viewer.