/[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.4 by edhill, Sat Dec 18 02:18:55 2004 UTC revision 1.5 by edhill, Wed Dec 29 02:13:38 2004 UTC
# Line 102  C       Update the record dimension by w Line 102  C       Update the record dimension by w
102    
103          CALL MNC_CW_ADD_GNAME('diag_levels', 1,          CALL MNC_CW_ADD_GNAME('diag_levels', 1,
104       &       dim, dn, ib, ie, myThid)       &       dim, dn, ib, ie, myThid)
105          CALL MNC_CW_ADD_VNAME('level_indicies', 'diag_levels',          CALL MNC_CW_ADD_VNAME(dn(1), 'diag_levels',
106       &       0,0, myThid)       &       0,0, myThid)
107          CALL MNC_CW_ADD_VATTR_TEXT('level_indicies','description',          CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
108       &    'Idicies of vertical levels within the data source arrays',       &    'Idicies of vertical levels within the data source arrays',
109       &       myThid)       &       myThid)
110                    
111          CALL MNC_CW_RL_W('I',diag_mnc_bn,0,0,          CALL MNC_CW_RL_W('I',diag_mnc_bn,0,0,
112       &       'level_indicies', levs(1,listnum), myThid)       &       dn(1), levs(1,listnum), myThid)
113    
114          CALL MNC_CW_DEL_VNAME('level_indicies', myThid)          CALL MNC_CW_DEL_VNAME(dn(1), myThid)
115          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)
116        ENDIF        ENDIF
117  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
# Line 230  C         Prepare for mdsio optionality Line 230  C         Prepare for mdsio optionality
230                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)
231              ENDDO              ENDDO
232    
233              dn(1)(1:2) = 'Xd'  C           XY dimensions
234              dim(1) = sNx + 2*OLx              dim(1)       = sNx + 2*OLx
235              ib(1)  = OLx + 1              dim(2)       = sNy + 2*OLy
236              ie(1)  = OLx + sNx              ib(1)        = OLx + 1
237              dn(2)(1:2) = 'Yd'              ib(2)        = OLy + 1
238              dim(2) = sNy + 2*OLy              IF (gdiag(n)(2:2) .EQ. 'M') THEN
239              ib(2)  = OLy + 1                dn(1)(1:2) = 'X'
240              ie(2)  = OLy + sNy                ie(1)      = OLx + sNx
241                  dn(2)(1:2) = 'Y'
242                  ie(2)      = OLy + sNy
243                ELSEIF (gdiag(n)(2:2) .EQ. 'U') THEN
244                  dn(1)(1:3) = 'Xp1'
245                  ie(1)      = OLx + sNx + 1
246                  dn(2)(1:2) = 'Y'
247                  ie(2)      = OLy + sNy
248                ELSEIF (gdiag(n)(2:2) .EQ. 'V') THEN
249                  dn(1)(1:2) = 'X'
250                  ie(1)      = OLx + sNx
251                  dn(2)(1:3) = 'Yp1'
252                  ie(2)      = OLy + sNy + 1
253                ELSEIF (gdiag(n)(2:2) .EQ. 'Z') THEN
254                  dn(1)(1:3) = 'Xp1'
255                  ie(1)      = OLx + sNx + 1
256                  dn(2)(1:3) = 'Yp1'
257                  ie(2)      = OLy + sNy + 1
258                ENDIF
259                
260  C           Z is special since it varies  C           Z is special since it varies
261              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)
262              dim(3) = Nr+Nrphys              dim(3) = Nr+Nrphys
263              ib(3)  = 1              ib(3)  = 1
264              ie(3)  = nlevels(listnum)              ie(3)  = nlevels(listnum)
265    
266              CALL MNC_CW_ADD_GNAME('diag_cw_temp', 3,  C           Time dimension
267                dn(4)(1:1) = 'T'
268                dim(4) = -1
269                ib(4)  = 1
270                ie(4)  = 1
271    
272                CALL MNC_CW_ADD_GNAME('diag_cw_temp', 4,
273       &             dim, dn, ib, ie, myThid)       &             dim, dn, ib, ie, myThid)
274              CALL MNC_CW_ADD_VNAME(cdiag(m), 'diag_cw_temp',              CALL MNC_CW_ADD_VNAME(cdiag(m), 'diag_cw_temp',
275       &             4,5, myThid)       &             4,5, myThid)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22