122 |
INTEGER js |
INTEGER js |
123 |
#else |
#else |
124 |
INTEGER nLoop |
INTEGER nLoop |
125 |
|
_RL tmpVar |
126 |
#endif |
#endif |
127 |
#ifdef TARGET_NEC_SX |
#ifdef TARGET_NEC_SX |
128 |
INTEGER ic, ii, icnt |
INTEGER ic, ii, icnt |
360 |
ENDDO |
ENDDO |
361 |
C # of pts = nyIn+2 ; # of interval = nyIn+1 ; evaluate nLoop as |
C # of pts = nyIn+2 ; # of interval = nyIn+1 ; evaluate nLoop as |
362 |
C 1 + truncated log2(# interval -1); add epsil=1.e-3 for safey |
C 1 + truncated log2(# interval -1); add epsil=1.e-3 for safey |
363 |
nLoop = 1 + INT( LOG(DFLOAT(nyIn)+1. _d -3)/LOG(2. _d 0) ) |
tmpVar = nyIn + 1. _d -3 |
364 |
|
nLoop = 1 + INT( LOG(tmpVar)/LOG(2. _d 0) ) |
365 |
DO l=1,nLoop |
DO l=1,nLoop |
366 |
DO j=1,sNy |
DO j=1,sNy |
367 |
DO i=1,sNx |
DO i=1,sNx |