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

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

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

revision 1.39 by mlosch, Tue May 27 08:37:19 2008 UTC revision 1.40 by jmc, Tue Nov 18 21:41:06 2008 UTC
# Line 38  C     myThid  :: my Thread Id number Line 38  C     myThid  :: my Thread Id number
38        INTEGER listId, myIter, myThid        INTEGER listId, myIter, myThid
39  CEOP  CEOP
40    
41    C     !FUNCTIONS:
42          INTEGER ILNBLNK
43          EXTERNAL ILNBLNK
44    #ifdef ALLOW_FIZHI
45          _RL   getcon
46          EXTERNAL getcon
47    #endif
48    
49  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
50  C     i,j,k :: loop indices  C     i,j,k :: loop indices
51  C     lm    :: loop index (averageCycle)  C     lm    :: loop index (averageCycle)
# Line 59  C              diagnostic storage qdiag Line 67  C              diagnostic storage qdiag
67        INTEGER md, ndId, ip, im        INTEGER md, ndId, ip, im
68        INTEGER mate, mVec        INTEGER mate, mVec
69        CHARACTER*10 gcode        CHARACTER*10 gcode
70        _RL undef, getcon        _RL undef
71        _RL tmpLev        _RL tmpLev
       EXTERNAL getcon  
       INTEGER ILNBLNK  
       EXTERNAL ILNBLNK  
72        INTEGER ilen        INTEGER ilen
73    
74        INTEGER ioUnit        INTEGER ioUnit
# Line 98  C              diagnostic storage qdiag Line 103  C              diagnostic storage qdiag
103  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
104    
105        ioUnit= standardMessageUnit        ioUnit= standardMessageUnit
106          undef = UNSET_RL
107    #ifdef ALLOW_FIZHI
108    c     IF ( useFIZHI ) undef = getcon('UNDEF')
109        undef = getcon('UNDEF')        undef = getcon('UNDEF')
110    #endif
111        WRITE(suff,'(I10.10)') myIter        WRITE(suff,'(I10.10)') myIter
112        ilen = ILNBLNK(fnames(listId))        ilen = ILNBLNK(fnames(listId))
113        WRITE( fn, '(A,A,A)' ) fnames(listId)(1:ilen),'.',suff(1:10)        WRITE( fn, '(A,A,A)' ) fnames(listId)(1:ilen),'.',suff(1:10)
# Line 156  C       beginning and ending times for e Line 165  C       beginning and ending times for e
165       &       'Idicies of vertical levels within the source arrays',       &       'Idicies of vertical levels within the source arrays',
166       &       myThid)       &       myThid)
167  C     suppress the missing value attribute (iflag = 0)  C     suppress the missing value attribute (iflag = 0)
168          IF (useMissingValue)          IF (useMissingValue)
169       &       CALL MNC_CW_VATTR_MISSING('diag_levels', 0,       &       CALL MNC_CW_VATTR_MISSING('diag_levels', 0,
170       I       misval_r8, misval_r4, misval_int,       I       misval_r8, misval_r4, misval_int,
171       I       myThid )       I       myThid )
172    
173          CALL MNC_CW_RL_W('D',diag_mnc_bn,0,0,          CALL MNC_CW_RL_W('D',diag_mnc_bn,0,0,
# Line 221  C         for averaged levels. Line 230  C         for averaged levels.
230  C     suppress the missing value attribute (iflag = 0)  C     suppress the missing value attribute (iflag = 0)
231            IF (useMissingValue)            IF (useMissingValue)
232       &         CALL MNC_CW_VATTR_MISSING(dn(1), 0,       &         CALL MNC_CW_VATTR_MISSING(dn(1), 0,
233       I         misval_r8, misval_r4, misval_int,       I         misval_r8, misval_r4, misval_int,
234       I         myThid )       I         myThid )
235            CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)            CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)
236            CALL MNC_CW_DEL_VNAME(dn(1), myThid)            CALL MNC_CW_DEL_VNAME(dn(1), myThid)
# Line 473  C     Missing values only for scalar dia Line 482  C     Missing values only for scalar dia
482              IF ( useMisValForThisDiag ) THEN              IF ( useMisValForThisDiag ) THEN
483  C     assign missing values and set flag for adding the netCDF atttibute  C     assign missing values and set flag for adding the netCDF atttibute
484               CALL MNC_CW_VATTR_MISSING(cdiag(ndId), 2,               CALL MNC_CW_VATTR_MISSING(cdiag(ndId), 2,
485       I            misval_r8, misval_r4, misval_int,       I            misval_r8, misval_r4, misval_int,
486       I            myThid )       I            myThid )
487  C     and now use the missing values for masking out the land points  C     and now use the missing values for masking out the land points
488               DO bj = myByLo(myThid), myByHi(myThid)               DO bj = myByLo(myThid), myByHi(myThid)
# Line 494  C     suppress the missing value attribu Line 503  C     suppress the missing value attribu
503  C     Note: We have to call the following subroutine for each mnc that has  C     Note: We have to call the following subroutine for each mnc that has
504  C     been created "on the fly" by mnc_cw_add_vname and will be deleted  C     been created "on the fly" by mnc_cw_add_vname and will be deleted
505  C     by mnc_cw_del_vname, because all of these variables use the same  C     by mnc_cw_del_vname, because all of these variables use the same
506  C     identifier so that mnc_cw_vfmv(indv) needs to be overwritten for  C     identifier so that mnc_cw_vfmv(indv) needs to be overwritten for
507  C     each of these variables  C     each of these variables
508               CALL MNC_CW_VATTR_MISSING(cdiag(ndId), 0,               CALL MNC_CW_VATTR_MISSING(cdiag(ndId), 0,
509       I            misval_r8, misval_r4, misval_int,       I            misval_r8, misval_r4, misval_int,
510       I            myThid )       I            myThid )
511              ENDIF              ENDIF
512    
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22