MITgcm_contrib/jmc_script/cmpnum.f

C $Header: $
C $Name:  $
      program cmpnum
      implicit none
      real*8 a,b,diff
      integer linnum,tmp,best
      open(17,file='cmpnum.log',status='unknown')
      best=-22
  99  read(*,*,end=70,err=60) linnum,a,b
      diff=0.5*(abs(a)+abs(b))
c     print *,a,b,diff,abs(a-b)/diff
c     if (diff.gt.1.e-12) then
      if (diff.gt.0.) then
        diff=abs(a-b)/diff
        if (diff.gt.0.) then
c         print *,int(log10(diff)),diff
          tmp=nint(log10(diff))
        else
          tmp=-16;
        endif
        best=max(best,tmp)
      else
        tmp =-22
      endif
      write(17,'(I4,1P2E22.14,A,I4)') linnum,a,b,' :',-tmp
      goto 99
  60  stop 'cmpnum: An error occured reading a,b'
  70  continue
      close(17)
      print *,-best
      end
1	C $Header: $
2	C $Name: $
3	program cmpnum
4	implicit none
5	real*8 a,b,diff
6	integer linnum,tmp,best
7	open(17,file='cmpnum.log',status='unknown')
8	best=-22
9	99 read(,,end=70,err=60) linnum,a,b
10	diff=0.5*(abs(a)+abs(b))
11	c print *,a,b,diff,abs(a-b)/diff
12	c if (diff.gt.1.e-12) then
13	if (diff.gt.0.) then
14	diff=abs(a-b)/diff
15	if (diff.gt.0.) then
16	c print *,int(log10(diff)),diff
17	tmp=nint(log10(diff))
18	else
19	tmp=-16;
20	endif
21	best=max(best,tmp)
22	else
23	tmp =-22
24	endif
25	write(17,'(I4,1P2E22.14,A,I4)') linnum,a,b,' :',-tmp
26	goto 99
27	60 stop 'cmpnum: An error occured reading a,b'
28	70 continue
29	close(17)
30	print *,-best
31	end