I have received some code that uses quadruple precision in some parts:
SUBROUTINE GAUSSQ(NS,C,B,NSD,A) IMPLICIT REAL*16 (A-H,O-Z) DIMENSION C(NS),B(NS),A(NSD,NS) IF (NS.EQ.2) THEN C(1)= .211324865405187117745425609749Q+00 C(2)= .788675134594812882254574390251Q+00 B(1) = 0.5Q+00 A(1,1)= .250000000000000000000000000000Q+00 A(1,2)= -.386751345948128822545743902510Q-01 A(2,1)= .538675134594812882254574390251Q+00 A(2,2)= .250000000000000000000000000000Q+00 END IF
and so on. However, I can’t compile this code with pgf77, pgf90, etc. I get an error
PGFTN-W-0031-Illegal data type length specifier for real
Is there a way to get this to compile and run?