135 |
timeRec(1) = myTime-deltaTClock*0.5 _d 0 |
timeRec(1) = myTime-deltaTClock*0.5 _d 0 |
136 |
timeRec(1) = (timeRec(1)-phase(listId))/freq(listId) |
timeRec(1) = (timeRec(1)-phase(listId))/freq(listId) |
137 |
i = INT( timeRec(1) ) |
i = INT( timeRec(1) ) |
138 |
IF ( timeRec(1).LT.0. ) i = i - 1 |
IF ( timeRec(1).LT.0. ) THEN |
139 |
|
tmpLev = FLOAT(i) |
140 |
|
IF ( timeRec(1).NE.tmpLev ) i = i - 1 |
141 |
|
ENDIF |
142 |
timeRec(1) = phase(listId) + freq(listId)*FLOAT(i) |
timeRec(1) = phase(listId) + freq(listId)*FLOAT(i) |
143 |
c if ( listId.eq.2 ) write(0,*) 'f',i,timeRec(1)/deltaTClock |
c if ( listId.eq.2 ) write(0,*) 'f',i,timeRec(1)/deltaTClock |
144 |
timeRec(1) = MAX( timeRec(1), startTime ) |
timeRec(1) = MAX( timeRec(1), startTime ) |
152 |
timeRec(1) = baseTime + deltaTClock*FLOAT(i) |
timeRec(1) = baseTime + deltaTClock*FLOAT(i) |
153 |
c if ( listId.eq.2 ) write(0,*) i,timeRec(1)/deltaTClock |
c if ( listId.eq.2 ) write(0,*) i,timeRec(1)/deltaTClock |
154 |
ENDIF |
ENDIF |
155 |
|
C-- Convert time to iteration number (debug) |
156 |
|
c DO i=1,nTimRec |
157 |
|
c timeRec(i) = timeRec(i)/deltaTClock |
158 |
|
c ENDDO |
159 |
|
|
160 |
#ifdef ALLOW_MNC |
#ifdef ALLOW_MNC |
161 |
C-- this is a trick to reverse the order of the loops on md (= field) |
C-- this is a trick to reverse the order of the loops on md (= field) |