newbie question
newbie question
(OP)
can this be correct:
...
DO 2 J=1,80
IF(Y.GT.1.0) Y=0.6
F= -A*RP(I)+(Y+Y*Y+Y**3-Y**4)/(1.-Y)**3-B*Y*Y+C*Y**D
IF(ABS(F)-1.E-6) 4,4,3
3 DFDY=(1.+4.*Y+4.*Y*Y-4.*Y**3+Y**4)/(1.-Y)**4-2.*B*Y+D*C*Y**(D-1.)
2 Y= Y-F/DFDY
...
(do statement that invokes line no. 2)
...
DO 2 J=1,80
IF(Y.GT.1.0) Y=0.6
F= -A*RP(I)+(Y+Y*Y+Y**3-Y**4)/(1.-Y)**3-B*Y*Y+C*Y**D
IF(ABS(F)-1.E-6) 4,4,3
3 DFDY=(1.+4.*Y+4.*Y*Y-4.*Y**3+Y**4)/(1.-Y)**4-2.*B*Y+D*C*Y**(D-1.)
2 Y= Y-F/DFDY
...
(do statement that invokes line no. 2)
RE: newbie question
You could change it such that :
2 CONTINUE is the last statement if you wish.
Many would argue that using CONTINUE is cleaner.
If you are using Fortran 90 or 95 ( rather than the older Fortran 77 ), you can always end your loops ( if you wish ) with :
END DO
Most Fortran 77 compilers also accept END DO as a non-standard extension to their compilers implementation of Fortran. Note that a Fortran 90/95 compiler should have no problems compiling standard conforming Fortran 77 codes ( with a few very minor exceptions that is ).
Dan
www.dtware.com