/[MITgcm]/MITgcm/pkg/diagnostics/diagstats_output.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagstats_output.F

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

revision 1.1 by jmc, Fri May 20 07:28:52 2005 UTC revision 1.2 by jmc, Sun Jun 26 16:51:49 2005 UTC
# Line 15  C     !INTERFACE: Line 15  C     !INTERFACE:
15    
16  C     !DESCRIPTION:  C     !DESCRIPTION:
17  C     Write output for diagnostics fields.  C     Write output for diagnostics fields.
18          
19  C     !USES:  C     !USES:
20        IMPLICIT NONE        IMPLICIT NONE
21  #include "SIZE.h"  #include "SIZE.h"
# Line 41  C     myThid  :: my Thread Id number Line 41  C     myThid  :: my Thread Id number
41  CEOP  CEOP
42    
43  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
44        INTEGER j, m, ndId        INTEGER j, m, ndId, iSp, iSm
45        CHARACTER*8 parms1        CHARACTER*8 parms1
46        CHARACTER*3 mate_index        CHARACTER*3 mate_index
47        INTEGER mate        INTEGER mate
# Line 63  C---+----1----+----2----+----3----+----4 Line 63  C---+----1----+----2----+----3----+----4
63        DO m = 1,diagSt_nbFlds(listId)        DO m = 1,diagSt_nbFlds(listId)
64         ndId = jSdiag(m,listId)         ndId = jSdiag(m,listId)
65         parms1 = gdiag(ndId)(1:8)         parms1 = gdiag(ndId)(1:8)
66         IF ( iSdiag(ndId).NE.0 .AND. parms1(5:5).NE.'D' ) THEN         IF ( iSdiag(m,listId).NE.0 .AND. parms1(5:5).NE.'D' ) THEN
67  C--    Start processing 1 Fld :  C--    Start processing 1 Fld :
68    
69           IF ( parms1(5:5).EQ.'C' ) THEN           IF ( parms1(5:5).EQ.'C' ) THEN
# Line 77  C          ----------------------------- Line 77  C          -----------------------------
77    
78           DO j=0,nRegions           DO j=0,nRegions
79            IF ( diagSt_region(j,listId).GT.0 ) THEN            IF ( diagSt_region(j,listId).GT.0 ) THEN
80                iSp = ABS(iSdiag(m,listId))
81                iSm = mSdiag(m,listId)
82              CALL DIAGSTATS_GLOBAL(              CALL DIAGSTATS_GLOBAL(
83       O                       statGlob(0,0,j), tmp_Glob,       O                       statGlob(0,0,j), tmp_Glob,
84       I                       undef, nLev,       I                       undef, nLev, j,
85       I                       j, ndId,mate,myThid)       I                       ndId, mate, iSp, iSm, myThid )
86    
87  C-          Check for empty Diag (= not filled or using empty mask)  C-          Check for empty Diag (= not filled or using empty mask)
88              IF ( statGlob(0,0,j).EQ.0. ) THEN              IF ( statGlob(0,0,j).EQ.0. ) THEN
# Line 108  C-          Check for empty Diag (= not Line 110  C-          Check for empty Diag (= not
110    
111              IF ( debugLevel .GE. debLevA ) THEN              IF ( debugLevel .GE. debLevA ) THEN
112               _BEGIN_MASTER( myThid )               _BEGIN_MASTER( myThid )
113                WRITE(ioUnit,'(A,I3,3A,I3,A,1PE10.3,2A)')                WRITE(ioUnit,'(A,I3,3A,I3,A,1PE10.3,2A)')
114       &         ' Compute Stats, Diag. # ',ndId, '  ', cdiag(ndId),       &         ' Compute Stats, Diag. # ',ndId, '  ', cdiag(ndId),
115       &         '  vol(',j,' ):', statGlob(0,0,j),'  Parms: ',gdiag(ndId)       &         '  vol(',j,' ):', statGlob(0,0,j),'  Parms: ',gdiag(ndId)
116               IF ( mate.GT.0 ) THEN               IF ( mate.GT.0 ) THEN
117                WRITE(ioUnit,'(A,I3,3A,I3,2(A,1PE10.3))')                WRITE(ioUnit,'(A,I3,3A,I3,2(A,1PE10.3))')
118       &         '    use Counter Mate  # ', mate,'  ',cdiag(mate),       &         '    use Counter Mate  # ', mate,'  ',cdiag(mate),
119       &         '  vol(',j,' ):',tmp_Glob(0,0), ' integral',tmp_Glob(1,0)       &         '  vol(',j,' ):',tmp_Glob(0,0), ' integral',tmp_Glob(1,0)
120               ENDIF               ENDIF
# Line 131  C--    end of Processing Fld # m Line 133  C--    end of Processing Fld # m
133         ENDIF         ENDIF
134        ENDDO        ENDDO
135    
136        RETURN        RETURN
137        END        END
138                                                                        
139  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22