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

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

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

revision 1.4 by molod, Fri Feb 18 19:44:11 2005 UTC revision 1.6 by jmc, Sun Jun 26 16:51:49 2005 UTC
# Line 74  C we are trying to fill a valid diagnost Line 74  C we are trying to fill a valid diagnost
74        ipointer = 0        ipointer = 0
75        DO n=1,nlists        DO n=1,nlists
76         DO m=1,nActive(n)         DO m=1,nActive(n)
77          IF ( chardiag.EQ.flds(m,n) ) THEN          IF ( chardiag.EQ.flds(m,n) .AND. idiag(m,n).GT.0 ) THEN
78           ndiagnum = jdiag(m,n)           ndiagnum = jdiag(m,n)
79           IF (ndiag(ndiagnum).GE.0) ipointer = idiag(ndiagnum)           ipointer = idiag(m,n)
80          ENDIF           IF ( ndiagnum.NE.0 .AND. ndiag(ipointer,1,1).GE.0 ) THEN
81         ENDDO  C--   do the filling: start here:
82        ENDDO  
83           IF ( (ABS(kLev).LE.1) .AND. (npeice.eq.1) ) THEN
84        print *,' in diag_vegtile_fill for ',chardiag,  C Increment the counter for the diagnostic
85       .     ' ndiagnum ',ndiagnum,' ipointer ',ipointer           ndiag(ipointer,bi,bj) = ndiag(ipointer,bi,bj) + 1
86           ENDIF
 C If-sequence to see if we are a valid and an active diagnostic  
   
       IF ( ndiagnum.NE.0 .AND. ipointer.NE.0 ) THEN  
   
 C Increment the counter for the diagnostic (if we are at bi=bj=myThid=1)  
        _BEGIN_MASTER(myThid)  
        IF((bi.EQ.1).AND.(bj.EQ.1).AND.(ABS(kLev).LE.1)  
      .         .AND.(npeice.eq.1))  
      .                     ndiag(ndiagnum) = ndiag(ndiagnum) + 1  
        _END_MASTER(myThid)  
87    
88         offset = ib*(npeice-1)         offset = ib*(npeice-1)
89         Lena    = min(ib,numpts-offset)         Lena    = min(ib,numpts-offset)
# Line 148  C-      Check for consistency with Nb of Line 138  C-      Check for consistency with Nb of
138    
139         DO k = kFirst,kLast         DO k = kFirst,kLast
140          kd = kd0 + ksgn*k          kd = kd0 + ksgn*k
         print *,' level slot in qdiag= ',kd, ' level from input= ',k  
141          if( check ) then          if( check ) then
142           undef = getcon('UNDEF')           undef = getcon('UNDEF')
143           do i= 1,Lena           do i= 1,Lena
# Line 172  C-      Check for consistency with Nb of Line 161  C-      Check for consistency with Nb of
161          endif          endif
162         ENDDO         ENDDO
163    
164        ELSE  C--   do the filling: ends here.
165  C     IF (myThid.EQ.1) WRITE(6,1000) chardiag           ENDIF
166        ENDIF          ENDIF
167           ENDDO
168          ENDDO
169    
170   1000 format(' ',' Warning: Trying to write to diagnostic ',a8,   1000 format(' ',' Warning: Trying to write to diagnostic ',a8,
171       .        ' But it is not a valid (or active) name ')       .        ' But it is not a valid (or active) name ')
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22