Использование ran1 FORTRAN 77 подпрограмма в FORTRAN 90

Question

Я пытаюсь использовать ran1 из Numerical Recipes в моем коде FORTRAN 90. Я думаю, что общий способ состоит в том, чтобы скомпилировать старую подпрограмму отдельно, а затем использовать объектный файл. Но здесь я хочу знать, какие изменения необходимо использовать непосредственно в моем коде.

FUNCTION ran1(idum) 
INTEGER idum,IA,IM,IQ,IR,NTAB,NDIV 
REAL ran1,AM,EPS,RNMX 
PARAMETER (IA=16807,IM=2147483647,AM=1./IM,IQ=127773,IR=2836, 
! NTAB=32,NDIV=1+(IM-1)/NTAB,EPS=1.2e-7,RNMX=1.-EPS) 
! “Minimal” random number generator of Park and Miller with Bays-Durham shuffle and 
! added safeguards. Returns a uniform random deviate between 0.0 and 1.0 (exclusive of 
! the endpoint values). Call with idum a negative integer to initialize; thereafter, do not 
! alter idum between successive deviates in a sequence. RNMX should approximate the largest 
! floating value that is less than 1. 
INTEGER j,k,iv(NTAB),iy 
SAVE iv,iy 
DATA iv /NTAB*0/, iy /0/ 
iy = 0 
if (idum.le.0.or.iy.eq.0) then !Initialize. 
idum=max(-idum,1) 
! Be sure to prevent idum = 0. 
do 11 j=NTAB+8,1,-1 
! Load the shuffle table (after 8 warm-ups). 
k=idum/IQ 
idum=IA*(idum-k*IQ)-IR*k 
if (idum.lt.0) idum=idum+IM 
if (j.le.NTAB) iv(j)=idum! Compute idum=mod(IA*idum,IM) without overflows by 
enddo 11 
iy=iv(1) 
endif 
k=idum/IQ 
idum=IA*(idum-k*IQ)-IR*k 
! Compute idum=mod(IA*idum,IM) without overflows by 
if (idum.lt.0) idum=idum+IM  ! Schrage’s method. 
j=1+iy/NDIV 
iy=iv(j)    ! Output previously stored value and refill the shuffle table. 

iv(j)=idum 
ran1=min(AM*iy,RNMX) ! Because users don’t expect endpoint values. 
return 
END 

Answer 1

Ваш код неисправен. Похоже, вы скопировали его вручную из книги, но не совсем. Вторая проблема действительно присутствует даже в книге.

Во-первых, должно быть продолжение строки и не комментарий в заявлении параметров на второй и третьей линии

PARAMETER (IA=16807,IM=2147483647,AM=1./IM,IQ=127773,IR=2836, & 
    NTAB=32,NDIV=1+(IM-1)/NTAB,EPS=1.2e-7,RNMX=1.-EPS) 

(преобразован в свободной форме, можно найти в книге для оригинала)

Во-вторых, цикл представляет собой странное сочетание петли делать с цифровой меткой и сделать петлю с end do:

do 11 j=NTAB+8,1,-1 
... 
enddo 11 

должен быть

do j=NTAB+8,1,-1 
... 
enddo 

или

do 11 j=NTAB+8,1,-1 
... 
11 continue 

Там может быть больше проблем нет.