natl_box_forward/output/comp_orig.F

C     Compare KPPhbl_01_01.00000_00720_240 to KPPhbl_orig.
C     for North Atlantic test domain.
C     Maximum difference between two files is reported.
C     On SGI compile using "f77 -bytereclen comp_jpl.f"

      program comp_orig

      implicit none

      integer nx, ny, nrec
      parameter (nx=20, ny=16, nrec=4)

      integer lbuffer1, lbuffer2, i, j, k, irec
      character fni*80, fno*80
      real*4 fld1(nx,ny), fld2(nx,ny), maxdiff

      fni='KPPhbl.001.001.data'
      fno='KPPhbl_orig'
      lBuffer1=nx*ny*4
      lBuffer2=nx*ny*4
      open(1,file=fni,access='direct',recl=lBuffer1)
      open(2,file=fno,access='direct',recl=lBuffer2)

      irec=0
      maxdiff=0.
      do k=1,nrec
         irec=irec+1
         read(1,rec=irec) ((fld1(i,j),i=1,nx),j=1,ny)
         read(2,rec=irec) ((fld2(i,j),i=1,nx),j=1,ny)
         do i=1,nx
            do j=1,ny
               maxdiff=max(maxdiff,abs(fld1(i,j)-fld2(i,j)))
            enddo
         enddo
      enddo
      close(1)
      close(2)

      if(maxdiff.lt.5e-4) then
         print*,'North Atlantic test passed.  KPPhbl file is unchanged.'
      else
         print*,'WARNING! WARNING! WARNING! WARNING! WARNING! WARNING!'
         print*,'North Atlantic test failed.  KPPhbl file has changed.'
      endif
      print*,'Maximum difference is',maxdiff,' m.'

      stop
      end
1	heimbach	1.1.2.1	C Compare KPPhbl_01_01.00000_00720_240 to KPPhbl_orig.
2			C for North Atlantic test domain.
3			C Maximum difference between two files is reported.
4			C On SGI compile using "f77 -bytereclen comp_jpl.f"
5
6			program comp_orig
7
8			implicit none
9
10			integer nx, ny, nrec
11			parameter (nx=20, ny=16, nrec=4)
12
13			integer lbuffer1, lbuffer2, i, j, k, irec
14			character fni80, fno80
15			real*4 fld1(nx,ny), fld2(nx,ny), maxdiff
16
17			fni='KPPhbl.001.001.data'
18			fno='KPPhbl_orig'
19			lBuffer1=nxny4
20			lBuffer2=nxny4
21			open(1,file=fni,access='direct',recl=lBuffer1)
22			open(2,file=fno,access='direct',recl=lBuffer2)
23
24			irec=0
25			maxdiff=0.
26			do k=1,nrec
27			irec=irec+1
28			read(1,rec=irec) ((fld1(i,j),i=1,nx),j=1,ny)
29			read(2,rec=irec) ((fld2(i,j),i=1,nx),j=1,ny)
30			do i=1,nx
31			do j=1,ny
32			maxdiff=max(maxdiff,abs(fld1(i,j)-fld2(i,j)))
33			enddo
34			enddo
35			enddo
36			close(1)
37			close(2)
38
39			if(maxdiff.lt.5e-4) then
40			print*,'North Atlantic test passed. KPPhbl file is unchanged.'
41			else
42			print*,'WARNING! WARNING! WARNING! WARNING! WARNING! WARNING!'
43			print*,'North Atlantic test failed. KPPhbl file has changed.'
44			endif
45			print*,'Maximum difference is',maxdiff,' m.'
46
47			stop
48			end