 | | call therm (tk, rho, a, cjth, cjtt, cp2, cv, dpdr, dpdt, g, h, pmpa, s, u) error = max (error, abs (cp - cp2)) write (*, '(3f12.6,2f14.8)') pbar, tc, rho, cp, cp2 end do write (*, '(a)') ' ' write (*, '(a,g14.6)') ' Maximum absolute error was ', error return end subroutine test02 !******************************************************************************* ! |