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 | |
2 | #include "DEBUG_OPTIONS.h" |
3 | |
4 | SUBROUTINE DEBUG_FLD_STATS_RS( |
5 | I myNr, arr, |
6 | O theMin,theMax,theMean,theSD, |
7 | I myThid ) |
8 | C /==========================================================\ |
9 | C | SUBROUTINE DEBUG_FLD_STATS_RS | |
10 | C | o Calculate bare statistics of global array "_RS arr" | |
11 | C |==========================================================| |
12 | C \==========================================================/ |
13 | IMPLICIT NONE |
14 | |
15 | C === Global data === |
16 | #include "SIZE.h" |
17 | #include "EEPARAMS.h" |
18 | |
19 | C === Routine arguments === |
20 | INTEGER myNr |
21 | _RS arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy) |
22 | _RL theMin |
23 | _RL theMax |
24 | _RL theMean |
25 | _RL theSD |
26 | INTEGER myThid |
27 | |
28 | C === Local variables ==== |
29 | INTEGER bi,bj,I,J,K |
30 | INTEGER numPnts |
31 | LOGICAL noPnts |
32 | _RL tmpVal,rNumPnts |
33 | _RL theVar |
34 | |
35 | theMin=0. |
36 | theMax=0. |
37 | theMean=0. |
38 | theSD=0. |
39 | theVar=0. |
40 | numPnts=0 |
41 | noPnts=.TRUE. |
42 | |
43 | DO bj=myByLo(myThid),myByHi(myThid) |
44 | DO bi=myBxLo(myThid),myBxHi(myThid) |
45 | DO K=1,myNr |
46 | DO J=1,sNy |
47 | DO I=1,sNx |
48 | tmpVal=arr(I,J,K,bi,bj) |
49 | IF (tmpVal.NE.0. .AND. noPnts) THEN |
50 | theMin=tmpVal |
51 | theMax=tmpVal |
52 | noPnts=.FALSE. |
53 | ENDIF |
54 | IF (tmpVal.NE.0.) THEN |
55 | theMin=min(theMin,tmpVal) |
56 | theMax=max(theMax,tmpVal) |
57 | theMean=theMean+tmpVal |
58 | theVar=theVar+tmpVal**2 |
59 | numPnts=numPnts+1 |
60 | ENDIF |
61 | ENDDO |
62 | ENDDO |
63 | ENDDO |
64 | ENDDO |
65 | ENDDO |
66 | |
67 | theMin=-theMin |
68 | _GLOBAL_MAX_R8(theMin,myThid) |
69 | theMin=-theMin |
70 | _GLOBAL_MAX_R8(theMax,myThid) |
71 | _GLOBAL_SUM_R8(theMean,myThid) |
72 | _GLOBAL_SUM_R8(theVar,myThid) |
73 | tmpVal=FLOAT(numPnts) |
74 | _GLOBAL_SUM_R8(tmpVal,myThid) |
75 | numPnts=INT(tmpVal+0.5) |
76 | |
77 | IF (tmpVal.GT.0.) THEN |
78 | rNumPnts=1./tmpVal |
79 | theMean=theMean*rNumPnts |
80 | theVar=theVar*rNumPnts |
81 | |
82 | DO bj=myByLo(myThid),myByHi(myThid) |
83 | DO bi=myBxLo(myThid),myBxHi(myThid) |
84 | DO K=1,myNr |
85 | DO J=1,sNy |
86 | DO I=1,sNx |
87 | tmpVal=arr(I,J,K,bi,bj) |
88 | IF (tmpVal.NE.0.) THEN |
89 | theSD=theSD+(tmpVal-theMean)**2 |
90 | ENDIF |
91 | ENDDO |
92 | ENDDO |
93 | ENDDO |
94 | ENDDO |
95 | ENDDO |
96 | |
97 | _GLOBAL_SUM_R8(theSD,myThid) |
98 | |
99 | theSD=sqrt(theSD*rNumPnts) |
100 | c theSD=sqrt(theVar-theMean**2) |
101 | ENDIF |
102 | |
103 | RETURN |
104 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |