파이프의 직경을 결정하는 문제입니다.

[신용욱]

조선소 기장설계팀에서 일하고 있는 직딩입니다.
현장에서 자주 설계문제에 대한 검증을 해야 하는 업무들을 접하게 되는데 아래와 같이 간단한 문제를 업무에 한번 적용해보고자 합니다. 그런데 포트란은 꽝이라... 고수님들 답변 바랍니다. 포멧 문에 문제가 있는거 같은데 이리 저리 고쳐 봐도 잘 안되네요..^^
내용을 간단히 설명하면 관내 유동에서 관의 직경을 결정하는 문제입니다. 아주 간단한 소스지만 현장에서 아주 유용하게 사용할 수 있을 거 같네요...^^

C EXAMPLE SOLUTION FOR CASE 3
C
REAL K,L,NU
WRITE (6,10)
G=9.80965
PI=3.14159
F=0.02
1 N=0
READ(5,11)I,HEAD,L,Q,NU,ESP
IF(I.EQ.999) GOTO 99
2 D=(F*L*Q*Q/(2.*G*HEAD*(PI/4.)**2))**.2
K=EPS/D
V=Q*4./(PI*D*D)
R=V*D/NU
IF((K.LE..00001).OR.(K.GE..04).OR.(R.LE.10000)) GOTO 5
A=.094*K**.225+.53*K
B=88.*K**.44
C=1.62*K**.134
F1=A+B/R**C
IF(ABS(F1-F).LE..0005) GOTO 3
N=N+1
IF(N.GE.8) GOTO 4
F=FI
GOTO 2
3 WRITE(6,12)L,Q,NU,HEAD,K,R,F,D
GOTO 1
4 WRITE(6,13)L,Q,MU,HEAD
GOTO 1
5 WRITE(6,14)L,Q,NU,HEAD,K,R
GOTO 1
10 FORMAT(1H1,///18X, 44HSOLUTION OF SIMPLE PIPE PROBLEM FOR
1 DIAMETER,
2 ///7X,1HL,8X,1HQ,9X,2HNU,8X,4HHEAD,7X,5HEPS/D,8X,HR 9X,
3 HF,8X,1HD,1)
11 FORMAT(13,F12.3,F15.2,F15.4,F15.9,F12.8)
12 FORMAT(F11.2,F9.4,F12.8,F9.3,F12.6,F12.1,F9.5,F9.5/)
13 FORMAT(11.2,F9.4,F12.8,F9.3,10X,24HANSWER DOES NOT CONVERG,/)
14 FORMAT(F11.2,F9.4,F12.8,F9.3,F12.6,F12.1,2X,16HINV ALID
1 SOLUTION,/)
99 STOP
END