- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/math_lib.F90
r5099 r5158 96 96 else 97 97 sumo = 0._wp 98 doj=i1,i298 DO j=i1,i2 99 99 deltah = abs(s(i1+1)-s(i1)) 100 100 sumo = sumo + f(j)*deltah … … 204 204 end if 205 205 206 doi = 2, ntab206 DO i = 2, ntab 207 207 if ( xtab(i) <= xtab(i-1) ) then 208 208 lerror = .true. … … 239 239 ihi = ntab 240 240 241 doi = 1, ntab241 DO i = 1, ntab 242 242 if ( a <= xtab(i) ) then 243 243 exit … … 249 249 ilo = min ( ilo, ntab - 1 ) 250 250 251 doi = 1, ntab251 DO i = 1, ntab 252 252 if ( xtab(i) <= b ) then 253 253 exit … … 263 263 !ds sum1 = 0.0D+00 264 264 265 doi = ilo, ihi265 DO i = ilo, ihi 266 266 267 267 x1 = xtab(i-1) … … 371 371 ga=1._wp 372 372 m1=x-1 373 dok=2,m1373 DO k=2,m1 374 374 ga=ga*k 375 375 enddo … … 382 382 m=int(z) 383 383 r=1._wp 384 dok=1,m384 DO k=1,m 385 385 r=r*(z-k) 386 386 enddo … … 390 390 endif 391 391 gr=g(26) 392 dok=25,1,-1392 DO k=25,1,-1 393 393 gr=gr*z+g(k) 394 394 enddo
Note: See TracChangeset
for help on using the changeset viewer.