Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/misc
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/misc
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_formcoord.f90
r5117 r5158 32 32 else 33 33 dxmin = abs(x(2) - x(1)) 34 doi = 2, n - 134 DO i = 2, n - 1 35 35 dx = abs(x(i + 1) - x(i)) 36 36 IF (dx<dxmin) dxmin = dx -
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 -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_libmath.f90
r5122 r5158 14 14 ismax = 1 15 15 sxmax = sx(1) 16 doi = 1, n - 116 DO i = 1, n - 1 17 17 ix = ix + incx 18 18 IF(sx(ix)>sxmax) THEN -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_q_sat.f90
r5117 r5158 45 45 ! ------------------------------------------------------------------ 46 46 47 doip = 1, np47 DO ip = 1, np 48 48 49 49 ! WRITE(*,*)'kelvin,millibar=',kelvin,millibar -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_slopes.f90
r5119 r5158 40 40 !------------------------------------------------------------------------------- 41 41 ! Local: 42 INTEGER :: n, i , j, sta(2), sto(2)42 INTEGER :: n, i 43 43 REAL :: xc(SIZE(f,1)) ! (n) cell centers 44 44 REAL :: h(2:SIZE(f,1)-1), delta_xc(2:SIZE(f,1)-1) ! (2:n-1) 45 REAL :: fm, ff, fp , dx45 REAL :: fm, ff, fp 46 46 !------------------------------------------------------------------------------- 47 47 n=SIZE(f,ix) -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_ssum_scopy.f90
r5120 r5158 32 32 ssum = 0. 33 33 ix = 1 34 doi = 1, n34 DO i = 1, n 35 35 ssum = ssum + sx(ix) 36 36 ix = ix + incx -
LMDZ6/branches/Amaury_dev/libf/misc/lmdz_write_field.f90
r5119 r5158 32 32 33 33 GetFieldIndex = -1 34 doi = 1, NbField34 DO i = 1, NbField 35 35 IF (TrueName==FieldName(i)) THEN 36 36 GetFieldIndex = i … … 152 152 Nb = 0 153 153 Pos = 2 154 doi = 1, Dim(1)154 DO i = 1, Dim(1) 155 155 nb = nb + 1 156 156 … … 195 195 spacing = '(t2,"' // repeat('-', ColumnSize * NbCol) // '")' 196 196 197 doi = 1, Dim(2)197 DO i = 1, Dim(2) 198 198 nb = 0 199 199 Pos = 2 200 doj = 1, Dim(1)200 DO j = 1, Dim(1) 201 201 nb = nb + 1 202 202 … … 255 255 spacing = '(t2,"' // repeat('-', ColumnSize * NbCol) // '")' 256 256 257 doi = 1, Dim(3)258 259 doj = 1, Dim(2)257 DO i = 1, Dim(3) 258 259 DO j = 1, Dim(2) 260 260 nb = 0 261 261 Pos = 2 262 262 263 dok = 1, Dim(1)263 DO k = 1, Dim(1) 264 264 nb = nb + 1 265 265
Note: See TracChangeset
for help on using the changeset viewer.