! { dg-do compile } ! { dg-options "-O2 -floop-nest-optimize" } SUBROUTINE STONG(IGAUSS) DIMENSION EXX(6) PARAMETER (MXSH=1000, MXGTOT=5000) COMMON /NSHEL / EX(MXGTOT),CS(MXGTOT),NSHELL 100 CONTINUE NSHELL = NSHELL+1 IF(NSHELL.GT.MXSH) THEN RETURN END IF DO 320 I = 1,IGAUSS K = K1+I-1 EX(K) = EXX(I)*SCALE 320 CONTINUE IF(TNORM.GT.TOLNRM) THEN STOP END IF DO 460 IG = K1,K2 CS(IG) = FACS*CS(IG) 460 CONTINUE GO TO 100 END