Parent Directory | Revision Log | Revision Graph
o first check-in for the "branch-genmake2" merge o verification suite as run on shelley (gcc 3.2.2): Wed Oct 8 23:42:29 EDT 2003 T S U V G D M c m s m s m s m s E p a R g m m e . m m e . m m e . m m e . N n k u 2 i a a d i a a d i a a d i a a d 2 d e n d n x n . n x n . n x n . n x n . OPTFILE=NONE Y Y Y Y 13 16 16 16 0 16 16 16 16 16 16 16 16 13 12 0 0 pass adjustment.128x64x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 0 0 16 16 0 0 pass adjustment.cs-32x32x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 22 0 16 16 22 0 pass adjust_nlfs.cs-32x32x1 Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O advect_cs Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O advect_xy Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O advect_xz Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_cs Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass aim.5l_Equatorial_Channel Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass aim.5l_LatLon Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass exp0 Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass exp1 Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass exp2 Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass exp4 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass global_ocean.90x40x15 Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL global_ocean.cs32x15 Y Y Y Y 6 11 12 13 13 12 13 16 13 9 9 9 9 10 9 9 11 FAIL global_ocean_pressure Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass global_with_exf Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16 9 16 pass hs94.128x64x5 Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass hs94.1x64x5 Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass hs94.cs-32x32x5 Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL ideal_2D_oce Y Y Y Y 8 16 16 16 16 16 16 16 16 13 13 8 16 16 16 16 16 FAIL internal_wave Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass inverted_barometer Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL lab_sea Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL natl_box Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass plume_on_slope Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass solid-body.cs-32x32x1
1 | edhill | 1.2 | C $Header: /u/u3/gcmpack/MITgcm/pkg/debug/debug_fld_stats_rl.F,v 1.1.12.1 2003/10/02 18:30:07 adcroft Exp $ |
2 | jmc | 1.1 | C $Name: $ |
3 | |||
4 | edhill | 1.2 | #include "DEBUG_OPTIONS.h" |
5 | jmc | 1.1 | |
6 | SUBROUTINE DEBUG_FLD_STATS_RL( | ||
7 | I myNr, arr, | ||
8 | O theMin,theMax,theMean,theSD, | ||
9 | I myThid ) | ||
10 | C /==========================================================\ | ||
11 | C | SUBROUTINE DEBUG_FLD_STATS_RL | | ||
12 | C | o Calculate bare statistics of global array "_RL arr" | | ||
13 | C |==========================================================| | ||
14 | C \==========================================================/ | ||
15 | IMPLICIT NONE | ||
16 | |||
17 | C === Global data === | ||
18 | #include "SIZE.h" | ||
19 | #include "EEPARAMS.h" | ||
20 | |||
21 | C === Routine arguments === | ||
22 | INTEGER myNr | ||
23 | _RL arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy) | ||
24 | _RL theMin | ||
25 | _RL theMax | ||
26 | _RL theMean | ||
27 | _RL theSD | ||
28 | INTEGER myThid | ||
29 | |||
30 | C === Local variables ==== | ||
31 | INTEGER bi,bj,I,J,K | ||
32 | INTEGER numPnts | ||
33 | LOGICAL noPnts | ||
34 | _RL tmpVal,rNumPnts | ||
35 | _RL theVar | ||
36 | |||
37 | theMin=0. | ||
38 | theMax=0. | ||
39 | theMean=0. | ||
40 | theSD=0. | ||
41 | theVar=0. | ||
42 | numPnts=0 | ||
43 | noPnts=.TRUE. | ||
44 | |||
45 | DO bj=myByLo(myThid),myByHi(myThid) | ||
46 | DO bi=myBxLo(myThid),myBxHi(myThid) | ||
47 | DO K=1,myNr | ||
48 | DO J=1,sNy | ||
49 | DO I=1,sNx | ||
50 | tmpVal=arr(I,J,K,bi,bj) | ||
51 | IF (tmpVal.NE.0. .AND. noPnts) THEN | ||
52 | theMin=tmpVal | ||
53 | theMax=tmpVal | ||
54 | noPnts=.FALSE. | ||
55 | ENDIF | ||
56 | IF (tmpVal.NE.0.) THEN | ||
57 | theMin=min(theMin,tmpVal) | ||
58 | theMax=max(theMax,tmpVal) | ||
59 | theMean=theMean+tmpVal | ||
60 | theVar=theVar+tmpVal**2 | ||
61 | numPnts=numPnts+1 | ||
62 | ENDIF | ||
63 | ENDDO | ||
64 | ENDDO | ||
65 | ENDDO | ||
66 | ENDDO | ||
67 | ENDDO | ||
68 | |||
69 | theMin=-theMin | ||
70 | _GLOBAL_MAX_R8(theMin,myThid) | ||
71 | theMin=-theMin | ||
72 | _GLOBAL_MAX_R8(theMax,myThid) | ||
73 | _GLOBAL_SUM_R8(theMean,myThid) | ||
74 | _GLOBAL_SUM_R8(theVar,myThid) | ||
75 | tmpVal=FLOAT(numPnts) | ||
76 | _GLOBAL_SUM_R8(tmpVal,myThid) | ||
77 | numPnts=INT(tmpVal+0.5) | ||
78 | |||
79 | IF (tmpVal.GT.0.) THEN | ||
80 | rNumPnts=1./tmpVal | ||
81 | theMean=theMean*rNumPnts | ||
82 | theVar=theVar*rNumPnts | ||
83 | |||
84 | DO bj=myByLo(myThid),myByHi(myThid) | ||
85 | DO bi=myBxLo(myThid),myBxHi(myThid) | ||
86 | DO K=1,myNr | ||
87 | DO J=1,sNy | ||
88 | DO I=1,sNx | ||
89 | tmpVal=arr(I,J,K,bi,bj) | ||
90 | IF (tmpVal.NE.0.) THEN | ||
91 | theSD=theSD+(tmpVal-theMean)**2 | ||
92 | ENDIF | ||
93 | ENDDO | ||
94 | ENDDO | ||
95 | ENDDO | ||
96 | ENDDO | ||
97 | ENDDO | ||
98 | |||
99 | _GLOBAL_SUM_R8(theSD,myThid) | ||
100 | |||
101 | theSD=sqrt(theSD*rNumPnts) | ||
102 | c theSD=sqrt(theVar-theMean**2) | ||
103 | ENDIF | ||
104 | |||
105 | RETURN | ||
106 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |