Fortranで正規分布に従う疑似乱数を発生させるプログラムが実行できない

Question

前提

プログラミング初心者です。
ネット上に載っていた「正規分布に従う疑似乱数を発生させるプログラム」を自分の環境で実行したのですが、READ文で読み込む乱数の数を入力しても、プログラムの中身が実行されずに終了してしまいました。
このプログラム自体が間違っているのか、私の操作が違うのか、教えていただけないでしょうか？

実現したいこと

以下のCLTプログラムの実行結果を正しく表示させたい

該当のソースコード

implicit none
! CTL method
real   hist(-60:60)
real*8 x, y
real*8 Ntot, Nin
real   z
integer i, j

open(1,file="hist.csv")

Nin = 0.0
Ntot = 0.0
write(6,*) 'Enter Ntot :'
read(5,*) Ntot

do i=-60,60
  hist(i) = 0.
enddo

do i=1,int(Ntot)
 call random_number(y)
 do j=1,11
   call random_number(x)
   y = y + x
 enddo 
 y = y - 6d0
 hist( int(y/0.1) ) = hist( int(y/0.1) ) + 1.
 
enddo

do i=-60,60 ! Normalize 0.1 ha kizamihaba
  hist(i) = hist(i) / Ntot / 0.1
enddo

do i=-60,60
  write(1,*) i*0.1, ',',  hist(i)
enddo


end

試したこと

コンパイルエラーがないことを確認して実行し、READ文でNtotに値を入力したが、何も表示されなかった。

Comments

[melian]

open(1,file="hist.csv")

としているので計算結果は hist.csv に出力されているのではないでしょうか。

[y_waiwai]

質問は編集できます
無駄に質問建てないで追記修正しましょう