/[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.1.2.1 - (show annotations) (download)
Wed Jan 24 17:08:17 2001 UTC (23 years, 4 months ago) by adcroft
Branch: branch-atmos-merge
CVS Tags: branch-atmos-merge-zonalfilt
Changes since 1.1: +42 -0 lines
This routine is not called by default and not enabled by any CPP options.
It useful to have it around though.

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

  ViewVC Help
Powered by ViewVC 1.1.22