/[MITgcm]/MITgcm/pkg/debug/chksum_tiled.F
ViewVC logotype

Contents of /MITgcm/pkg/debug/chksum_tiled.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3.20.1 - (show annotations) (download)
Thu Oct 2 18:30:07 2003 UTC (20 years, 8 months ago) by adcroft
Branch: branch-genmake2
Changes since 1.3: +3 -3 lines
Mis-type CONFIG as CONF !!!

1 C $Header: /u/gcmpack/MITgcm/pkg/debug/chksum_tiled.F,v 1.3 2001/04/10 22:35:27 heimbach Exp $
2 C $Name: $
3
4 #include "DEBUG_OPTIONS.h"
5
6 SUBROUTINE CHKSUM_TILED( str,fld,NN,bi,bj,myThid )
7 IMPLICIT NONE
8 C Global/common
9 #include "SIZE.h"
10 #include "EEPARAMS.h"
11 C Routines arguments
12 CHARACTER*(*) str
13 INTEGER NN,bi,bj,myThid
14 _RL fld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,NN,nSx,nSy)
15 C Local
16 INTEGER I,J,K
17 INTEGER cnt
18 _RL Fmn,Fmin,Fmax,Frms
19
20 Fmn=0.
21 Frms=0.
22 Fmin=fld(1,1,1,bi,bj)
23 Fmax=fld(1,1,1,bi,bj)
24 cnt=0
25 DO K=1,NN
26 DO J=1,sNy
27 DO I=1,sNx
28 Fmin=min(Fmin,fld(I,J,K,bi,bj))
29 Fmax=max(Fmax,fld(I,J,K,bi,bj))
30 Fmn=Fmn+fld(I,J,K,bi,bj)
31 cnt=cnt+1
32 ENDDO
33 ENDDO
34 ENDDO
35 Fmn=Fmn/float(cnt)
36 DO K=1,NN
37 DO J=1,sNy
38 DO I=1,sNx
39 Frms=Frms+(fld(I,J,K,bi,bj)-Fmn)**2
40 ENDDO
41 ENDDO
42 ENDDO
43 Frms=sqrt(Frms/float(cnt))
44
45 write(*,'(a,1p4e12.5,1x,a)')
46 & 'CHKSUM_TILED: ',Fmin,Fmax,Fmn,Frms,str
47
48 RETURN
49 END

  ViewVC Help
Powered by ViewVC 1.1.22