- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/ecsort_shared.h
r3908 r5158 421 421 if (imod == 0) iadd = 0 422 422 ista(1) = 1 423 doj=2,inumt423 DO j=2,inumt 424 424 ista(j) = ista(j-1) + inc + iadd 425 425 if (iadd > 0 .and. j > imod) iadd = 0 426 426 enddo 427 427 ista(inumt+1) = n + 1 428 doj=1,inumt428 DO j=1,inumt 429 429 ichunk(j) = ista(j+1) - ista(j) 430 430 enddo … … 452 452 ITID = OML_MY_THREAD() 453 453 !$OMP DO SCHEDULE(DYNAMIC,1) 454 doj=1,inumt454 DO j=1,inumt 455 455 j1 = ista(j) 456 456 inum = ichunk(j) … … 490 490 inc = 2 491 491 imax = (inumt+inc-1)/inc 492 dojj=1,imax492 DO jj=1,imax 493 493 if (LLdebug) write(0,1001) jj,'<before_merge> jj,inc,imax,inumt=',jj,inc,imax,inumt 494 494 !$OMP PARALLEL PRIVATE(j,j1,j2,inum,iamax,ibmax,jmid,iret,ZHOOK_SUBHANDLE3,ITID) … … 496 496 ITID = OML_MY_THREAD() 497 497 !$OMP DO SCHEDULE(DYNAMIC,1) 498 doj=1,inumt,inc498 DO j=1,inumt,inc 499 499 j1 = j 500 500 j2 = j + inc - 1 … … 568 568 allocate(data(n)) 569 569 !$OMP DO SCHEDULE(DYNAMIC,1) 570 doj=1,sda570 DO j=1,sda 571 571 data(1:n) = a(iindex(1:n),j) 572 572 a(1:n,j) = data(1:n) … … 580 580 allocate(data(n)) 581 581 !$OMP DO SCHEDULE(DYNAMIC,1) 582 doi=1,lda582 DO i=1,lda 583 583 data(1:n) = a(i,iindex(1:n)) 584 584 a(i,1:n) = data(1:n) … … 613 613 WRITE(0,1000) jj,cdstr//': n,key,k1,k2,kadd,a(index(k1:k2)+kadd,:)=',& 614 614 & n,key,k1,k2,kadd 615 doj=k1,k2615 DO j=k1,k2 616 616 i = index(j)+kadd 617 617 WRITE(0,'(2i6)',advance='no') j,i-kadd … … 649 649 i = left 650 650 j = 2*left 651 dowhile (j <= right)651 DO while (j <= right) 652 652 if (j < right) then 653 653 if (a(index(j)) < a(index(j+1))) j = j + 1
Note: See TracChangeset
for help on using the changeset viewer.