93 |
C assume Interface at middle between 2 Center |
C assume Interface at middle between 2 Center |
94 |
drF(1) = delRc(1) |
drF(1) = delRc(1) |
95 |
DO k=2,Nr |
DO k=2,Nr |
96 |
drF(k-1) = 0.5 _d 0 *delRc(k) + drF(k-1) |
c drF(k-1) = 0.5 _d 0 *delRc(k) + drF(k-1) |
97 |
|
c drF( k ) = 0.5 _d 0 *delRc(k) |
98 |
|
C- note: change the order to prevent some compilers to produce wrong code |
99 |
|
C when trying to optimise this loop : |
100 |
drF( k ) = 0.5 _d 0 *delRc(k) |
drF( k ) = 0.5 _d 0 *delRc(k) |
101 |
|
drF(k-1) = 0.5 _d 0 *delRc(k) + drF(k-1) |
102 |
ENDDO |
ENDDO |
103 |
drF(Nr) = delRc(Nr+1) + drF(Nr) |
drF(Nr) = delRc(Nr+1) + drF(Nr) |
104 |
ENDIF |
ENDIF |