--- MITgcm/pkg/debug/chksum_tiled.F 2001/01/24 17:08:17 1.1 +++ MITgcm/pkg/debug/chksum_tiled.F 2001/02/02 21:36:29 1.2 @@ -0,0 +1,49 @@ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/debug/chksum_tiled.F,v 1.2 2001/02/02 21:36:29 adcroft Exp $ +C $Name: $ + +#include "CPP_OPTIONS.h" + + SUBROUTINE CHKSUM_TILED( str,fld,NN,bi,bj,myThid ) + IMPLICIT NONE +C Global/common +#include "SIZE.h" +#include "EEPARAMS.h" +C Routines arguments + CHARACTER*(*) str + INTEGER NN,bi,bj,myThid + _RL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,NN,nSx,nSy) +C Local + INTEGER I,J,K + INTEGER cnt + _RL Fmn,Fmin,Fmax,Frms + + Fmn=0. + Frms=0. + Fmin=fld(1,1,1,bi,bj) + Fmax=fld(1,1,1,bi,bj) + cnt=0 + DO K=1,NN + DO J=1,sNy + DO I=1,sNx + Fmin=min(Fmin,fld(I,J,K,bi,bj)) + Fmax=max(Fmax,fld(I,J,K,bi,bj)) + Fmn=Fmn+fld(I,J,K,bi,bj) + cnt=cnt+1 + ENDDO + ENDDO + ENDDO + Fmn=Fmn/float(cnt) + DO K=1,NN + DO J=1,sNy + DO I=1,sNx + Frms=Frms+(fld(I,J,K,bi,bj)-Fmn)**2 + ENDDO + ENDDO + ENDDO + Frms=sqrt(Frms/float(cnt)) + + write(0,'(a,1p4e12.5,1x,a)') + & 'CHKSUM_TILED: ',Fmin,Fmax,Fmn,Frms,str + + RETURN + END