hdiff output

r31323/mc_gbh.F90 2016-10-16 19:30:12.308481578 +0100 r31322/mc_gbh.F90 2016-10-16 19:30:12.556484871 +0100
185:   USE COMMONS, ONLY : TEMP, QALCSV, MYUNIT185:   USE COMMONS, ONLY : TEMP, QALCSV, MYUNIT
186:   !186:   !
187:   IMPLICIT NONE187:   IMPLICIT NONE
188:   !188:   !
189:   INTEGER, INTENT(IN) :: N189:   INTEGER, INTENT(IN) :: N
190:   DOUBLE PRECISION, INTENT(IN) :: VALUES(N)190:   DOUBLE PRECISION, INTENT(IN) :: VALUES(N)
191:   INTEGER, INTENT(OUT) :: I191:   INTEGER, INTENT(OUT) :: I
192:   !192:   !
193:   DOUBLE PRECISION :: PSUM(N), X, Y, DPRAND, ELOWEST193:   DOUBLE PRECISION :: PSUM(N), X, Y, DPRAND, ELOWEST
194:   !194:   !
195:   !X=0.0D0 ! initialise total sum 
196:   !DO I=1,N-1 
197:   !   Y = MIN(DEXP(-(VALUES(I)-VALUES(N))/TEMP(1)), 1.0D0) 
198:   !   X = X + Y 
199:   !   PSUM(I) = X ! store partial sum 
200:   !ENDDO 
201:   !PSUM(N) = DBLE(N-1) 
202:   ! 
203:   ELOWEST=MINVAL(VALUES) 
204:   X=0.0D0 ! initialise total sum195:   X=0.0D0 ! initialise total sum
205:   DO I=1,N-1196:   DO I=1,N-1
206:      IF(VALUES(I) > VALUES(N) .AND. VALUES(N) < ELOWEST) THEN197:      Y = MIN(DEXP(-(VALUES(I)-VALUES(N))/TEMP(1)), 1.0D0)
207:         Y=0.0D0 
208:      ELSE 
209:         Y=DEXP(-(VALUES(I)-ELOWEST)/TEMP(1)) 
210:      ENDIF 
211:      X = X + Y198:      X = X + Y
212:      PSUM(I) = X ! store partial sum199:      PSUM(I) = X ! store partial sum
213:   ENDDO200:   ENDDO
214:   IF(VALUES(N) < ELOWEST) THEN201:   PSUM(N) = DBLE(N-1)
215:      PSUM(N) = PSUM(N-1) ! Disallow rejection 
216:   ELSE 
217:      PSUM(N) = DBLE(N-1) ! Allow rejection 
218:   ENDIF 
219:   !202:   !
220:   X=PSUM(N)*DPRAND()203:   X=PSUM(N)*DPRAND()
221:   !204:   !
222:   I=1205:   I=1
223:   DO WHILE (X > PSUM(I))206:   DO WHILE (X > PSUM(I))
224:      I=I+1207:      I=I+1
225:   ENDDO208:   ENDDO
226:   !209:   !
227:   !IF(QALCSV) THEN210:   !IF(QALCSV) THEN
228:   !   WRITE(MYUNIT, *) 'choose_from_list> vals=', VALUES211:   !   WRITE(MYUNIT, *) 'choose_from_list> vals=', VALUES


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0