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

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

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

revision 1.4 by jmc, Fri May 20 07:17:07 2005 UTC revision 1.5 by jmc, Tue Feb 5 15:31:19 2008 UTC
# Line 5  C $Name$ Line 5  C $Name$
5    
6  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7  CBOP  CBOP
8  C     !ROUTINE: DIAGNOSTICS_READ_PICKUP  C     !ROUTINE: DIAGNOSTICS_READ_PICKUP
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE DIAGNOSTICS_READ_PICKUP(        SUBROUTINE DIAGNOSTICS_READ_PICKUP(
11       I     myThid )       I     myThid )
# Line 72  C     Add pickup capability Line 72  C     Add pickup capability
72  C         Update the record dimension by writing the iteration number  C         Update the record dimension by writing the iteration number
73            CALL MNC_FILE_CLOSE_ALL_MATCHING(diag_mnc_bn, myThid)            CALL MNC_FILE_CLOSE_ALL_MATCHING(diag_mnc_bn, myThid)
74            CALL MNC_CW_SET_UDIM(diag_mnc_bn, 1, myThid)            CALL MNC_CW_SET_UDIM(diag_mnc_bn, 1, myThid)
75              
76  C         Read the qdiag() array  C         Read the qdiag() array
77            d_cw_name(1:NLEN) = dn_blnk(1:NLEN)            d_cw_name(1:NLEN) = dn_blnk(1:NLEN)
78            DO ii = 1,CW_DIMS            DO ii = 1,CW_DIMS
# Line 88  C         Read the qdiag() array Line 88  C         Read the qdiag() array
88            ib(2)      = OLy + 1            ib(2)      = OLy + 1
89            ie(2)      = OLy + sNy + 1            ie(2)      = OLy + sNy + 1
90            dn(3)(1:2) = 'Zd'            dn(3)(1:2) = 'Zd'
91            dim(3)     = numdiags            dim(3)     = numDiags
92            ib(3)      = 1            ib(3)      = 1
93            ie(3)      = numdiags            ie(3)      = numDiags
94            dn(4)(1:1) = 'T'            dn(4)(1:1) = 'T'
95            dim(4)     = -1            dim(4)     = -1
96            ib(4)      = 1            ib(4)      = 1
97            ie(4)      = 1            ie(4)      = 1
98            CALL MNC_CW_ADD_GNAME(d_cw_name, 4,            CALL MNC_CW_ADD_GNAME(d_cw_name, 4,
99       &         dim, dn, ib, ie, myThid)       &         dim, dn, ib, ie, myThid)
100            CALL MNC_CW_ADD_VNAME(d_cw_name, d_cw_name,            CALL MNC_CW_ADD_VNAME(d_cw_name, d_cw_name,
101       &         4,5, myThid)       &         4,5, myThid)
102            CALL MNC_CW_RL_R('D',diag_mnc_bn,0,0,            CALL MNC_CW_RL_R('D',diag_mnc_bn,0,0,
103       &         d_cw_name, qdiag, myThid)       &         d_cw_name, qdiag, myThid)
104            CALL MNC_CW_DEL_VNAME(d_cw_name, myThid)            CALL MNC_CW_DEL_VNAME(d_cw_name, myThid)
105            CALL MNC_CW_DEL_GNAME(d_cw_name, myThid)            CALL MNC_CW_DEL_GNAME(d_cw_name, myThid)
106              
107  C         Read the ndiag() array  C         Read the ndiag() array
108            d_cw_name(1:NLEN) = dn_blnk(1:NLEN)            d_cw_name(1:NLEN) = dn_blnk(1:NLEN)
109            DO ii = 1,CW_DIMS            DO ii = 1,CW_DIMS
# Line 111  C         Read the ndiag() array Line 111  C         Read the ndiag() array
111            ENDDO            ENDDO
112            d_cw_name(1:10) = 'diag_count'            d_cw_name(1:10) = 'diag_count'
113            dn(1)(1:2) = 'Nd'            dn(1)(1:2) = 'Nd'
114            dim(1)     = numdiags            dim(1)     = numDiags
115            ib(1)      = 1            ib(1)      = 1
116            ie(1)      = numdiags            ie(1)      = numDiags
117            dn(2)(1:1) = 'T'            dn(2)(1:1) = 'T'
118            dim(2)     = -1            dim(2)     = -1
119            ib(2)      = 1            ib(2)      = 1
120            ie(2)      = 1            ie(2)      = 1
121            CALL MNC_CW_ADD_GNAME(d_cw_name, 2,            CALL MNC_CW_ADD_GNAME(d_cw_name, 2,
122       &         dim, dn, ib, ie, myThid)       &         dim, dn, ib, ie, myThid)
123            CALL MNC_CW_ADD_VNAME(d_cw_name, d_cw_name,            CALL MNC_CW_ADD_VNAME(d_cw_name, d_cw_name,
124       &         4,5, myThid)       &         4,5, myThid)
125            CALL MNC_CW_ADD_VATTR_TEXT(d_cw_name,'description',            CALL MNC_CW_ADD_VATTR_TEXT(d_cw_name,'description',
126       &         'diagnostics state',myThid)       &         'diagnostics state',myThid)
127  C- jmc: get warnings when I compile this S/R because something is not right  C- jmc: get warnings when I compile this S/R because something is not right
128  C       in the type or one or more arguments. commented out for now  C       in the type or one or more arguments. commented out for now
129  c         CALL MNC_CW_RL_R('I',diag_mnc_bn,0,0,  c         CALL MNC_CW_RL_R('I',diag_mnc_bn,0,0,
130  c    &         d_cw_name, ndiag, myThid)  c    &         d_cw_name, ndiag, myThid)
131            CALL MNC_CW_DEL_VNAME(d_cw_name, myThid)            CALL MNC_CW_DEL_VNAME(d_cw_name, myThid)
# Line 133  c    &         d_cw_name, ndiag, myThid) Line 133  c    &         d_cw_name, ndiag, myThid)
133    
134          ENDIF          ENDIF
135  #endif /* ALLOW_MNC */  #endif /* ALLOW_MNC */
136              
137  #ifdef ALLOW_MDSIO  #ifdef ALLOW_MDSIO
138          IF (diag_pickup_read_mdsio) THEN          IF (diag_pickup_read_mdsio) THEN
139            _BEGIN_MASTER(myThid)            _BEGIN_MASTER(myThid)
# Line 144  C         Read qdiag() Line 144  C         Read qdiag()
144            ENDDO            ENDDO
145            write(fn,'(A,I10.10)') 'pickup_qdiag', nIter0            write(fn,'(A,I10.10)') 'pickup_qdiag', nIter0
146            CALL MDSREADFIELD(fn,readBinaryPrec,'RL',            CALL MDSREADFIELD(fn,readBinaryPrec,'RL',
147       &         numdiags,qdiag,1,myThid)       &         numDiags,qdiag,1,myThid)
148    
149  C         Read ndiag()  C         Read ndiag()
150            DO i = 1,80            DO i = 1,80
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22