FUNCTION TABA (NPC,SF,T,PA) TABS 915 DIMENSION PA(2,NPC) TABS 916 C TABS 917 C SPECTRUM INTERPOLATION TABS 918 C TABS 919 DO 100 I=2,NPC TABS 920 T1=PA(1,I-1) TABS 921 T2=PA(1,I ) TABS 922 IF(T.LE.T2) GO TO 200 TABS 923 100 CONTINUE TABS 924 200 R1=(T2-T)/(T2-T1) TABS 925 R2=(T-T1)/(T2-T1) TABS 926 TABA=SF*(PA(2,I-1)*R1+PA(2,I)*R2) TABS 927 RETURN TABS 928 END

In modern Fortran I is equal to five as NPC is 4, any one got any ideas what value I would have had in Fortran 66 at the end of a do loop?

I will need to look at the original formula as well.