Line data Source code
1 : MODULE m_fitchk 2 : CONTAINS 3 24 : SUBROUTINE fitchk(f1,f2,av,rms,dmx) 4 : ! ************************************************ 5 : ! compare functions f1 and f2 6 : ! ************************************************ 7 : IMPLICIT NONE 8 : REAL,INTENT(OUT):: av,dmx,rms 9 : REAL,INTENT(IN):: f1(:),f2(:) 10 : ! .. Local Scalars .. 11 : REAL d 12 : INTEGER i 13 : 14 24 : av = 0. 15 24 : rms = 0. 16 24 : dmx = 0. 17 8424 : DO i = 1,SIZE(f1) 18 8400 : av = av + f1(i) 19 8400 : d = (f1(i)-f2(i))**2 20 8400 : dmx = MAX(d,dmx) 21 8424 : rms = rms + d 22 : ENDDO 23 24 : av = av/size(f1) 24 24 : IF (abs(av).LT.1.e-30) THEN 25 0 : rms = 0. 26 0 : dmx = 0. 27 0 : RETURN 28 : END IF 29 24 : rms = sqrt(rms/size(f1))/av*100. 30 24 : dmx = sqrt(dmx)/av*100. 31 : END SUBROUTINE 32 : END